Search NASA⌕ Search

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

An explicit mixed numerical method for mesoscale model

A mixed numerical method has been developed for mesoscale models. The technique consists of a forward difference scheme for time tendency terms, an upstream scheme for advective terms, and a central scheme for the other terms in a physical system. It is shown that the mixed method is conditionally stable and highly accurate for approximating the system of either shallow-water equations in one dimension or primitive equations in three dimensions. Since the technique is explicit and two time level, it conserves computer and programming resources.

Hsu, H.-M.↗
NASA NTRS1995

Numerical mixing calculations of confined reacting jet flows in a cylindrical duct

The results reported in this paper describe some of the main flow characteristics and NOx production results which develop in the mixing process in a constant cross-sectional cylindrical duct. A 3-dimensional numerical model has been used to predict the mixing flow field and NOx characteristics in a mixing section of an RQL combustor. Eighteen configurations have been analyzed in a circular geometry in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying three parameter: (1) jet-to-mainstream momentum-flux ration (J), (2) orifice shape or orifice aspect ratio, and (3) slot slant angle. The results indicate that the mixing flow field and NOx production significantly vary with the value of the jet penetration and subsequently, slanting elongated slots generally improve the NOx production at high J conditions. Round orifices produce low NOx at low J due to the strong jet penetration. The NOx production trends do not correlate with the mixing non-uniformity parameters described herein.

Oechsle, Victor L.↗
NASA NTRS1990

Numerical simulation of low Prandtl number turbulent mixing

Numerical simulations of turbulent mixing of strongly diffusive scalar fields were carried out with and without subgrid-scale modeling of the small-scale strain field. For low Reynolds number flows, when the rate of strain field (determined primarily by the small scales) is fully resolved, the scalar microstructure was found to collapse under Batchelor rate-of-strain scaling even for small Prandtl numbers, in agreement with Kerr. For high Reynolds number flows, when small-scale straining is modeled with a subgrid-scale model, the scalar microstructure follows the Batchelor, Howells, and Townsend predictions that the small-scale rate-of-strain is irrelevant.

Gibson, C.↗
NASA NTRS1985

Resistive wave dissipation on magnetic inhomogeneities Normal modes and phase mixing

Numerical solutions of the linearized, resistive MHD equations indicate the presence of two characteristic forms of wave decay on magnetic inhomogeneities. Solutions for relatively long-wavelength disturbances and large values of the resistivity have the behavior of decaying normal modes. In the opposite limits, phase mixing becomes important, with the eventual build-up of large spatial gradients in which the resistive terms dominate. This composite behavior is shown to be conceptually consistent with analytic results, which predict that the complete solution consists of the sum of collective and noncollective contributions. The numerical simulations go beyond analytic theory by defining where each contribution prevails and computing wave decay when phase mixing is effective. The dispersion relation in the normal-mode regime is also determined analytically, with some approximations, and is in good agreement with the numerical solution. The decay time of the normal modes varies with Lundquist number at S1/6, while the phase-mixed decay time scales as S1/3.

Steinolfson, R. S.↗
NASA NTRS1995

Analysis and modeling of subgrid scalar mixing using numerical data

Direct numerical simulations (DNS) of passive scalar mixing in isotropic turbulence is used to study, analyze and, subsequently, model the role of small (subgrid) scales in the mixing process. In particular, we attempt to model the dissipation of the large scale (supergrid) scalar fluctuations caused by the subgrid scales by decomposing it into two parts: (1) the effect due to the interaction among the subgrid scales; and (2) the effect due to interaction between the supergrid and the subgrid scales. Model comparisons with DNS data show good agreement. This model is expected to be useful in the large eddy simulations of scalar mixing and reaction.

Girimaji, Sharath S.↗
NASA NTRS1990

Numerical simulations of compressible mixing layers

Direct numerical simulations of two-dimensional temporally growing compressible mixing layers are presented. The Kelvin-Helmholtz instability is initially excited by a white-noise perturbation superimposed onto a hyperbolic tangent meanflow profile. The linear regime is studied at low resolution in the case of two flows of equal temperatures, for convective Mach numbers from 0.1 to 1 and for different values of the Reynolds number. At higher resolution, the complete evolution of a two-eddy mixing layer between two flows of different temperatures is simulated at moderate Reynolds number. Similarities and differences between flows of equal convective Mach numbers are discussed.

Normand, Xavier↗
NASA NTRS1992

Numerical study on mixing of sprayed liquid in an LNG storage tank

This paper presents a numerical method to simulate the mixing of heavier LNG sprayed on lighter layer. Numerical results for evolutions of flow field and density field are obtained in a rectangular computational domain which includes the vicinity of the liquid surface. At the surface boundary, uniform distributions of the fluid velocity and the density are assumed. Detail structure of flow caused by impingements of liquid drops are neglected. But, to trigger a realistic motion, a series of random numbers is employed. It is used as an initial distribution of the density near the surface. This method successfully gives a realistic simulation of the mixing process. Numerical results for mixing velocity shows good agreement with experimental data.

Uchida, Hiroyuki↗
NASA NTRS1992

Bleed-boundary conditions for numerically simulated mixed-compression supersonic inlet flow

A numerical study was performed to investigate boundary conditions for bleed openings that are needed in the numerical simulation of critical flow through turbojet inlets with bleed. Nine different boundary conditions (BCs) were devised and examined. Of the BCs examined, three were able to stabilize the terminal shock near the inlet's throat. This numerical study was based on the ensemble-averaged, unsteady, three-dimensional, compressible Navier-Stokes equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the Navier-Stokes equations were obtained by using a hybrid implicit-explicit method (Beam and Warming/MacCormack). The multiblock grid system used in the solution procedure was generated by solving a set of elliptic partial differential equations. Solutions obtained were compared with available experimental data.

Chyu, W. J.↗
NASA NTRS1992

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Averaged Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗
NASA NTRS1993

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Average Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗
NASA NTRS1987

Numerical simulation of excited jet mixing layers

A numerical simulation of unsteady flow in jet mixing layers, both with and without external excitation, has been performed by solving the time-dependent compressible Navier-Stokes equations. Computations were performed on a CRAY X-MP computer using MacCormick's explicit finite difference algorithm. Different excitation methods were investigated and were shown to be very effective in controlling the well organized periodic production, shedding and pairing of large scale vortex structures. It is found that pressure excitation was generally more effective than temperature excitation, and that grid refinement results in substantial improvement in the resolution of unsteady features. The location and orientation, in addition to the frequency, of the excitation source are shown to have a significant influence on the production and interaction of large scale vortex structures in the jet mixing layer.

Scott, J. N.↗
NASA NTRS1993

Nonreflecting boundary conditions for linearized unsteady aerodynamic calculations

The present method for the implementation of nonreflecting boundary conditions in 2D and 3D linearized unsteady flow computations is applied to cases of unsteady flows in turbomachine blade rows. The eigenmodes of a discrete representation of the governing equations are computed and used to construct nonreflecting boundary conditions. In 3D, a mixed numerical method is used; in 2D, the discrete representation of the governing equations is obtained from the discretized equations used by the flow solver itself. Wavenumbers and radial mode shapes are computed.

Hall, Kenneth C.↗
NASA NTRS2010

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗
NASA NTRS1982

Numerical study of terrain-induced mesoscale motions in a mixed layer

Numerical integrations using a potential enstrophy-conserving scheme are presented for the flow within a mixed layer over hilly terrain using the hydrostatic shallow-water equations with a quadratic drag law. The mesoscale area treated is 150 km on a side; cyclic lateral boundary conditions are used. It is found that for the idealized conditions treated (no surface heating, no entrainment and no pressure adjustments aloft), the topography quickly induces a steady state flow pattern by means of surface friction. Unsteadiness does not occur unless a surface-friction Reynolds number is greater than approximately 100. Effects of varying the Rossby number, Froude number and terrain-height parameter are examined.

Han, Y.-J.↗
NASA NTRS1990

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗
NASA NTRS1992

Numerical Study of Mixing of Two Fluids Under Low Gravity

The mixing characteristics of two fluids inside a cavity due to buoyancy driven flow fields for low gravity conditions is investigated via numerical experiments. The buoyancy driven flow, depending on the parametric region, stretches and deforms the material interface into a wave morphological pattern. The morphological pattern affects the resulting stratification thickness of the mixed region. Three basic mixing regimes occur: convective, diffusive, and chaotic. In the convective regime, an overturning motion occurs which gives rise to a stable wave formation. This wave oscillates and its decay leads to a stable stratification. Whereas, in the diffusive regime, the length of the interface remains constant while mixing occurs. This limiting behavior is very important to materials processing in space, and it admits a closed form solution corresponding to vanishing convective terms which agrees with computational results. Finally, in the chaotic regime, the material interface continuously stretches and folds on itself similar to a horseshoe map. The length of stretch of the interface increases exponentially. Internal wavebreaking occurs for this case. This wavebreaking generates local turbulence, and provides an effective mechanism for mixing.

Duval, Walter M. B.↗