Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Eulerian model”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Transport of Passive Tracers in Baroclinic Wave Life Cycles

The transport of passive tracers in idealized baroclinic wave life cycles is studied using output from the National Center for Atmospheric Research Community Climate Model (CCM2). Two life cycles, LCn and LCs, are simulated, starting with baroclinically unstable initial conditions similar to those used by Thorncroft et al. in their study of two life cycle paradigms. The two life cycles LCn and LCs have different initial horizontal wind shear structures that result in distinctive nonlinear development. In terms of potential vorticity-potential temperature (PV-theta) diagnostics, the LCn case is characterized by thinning troughs that are advected anti-cyclonically and equatorward, while the LCs case has broadening troughs that wrap up cyclonically and poleward. Four idealized passive tracers are included in the model to be advected by the semi-Lagrangian transport scheme of the CCM2, and their evolutions are investigated throughout the life cycles. Tracer budgets are analyzed in terms of the transformed Eulerian mean constituent transport formalism in pressure coordinates and also in isentropic coordinates. Results for both LCn and LCs show transport that is downgradient with respect to the background structure of the tracer field, but with a characteristic spatial structure that maximizes in the middle to high latitudes. For the idealized tropospheric tracers in this study, this represents a net upward and poleward transport that enhances concentrations at high latitudes. These results vary little with the initial distribution of the constituent field. The time tendency of the tracer is influenced most strongly by the eddy flux term. with the largest transport occurring during the nonlinear growth stage of the life cycle. The authors also study the transport of a lower-stratospheric tracer, to examine stratosphere-troposphere exchange for baroclinic waves.

Stone, Elizabeth M.↗
NASA NTRS2019

Antarctic Ice Shelf Thickness Change from Multimission Lidar Mapping

We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.

Sutterley, Tyler C.↗
NASA NTRS1980

Numerical calculation of shock-induced initiation of detonations

Results of numerical calculations of the impact of steel cylinders and spheres on the plastic bonded high explosive PBX 9501 are described. The calculations were carried out by a reactive, multicomponent, two dimensional, Eulerian hydrodynamic computer code, 2DE. The 2DE computer code is a finite difference code that uses the donor acceptor cell method to compute mixed cell fluxes. The parameters in the Forest Fire burn model are developed from experiments where the induced shock approximates a plane wave and are applied, in this case, to a situation where the induced shock is a divergent wave with curvature that depends on the size and shape of the projectile. The calculated results are compared with results from experiments involving instrumented mock and live high explosives, with projectiles of varying size, shapes, and velocities.

Cort, G. E.↗
NASA NTRS2000

Unsteady Cascade Aerodynamic Response Using a Multiphysics Simulation Code

The multiphysics code Spectrum(TM) is applied to calculate the unsteady aerodynamic pressures of oscillating cascade of airfoils representing a blade row of a turbomachinery component. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena, in the present case being between fluids and structures. Interaction constraints are enforced in a fully coupled manner using the augmented-Lagrangian method. The arbitrary Lagrangian-Eulerian method is utilized to account for deformable fluid domains resulting from blade motions. Unsteady pressures are calculated for a cascade designated as the tenth standard, and undergoing plunging and pitching oscillations. The predicted unsteady pressures are compared with those obtained from an unsteady Euler co-de refer-red in the literature. The Spectrum(TM) code predictions showed good correlation for the cases considered.

Lawrence, C.↗
NASA NTRS1976

Impact ejecta on the moon

The response of a lunar-sized object to the impact of meteoroids no more than about 100 km in radius is studied by means of a numerical model. The partitioning of impact energy into the kinetic and internal energy of the ejecta is obtained by using the conservation of mass, momentum, and energy conservation equations in finite-difference form within an Eulerian framework with approximate equations of state. The calculations are performed for a 15 km/sec impact of an iron object 5 cm in radius on a gabbroic anorthosite surface. Ejecta ballistic analysis is then performed. Most of the material lost escaping the moon is lunar crust material. Only 0.2% of the meteoroid escapes, all in the vapor phase.

Okeefe, J. D.↗
NASA NTRS1984

Error analysis of analytic solutions for self-excited near-symmetric rigid bodies - A numerical study

Analytic error bounds are presented for the solutions of approximate models for self-excited near-symmetric rigid bodies. The error bounds are developed for analytic solutions to Euler's equations of motion. The results are applied to obtain a simplified analytic solution for Eulerian rates and angles. The results of a sample application of the range and error bound expressions for the case of the Galileo spacecraft experiencing transverse torques demonstrate the use of the bounds in analyses of rigid body spin change maneuvers.

Kia, T.↗
NASA NTRS1983

Effect of particle presence on the incompressible inviscid flow through a two dimensional compressor cascade

A computational model for calculating the flow field, in the presence of solid particles, through a two dimensional compressor cascade is presented. The method treats the particle phase in the Lagrangian system and the fluid phase in the Eulerian system. The equations of momentum of the fluid phase are modified to account for the momentum exchange between the two phases. The resulting modified momentum and continuity equations are reduced to the conventional stream function-vorticity formulation. The analysis yields the change in the blade surface pressure distribution, the total pressure and velocity in the flow field. Based on this, it is possible to predict the decrease in compressor performance under the presence of particles.

Balan, C.↗
NASA NTRS2020

Description and Evaluation of the Specified-Dynamics Experiment in the Chemistry-Climate Model Initiative

We provide an overview of the REF-C1SD specified-dynamics experiment that was conducted as part of phase 1 of the Chemistry-Climate Model Initiative (CCMI). The REF-C1SD experiment, which consisted of mainly nudged general circulation models (GCMs) constrained with (re)analysis fields, was designed to examine the influence of the large-scale circulation on past trends in atmospheric composition. The REF-C1SD simulations were produced across various model frameworks and are evaluated in terms of how well they represent different measures of the dynamical and transport circulations. In the troposphere there are large (~40 %) differences in the climatological mean distributions, seasonal cycle amplitude, and trends of the meridional and vertical winds. In the stratosphere there are similarly large (~50 %) differences in the magnitude, trends and seasonal cycle amplitude of the transformed Eulerian mean circulation and among various chemical and idealized tracers. At the same time, interannual variations in nearly all quantities are very well represented, compared to the underlying reanalyses. We show that the differences in magnitude, trends and seasonal cycle are not related to the use of different reanalysis products; rather, we show they are associated with how the simulations were implemented, by which we refer both to how the large-scale flow was prescribed and to biases in the underlying free-running models. In most cases these differences are shown to be as large or even larger than the differences exhibited by free-running simulations produced using the exact same models, which are also shown to be more dynamically consistent. Overall, our results suggest that care must be taken when using specified-dynamics simulations to examine the influence of large-scale dynamics on composition.

Clara Orbe↗
NASA NTRS1992

A hydrodynamic treatment of the cold dark matter cosmological scenario

The evolution of structure in a postrecombination Friedmann-Robertson-Walker universe containing both gaseous baryons and cold dark matter (CDM) is studied by means of an Eulerian code coupled with a standard particle-mesh code. Ionization state and radiative opacity are calculated in detail, and the hydrodynamic simulations make it possible to compute properties of gas distribution on scales larger than three cell sizes. The model yields a soft X-ray background consistent with the latest cosmic nucleosynthesis values, and can accurately reproduce the galaxy-galaxy two-point correlation. The rate of galaxy formation peaks at a relatively late epoch. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum peaks in the vicinity of m(b) = 10 exp 9.2 solar masses with a reasonable fit to the Schecter luminosity function if the baryon mass to blue light ratio is approximately 4. Overall, the simulations provide strong support for the CMD scenario. Of particular interest is that, while the baryons are not biased on scales greater than 1/h Mpc, the galaxies are, and that the 'galaxies' have a correlation function of the required slope and the correct amplitude.

Cen, Renyue↗
NASA NTRS2001

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗
NASA NTRS2014

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Yang, H. Q.↗
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗
NASA NTRS1985

Tracer transport by planetary waves - A comparison of explicit and parameterized models

A two-dimensional stratospheric tracer transport parameterization is developed using the theory of small-amplitude generalized Lagrangian means and tested by two experiments. The first tests the parameterization for mean climatological winter conditions without wave transience; the second tests it in the presence of transience and critical layers, when nonacceleration conditions are strongly violated. Tracer transport by Eulerian velocity, the residual velocity, and the nondivergent mean meridional Lagrangian velocity are compared.

Hess, P. G.↗
NASA NTRS1991

The meandering Gulf Stream as seen by the Geosat altimeter - Surface transport, position, and velocity variance from 73 deg to 46 deg W

Results are presented of an analysis of the surface geostrophic velocity field for the Gulf Stream region for the position, structure, and surface transport of the Gulf Stream for 2.5 yr of the Geosat altimeter Exact Repeat Mission. Synthetic data using a Gaussian velocity profile were generated and fit to the sea surface residual heights to create a synthetic mean sea surface height field and profiles of absolute geostrophic currents. An analysis of the model parameters and the actual geostrophic velocity profiles revealed two different flow regimes for the Gulf Stream connected by a narrow transition region coincident with the New England Seamount Chain. The upstream region was found to exhibit relatively straight Gulf Stream paths, long Eulerian time scales, and eastward propagating meanders. The downstream region had more large meanders, no consistent propagation direction, and shorter Eulerian time scales. A 25-percent reduction in surface transport occurred in the transition region, with a corresponding reduction in current speed and no change in Gulf Stream width.

Kelly, Kathryn A.↗
NASA NTRS1991

Three-Dimensional Analysis and Modeling of a Wankel Engine

A new computer code, AGNI-3D, has been developed for the modeling of combustion, spray, and flow properties in a stratified-charge rotary engine (SCRE). The mathematical and numerical details of the new code are described by the first author in a separate NASA publication. The solution procedure is based on an Eulerian-Lagrangian approach where the unsteady, three-dimensional Navier-Stokes equations for a perfect gas-mixture with variable properties are solved in generalized, Eulerian coordinates on a moving grid by making use of an implicit finite-volume, Steger-Warming flux vector splitting scheme. The liquid-phase equations are solved in Lagrangian coordinates. The engine configuration studied was similar to existing rotary engine flow-visualization and hot-firing test rigs. The results of limited test cases indicate a good degree of qualitative agreement between the predicted and measured pressures. It is conjectured that the impulsive nature of the torque generated by the observed pressure nonuniformity may be one of the mechanisms responsible for the excessive wear of the timing gears observed during the early stages of the rotary combustion engine (RCE) development. It was identified that the turbulence intensities near top-dead-center were dominated by the compression process and only slightly influenced by the intake and exhaust processes. Slow mixing resulting from small turbulence intensities within the rotor pocket and also from a lack of formation of any significant recirculation regions within the rotor pocket were identified as the major factors leading to incomplete combustion. Detailed flowfield results during exhaust and intake, fuel injection, fuel vaporization, combustion, mixing and expansion processes are also presented. The numerical procedure is very efficient as it takes 7 to 10 CPU hours on a CRAY Y-MP for one entire engine cycle when the computations are performed over a 31 x16 x 20 grid.

Raju, M. S.↗
NASA NTRS2000

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗
NASA NTRS1985

Vector diffraction analysis of reflector antennas with mesh surfaces

Reflector antennas with mesh surfaces are used extensively in many satellite and ground antenna systems. A strip-aperture modeling of commonly used mesh surfaces is presented which provides considerable versatility in characterizing the mesh cells. The mesh transmission coefficients are constructed using a Floquet-modal expansion in conjuction with two dominant aperture modes. To account for the mesh local coordinates, the Eulerian angle transformation is invoked to obtain the total induced current on the curved reflector surface. General formulas are presented to show how the solid surface induced current is modified due to the transmission through the mesh. The effects of a variety of mesh configurations on both the co-polar and cross-polar patterns of reflector antennas are studied by numerically evaluating the vector diffraction integral using the Jacobi-Bessel expansion. For some special cases, a comparison is made with the results of the commonly used wire-grid formulation. Many of the numerical data are tailored to the dimensions of a conceptually designed mesh deployable offset reflector of the land mobile satellite system (LMSS).

Rahmat-Samii, Y.↗
NASA NTRS1987

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗