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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A new model of resonance in the winter stratosphere

It is generally accepted that the planetary waves observed in the winter stratosphere are primarily a response to dynamical forcing from the troposphere. Nevertheless, the mechanism by which wave amplitudes sometimes become larger remains uncertain. It is possible that anomalously large waves in the stratosphere might simply be the result of anomalously large tropospheric forcing. However, it has also been suggested that they are a response to a stratospheric-tropospheric cavity being in a near-resonant configuration. It has been suggested that nonlinear self-tuning effects could play an important role in the behavior of such a cavity. Self-tuning may occur when a system starts to one side of resonance, such that the mean-state change induced by growing waves brings the system closer to resonance. A new model of the stratospheric cavity is introduced and is then used to re-examine the possibility of wave growth in the real atmosphere and in atmospheric models due to self-tuning effects. The new model is based on the picture of the winter-time stratosphere which has been revealed by the observations of Ertel's potential vorticity, Q. Isentropic maps of Q show two rather distinct regions, the first containing the circumpolar vortex, where gradients of Q are large and Rossby waves may propagate easily. Surrounding this is a second, low-latitude region where the gradients are generally weak and where, because the Eulerian-mean flow is comparable with their phase speed, Rossby waves must be continually breaking. As the waves are observed to grow the relative sizes of these two regions change in time. This leads to the interesting possibility that self-tuning, mainly due to irreversible changes in the size of the polar vortex, is taking place.

Haynes, P. H.↗
NASA NTRS1993

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗
NASA NTRS1983

Dynamics of the middle atmosphere

It is demonstrated by means of model calculations that while the general configuration of middle atmosphere dynamics (i.e., winter westerlies and summer easterlies) is determined by differential solar heating, the deviations of the zonally-averaged temperature field from radiative equilibrium and the closure of the jet structures with increasing altitude result from the action of zonal mean momentum dissipation process. The apparent heating and acceleration of the mean zonal state by planetary wave heat and momentum fluxes are examined using an Eulerian framework, and it is demonstrated that these are overestimates of their net effect. It is argued that since decelerations of the mean zonal flow are required in both winter and summer, and planetary waves are known to be very weak in the summer middle atmosphere, gravity waves are probably responsible for most of the middle atmosphere momentum dissipation as a result of their attenuation with height above their wave breaking altitude. Radar observations of middle atmosphere dynamics together with some theoretical work indicates that the effects of breaking gravity waves should be included in the thermodynamic equation as well as the momentum equation.

Geller, M. A.↗
NASA NTRS1993

Tracking of raindrops in flow over an airfoil

The splashback that occurs when raindrops impact an airfoil results in an 'ejecta fog' of small droplets around the leading edge. Acceleration of these droplets by the air flow field is a momentum sink for the air flow and has been hypothesized to contribute to the degradation of airfoil performance in heavy rain. Presented here is a one-way coupled Eulerian-Lagrangian particle tracking scheme to evaluate droplet number densities and momentum sink terms around a NACA 64-210 airfoil section for three rainfall rates. A laminar air flow field is determined with a standard CFD code and input to the particle tracking algorithm. Raindrops are assumed to be non-interacting, non-deforming, non-evaporating, and non-spinning spheres and are tracked through the curvilinear grid used by the air flow code. A simple model is used to simulate impacts and the resulting splashback on the airfoil surface.

Valentine, James R.↗
NASA NTRS1998

Three-Dimensional Hydrodynamic Simulations of the Postimpact Proto-Earth

Astrophysical fluid configurations are susceptible to a variety of nonaxisymmetric instabilities under the combined effects of rotation, self-gravity, and thermal pressure. When strong enough, they can induce rapid transport of mass and angular momentum. Our own previous studies of nonaxisymmetric instabilities in model protostars and protostellar disks show that significant transport can occur on orbital timescales and that material can be ejected to large distances. In this contribution, we present three-dimensional simulations of the circumterrestrial debris belt that may have resulted from a giant impact. Our three-dimensional hydrodynamics code with self-gravity and artificial viscosity is fully second-order in space and time; the equations of hydrodynamics and the Poisson equation are solved on an Eulerian cylindrical grid. In the preliminary calculations presented here, we use a simplified EOS where the central proto-Earth is treated as an n = 1/2 polytropic fluid, surrounded by a more compressible, rapidly rotating, fluid disk that represents silicate vapor. Our initial disk parameters are generated from the endstate data of recent smoothed particle hydrodynamics giant-impact calculations. Ultimately, we wish to detennine under what conditions nonaxisymmetric instabilities grow in the postimpact disk and whether they facilitate the transport of material outside the proto-Earth's Roche Limit, leading to the formation of the Moon. In future work, we hope to include a more realistic EOS and the consequences of heating, cooling, and phase transitions.

Pickett, B. K.↗
NASA NTRS1975

Three-dimensional finite element analysis for high velocity impact

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Chan, S. T. K.↗
NASA NTRS2003

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

Lin, Shian-Jiann↗
NASA NTRS1991

The formation of terrains antipodal to major impacts

The SALE, or Simplified Arbitrary Lagrangian Eulerian code is presently used to calculate the antipodal pressures generated by large impacts on icy satellites and planets of various compositions and core sizes. While in a high density core/lower density mantle planet there is a refraction of the incident pressure wave away from the antipode by the core, resulting in a low antipodal pressure, a planet with a low velocity core with a higher-wave speed mantle will focus the waveform to maximize antipodal pressure. Specific impact basins on the moon, Mercury, and several icy satellites are modeled to ascertain whether the impacts were capable of producing the observed antipodal features.

Watts, Alison W.↗
NASA NTRS1998

EUPDF: An Eulerian-Based Monte Carlo Probability Density Function (PDF) Solver. User's Manual

EUPDF is an Eulerian-based Monte Carlo PDF solver developed for application with sprays, combustion, parallel computing and unstructured grids. It is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and spray solvers. The solver accommodates the use of an unstructured mesh with mixed elements of either triangular, quadrilateral, and/or tetrahedral type. The manual provides the user with the coding required to couple the PDF code to any given flow code and a basic understanding of the EUPDF code structure as well as the models involved in the PDF formulation. The source code of EUPDF will be available with the release of the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗
NASA NTRS1981

The Lagrangian-mean motions forced by steady, dissipating equatorial waves. I

Waves are treated with a normal mode structure in order to determine the steady mean motion of the atmosphere that can be induced by dissipating equatorial waves. A model is developed which comprises a continuously stratified atmosphere at rest on the equatorial beta-plane. It is assumed that waves are excited by the corrugated bottom and are in a steady state, that dissipation is due to Newtonian cooling and Rayleigh friction, steadiness in wave magnitude is up to the second order, the waves have a long wave length, wave induced mean flows do not affect the waves, mean flows are steady, and dissipation mechanisms for the mean flows are the same as for the waves. Disturbance equations are formulated, along with Eulerian- and Lagrangian-mean flows, and the nonexistence of cross equatorial mean flows is demonstrated. Kelvin waves are shown to possess a Lagrangian-mean meridional circulation which is the same as the Eulerian-mean circulation. In the Boussinesq limit, however, neither the Eulerian- nor the Lagrangian-mean meridional circulations are caused by Kelvin waves. Further examination is made of Rossby-gravity waves and n = 1 westward propagating inertio-gravity waves.

Takahashi, M.↗
NASA NTRS1992

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗
NASA NTRS1988

Rapid slewing of the orbiting Spacecraft Control Laboratory Experiment (SCOLE) using LQR techniques

The rotational equations of motion, describing the dynamics of the (rigidized) proposed orbiting Spacecraft Control Laboratory Experiment (SCOLE) during the station keeping phase, are derived using the Eulerian formulation. When the attitude angles (roll, pitch, and yaw) are assumed small, a stability analysis is conducted for the system. It is seen that the pitch equation decouples from the roll and yaw equations when the interface between the mast on the reflector is not offset or the offset is only along the Shuttle roll axis. When a second offset is introduced along the pitch axis of the system and when the gravity-gradient torques are present in the dynamics, the system assumes a new equilibrium position. The linear regulator theory is used to derive a control law for the linear model of the rigidized SCOLE. This law is applied to the nonlinear model of the same configuration of the system and preliminary single axis slewing maneuvers (20 deg amplitude) are simulated.

Diarra, Cheick M.↗
NASA NTRS2007

Chemical Data Assimilation Estimates of Continental US Ozone and Nitrogen Budgets during INTEX-A

Global ozone analyses, based on assimilation of stratospheric profile and ozone column measurements, and NOy predictions from the Real-time Air Quality Modeling System (RAQMS) are used to estimate the ozone and NOy budget over the Continental US during the July-August 2004 Intercontinental Chemical Transport Experiment-North America (INTEX-A). Comparison with aircraft, satellite, surface, and ozonesonde measurements collected during the INTEX-A show that RAQMS captures the main features of the global and Continental US distribution of tropospheric ozone, carbon monoxide, and NOy with reasonable fidelity. Assimilation of stratospheric profile and column ozone measurements is shown to have a positive impact on the RAQMS upper tropospheric/lower stratosphere ozone analyses, particularly during the period when SAGE III limb scattering measurements were available. Eulerian ozone and NOy budgets during INTEX-A show that the majority of the Continental US export occurs in the upper troposphere/lower stratosphere poleward of the tropopause break, a consequence of convergence of tropospheric and stratospheric air in this region. Continental US photochemically produced ozone was found to be a minor component of the total ozone export, which was dominated by stratospheric ozone during INTEX-A. The unusually low photochemical ozone export is attributed to anomalously cold surface temperatures during the latter half of the INTEX-A mission, which resulted in net ozone loss during the first 2 weeks of August. Eulerian NOy budgets are shown to be very consistent with previously published estimates. The NOy export efficiency was estimated to be 24 percent, with NOx+PAN accounting for 54 percent of the total NOy export during INTEX-A.

Pierce, Robert B.↗
NASA NTRS1995

Interstellar clouds in high-speed, supersonic flows: Two-dimensional simulations

We present a series of gasdynamical simulations of the interaction of a dense, cool interstellar cloud with a high-speed, supersonic wind that confines and accelerates the embedded cloud. Our goal is to attempt to determine if such clouds can survive various potentially disruptive instabilities, that occur at their peripheries, long enough to be accelerated to speeds which are comparable to the wind velocity. These simulations are performed using two-dimensional, Eulerian gas dynamics on both an axisymmetric (about the cloud axis) and 'slab' geometric grid. The spatial and temporal resolutions of the simulations are varied over a wide range to investigate the effects of small-scale instabilities on the overall acceleration of clouds and the development of large-scale, disruptive instabilities. Also, we study the effects of wind/cloud Mach number variations by changing the wind speed constant at about 12 km/s (which corresponds to a cloud temperature of 10,000 K). The current simulations track the evolution of clouds as they are accelerated to speeds approximately 4-5 times greater than their internal sound speeds. Furthermore, the models with the highest resolution were extended far beyond quasi-linear Rayleigh-Taylor growth times reaching 6-7 Rayleigh-Taylor growth times for the largest scale instabilities before being terminated because of the accumulation of errors at the rear grid boundary.

Schiano, A. V. R.↗
NASA NTRS1994

Comparison of dynamical approximation schemes for non-linear gravitational clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.

Melott, Adrian L.↗
NASA NTRS1985

The dynamics and control of large flexible space structures, 8

A development of the in plane open loop rotational equations of motion for the proposed Spacecraft Control Laboratory Experiment (SCOLE) in orbit configuration is presented based on an Eulerian formulation. The mast is considered to be a flexible beam connected to the (rigid) shuttle and the reflector. Frequencies and mode shapes are obtained for the mast vibrational appendage modes (assumed to be decoupled) for different boundary conditions based on continuum approaches and also preliminary results are obtained using a finite element representation of the mast reflector system. The linearized rotational in plane equation is characterized by periodic coefficients and open loop system stability can be examined with an application of the Floquet theorem. Numerical results are presented to illustrate the potential instability associated with actuator time delays even for delays which represent only a small fraction of the natural period of oscillation of the modes contained in the open loop model of the system. When plant and measurement noise effects are added to the previously designed deterministic model of the hoop column system, it is seen that both the system transient and steady state performance are degraded. Mission requirements can be satisfied by appropriate assignment of cost function weighting elements and changes in the ratio of plant noise to measurement noise.

Bainum, P. M.↗