Characteristics of the turbulent/nonturbulent interface in boundary layers, jets and shear-free turbulence

The characteristics of turbulent/nonturbulent interfaces (TNTI) from boundary layers, jets and shear-free turbulence are compared using direct numerical simulations. The TNTI location is detected by assessing the volume of turbulent flow as function of the vorticity magnitude and is shown to be equivalent to other procedures using a scalar field. Vorticity maps show that the boundary layer contains a larger range of scales at the interface than in jets and shear-free turbulence where the change in vorticity characteristics across the TNTI is much more dramatic. The intermittency parameter shows that the extent of the intermittency region for jets and boundary layers is similar and is much bigger than in shear-free turbulence, and can be used to compute the vorticity threshold defining the TNTI location. The statistics of the vorticity jump across the TNTI exhibit the imprint of a large range of scales, from the Kolmogorov micro-scale to scales much bigger than the Taylor scale. Finally, it is shown that contrary to the classical view, the low-vorticity spots inside the jet are statistically similar to isotropic turbulence, suggesting that engulfing pockets simply do not exist in jets

Electric current of an electrified jet issuing from a long metallic tube.

This paper presents an analysis of the transport of electric current in a jet of an electrically conducting liquid discharging from a metallic tube into a gas or a vacuum, and subject to an electric field due to a high voltage applied between the tube and a far electrode. The flow, the surface charge and the electric field are computed in the current transfer region of the jet, where conduction current in the liquid becomes surface current due to the convection of electric charge accumulated at its surface. The electric current computed as a function of the flow rate of the liquid injected through the tube increases first as the square root of this flow rate, levels to a nearly constant value when the flow rate is increased and finally sets to a linear increase when the flow rate is further increased. The current increases linearly with the applied voltage at small and moderate values of this variable, and faster than linearly at high voltages. The characteristic length and structure of the current transfer region are determined. Order-of-magnitude estimates for jets which are only weakly stretched by the electric stresses are worked out that qualitatively account for some of the numerical results.

Second-order Lagrangians admitting a first-order Hamiltonian formalism

Let p: E —» JV be an arbitrary fibred manifold over a connected n-dimensional manifold N oriented by a volume form v = dx1^-...^dxn, and let pk: JkE → N be the bundle of K-jets of local sections of p, with projections Plk : JkE → JlE for every k ≥ 1

A new generation of real-time systems in the JET tokamak

Recently a new recipe for developing and deploying real-time systems has become increasingly adopted in the JET tokamak. Powered by the advent of x86 multi-core technology and the reliability of the JET’s well established Real-Time Data Network (RTDN) to handle all real-time I/O, an official Linux vanilla kernel has been demonstrated to be able to provide realtime performance to user-space applications that are required to meet stringent timing constraints. In particular, a careful rearrangement of the Interrupt ReQuests’ (IRQs) affinities together with the kernel’s CPU isolation mechanism allows to obtain either soft or hard real-time behavior depending on the synchronization mechanism adopted. Finally, the Multithreaded Application Real-Time executor (MARTe) framework is used for building applications particularly optimised for exploring multicore architectures. In the past year, four new systems based on this philosophy have been installed and are now part of the JET’s routine operation. The focus of the present work is on the configuration and interconnection of the ingredients that enable these new systems’ real-time capability and on the impact that JET’s distributed real-time architecture has on system engineering requirements, such as algorithm testing and plant commissioning. Details are given about the common real-time configuration and development path of these systems, followed by a brief description of each system together with results regarding their real-time performance. A cycle time jitter analysis of a user-space MARTe based application synchronising over a network is also presented. The goal is to compare its deterministic performance while running on a vanilla and on a Messaging Real time Grid (MRG) Linux kernel.

Advances in global instability computations: from incompressible to hypersonic flow

Esta tesis constituye un gran avance en el conocimiento del estudio y análisis de inestabilidades hidrodinámicas desde un punto de vista físico y teórico, como consecuencia de haber desarrollado innovadoras técnicas para la resolución computacional eficiente y precisa de la parte principal del espectro correspondiente a los problemas de autovalores (EVP) multidimensionales que gobiernan la inestabilidad de flujos con dos o tres direcciones espaciales inhomogéneas, denominados problemas de estabilidad global lineal. En el contexto del trabajo de desarrollo de herramientas computacionales presentado en la tesis, la discretización mediante métodos de diferencias finitas estables de alto orden de los EVP bidimensionales y tridimensionales que se derivan de las ecuaciones de Navier-Stokes linealizadas sobre flujos con dos o tres direcciones espaciales inhomogéneas, ha permitido una aceleración de cuatro órdenes de magnitud en su resolución. Esta mejora de eficiencia numérica se ha conseguido gracias al hecho de que usando estos esquemas de diferencias finitas, técnicas eficientes de resolución de problemas lineales son utilizables, explotando el alto nivel de dispersión o alto número de elementos nulos en las matrices involucradas en los problemas tratados. Como más notable consecuencia cabe destacar que la resolución de EVPs multidimensionales de inestabilidad global, que hasta la fecha necesitaban de superordenadores, se ha podido realizar en ordenadores de sobremesa. Además de la solución de problemas de estabilidad global lineal, el mencionado desarrollo numérico facilitó la extensión de las ecuaciones de estabilidad parabolizadas (PSE) lineales y no lineales para analizar la inestabilidad de flujos que dependen fuertemente en dos direcciones espaciales y suavemente en la tercera con las ecuaciones de estabilidad parabolizadas tridimensionales (PSE-3D). Precisamente la capacidad de extensión del novedoso algoritmo PSE-3D para el estudio de interacciones no lineales de los modos de estabilidad, desarrollado íntegramente en esta tesis, permite la predicción de transición en flujos complejos de gran interés industrial y por lo tanto extiende el concepto clásico de PSE, el cuál ha sido empleado exitosamente durante las pasadas tres décadas en el mismo contexto para problemas de capa límite bidimensional. Típicos ejemplos de flujos incompresibles se han analizado en este trabajo sin la necesidad de recurrir a restrictivas presuposiciones usadas en el pasado. Se han estudiado problemas vorticales como es el caso de un vórtice aislado o sistemas de vórtices simulando la estela de alas, en los que la homogeneidad axial no se impone y así se puede considerar la difusión viscosa del flujo. Además, se ha estudiado el chorro giratorio turbulento, cuya inestabilidad se utiliza para mejorar las características de funcionamiento de combustores. En la tesis se abarcan adicionalmente problemas de flujos compresibles. Se presenta el estudio de inestabilidad de flujos de borde de ataque a diferentes velocidades de vuelo. También se analiza la estela formada por un elemento rugoso aislado en capa límite supersónica e hipersónica, mostrando excelentes comparaciones con resultados obtenidos mediante simulación numérica directa. Finalmente, nuevas inestabilidades se han identificado en el flujo hipersónico a Mach 7 alrededor de un cono elíptico que modela el vehículo de pruebas en vuelo HIFiRE-5. Los resultados comparan favorablemente con experimentos en vuelo, lo que subraya aún más el potencial de las metodologías de análisis de estabilidad desarrolladas en esta tesis. ABSTRACT The present thesis constitutes a step forward in advancing the frontiers of knowledge of fluid flow instability from a physical point of view, as a consequence of having been successful in developing groundbreaking methodologies for the efficient and accurate computation of the leading part of the spectrum pertinent to multi-dimensional eigenvalue problems (EVP) governing instability of flows with two or three inhomogeneous spatial directions. In the context of the numerical work presented in this thesis, the discretization of the spatial operator resulting from linearization of the Navier-Stokes equations around flows with two or three inhomogeneous spatial directions by variable-high-order stable finite-difference methods has permitted a speedup of four orders of magnitude in the solution of the corresponding two- and three-dimensional EVPs. This improvement of numerical performance has been achieved thanks to the high-sparsity level offered by the high-order finite-difference schemes employed for the discretization of the operators. This permitted use of efficient sparse linear algebra techniques without sacrificing accuracy and, consequently, solutions being obtained on typical workstations, as opposed to the previously employed supercomputers. Besides solution of the two- and three-dimensional EVPs of global linear instability, this development paved the way for the extension of the (linear and nonlinear) Parabolized Stability Equations (PSE) to analyze instability of flows which depend in a strongly-coupled inhomogeneous manner on two spatial directions and weakly on the third. Precisely the extensibility of the novel PSE-3D algorithm developed in the framework of the present thesis to study nonlinear flow instability permits transition prediction in flows of industrial interest, thus extending the classic PSE concept which has been successfully employed in the same context to boundary-layer type of flows over the last three decades. Typical examples of incompressible flows, the instability of which was analyzed in the present thesis without the need to resort to the restrictive assumptions used in the past, range from isolated vortices, and systems thereof, in which axial homogeneity is relaxed to consider viscous diffusion, as well as turbulent swirling jets, the instability of which is exploited in order to improve flame-holding properties of combustors. The instability of compressible subsonic and supersonic leading edge flows has been solved, and the wake of an isolated roughness element in a supersonic and hypersonic boundary-layer has also been analyzed with respect to its instability: excellent agreement with direct numerical simulation results has been obtained in all cases. Finally, instability analysis of Mach number 7 ow around an elliptic cone modeling the HIFiRE-5 flight test vehicle has unraveled flow instabilities near the minor-axis centerline, results comparing favorably with flight test predictions.

Numerical analysis of deflagration initation by a hot jet

Numerical simulations of axisymmetric reactive jets with one-step Arrhenius kinetics are used to investigate the problem of deflagration initiation in a premixed fuel–air mixture by the sudden discharge of a hot jet of its adiabatic reaction products. For the moderately large values of the jet Reynolds number considered in the computations, chemical reaction is seen to occur initially in the thin mixing layer that separates the hot products from the cold reactants. This mixing layer is wrapped around by the starting vortex, thereby enhancing mixing at the jet head, which is followed by an annular mixing layer that trails behind, connecting the leading vortex with the orifice rim. A successful deflagration is seen to develop for values of the orifice radius larger than a critical value a c in the order of the flame thickness of the planar deflagration δL. Introduction of appropriate scales provides the dimensionless formulation of the problem, with flame initiation characterised in terms of a critical Damköhler number Δc=(a d/δL)2, whose parametric dependence is investigated. The numerical computations reveal that, while the jet Reynolds number exerts a limited influence on the criticality conditions, the effect of the reactant diffusivity on ignition is much more pronounced, with the value of Δc increasing significantly with increasing Lewis numbers. The reactant diffusivity affects also the way ignition takes place, so that for reactants with the flame develops as a result of ignition in the annular mixing layer surrounding the developing jet stem, whereas for highly diffusive reactants with Lewis numbers sufficiently smaller than unity combustion is initiated in the mixed core formed around the starting vortex. The analysis provides increased understanding of deflagration initiation processes, including the effects of differential diffusion, and points to the need for further investigations corporating detailed chemistry models for specific fuel–air mixtures.

Variable-density jet flows induced by concentrated sources of momentum and energy

The planar and axisymmetric variable-density flows induced in a quiescent gas by a concentrated source of momentum that is simultaneously either a source or a sink of energy are investigated for application to the description of the velocity and temperature far fields in laminar gaseous jets with either large or small values of the initial jet-to-ambient temperature ratio. The source fluxes of momentum and heat are used to construct the characteristic scales of velocity and length in the region where the density differences are of the order of the ambient density, which is slender for the large values of the Reynolds number considered herein. The problem reduces to the integration of the dimensionless boundary-layer conservation equations, giving a solution that depends on the gas transport properties but is otherwise free of parameters. The boundary conditions at the jet exit for integration are obtained by analysing the self-similar flow that appears near the heat source in planar and axisymmetric configurations and also near the heat sink in the planar case. Numerical integrations of the boundary-layer equations with these conditions give solutions that describe accurately the velocity and temperature fields of very hot planar and round jets and also of very cold plane jets in the far field region where the density and temperature differences are comparable to the ambient values. Simple scaling arguments indicate that the point source description does not apply, however, to cold round jets, whose far field region is not large compared with the jet development region, as verified by numerical integrations

One-dimensional linear analysis of the compound jet

The stability of an infinitely long compound liquid column is analysed by using a one-dimensional inviscid slice model. Results obtained from this one-dimensional linear analysis are applicable to the study of compound capillary jets, which are used in the ink-jet printing technique. Stability limits and the breaking regimes of such fluid configurations are established, and, whenever possible, theoretical results are compared with experimental ones.

Desarrollo de un método numérico para flujos multimaterial en 2D en ARWEN

En esta tesis se presenta un método numérico para resolver las ecuaciones de Euler para flujos multimaterial en malla euleriana. Este solver se ha acoplado en el código hidrodinámico en dos dimensiones con transporte de radiación desarrollado en el Instituto de Fusión Nuclear de la UPM bajo la dirección del profesor Pedro Velarde, ARWEN. Los objetivos de este trabajo son: Desarrollo e implementación de un método de Godunov unsplit de alto orden multimaterial en 2D para malla euleriana en geometría cartesiana y geometría cilíndrica. Se presenta una extensión del trabajo realizado por Miller y Puckett (36) a una formulación unsplit. Además, se ha prestado especial atención al acoplamiento con el transporte de radiación y la conducción de calor. El método presentado se ha probado en una gran cantidad de problemas. Aplicación del código multimaterial al estudio de experimentos reales: • Simulación de una propuesta de experimento de laboratorio para reproducir la etapa de arrancamiento de material de la interacción entre el gas proveniente de la explosión de una supernova y la estrella secundaria en un escenario degenarado (SD). • Formación de jets en el laboratorio producidos por la colisión de dos plasmas. ABSTRACT We present a solver for the Euler equations for multimaterial flows in eulerian mesh. This solver has been coupled in the 2D AMR radiation transport code developed at Instituto de Fusión Nuclear (UPM) under the direction of professor Pedro Velarde, ARWEN. The main goals of this thesis are: Development and implementation of an 2D unsplit high-order Godunov method for multimaterial flows in eulerian mesh for cartesian and axialsimetry geometry. We present an extension of the work of Miller and Puckett (36) to an unsplit formulation. Also, we have paid special attention to the coupling with radiation transport and heat conduction. The method has been tested in a wide variety of problems. Application of the multimaterial solver to the study of real experiments: • Simulation of a proposal of a laboratory experiment aimed to reproducing the stripping stage of the interaction between the gas ejected during a supernova explosion and the secondary star in the Single Degenerate scenario. • Experiments of plasma jets in the laboratory obtained by the collission of two hot plasmas.

Energetic performance of landfill and digester biogas in a domestic cooker

The energetic performance of landfill biogas (LB) and biodigester biogas (BB) from municipal waste was examined in consumption tests. These tests were performed in situ at a gas generation plant associated with a landfill facility in Madrid (Spain) and following the standard UNE-EN 30-2-1 (1999). The jets of a domestic cooker commonly used for natural gas (NG) or liquefied petroleum gas (LPG) were modified to operate with the biogases produced at the facility. The working pressures best suited to the tested gases, i.e., to avoid flashback and flame lift, and to ensure the stability and correct functioning of the flame during combustion, were determined by trial and error. Both biogases returned optimum energetic performance for the transfer of heat to water in a metallic recipient (as required by the above standard) at a supply pressure of 10 mbar. Domestic cookers are normally supplied with NG at a pressure of 20 mbar, at which pressure the energetic performance of G20 reference gas was higher than that of both biogases (52.84% compared to 38.06% and 49.77% respectively). Data concerning these issues involving also unexplored feedstock are required for the correct conversions of domestic cookers in order to avoid risks of serious personal injuries or property damages.

Strategies for Minimum Distance in Simulated Evacuation of Transport Airplanes

The present paper describes the preliminary stages of the development of a new, comprehensive model conceived to simulate the evacuation of transport airplanes in certification studies. Two previous steps were devoted to implementing an efficient procedure to define the whole geometry of the cabin, and setting up an algorithm for assigning seats to available exits. Now, to clarify the role of the cabin arrangement in the evacuation process, the paper addresses the influence of several restrictions on the seat-to-exit assignment algorithm, maintaining a purely geometrical approach for consistency. Four situations are considered: first, an assignment method without limitations to search the minimum for the total distance run by all passengers along their escaping paths; second, a protocol that restricts the number of evacuees through each exit according to updated FAR 25 capacity; third, a procedure which tends to the best proportional sharing among exits but obliges to each passenger to egress through the nearest fore or rear exits; and fourth, a scenario which includes both restrictions. The four assignment strategies are applied to turboprops, and narrow body and wide body jets. Seat to exit distance and number of evacuees per exit are the main output variables. The results show the influence of airplane size and the impact of non-symmetries and inappropriate matching between size and longitudinal location of exits.