Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian modeling”

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 NTRS2005

3 Lectures: "Lagrangian Models", "Numerical Transport Schemes", and "Chemical and Transport Models"

The topics for the three lectures for the Canadian Summer School are Lagrangian Models, numerical transport schemes, and chemical and transport models. In the first lecture I will explain the basic components of the Lagrangian model (a trajectory code and a photochemical code), the difficulties in using such a model (initialization) and show some applications in interpretation of aircraft and satellite data. If time permits I will show some results concerning inverse modeling which is being used to evaluate sources of tropospheric pollutants. In the second lecture I will discuss one of the core components of any grid point model, the numerical transport scheme. I will explain the basics of shock capturing schemes, and performance criteria. I will include an example of the importance of horizontal resolution to polar processes. We have learned from NASA's global modeling initiative that horizontal resolution matters for predictions of the future evolution of the ozone hole. The numerical scheme will be evaluated using performance metrics based on satellite observations of long-lived tracers. The final lecture will discuss the evolution of chemical transport models over the last decade. Some of the problems with assimilated winds will be demonstrated, using satellite data to evaluate the simulations.

Anne Douglass↗
NASA NTRSDate not supplied

Modeling of Cavity Residence Time via Modeled Lagrangian Particle Tracking

One of the critical components in a scramjet is the cavity flame holder. In evaluating the sizing of this component, the residence time of the fuel in the cavity is a common parameter used. However, to obtain the residence time via CFD is an expensive undertaking. One method is to perform an LES simulation of the cavity and use Lagrangian particle tracking to obtain both the mean and probability density function of the time that the particles remain in the cavity. While this method affords a significant amount of information, it is very costly in CPU run time. A simpler method involves performing a time accurate RANS simulation of the cavity, and when converged, overlay a traceable scalar in the cavity region, and then record the decay of this scalar’s flux at a downstream location. The time scale of the decay rate can then be estimated and a residence time inferred. This method is less computationally expensive than the LES approach, but still requires a time-accurate solution while only providing a mean cavity residence time evaluation. In this paper a method of determining the mean cavity residence time as well as the Probability Density Function (PDF) of the residence time is proposed and demonstrated. This method is based on Lagrangian particle tracking modeled by the Generalized Langevin Equation and is implemented as a post-processing step for a steady RANS solution. Computational cost is on the order of minutes and the results are in agreement with values obtained by the LES and scalar decay methods.

Langevin↗
NASA NTRS2020

Lagrangian Modeling and Simulation of the Free Surface-Affected Dynamics of Underwater Vehicles

The paper deals with the development of a novelsimulator for underwater vehicles that takes into account theinteraction between the submerged body and the free watersurface. By using a Lagrangian approach, the potential flowtheory is combined with the 6-DoF equations of motions, inorder to derive a mathematical model of the system dynamics tobe solved in time-domain. A numerical model is then specificallydeveloped and an extensive simulation campaign was carried out,leading to the Submerged Bodies Simulator (SubBoS). The resultshighlight an extension of the prediction validity with respect tothe state-of-art modeling, typically based on stringent hypotheseson body motion (e.g. small displacements, pure surge motion).Index Terms—Lagrangian dynamics, modeling, simulation,potential theory, naval architecture.

Quadrelli, Marco B.↗
NASA NTRS2018

Southern California Megacity CO2, CH4, and CO Flux Estimates Using Ground- and Space-Based Remote Sensing and a Lagrangian Model

We estimate the overall CO2, CH4, and CO flux from the South Coast Air Basin using an inversion that couples Total Carbon Column Observing Network (TCCON) and Orbiting Carbon Observatory-2 (OCO-2) observations, with the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model and the Open-source Data Inventory for Anthropogenic CO2 (ODIAC). Using TCCON data we estimate the direct net CO2 flux from the So-CAB to be 104±26 Tg CO2 yr(exp -1) for the study period of July 2013–August 2016. We obtain a slightly higher estimate of 120±30 Tg CO2 yr(exp -1) using OCO-2 data. These CO2 emission estimates are on the low end of previous work. Our net CH4 (360±90 Gg CH4 y(exp -1)) flux estimate is in agreement with central values from previous top-down studies going back to 2010 (342–440 Gg CH4 yr(exp -1)). CO emissions are estimated at 487±122 Gg CO yr(exp -1), much lower than previous top-down estimates (1440 Gg CO yr(exp -1)). Given the decreasing emissions of CO, this finding is not unexpected. We perform sensitivity tests to estimate how much errors in the prior, errors in the covariance, different inversion schemes, or a coarser dynamical model influence the emission estimates. Overall, the uncertainty is estimated to be 25%, with the largest contribution from the dynamical model. Lessons learned here may help in future inversions of satellite data over urban areas.

Total Carbon Column Observing Network (TCCON)↗
NASA NTRS2003

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗
NASA NTRS1981

Quasi-Lagrangian models of nascent thermals

Slab-symmetric and axisymmetric geometries are considered in the modeling of thermals by considering their boundaries as self-advecting vortex sheets, for the case of a neutral environment without turbulent processes. The initial accelerations for circular cross section thermals are obtained analytically, and a numerical method is developed to simulate the evolution of thermals by time-marching in which the vortex sheet is divided into finite segments whose positions are tracked in a quasi-Lagrangian fashion. Self-advection is considered, and a redistribution procedure is adopted to prevent the segments from bunching unmanageably. The induced field of motion is determined both inside and outside the thermal. It is shown that (1) the axisymmetric thermal rises more quickly than the slab-symmetric thermal, (2) vorticity maxima occur both within the concavities and on the trailing edges, and (3) the leading edges are very smooth.

Rambaldi, S.↗
NASA NTRS2013

CFD Lagrangian Modeling of Water Droplet Transport for ISS Hygiene Activity Application

The goal of this study was to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC) installed in Node 3. Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow enable identifying the paths of water transport. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 5-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain.

Son, Chang H.↗
NASA NTRS1983

A mechanistic model of Eulerian, Lagrangian mean, and Lagrangian ozone transport by steady planetary waves

Ozone transport is calculated for steady, dissipative planetary waves using the Eulerian, Lagrangian mean, and residual circulation. A Lagrangian model of parcel dynamics is used to interpret planetary wave-photochemistry interaction. In chemically active regions the mean field ozone changes are found to be significant only where there are large gradients in chemical sources and sinks along parcel trajectories. The largest changes in the mean field are found in the lower stratosphere and are due to the Lagrangian mean advection. When the Lagrangian mean advection is approximated by the residual circulation, errors in the transport velocities as large as 30 pct may occur.

Rood, R. B.↗
NASA NTRS2004

Lagrangian Assimilation of Satellite Data for Climate Studies in the Arctic

Under this grant we have developed and tested a new Lagrangian model of sea ice. A Lagrangian model keeps track of material parcels as they drift in the model domain. Besides providing a natural framework for the assimilation of Lagrangian data, it has other advantages: 1) a model that follows material elements is well suited for a medium such as sea ice in which an element retains its identity for a long period of time; 2) model cells can be added or dropped as needed, allowing the spatial resolution to be increased in areas of high variability or dense observations; 3) ice from particular regions, such as the marginal seas, can be marked and traced for a long time; and 4) slip lines in the ice motion are accommodated more naturally because there is no internal grid. Our work makes use of these strengths of the Lagrangian formulation.

Lindsay, Ronald W.↗
NASA NTRS1989

Lagrangian photochemical modeling studies of the 1987 Antarctic spring vortex. I - Comparison with AAOE observations

Results from the Lagrangian photochemical model integrated along computed air parcel trajectories intersected by the ER-2 aircraft are presented and compared with AAOE observations. According to the model, the BrO observations made from the ER-2 within the dehydrated denitrified region are consistent with there being approximately 5 parts per trillion by volume of BrO(y) at 428 K in spring. Within the high ClO region, ozone destruction rates are expected to exceed 2 percent/d with approximately 80 percent due to the ClO dimer mechanism.

Jones, R. L.↗
NASA NTRS1990

Overshooting thunderstorm cloud top dynamics as approximated by a linear Lagrangian parcel model with analytic exact solutions

Results are presented from a linear Lagrangian entraining parcel model of an overshooting thunderstorm cloud top. The model, which is similar to that of Adler and Mack (1986), gives analytic exact solutions for vertical velocity and temperature by representing mixing with Rayleigh damping instead of nonlinearly. Model results are presented for various combinations of stratospheric lapse rate, drag intensity, and mixing strength. The results are compared to those of Adler and Mack.

Schlesinger, Robert E.↗
NASA NTRS1987

Comprehensive modeling of turbulent particulate flows using Eulerian and Lagrangian schemes

This paper addresses turbulent particle dispersion and modulation effects in dilute gas-particle turbulent flows using the Eulerian and Lagrangian modeling approaches. Gradient diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangian formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. For the particle size and loading considered, the turbulence transport equations must be modified to account for the modulation effects. Models are proposed for both Eulerian and Lagrangian schemes. Comparisons and predictions are made in fully developed gas-solid pipe flow and confined coaxial jets laden with particles. For the monodispersed system investigated, Eulerian approach is less expensive and gives more consistent results than the Lagrangian approach. The Lagrangian technique should be further developed to eliminate current inherent inconsistencies especially with regard to the symmetry boundary condition.

Adeniji-Fashola, A.↗
NASA NTRS2009

Spatial and Temporal Extrapolation of Disdrometer Size Distributions Based on a Lagrangian Trajectory Model of Falling Rain

Methodologies to improve disdrometer processing, loosely based on mathematical techniques common to the field of particle flow and fluid mechanics, are examined and tested. The inclusion of advection and vertical wind field estimates appear to produce significantly improved results in a Lagrangian hydrometeor trajectory model, in spite of very strict assumptions of noninteracting hydrometeors, constant vertical air velocity, and time independent advection during the scan time interval. Wind field data can be extracted from each radar elevation scan by plotting and analyzing reflectivity contours over the disdrometer site and by collecting the radar radial velocity data to obtain estimates of advection. Specific regions of disdrometer spectra (drop size versus time) often exhibit strong gravitational sorting signatures, from which estimates of vertical velocity can be extracted. These independent wind field estimates become inputs and initial conditions to the Lagrangian trajectory simulation of falling hydrometeors.

Lane, John E.↗
NASA NTRS2006

Establishing Lagrangian Connections between Observations within Air Masses Crossing the Atlantic during the ICARTT Experiment

The International Consortium for Atmospheric Research on Transport and Transformation (ICARTT)-Lagrangian experiment was conceived with an aim to quantify the effects of photochemistry and mixing on the transformation of air masses in the free troposphere away from emissions. To this end attempts were made to intercept and sample air masses several times during their journey across the North Atlantic using four aircraft based in New Hampshire (USA), Faial (Azores) and Creil (France). This article begins by describing forecasts using two Lagrangian models that were used to direct the aircraft into target air masses. A novel technique is then used to identify Lagrangian matches between flight segments. Two independent searches are conducted: for Lagrangian model matches and for pairs of whole air samples with matching hydrocarbon fingerprints. The information is filtered further by searching for matching hydrocarbon samples that are linked by matching trajectories. The quality of these coincident matches is assessed using temperature, humidity and tracer observations. The technique pulls out five clear Lagrangian cases covering a variety of situations and these are examined in detail. The matching trajectories and hydrocarbon fingerprints are shown and the downwind minus upwind differences in tracers are discussed.

Methven, J.↗
NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗
NASA NTRS1989

Lagrangian photochemical modeling studies of the 1987 Antarctic spring vortex. II - Seasonal trends in ozone

A photochemical model consisting of 40 species and 107 reactions is integrated along 80-day air parcel trajectories calculated in the lower stratosphere for the springtime Antarctic. For the trajectory starting at 58 deg S, which may be regarded as outside the circumpolar vortex, only a small change in O3 occurs in the model. In contrast, for the air parcel starting in the vortex at 74 deg S, the O3 concentration is reduced by 93 percent during the 80 days from the beginning of August to late October. The model results for several species are compared with measurements from the Airborne Antarctic Ozone Experiment and, in general, good agreement is obtained. In the model, the dentrification of the air parcels in polar stratospheric clouds increases the amount of chlorine present in active form. Heterogeneous reactions maintain high active chlorine which destroys O3 via the formation of the ClO dimer. Results of calculations with reduced concentrations of inorganic chlorine show considerably reduced O3 destruction rates and compare favorably with the behavior of total O3 since the late 1970s. The remaining major uncertainties in the photochemical aspects of the Antarctic ozone hole are highlighted.

Austin, J.↗
NASA NTRS2008

A Modification and Analysis of Lagrangian Trajectory Modeling and Granular Dynamics of Lunar Dust Particles

A previously developed mathematical model is amended to more accurately incorporate the effects of lift and drag on single dust particles in order to predict their behavior in the wake of high velocity gas flow. The model utilizes output from a CFD or DSMC simulation of exhaust from a rocket nozzle hot gas jet. An extension of the Saffman equation for lift based on the research of McLaughlin (1991) and Mei (1992) is used, while an equation for the Magnus force modeled after the work of Oesterle (1994) and Tsuji et al (1985) is applied. A relationship for drag utilizing a particle shape factor (phi = 0.8) is taken from the work of Haider and Levenspiel (1989) for application to non-spherical particle dynamics. The drag equation is further adjusted to account for rarefaction and compressibility effects in rarefied and high Mach number flows according to the work of Davies (1945) and Loth (2007) respectively. Simulations using a more accurate model with the correction factor (Epsilon = 0.8 in a 20% particle concentration gas flow) given by Richardson and Zaki (1954) and Rowe (1961) show that particles have lower ejection angles than those that were previously calculated. This is more prevalent in smaller particles, which are shown through velocity and trajectory comparison to be more influenced by the flow of the surrounding gas. It is shown that particles are more affected by minor changes to drag forces than larger adjustments to lift forces, demanding a closer analysis of the shape and behavior of lunar dust particles and the composition of the surrounding gas flow.

Long, Jason M.↗
NASA NTRS2008

Lagrangian Trajectory Modeling of Lunar Dust Particles

Apollo landing videos shot from inside the right LEM window, provide a quantitative measure of the characteristics and dynamics of the ejecta spray of lunar regolith particles beneath the Lander during the final 10 [m] or so of descent. Photogrammetry analysis gives an estimate of the thickness of the dust layer and angle of trajectory. In addition, Apollo landing video analysis divulges valuable information on the regolith ejecta interactions with lunar surface topography. For example, dense dust streaks are seen to originate at the outer rims of craters within a critical radius of the Lander during descent. The primary intent of this work was to develop a mathematical model and software implementation for the trajectory simulation of lunar dust particles acted on by gas jets originating from the nozzle of a lunar Lander, where the particle sizes typically range from 10 micron to 500 micron. The high temperature, supersonic jet of gas that is exhausted from a rocket engine can propel dust, soil, gravel, as well as small rocks to high velocities. The lunar vacuum allows ejected particles to travel great distances unimpeded, and in the case of smaller particles, escape velocities may be reached. The particle size distributions and kinetic energies of ejected particles can lead to damage to the landing spacecraft or to other hardware that has previously been deployed in the vicinity. Thus the primary motivation behind this work is to seek a better understanding for the purpose of modeling and predicting the behavior of regolith dust particle trajectories during powered rocket descent and ascent.

Lane, John E.↗