Parameterization of the atmospheric boundary layer for offshore wind resource assessment with a limited length-scale k-ε model

The structure of the atmospheric boundary layer (ABL) is modelled with the limited- length-scale k-ε model of Apsley and Castro. Contrary to the standard k-ε model, the limited-length-scale k-ε model imposes a maximum mixing length which is derived from the boundary layer height, for neutral and unstable atmospheric situations, or by Monin-Obukhov length when the atmosphere is stably stratified. The model is first verified reproducing the famous Leipzig wind profile. Then the performance of the model is tested with measurements from FINO-1 platform using sonic anemometers to derive the appropriate maximum mixing length.

Inverse bremsstrahlung absorption in spherical laser targets

Inverse bremsstrahlung has been incorporated into an analytical model of the expanding corona of a laser-irradiated spherical target. Absorption decreases slowly with increasing intensity, in agreement with some numerical simulations, and contrary to estimates from simple models in use up to now, which are optimistic at low values of intensity and very pessimistic at high values. Present results agree well with experimental data from many laboratories; substantial absorption is found up to moderate intensities,say below IOl5 W cm-2 for 1.06 pm light. Anomalous absorption, wher, included in the analysis, leaves practically unaffected the ablation pressure and mass ablation rate, for given absorbed intensity. Universal results are given in dimensionless fom.

Current collection by an active spherical electrode in an unmagnetized plasma

A theoretical model for the steady-state response of anodic contactors that emit a plasma current Ii and collect electrons from a collisionless, unmagnetized plasma is presented. The use of a (kinetic) monoenergetic population for the attracted species, well known in passive probe theory, gives both accuracy and tractability to the theory. The monoenergetic population is proved to behave like an isentropic fluid with radial plus centripetal motion, allowing direct comparisons with ad hoc fluid models. Also, a modification of the original monoenergetic equations permits analysis of contactors operating in orbit-limited conditions. Besides that, the theory predicts that, only for plasma emissions above certain threshold current a presheath/double layer/core structure for the potential is formed (the core mode), while for emissions below that threshold, a plasma contactor behaves exactly as a positive-ion emitter with a presheath/sheath structure (the no-core mode). Ion emitters are studied as a particular case. Emphasis is placed on obtaining dimensionless charts and approximate asymptotic laws of the current-voltage characteristic.

Non-periodic driving of coupled oscillators:a spherical swing

Nonlinearly coupled, damped oscillators at 1:1 frequency ratio, one oscillator being driven coherently for efficient excitation, are exemplified by a spherical swing with some phase-mismatch between drive and response. For certain damping range, excitation is found to succeed if it lags behind, but to produce a chaotic attractor if it leads the response. Although a period-doubhng sequence, for damping increasing, leads to the attractor, this is actually born as a hard (as regards amplitude) bifurcation at a zero growth-rate parametric line; as damping decreases, an unstable fixed point crosses an invariant plane to enter as saddle-focus a phase-space domain of physical solutions. A second hard bifurcation occurs at the zero mismatch line, the saddle-focus leaving that domain. Times on the attractor diverge when approaching either fine, leading to exactly one-dimensional and noninvertible limit maps, which are analytically determined.

Spherical collectors versus bare tethers for drag, thrust, and power generation

Deorbit, power generation, and thrusting performances of a bare thin-tape tether and an insulated tether with a spherical electron collector are compared for typical conditions in low-Earth orbit and common values of length L = 4−20 km and cross-sectional area of the tether A = 1−5 mm2. The relative performance of moderately large spheres, as compared with bare tapes, improves but still lags as one moves from deorbiting to power generation and to thrusting: Maximum drag in deorbiting requires maximum current and, thus, fully reflects on anodic collection capability, whereas extracting power at a load or using a supply to push current against the motional field requires reduced currents. The relative performance also improves as one moves to smaller A, which makes the sphere approach the limiting short-circuit current, and at greater L, with the higher bias only affecting moderately the already large bare-tape current. For a 4-m-diameter sphere, relative performances range from 0.09 sphere-to-bare tether drag ratio for L = 4 km and A = 5 mm2 to 0.82 thrust–efficiency ratio for L = 20 km and A = 1 mm2. Extremely large spheres collecting the short-circuit current at zero bias at daytime (diameters being about 14 m for A = 1 mm2 and 31 m for A = 5 mm2) barely outperform the bare tape for L = 4 km and are still outperformed by the bare tape for L = 20 km in both deorbiting and power generation; these large spheres perform like the bare tape in thrusting. In no case was sphere or sphere-related hardware taken into account in evaluating system mass, which would have reduced the sphere performances even further.

Spherical collectors versus bare tethers for drag, thrust, and power generation

Performances of ED-tethers using either spherical collectors or bare tethers for drag, thrust, or power generation, are compared. The standard Parker-Murphy model of current to a full sphere, with neither space-charge nor plasmamotion effects considered, but modified to best fit TSS1R results, is used (the Lam, Al'pert/Gurevich space-charge limited model will be used elsewhere) In the analysis, the spherical collector is assumed to collect current well beyond its random-current value (thick-heath). Both average current in the bare-tether and current to the sphere are normalized with the short-circuit current in the absence of applied power, allowing a comparison of performances for all three applications in terms of characteristic dimensionless numbers. The sphere is always substantially outperformed by the bare-tether if ohmic effects are weak, though its performance improves as such effects increase.

Experimental evidence for super-resolution in the spherical geodesic waveguide

The previous publications (Miñano et al, 2011 and Gonzalez et al, 2012) have shown that using a Spherical Geodesic Waveguide (SGW) it can be achieved the super-resolution up to λ/3000, which is far below the classic Abbe diffraction limit, close to a set of discrete microwave frequencies. The SGW was designed and simulated in COMSOL as a thin geodesic waveguide bounded by an ideal and lossless metal. Herein we present the experimental results for a manufactured SGW, slightly modified due to fabrication requirements, showing the super-resolution up to λ/105.

Current-voltage response of a spherical plasma contactor

A theoretical model for a contactor, collecting electrons from an ambient, unmagnetized plasma and emitting a current Iiis discussed. The relation between Ii and the potential bias of the contactor is found to be crucial for the formation of a quasineutral core around the anode and, consequently, for the current colleted. Approximate analytical laws and charts for the current-voltage response are provided.

Inverse kinematics of a humanoid robot with non-spherical hip: a hybrid algorithm approach

This paper describes an approach to solve the inverse kinematics problem of humanoid robots whose construction shows a small but non negligible offset at the hip which prevents any purely analytical solution to be developed. Knowing that a purely numerical solution is not feasible due to variable efficiency problems, the proposed one first neglects the offset presence in order to obtain an approximate “solution” by means of an analytical algorithm based on screw theory, and then uses it as the initial condition of a numerical refining procedure based on the Levenberg‐Marquardt algorithm. In this way, few iterations are needed for any specified attitude, making it possible to implement the algorithm for real‐time applications. As a way to show the algorithm’s implementation, one case of study is considered throughout the paper, represented by the SILO2 humanoid robot.

NURBS-Based Geometry Parameterization for Aerodynamic Shape Optimization

La motivación de esta tesis es el desarrollo de una herramienta de optimización automática para la mejora del rendimiento de formas aerodinámicas enfocado en la industria aeronáutica. Este trabajo cubre varios aspectos esenciales, desde el empleo de Non-Uniform Rational B-Splines (NURBS), al cálculo de gradientes utilizando la metodología del adjunto continuo, el uso de b-splines volumétricas como parámetros de diseño, el tratamiento de la malla en las intersecciones, y no menos importante, la adaptación de los algoritmos de la dinámica de fluidos computacional (CFD) en arquitecturas hardware de alto paralelismo, como las tarjetas gráficas, para acelerar el proceso de optimización. La metodología adjunta ha posibilitado que los métodos de optimización basados en gradientes sean una alternativa prometedora para la mejora de la eficiencia aerodinámica de los aviones. La formulación del adjunto permite calcular los gradientes de una función de coste, como la resistencia aerodinámica o la sustentación, independientemente del número de variables de diseño, a un coste computacional equivalente a una simulación CFD. Sin embargo, existen problemas prácticos que han imposibilitado su aplicación en la industria, que se pueden resumir en: integrabilidad, rendimiento computacional y robustez de la solución adjunta. Este trabajo aborda estas contrariedades y las analiza en casos prácticos. Como resumen, las contribuciones de esta tesis son: • El uso de NURBS como variables de diseño en un bucle de automático de optimización, aplicado a la mejora del rendimiento aerodinámico de alas en régimen transónico. • El desarrollo de algoritmos de inversión de punto, para calcular las coordenadas paramétricas de las coordenadas espaciales, para ligar los vértices de malla a las NURBS. • El uso y validación de la formulación adjunta para el calculo de los gradientes, a partir de las sensibilidades de la solución adjunta, comparado con diferencias finitas. • Se ofrece una estrategia para utilizar la geometría CAD, en forma de parches NURBS, para tratar las intersecciones, como el ala-fuselaje. • No existen muchas alternativas de librerías NURBS viables. En este trabajo se ha desarrollado una librería, DOMINO NURBS, y se ofrece a la comunidad como código libre y abierto. • También se ha implementado un código CFD en tarjeta gráfica, para realizar una valoración de cómo se puede adaptar un código sobre malla no estructurada a arquitecturas paralelas. • Finalmente, se propone una metodología, basada en la función de Green, como una forma eficiente de paralelizar simulaciones numéricas. Esta tesis ha sido apoyada por las actividades realizadas por el Área de Dinámica da Fluidos del Instituto Nacional de Técnica Aeroespacial (INTA), a través de numerosos proyectos de financiación nacional: DOMINO, SIMUMAT, y CORESFMULAERO. También ha estado en consonancia con las actividades realizadas por el departamento de Métodos y Herramientas de Airbus España y con el grupo Investigación y Tecnología Aeronáutica Europeo (GARTEUR), AG/52. ABSTRACT The motivation of this work is the development of an automatic optimization strategy for large scale shape optimization problems that arise in the aeronautics industry to improve the aerodynamic performance; covering several aspects from the use of Non-Uniform Rational B-Splines (NURBS), the calculation of the gradients with the continuous adjoint formulation, the development of volumetric b-splines parameterization, mesh adaptation and intersection handling, to the adaptation of Computational Fluid Dynamics (CFD) algorithms to take advantage of highly parallel architectures in order to speed up the optimization process. With the development of the adjoint formulation, gradient-based methods for aerodynamic optimization become a promising approach to improve the aerodynamic performance of aircraft designs. The adjoint methodology allows the evaluation the gradients to all design variables of a cost function, such as drag or lift, at the equivalent cost of more or less one CFD simulation. However, some practical problems have been delaying its full implementation to the industry, which can be summarized as: integrability, computer performance, and adjoint robustness. This work tackles some of these issues and analyse them in well-known test cases. As summary, the contributions comprises: • The employment of NURBS as design variables in an automatic optimization loop for the improvement of the aerodynamic performance of aircraft wings in transonic regimen. • The development of point inversion algorithms to calculate the NURBS parametric coordinates from the space coordinates, to link with the computational grid vertex. • The use and validation of the adjoint formulation to calculate the gradients from the surface sensitivities in an automatic optimization loop and evaluate its reliability, compared with finite differences. • This work proposes some algorithms that take advantage of the underlying CAD geometry description, in the form of NURBS patches, to handle intersections and mesh adaptations. • There are not many usable libraries for NURBS available. In this work an open source library DOMINO NURBS has been developed and is offered to the community as free, open source code. • The implementation of a transonic CFD solver from scratch in a graphic card, for an assessment of the implementability of conventional CFD solvers for unstructured grids to highly parallel architectures. • Finally, this research proposes the use of the Green's function as an efficient paralellization scheme of numerical solvers. The presented work has been supported by the activities carried out at the Fluid Dynamics branch of the National Institute for Aerospace Technology (INTA) through national founding research projects: DOMINO, SIMUMAT, and CORESIMULAERO; in line with the activities carried out by the Methods and Tools and Flight Physics department at Airbus and the Group for Aeronautical Research and Technology in Europe (GARTEUR) action group AG/52.