Direct Numerical Simulation of Multi-Injection Mixing and Ignition in Diesel Engine Environments.
Abstract not provided.
Rieth, Martin↗
SEARCH · Search NASA
Results for “numerical mixing”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A detailed numerical model of a supersonic reacting mixing layer
A current research effort is underway at the NASA Langley Research Center to achieve a detailed understanding of important phenomena present when a supersonic flow undergoes chemical reaction. A computer program has been developed to study the details of such flows. The program has been constructed to consider the multicomponent diffusion and convection of important species, the finite-rate reaction of these species, and the resulting interaction between the fluid mechanics and chemistry. Code results from the analysis of a spatially developing and reacting mixing layer are presented, and conclusions are drawn regarding the structure of the evolving layer and its associated flame.
Drummond, J. P.↗
Aspects of supercritical turbulence: direct numerical simulation of an Osub2/Hsub2 temporal mixing layer
Explore the source record for details and available documents.
supercritical turbulence↗
Consistent large eddy simulation of a temporal mixing layer laden with evaporating drops; part 1: direct numerical simulation and formulation
Direct Numerical Simulation is conducted of a three-dimensional temporal mixing layer whose lower stream is initially laden with liquid drops which may evaporate during the simulation.
Large Eddy Simulation LES Direct Numerical Simulat↗
Geometry and dynamics of passive scalar structures in compressible turbulent mixing
We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 , and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Use of passive scalar tagging for the study of coherent structures in the plane mixing layer
Data obtained from the numerical simulation of a 2-D mixing layer were used to study the feasibility of using the instantaneous concentration of a passive scalar for detecting the typical coherent structures in the flow. The study showed that this technique works quite satisfactorily and yields results similar to those that can be obtained by using the instantaneous vorticity for structure detection. Using the coherent events educed by the scalar conditioning technique, the contribution of the coherent events to the total turbulent momentum and scalar transport was estimated. It is found that the contribution from the typical coherent events is of the same order as that of the time-mean value. However, the individual contributions become very large during the pairing of these structures. The increase is particularly spectacular in the case of the Reynolds shear stress.
Ramaprian, B. R.↗
NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE
In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.
detailed chemical kinetics↗
Numerical solution of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations
A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.
Hirsh, R. S.↗
A numerical study of the effects of heat diffusion through the base of the mixed layer
In the present paper, a simple numerical model is used to study the warming of the mixed layer during the early summer. It is shown that the springtime temperature increase in the layer below the mixed layer (for example, in the cold pool on a continental shelf) has a maximum value which occurs for a limiting value of the surface heat flux. This is a result of the positive feedback at large Richardson numbers between stability and vertical diffusion of heat. The springtime temperature increase in the mixed layer increases nonlinearly with surface heat flux, because of the same positive feedback. The effects of interseasonal fluctuations of the surface heat flux on the spring and summer mixed layer and deeper temperature increases can be as great as the effect of interseasonal fluctuations of the average heat flux.
Posmentier, E. S.↗
Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method
Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.
42 ENGINEERING↗
Parametric Study of a Mixer/Ejector Nozzle with Mixing Enhancement Devices
A numerical study employing a simplified model of the High Speed Civil Transport mixer/ejector nozzle has been conducted to investigate the effect of tabs (vortex generators) on the mixing process. More complete mixing of the primary and secondary flows within the confined ejector lowers peak exit velocity resulting in reduced jet noise. Tabs were modeled as vortex pairs and inserted into the computational model. The location, size, and number of tabs were varied and its effect on the mixing process is presented here both quantitatively and qualitatively. A baseline case (no tabs) along with six other cases involving two different vortex strengths at three different orientations have been computed and analyzed. The case with the highest vorticity (six vortices representing large tabs) gives the best mixing. It is shown that the influence of the vorticity acts primarily in the forward or middle portions of the duct, significantly alters the flow structure, and promotes some mixing in the lateral direction. Unmixed pockets were found at the top and bottom of the lobe, and more clever placement of tabs improved mixing in the vertical direction. The technique of replacing tabs with vortices shows promise as an efficient tool for quickly optimizing tab placement in lobed mixers.
DalBello, T.↗
Evolution of solar wind fluctuations and the influence of turbulent 'mixing'
We present various numerical and analytical solutions for the transport of solar wind turbulence. The model used takes into account the effects of convection, expansion, and wave propagation, as well as the recently illuminated effects of (non-WKB) 'mixing' terms. The radial evolution of the fluctuating kinetic energy, magnetic energy, and normalized cross helicity is computed, and, it is demonstrated that in appropriate limits the solutions converge to the WKB forms. In the general case, solutions which differ substantially from those predicted by WKB theory are obtained. The degree of turbulent 'mixing' shows considerable dependence on the nature of the turbulence, giving rise to varying levels, at 1 AU, of the ratio of 'inward' and 'outward' fluctuation energies and the ratio of kinetic to magnetic energies in the fluctuations. The transport properties described here may provide at least a partial explanation for the observed mixing of cross helicities with increasing heliocentric distance in the solar wind.
Oughton, S.↗
A numerical study of wave propagation in a confined mixing layer by eigenfunction expansions
It is well known that the growth rate of instability waves of a two-dimensional free shear layer is reduced greatly at supersonic convective Mach numbers. In previous works, it has been shown that new wave modes exist when the shear layers are bounded by a channel due to the coupling effect between the acoustic wave modes and the motion of the mixing layer. The present work studies the simultaneous propagation of multiple stability waves using numerical simulation. It is shown here that the coexistence of two wave modes in the flow field can lead to an oscillatory growth of disturbance energy with each individual wave mode propagating linearly. This is particularly important when the growth rates of the unstable waves are small. It is also shown here that the propagation of two neutrally stable wave modes can lead to a stationary periodic structure of rms fluctuations. In the numerical simulations presented here the forced wave modes are propagating at same frequency, but with different phase velocities. In order to track the growth of each wave mode as it propagates downstream, a numerical method that can effectively detect and separate the contribution of the individual wave is given. It is demonstrated that by a least square fitting of the disturbance field with eigenfunctions the amplitude of each wave mode can be found. Satisfactory results as compared to linear theory are obtained.
Hu, Fang Q.↗
A two-dimensional numerical simulation of a supersonic, chemically reacting mixing layer
Research has been undertaken to achieve an improved understanding of physical phenomena present when a supersonic flow undergoes chemical reaction. A detailed understanding of supersonic reacting flows is necessary to successfully develop advanced propulsion systems now planned for use late in this century and beyond. In order to explore such flows, a study was begun to create appropriate physical models for describing supersonic combustion, and to develop accurate and efficient numerical techniques for solving the governing equations that result from these models. From this work, two computer programs were written to study reacting flows. Both programs were constructed to consider the multicomponent diffusion and convection of important chemical species, the finite rate reaction of these species, and the resulting interaction of the fluid mechanics and the chemistry. The first program employed a finite difference scheme for integrating the governing equations, whereas the second used a hybrid Chebyshev pseudospectral technique for improved accuracy.
Drummond, J. Philip↗
Co-Optimization of Velocity and Charge-Depletion for Plug-In Hybrid Electric Vehicles: Accounting for Acceleration and Jerk Constraints
Abstract Recent advances in vehicle connectivity and automation technologies promote advanced control algorithms that co-optimize the longitudinal dynamics and powertrain operation of hybrid electric vehicles. Typically, a sequential optimization with the vehicle dynamics optimized followed by powertrain optimization is adopted to manage a number of complexities such as the inherent mixed-integer nature of the hybrid powertrain, the numerous state and control variables, the differing time scales of vehicle and powertrain subsystems, time-varying state constraints, and large horizon lengths. Instead, we solve the offline optimization problem in a centralize manner assuming exact knowledge of the lead vehicle's position over the entire trip by applying a discrete-time single shooting-based numerical approach, Discrete Mixed-Integer Shooting (DMIS), including a linearly increasing computational complexity to the problem horizon. In particular, the hierarchical problem structure is exploited to decompose the computationally intensive Hamiltonian minimization step into a set of low-dimensional optimizations. DMIS allows us to compute the direct fuel minimization problem including the vehicle and powertrain dynamics in a centralized manner to its full horizon while systematically tuning weighting factors that penalize passenger discomfort. For the first time, this study reveals that practically implemented sequential optimization exhibits similar fuel optimality as co-optimization when a certain level of passenger comfort is required.
Automation & Control Systems↗
Evaluation of a hybrid kinetics/mixing-controlled combustion model for turbulent premixed and diffusion combustion using KIVA-2
Two dimensional calculations were made of spark ignited premixed-charge combustion and direct injection stratified-charge combustion in gasoline fueled piston engines. Results are obtained using kinetic-controlled combustion submodel governed by a four-step global chemical reaction or a hybrid laminar kinetics/mixing-controlled combustion submodel that accounts for laminar kinetics and turbulent mixing effects. The numerical solutions are obtained by using KIVA-2 computer code which uses a kinetic-controlled combustion submodel governed by a four-step global chemical reaction (i.e., it assumes that the mixing time is smaller than the chemistry). A hybrid laminar/mixing-controlled combustion submodel was implemented into KIVA-2. In this model, chemical species approach their thermodynamics equilibrium with a rate that is a combination of the turbulent-mixing time and the chemical-kinetics time. The combination is formed in such a way that the longer of the two times has more influence on the conversion rate and the energy release. An additional element of the model is that the laminar-flame kinetics strongly influence the early flame development following ignition.
Nguyen, H. Lee↗
Evaluation of a hybrid kinetics/mixing-controlled combustion model for turbulent premixed and diffusion combustion using KIVA-II
Two-dimensional calculations were made of spark ignited premixed-charge combustion and direct injection stratified-charge combustion in gasoline fueled piston engines. Results are obtained using kinetic-controlled combustion submodel governed by a four-step global chemical reaction or a hybrid laminar kinetics/mixing-controlled combustion submodel that accounts for laminar kinetics and turbulent mixing effects. The numerical solutions are obtained by using KIVA-2 computer code which uses a kinetic-controlled combustion submodel governed by a four-step global chemical reaction (i.e., it assumes that the mixing time is smaller than the chemistry). A hybrid laminar/mixing-controlled combustion submodel was implemented into KIVA-2. In this model, chemical species approach their thermodynamics equilibrium with a rate that is a combination of the turbulent-mixing time and the chemical-kinetics time. The combination is formed in such a way that the longer of the two times has more influence on the conversion rate and the energy release. An additional element of the model is that the laminar-flame kinetics strongly influence the early flame development following ignition.
Nguyen, H. Lee↗
Shock-driven three-fluid mixing with various chevron interface configurations
When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗