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 163 records · Page 9

Electromagnetic wave extinction within a forested canopy

A forested canopy is modeled by a collection of randomly oriented finite-length cylinders shaded by randomly oriented and distributed disk- or needle-shaped leaves. For a plane wave exciting the forested canopy, the extinction coefficient is formulated in terms of the extinction cross sections (ECSs) in the local frame of each forest component and the Eulerian angles of orientation (used to describe the orientation of each component). The ECSs in the local frame for the finite-length cylinders used to model the branches are obtained by using the forward-scattering theorem. ECSs in the local frame for the disk- and needle-shaped leaves are obtained by the summation of the absorption and scattering cross-sections. The behavior of the extinction coefficients with the incidence angle is investigated numerically for both deciduous and coniferous forest. The dependencies of the extinction coefficients on the orientation of the leaves are illustrated numerically.

Karam, M. A.↗

Role of shock dynamics in transonic flutter

A computational study of the influence of shock motion on flutter and divergence in transonic flow is presented. The numerical scheme models the entire fluid-structure system as a single continuum dynamics problem, by using a mixed Eulerian-Lagrangian formulation. No assumptions of small displacements are made, but the effect of viscosity is neglected. The results from this study indicate that the shock dynamics gives rise to limit cycles and highly nonlinear aeroelastic phenomena, such as weak divergence and flutter-divergence interactions. Although the shocks typically are destabilizing at the linear flutter boundary, they often have a strongly stabilizing effect for moderate-amplitude motions. The shocks are thus capable of quenching an emerging bending-torsion flutter motion and turning it into limit cycle flutter. The usefulness of classical flutter and divergence boundary diagrams is severely limited in transonic flow, because much of the global dynamic stability information is lost in such a presentation.

Bendiksen, Oddvar O.↗

Consistent Large-Eddy Simulation of a Temporal Mixing Layer Laden with Evaporating Drops: A Posteriori Modelling - Part 2

Large-eddy simulation (LES) is conducted of a three-dimensional temporal mixing layer whose lower stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in an Eulerian frame for two perfect gas species (carrier gas and vapour emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame. The effect of drop evaporation on the gas phase is considered through mass, species, momentum and energy source terms. The drop evolution is modelled using physical drops, or using computational drops to represent the physical drops. Simulations are performed using various LES models previously assessed on a database obtained from direct numerical simulations (DNS). These LES models are for: (i) the subgrid-scale (SGS) fluxes and (ii) the filtered source terms (FSTs) based on computational drops. The LES, which are compared to filtered-and-coarsened (FC) DNS results at the coarser LES grid, are conducted with 64 times fewer grid points than the DNS, and up to 64 times fewer computational than physical drops. It is found that both constant-coefficient and dynamic Smagorinsky SGS-flux models, though numerically stable, are overly dissipative and damp generated small-resolved-scale (SRS) turbulent structures. Although the global growth and mixing predictions of LES using Smagorinsky models are in good agreement with the FC-DNS, the spatial distributions of the drops differ significantly. In contrast, the constant-coefficient scale-similarity model and the dynamic gradient model perform well in predicting most flow features, with the latter model having the advantage of not requiring a priori calibration of the model coefficient. The ability of the dynamic models to determine the model coefficient during LES is found to be essential since the constant-coefficient gradient model, although more accurate than the Smagorinsky model, is not consistently numerically stable despite using DNS-calibrated coefficients. With accurate SGS-flux models, namely scale-similarity and dynamic gradient, the FST model allows up to a 32-fold reduction in computational drops compared to the number of physical drops, without degradation of accuracy; a 64-fold reduction leads to a slight decrease in accuracy.

two-phase flow↗

(abstract) A Polarimetric Model for Effects of Brine Infiltrated Snow Cover and Frost Flowers on Sea Ice Backscatter

A polarimetric scattering model is developed to study effects of snow cover and frost flowers with brine infiltration on thin sea ice. Leads containing thin sea ice in the Artic icepack are important to heat exchange with the atmosphere and salt flux into the upper ocean. Surface characteristics of thin sea ice in leads are dominated by the formation of frost flowers with high salinity. In many cases, the thin sea ice layer is covered by snow, which wicks up brine from sea ice due to capillary force. Snow and frost flowers have a significant impact on polarimetric signatures of thin ice, which needs to be studied for accessing the retrieval of geophysical parameters such as ice thickness. Frost flowers or snow layer is modeled with a heterogeneous mixture consisting of randomly oriented ellipsoids and brine infiltration in an air background. Ice crystals are characterized with three different axial lengths to depict the nonspherical shape. Under the covering multispecies medium, the columinar sea-ice layer is an inhomogeneous anisotropic medium composed of ellipsoidal brine inclusions preferentially oriented in the vertical direction in an ice background. The underlying medium is homogeneous sea water. This configuration is described with layered inhomogeneous media containing multiple species of scatterers. The species are allowed to have different size, shape, and permittivity. The strong permittivity fluctuation theory is extended to account for the multispecies in the derivation of effective permittivities with distributions of scatterer orientations characterized by Eulerian rotation angles. Polarimetric backscattering coefficients are obtained consistently with the same physical description used in the effective permittivity calculation. The mulitspecies model allows the inclusion of high-permittivity species to study effects of brine infiltrated snow cover and frost flowers on thin ice. The results suggest that the frost cover with a rough interface significantly increases the backscatter from thin saline ice and the polarimetric signature becomes closer to the isotropic characteristics. The snow cover also modifies polarimetric signatures of thin sea ice depending on the snow mixture and the interface condition.

thin sea ice radar backscatter frost flowers snow ↗

Atomization simulations using an Eulerian-VOF-Lagrangian method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservations are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present innovative approach by simulating benchmark problems including the coaxial jet atomization.

Chen, Yen-Sen↗

Thermodynamic balance of three-dimensional stratospheric winds derived from a data assimilation procedure

The NASA/Goddard three-dimensional chemistry and transport model is driven by winds from a stratospheric data assimilation system. Synoptic- and planetary-scale patterns, apparent in satellite observations of trace constituents, are successfully reproduced for seasonal integrations. As model integrations proceed, however, the quality of simulations decreases, and systematic differences between calculation and measurement appear. The differences are explained by examining the zonal-mean residual circulation. The vertical velocity w-bar (sup star) is calculated two ways: (1) from the diabatic heating rates and temperature tendency and (2) from the Eulerian vertical velocity and the horizontal eddy heat flux convergence. The results from these calculations differ substantially. Periodic insertion of observational data during the assimilation process continually shocks the general circulation model and produces these differences, which leads to an overestimate of the mean vertical heat and cconstituent transport. Such differences are expected to be general to all data assimilation products. This interpretation is corroborated by two-dimensional (2D) model calculations. When w-bar(sup star) is calculated from (2), the 2D ozone evolution is unrealistic and qualitatively similar to the 3D model simulation. The 2D ozone evolution is reasonable when w-bar (sup star) from (1).

Weaver, Clark J.↗

A Priori Subgrid Scale Modeling for a Droplet Laden Temporal Mixing Layer

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using a direct numerical simulation (DNS) database. The DNS is for a Reynolds number (based on initial vorticity thickness) of 600, with droplet mass loading of 0.2. The gas phase is computed using a Eulerian formulation, with Lagrangian droplet tracking. Since Large Eddy Simulation (LES) of this flow requires the computation of unfiltered gas-phase variables at droplet locations from filtered gas-phase variables at the grid points, it is proposed to model these by assuming the gas-phase variables to be given by the filtered variables plus a correction based on the filtered standard deviation, which can be computed from the sub-grid scale (SGS) standard deviation. This model predicts unfiltered variables at droplet locations better than simply interpolating the filtered variables. Three methods are investigated for modeling the SGS standard deviation: Smagorinsky, gradient and scale-similarity. When properly calibrated, the gradient and scale-similarity methods give results in excellent agreement with the DNS.

Okongo, Nora↗

How Very Massive Metal-Free Stars Start Cosmological Reionization

The initial conditions and relevant physics for the formation of the earliest galaxies are well specified in the concordance cosmology. Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include nonequilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the case B approximation using adaptively ray-traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of approx.10(exp 6). These first sources of reionization are highly intermittent and anisotropic and first photoionize the small-scale voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf-sized galaxies, the escape fraction of UV radiation decreases and the H II regions only break out on some sides of the galaxies, making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately five ionizing photons are needed per sustained ionization when star formation in 10(exp 6) stellar Mass halos is dominant in the calculation. As the halos become larger than approx.10(exp 7) Stellar Mass, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Radiative feedback decreases clumping factors by 25% when compared to simulations without star formation and increases the average temperature of ionized gas to values between 3000 and 10,000 K.

Wise, John H.↗

Simulations of Turbulent Momentum and Scalar Transport in Confined Swirling Coaxial Jets

This paper presents the numerical simulations of confined three-dimensional coaxial water jets. The objectives are to validate the newly proposed nonlinear turbulence models of momentum and scalar transport, and to evaluate the newly introduced scalar APDF and DWFDF equation along with its Eulerian implementation in the National Combustion Code(NCC). Simulations conducted include the steady RANS, the unsteady RANS (URANS), and the time-filtered Navier-Stokes (TFNS); both without and with invoking the APDF or DWFDF equation.

combustors↗

A LaGrangian Perspective of the Surface Energy Budget Contributions to Sea Ice Evolution in the ‘New Arctic’

Recent Arctic warming outpaces global warming and is characterized by dramatic changes in sea ice (record sea ice loss, thinner and more seasonal ice cover), warming SSTs, and a lengthening melt season. These changes impact the surface energy budget and the sensitivity of sea ice to radiative forcing, namely an increased sensitivity of thinner ice to local radiative forcings and a weaker buffering of surface-atmosphere exchanges of heat and moisture. However, these processes are a large source of uncertainty in global climate models with prior studies finding that poor representations of downwelling longwave radiation and sea ice albedo lead to model spread in sea ice conditions. A limitation of prior work is examining sea ice-atmosphere interactions through an Eulerian framework. This study utilizes a sea ice parcel dataset in a LaGrangian framework, allowing for tracking of surface energy budget-sea ice interactions over the lifespan of individual sea ice parcels; this framework also allows for the quantification of the sea ice response to episodic weather events along with seasonal drivers, including the drivers of rapid growth and melt events. First, lag-lead relationships between relevant surface energy budget terms and sea ice thickness changes are explored to understand the atmospheric conditions preceding rapid growth and melt. The impact of episodic events- moisture intrusion, warm air outbreaks, cyclones- on sea ice growth are also investigated by stratifying first-year parcels by sea ice thickness at the end of the growth season and analyzing the effect of atmospheric variability (number and intensity of episodic events) on parcel thickness change. Lastly, parcel growth and melt and the associated radiative anomalies are compared to snowfall events/snow thickness to determine how snowfall modifies the interaction between the local surface energy budget and sea ice mass balance.

Patrick Taylor↗

A posterirori study of models for large eddy simulations of drop-laden flows

Large Eddy Simulation (LES) is conducted of a three-dimensional temporal mixing layer whose stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in Eulerian frame for two perfect gas species (carrier gas and vapor emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame.

spray modeling↗

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Our work explores the use of extended reality (XR) to im- prove scientific discovery with numerical weather/climate models that inform Earth science digital twins, specifically the NASA Goddard Earth Observing System (GEOS) global atmospheric model. The overall project is named the Vi- sualization And Lagrangian dynamics Immersive eXtended Reality Toolkit (VALIXR), which has two main areas of focus: (1) enhancing the understanding of and interaction with model output data through advanced visualizations in the XR environment, and (2) the integration of Lagrangian dynamics into the GEOS model, which allows a natural, feature-specific analysis of Earth science phenomena as op- posed to traditional, fixed-point Eulerian dynamics. Here, we report initial work on these focus areas.

Thomas Grubb↗

Orientation of the earth by numerical integration

A fundamental problem is the determination of the orientation of the earth in the celestial coordinate system. Classical reductions for precession and nutation can be expected to be consistent with the present-day observations, however, corrections to the classical theory are difficult to model because of the large number of coefficients involved. Consequently, a portion of the research has been devoted to numerically integrating the Eulerian equations of motion for a rigid earth and considering the six initial conditions of the integration as unknowns. Comparison of the three adjusted Eulerian angles from the numerical integration over 1000 days indicates agreement with classical theory to within 0.003 seconds of arc.

Fajemirokun, F. A.↗

Gas dynamics for accretion disk simulations

The behavior of accretion disks can largely be understood in terms of the basic physical processes of mass, energy, and momentum conservation. Despite this, detailed modeling of these systems using modern computational techniques is challenging and controversial. Disturbing differences exist between methods used widely in astrophysics, namely Eulerian finite-difference techniques and particle codes such as SPH. Therefore neither technique is fully satisfactory for accretion disk simulations. This paper describes a new fully Lagrangian method designed to resolve these difficulties.

Whitehurst, R.↗

Relativistic tidal interaction of a white dwarf with a massive black hole

We compute encounters of a realistic white dwarf model with a massive black hole in the regime where relativistic effects are important, using a three-dimensional, finite-difference, Eulerian, piecewise parabolic method (PPM) hydrodynamical code. Both disruptive and nondisruptive encounters are considered. We identify and discuss relativistic effects important for the problem: relativistic shift of the pericenter distance, time delay, relativistic precession, and the tensorial structure of the tidal forces. In the nondisruptive case, stripping of matter takes place. In the surface layers of the surviving core, complicated hydrodynamical phenomena are revealed. In both disruptive and nondispruptive encounters, material flows out in the form of two thin, S-shaped, supersonic jets. Our results provide realistic initial conditions for the subsequent investigation of the dynamics of the debris in the field of the black hole. We evaluate the critical conditions for complete disruption of the white dwarf, and compare our results with the corresponding results for nonrelativistic encounters.

Frolov, V. P.↗

Comprehensive modeling of a liquid rocket combustion chamber

An analytical model for the simulation of detailed three-phase combustion flows inside a liquid rocket combustion chamber is presented. The three phases involved are: a multispecies gaseous phase, an incompressible liquid phase, and a particulate droplet phase. The gas and liquid phases are continuum described in an Eulerian fashion. A two-phase solution capability for these continuum media is obtained through a marriage of the Implicit Continuous Eulerian (ICE) technique and the fractional Volume of Fluid (VOF) free surface description method. On the other hand, the particulate phase is given a discrete treatment and described in a Lagrangian fashion. All three phases are hence treated rigorously. Semi-empirical physical models are used to describe all interphase coupling terms as well as the chemistry among gaseous components. Sample calculations using the model are given. The results show promising application to truly comprehensive modeling of complex liquid-fueled engine systems.

Liang, P.-Y.↗

Eddy mixing coefficient values derived from simulations with the NASA Ames Mars GCM

Values of eddy mixing coefficients, especially that for vertical mixing, are of particular importance for various photochemical models of the Mars atmosphere, which are either globally or zonally averaged. These models represent atmospheric transport processes in terms of eddy diffusion. While this is not appropriate for the advective transport by the mean meridional circulation (in general, a part of the eddy diffusion tensor actually does represent advection--a 'correction' to the Eulerian-mean circulation), it can be applied to simulate both the small- and larger-scale eddy transports of an atmospheric constituent. An effort is underway to estimate values of the eddy mixing coefficients for the Mars atmosphere using circulation data generated with the NASA Ames Mars GCM. This model simulates the three-dimensional winds in the atmosphere, and these are then used as inputs for a tracer transport model. The latter model has previously been used in interactive dust transport experiments with the Ames GCM. Idealized tracer transport experiments, with a conservative or nearly conservative hypothetical tracer, are performed to generate data from which the eddy diffusion coefficients can be estimated. Carrying out matched pairs of tracer experiments, using two very different initial states, permits all four components of the diffusion tensor to be determined. A large number of GCM experiments have been conducted, spanning a range of seasons and atmospheric dust loading, allowing the eddy mixing coefficients to be estimated for a variety of atmospheric conditions.

Barnes, J. R.↗