Direct Numerical Simulation of Conventional and Alternative Jet Fuel Mixing and Ignition in Compression Ignition 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.
Direct Numerical Simulation of Conventional and Alternative Jet Fuel Mixing and Ignition in Compression Ignition Engine Environments.
Abstract not provided.
Temme, Jacob↗
Sensitivity of the numerical analysis of the three-fluid plasma mixed initial-boundary value problem.
Explore the source record for details and available documents.
Johnson, B. H.↗
A numerical method for systems of conservation laws of mixed type admitting hyperbolic flux splitting
The present treatment of elliptic regions via hyperbolic flux-splitting and high order methods proposes a flux splitting in which the corresponding Jacobians have real and positive/negative eigenvalues. While resembling the flux splitting used in hyperbolic systems, the present generalization of such splitting to elliptic regions allows the handling of mixed-type systems in a unified and heuristically stable fashion. The van der Waals fluid-dynamics equation is used. Convergence with good resolution to weak solutions for various Riemann problems are observed.
Shu, Chi-Wang↗
Applicability of the Effective-Medium Approximation to Heterogeneous Aerosol Particles.
The effective-medium approximation (EMA) is based on the assumption that a heterogeneous particle can have a homogeneous counterpart possessing similar scattering and absorption properties. We analyze the numerical accuracy of the EMA by comparing superposition T-matrix computations for spherical aerosol particles filled with numerous randomly distributed small inclusions and Lorenz-Mie computations based on the Maxwell-Garnett mixing rule. We verify numerically that the EMA can indeed be realized for inclusion size parameters smaller than a threshold value. The threshold size parameter depends on the refractive-index contrast between the host and inclusion materials and quite often does not exceed several tenths, especially in calculations of the scattering matrix and the absorption cross section. As the inclusion size parameter approaches the threshold value, the scattering-matrix errors of the EMA start to grow with increasing the host size parameter and or the number of inclusions. We confirm, in particular, the existence of the effective-medium regime in the important case of dust aerosols with hematite or air-bubble inclusions, but then the large refractive-index contrast necessitates inclusion size parameters of the order of a few tenths. Irrespective of the highly restricted conditions of applicability of the EMA, our results provide further evidence that the effective-medium regime must be a direct corollary of the macroscopic Maxwell equations under specific assumptions.
aerosols↗
Disintegration of fluids under supercritical conditions from mixing layer studies
Databases of transitional states obtained from Direct Numerical simulations (DNS) of temporal, supercritical mixing layers for two species systems, O2/H2 and C7H16/N2, are analyzed to elucidate species-specific turbulence aspects and features of fluid disintegration.
mixing layers supercritical combustion↗
Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 2: A semi-discrete error analysis framework for assessing coupling schemes
Abstract. Part 1 (Wan et al., 2024) of this study discusses the motivation and empirical evaluation of a revision to the aerosol-related numerical process coupling in the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) to address the previously reported issue of strong sensitivity of the simulated dust aerosol lifetime and dry removal rate to the model's vertical resolution. This paper complements that empirical justification of the revised scheme with a mathematical justification leveraging a semi-discrete analysis framework for assessing the splitting error of process coupling methods. The framework distinguishes the error due to numerical splitting from the error due to the time integration method(s) used for each individual process. Such a distinction results in a framework that provides an intuitive understanding of the causes of the splitting error. The application of this framework to the dust life cycle in EAMv1 confirms (i) that the original EAMv1 scheme artificially strengthens the effect of dry removal processes and (ii) that the revised splitting reduces that artificial strengthening. While the error analysis framework is presented in the context of the dust life cycle in EAMv1, the framework can be broadly leveraged to evaluate process coupling schemes, both in other physical problems and for any number of processes. This framework will be particularly powerful when the various process implementations support a variety of time integration approaches. Whereas traditional local truncation error approaches require separate consideration of each combination of time integration methods, this framework enables evaluation of coupling schemes independent of particular time integration approaches for each process while still allowing for the incorporation of these specific time integration errors if so desired. The framework also explains how the splitting error terms result from (i) the integration of individual processes in isolation from other processes and (ii) the choices of input state and time step size for the isolated integration of processes. Such a perspective has the potential for the rapid development of alternative coupling approaches that utilize knowledge both about the desired accuracy and about the computational costs of individual processes.
58 GEOSCIENCES↗
Numerical prediction of wakes in cascades and compressor rotors including the effects of mixing. I - Cascade wakes including the effects of incidence and freestream turbulence
The results of a numerical investigation to predict the flow in the wake regions of compressor cascades, and wakes and mixing in rotors are presented. Attention is given to the flow in compressor cascades including the effects of change in loading (incidence) and the inlet freestream turbulence intensity. The numerical analysis shows a slight increase in the total pressure loss coefficient through the cascade with increasing turbulence levels.
Lakshminarayana, B.↗
Mixing-length, shears, and differential rotation
Three aspects of convection in the solar envelope, the mixing length, the shears, and the differential rotation, are considered. Numerical experiments show that the mixing length is scaled by the pressure scale height rather than by the density scale height. Results indicate that the effective viscosity of convective turbulence on large scale shears may be approximated as 1/3 the product of the rms vertical velocity and the pressure scale height. It is suggested that the differential rotation of the solar convection zone may be viewed as an unstable axisymmetric mode which produces slower rotation at the surface.
Chan, K. L.↗
Phase Shift Interferometer and Growth Set Up to Step Pattern Formation During Growth From Solutions. Influence of the Oscillatory solution Flow on Stability
We have assembled an experimental setup based on Michelson interferometry with the growing crystal surface as one of the reflective surfaces. The crystallization part of the device allows optical monitoring of a face of a crystal growing at temperature stable within 0.05 C in a flow of solution of controlled direction and speed. The reference arm of the interferometer contains a liquid crystal element that allows controlled shifts of the phase of the interferograms. We employ an image-processing algorithm, which combines five images with a pi/2 phase difference between each pair of images. The images are transferred to a computer by a camera capable of capturing 60 frames per second. The device allows data collection on surface morphology and kinetics during the face layers growth over a relatively large area (approximately 4 sq. mm) in situ and in real time during growth. The estimated depth resolution of the phase shifting interferometry is approximately 50 Angstroms. The data will be analyzed in order to reveal and monitor step bunching during the growth process. The crystal chosen as a model for study in this work is KH2PO4 (KDP). This optically non-linear material is widely used in frequency doubling applications. There have been a number of studies of the kinetics of KDP crystallization that can serve as a benchmark for our investigations. However, so far, systematic quantitative characteristics of step interaction and bunching are missing. We intend to present our first quantitative results on the onset, initial stages and development of instabilities in moving step trains on vicinal crystal surfaces at varying supersaturation, flow rate, and flow direction. Behavior of a vicinal face growing from solution flowing normal to the steps and periodically changing its direction in time was considered theoretically. It was found that this oscillating flow reduces both stabilization and destabilization effects resulted from the unidirectional solution flow directed up the step stream and down the step stream. This reduction of stabilization and destabilization comes from effective mixing which entangles the phase shifts between the spatially periodic interface perturbation and the concentration wave induced by this perturbation. Numerical results and simplified mixing criterion will be discussed.
Chernov, Alex A.↗
Comparison of Mixing Characteristics for Several Fuel Injectors on an Open Plate and in a Ducted Flowpath Configuration at Hypervelocity Flow Conditions
In order to reduce the cost and complexity associated with fuel injection and mixing experiments for high-speed flows, and to further enable optical access to the test section for nonintrusive diagnostics, the Enhanced Injection and Mixing Project (EIMP) utilizes an open flat plate configuration to characterize inert mixing properties of various fuel injectors for hypervelocity applications. The experiments also utilize reduced total temperature conditions to alleviate the need for hardware cooling. The use of "cold" flows and non-reacting mixtures for mixing experiments is not new, and has been extensively utilized as a screening technique for scramjet fuel injectors. The impact of reduced facility-air total temperature, and the use of inert fuel simulants, such as helium, on the mixing character of the flow has been assessed in previous numerical studies by the authors. Mixing performance was characterized for three different injectors: a strut, a ramp, and a flushwall. The present study focuses on the impact of using an open plate to approximate mixing in the duct. Toward this end, Reynolds-averaged simulations (RAS) were performed for the three fuel injectors in an open plate configuration and in a duct. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared for the two configurations. In addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also computed for the reacting simulations.
Drozda, Tomasz G.↗
Observational and Model Studies of Large-Scale Mixing Processes in the Stratosphere
The following is the final technical report for grant NAGW-3442, 'Observational and Model Studies of Large-Scale Mixing Processes in the Stratosphere'. Research efforts in the first year concentrated on transport and mixing processes in the polar vortices. Three papers on mixing in the Antarctic were published. The first was a numerical modeling study of wavebreaking and mixing and their relationship to the period of observed stratospheric waves (Bowman). The second paper presented evidence from TOMS for wavebreaking in the Antarctic (Bowman and Mangus 1993). The third paper used Lagrangian trajectory calculations from analyzed winds to show that there is very little transport into the Antarctic polar vortex prior to the vortex breakdown (Bowman). Mixing is significantly greater at lower levels. This research helped to confirm theoretical arguments for vortex isolation and data from the Antarctic field experiments that were interpreted as indicating isolation. A Ph.D. student, Steve Dahlberg, used the trajectory approach to investigate mixing and transport in the Arctic. While the Arctic vortex is much more disturbed than the Antarctic, there still appears to be relatively little transport across the vortex boundary at 450 K prior to the vortex breakdown. The primary reason for the absence of an ozone hole in the Arctic is the earlier warming and breakdown of the vortex compared to the Antarctic, not replenishment of ozone by greater transport. Two papers describing these results have appeared (Dahlberg and Bowman; Dahlberg and Bowman). Steve Dahlberg completed his Ph.D. thesis (Dahlberg and Bowman) and is now teaching in the Physics Department at Concordia College. We also prepared an analysis of the QBO in SBUV ozone data (Hollandsworth et al.). A numerical study in collaboration with Dr. Ping Chen investigated mixing by barotropic instability, which is the probable origin of the 4-day wave in the upper stratosphere (Bowman and Chen). The important result from this paper is that even in the presence of growing, unstable waves, the mixing barriers around
Bowman, Kenneth P.↗
Numerical simulation of non-Newtonian free shear flows
Free shear flows, like those of mixing layers, are encountered in aerodynamics, in the atmosphere, and in the ocean as well as in many industrial applications such as flow reactors or combustion chambers. It is, therefore, crucial to understand the mechanisms governing the process of transition to turbulence in order to predict and control the evolution of the flow. Delaying transition to turbulence as far downstream as possible allows a gain in energy expenditure while accelerating the transition can be of interest in processes where high mixing is desired. Various methods, including the use of polymer additives, can be effective in controlling fluid flows. The drag reduction obtained by the addition of small amounts of high polymers has been an active area of research for the last three decades. It is now widely believed that polymer additives can affect the stability of a large variety of flows and that dilute solutions of these polymers have been shown to produce drag reductions of over 80 percent in internal flows and over 60 percent in external flows under a wide range of conditions. The major thrust of this work is to study the effects of polymer additives on the stability of the incompressible mixing layer through large scale numerical simulations. In particular, we focus on the two dimensional flow and examine how the presence of viscoelasticity may affect the typical structures of the flow, namely roll-up and pairing of vortices.
Homsy, G. M.↗
Numerical simulation of the flow and fuel-air mixing in an axisymmetric piston-cylinder arrangement
The implicit factored method of Beam and Warming was employed to describe the flow and the fuel-air mixing in an axisymmetric piston-cylinder configuration during the intake and compression strokes. The governing equations were established on the basis of laminar flow. The increased mixing due to turbulence was simulated by appropriately chosen effective transport properties. Calculations were performed for single-component gases and for two-component gases and for two-component gas mixtures. The flow field was calculated as functions of time and position for different geometries, piston speeds, intake-charge-to-residual-gas-pressure ratios, and species mass fractions of the intake charge. Results are presented in graphical form which show the formation, growth, and break-up of those vortices which form during the intake stroke and the mixing of fuel and air throughout the intake and compression strokes. It is shown that at bore-to-stroke ratio of less than unity, the vortices may break-up during the intake stroke. It is also shown that vortices which do not break-up during the intake stroke coalesce during the compression stroke. The results generated were compared to existing numerical solutions and to available experimental data.
Shih, T. I. P.↗
Scalar entrainment in the mixing layer
New definitions of entrainment and mixing based on the passive scalar field in the plane mixing layer are proposed. The definitions distinguish clearly between three fluid states: (1) unmixed fluid, (2) fluid engulfed in the mixing layer, trapped between two scalar contours, and (3) mixed fluid. The difference betwen (2) and (3) is the amount of fluid which has been engulfed during the pairing process, but has not yet mixed. Trends are identified from direct numerical simulations and extensions to high Reynolds number mixing layers are made in terms of the Broadwell-Breidenthal mixing model. In the limit of high Peclet number (Pe = ReSc) it is speculated that engulfed fluid rises in steps associated with pairings, introducing unmixed fluid into the large scale structures, where it is eventually mixed at the Kolmogorov scale. From this viewpoint, pairing is a prerequisite for mixing in the turbulent plane mixing layer.
Sandham, N. D.↗
Numerical study of terrain-induced mesoscale motions and hydrostatic form drag in a heated, growing mixed layer
Han et al. (1982) have found in a previous numerical study of terrain-induced mesoscale motions that the orography caused a steady-state flow pattern to occur. The study was concerned with a simplified case in which no surface heating occurred. The present investigation considers an extension of this study to the more realistic case of a heated, growing daytime mixed layer containing horizontal variations of potential temperature as well as velocity. The model is also extended to include three layers above the mixed layer. It is found for a heated, growing mixed layer, that the mesoscale form drag is a thermal-anomaly or buoyancy effect associated with horizontal variations of potential temperature within the layer.
Deardorff, J. W.↗
Results from computational analysis of a mixed compression supersonic inlet
A numerical study was performed to simulate the critical flow through a supersonic inlet. This flowfield has many phenomena such as shock waves, strong viscous effects, turbulent boundary layer development, boundary layer separations and mass flow suction through the walls (bleed). The computational tools used in this study were two full Navier-Stokes (FNS) codes. The supersonic inlet that was analyzed in this study is the Variable Diameter Centerbody (VDC), inlet. This inlet is a candidate concept for the next generation supersonic transport. Application of the code to the inlet geometry involved effort in generating an efficient grid geometry and specifying boundary conditions, particularly in the bleed region and at the outflow boundary. Results for a critical inlet operation compare favorably to Method of Characteristics predictions and experimental data.
Saunders, J. D.↗
Results from computational analysis of a mixed compression supersonic inlet
A numerical study was performed to simulate the critical flow through a supersonic inlet. This flow field has many phenomena such as shock waves, strong viscous effects, turbulent boundary layer development, boundary layer separations, and mass flow suction through the walls, (bleed). The computational tools used were two full Navier-Stokes (FNS) codes. The supersonic inlet that was analyzed is the Variable Diameter Centerbody, (VDC), inlet. This inlet is a candidate concept for the next generation supersonic involved effort in generating an efficient grid geometry and specifying boundary conditions, particularly in the bleed region and at the outflow boundary. Results for a critical inlet operation compare favorably to Method of Characteristics predictions and experimental data.
Saunders, J. D.↗