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

Development of the Joint NASA/NCAR General Circulation Model

The Data Assimilation Office at NASA/Goddard Space Flight Center is collaborating with NCAR/CGD in an ambitious proposal for the development of a unified climate, numerical weather prediction, and chemistry transport model which is suitable for global data assimilation of the physical and chemical state of the Earth's atmosphere. A prototype model based on the NCAR CCM3 physics and the NASA finite-volume dynamical core has been built. A unique feature of the NASA finite-volume dynamical core is its advanced tracer transport algorithm on the floating Lagrangian control-volume coordinate. The model currently has a highly idealized ozone production/loss chemistry derived from the observed 2D (latitude-height) climatology of the recent decades. Nevertheless, the simulated horizontal wave structure of the total ozone is in good qualitative agreement with the observed (TOMS). Long term climate simulations and NWP experiments have been carried out. Current up to date status and futur! e plan will be discussed in the conference.

Lin, S.-J.↗
DOE OSTI2020

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1985

Liquid rocket combustor computer code development

The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.

Liang, P. Y.↗
NASA NTRS1991

The formation of antipodal-impact terrains on Mars

The regions antipodal to Mars' three largest impact basins, Hellas, Isidis, and Argyre, were assessed for evidence of impact induced seismically disrupted terrains. Photogeology and preliminary computer modeling using the Simplified Arbitrary Lagrangian-Eulerian computer program suggest such terrain could have been formed by the Hellas and Isidis impacts. Previously, the only assessment of potential surface disruptions on Mars antipodal to large impacts was by Peterson (1978). Herein, the results give credence to Peterson's assertions. For Hellas, the computer model clearly suggests that antipodal pressures were probably strong enough to have fractured the crust and disrupted the surface, and may account for aspects of volcanism at Alba Patera. It is suggested that the Hellas impact produced sufficient focused energy at its antipode to produce deep fracture in the Martian crust centered below the current caldera for Alba Petera. With the generation and evolution of magma in the Tharsis area, the fracture system provided a ready conduit for early stage eruptions of Martians lavas derived from deep in the crust and/or upper mantle, perhaps of komatiitic composition.

Williams, D. A.↗
NASA NTRS2010

The Monotonic Lagrangian Grid for Rapid Air-Traffic Evaluation

The Air Traffic Monotonic Lagrangian Grid (ATMLG) is presented as a tool to evaluate new air traffic system concepts. The model, based on an algorithm called the Monotonic Lagrangian Grid (MLG), can quickly sort, track, and update positions of many aircraft, both on the ground (at airports) and in the air. The underlying data structure is based on the MLG, which is used for sorting and ordering positions and other data needed to describe N moving bodies and their interactions. Aircraft that are close to each other in physical space are always near neighbors in the MLG data arrays, resulting in a fast nearest-neighbor interaction algorithm that scales as N. Recent upgrades to ATMLG include adding blank place-holders within the MLG data structure, which makes it possible to dynamically change the MLG size and also improves the quality of the MLG grid. Additional upgrades include adding FAA flight plan data, such as way-points and arrival and departure times from the Enhanced Traffic Management System (ETMS), and combining the MLG with the state-of-the-art strategic and tactical conflict detection and resolution algorithms from the NASA-developed Stratway software. In this paper, we present results from our early efforts to couple ATMLG with the Stratway software, and we demonstrate that it can be used to quickly simulate air traffic flow for a very large ETMS dataset.

Kaplan, Carolyn↗
NASA NTRSDate not supplied

Enhanced Monitoring and Forecasting of South Asian Air Quality Episodes with Multi-Sensor Satellite Products and Dispersion Modeling

Air pollution and particulates represent a serious environmental, public health, and overarching societal concern in the Hindu Kush Himalaya (HKH) region of south-central Asia, especially during the boreal winter and spring months. Frequent contributors to poor air quality of the region include dust and particulate matter transport from the Middle East region and western India, persistent nocturnal fog and smog during the stable dry monsoon months, and biomass burning during the pre-monsoon months of boreal spring. Nocturnal and multi-day persistent stable fog and smog events can occur during the stable conditions of the dry monsoon months that lead to hazardous visibility reductions (e.g., at major airports such as Delhi, India)and poor air quality, particularly in urban corridors.Our project goal is therefore to design and implement a robust air quality and chemistry observation and modeling system using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) and the Weather Research and Forecasting coupled with Chemistry (WRF-Chem) models. This system effectively assimilates aerosol and trace gas retrievals from geostationary and polar-orbiting satellites to advance the monitoring and prediction capabilities of air quality and visibility reductions throughout the HKH region. For this presentation, we highlight several[multi-spectral] geostationary satellite products and dispersion modeling capabilities being developed and implemented for improved situational awareness and forecasting of air quality episodes.Launched in December 2018, the Advanced Meteorological Imager (AMI) aboard the Geostationary-Korean Multi-Purpose Satellite-2A (GEO-KOMPSAT-2A, or GK2A) provides 16 channels of multi-spectral information similar to the U.S. GOES-16/17 satellites. Proven community recipes for generating multi-spectral Red-Green-Blue (RGB) composite products to detect specific meteorological phenomena are applied to the AMI data for the HKH region. We implemented four such RGB products consisting of truecolor with Rayleigh correction to aid in smoke detection; dust RGB for monitoring larger particulate matter; nighttime microphysics for fog/smog detection;and natural color fire RGB and the shortwave 3.8 micron channel for fire hot-spot detection.We additionally configured HYSPLIT for dust/sand dispersion and concentration forecasts over HKH. The HYSPLIT dispersion simulations are initially being conducted with two different methodologies: (1) dust/sand lofting based on initializing plume releases from visual inspection of the dust RGB imagery; and (2) dust/sand release and dispersion using the internal HYSPLIT algorithm that determines lofting from input wind speeds and prescribed land use. The goal is to develop a near real-time HYSPLIT modeling solution for routinely simulating the transport and dispersion of dust plumes across the HKH region. This presentation illustrates these products and capabilities for various use cases from the 2019 to 2021 period, highlighting a particularly unhealthy episode during late March 2021 across Nepal.

Remote Sensing↗
NASA NTRS1981

Quasi-Langrangian models of nascent thermals

The motions in and around an isolated thermal were studied and rising motion in the core, and sinking motion on the outside were found; while the circulation resembled that of a vortex ring. In an entity cloud model, cloudy thermal is tracked, in a Lagrangian fashion, as a discrete entity; the field of motion in and around the thermal is not explicitly simulated. Field of motion cloud models, in which the equations of motion are numerically integrated on an Eulerian grid were developed. It is shown that the great potential of a hybrid cloud model can combine the simplicity of the entity models with the generality and flexibility of the field-of-motion models. A key problem to be overcome in the development of a hybrid model is the formulation of a mathematical framework within which the cloud dynamics can be represented.

Rambaldi, S.↗
DOE OSTI2021

Modeling iterative reconstruction and displacement field in the large scale structure

The next generation of galaxy surveys like the Dark Energy Spectroscopic Instrument and Euclid will provide datasets orders of magnitude larger than anything available to date. Our ability to model nonlinear effects in late time matter perturbations will be a key to unlock the full potential of these datasets, and the area of initial condition reconstruction is attracting growing attention. Iterative reconstruction developed by Marcel Schmittfull et al. is a technique designed to reconstruct the displacement field from the observed galaxy distribution. The nonlinear displacement field and initial linear density field are highly correlated. Therefore, reconstructing the nonlinear displacement field enables us to extract the primordial cosmological information better than from the late time density field at the level of the two-point statistics. This paper will test to what extent the iterative reconstruction can recover the true displacement field and construct a perturbation theory model for the postreconstructed field. We model the iterative reconstruction process with Lagrangian perturbation theory (LPT) up to third order for dark matter in real space and compare it with N-body simulations. Here we find that the simulated iterative reconstruction does not converge to the nonlinear displacement field, and the discrepancy mainly appears in the shift term, i.e., the term correlated directly with the linear density field. On the contrary, our 3LPT model predicts that the iterative reconstruction should converge to the nonlinear displacement field. We discuss the sources of discrepancy, including numerical noise/artifacts on small scales, and present an ad hoc phenomenological model that improves the agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Modeling of Particle Transport and Impact in Transient, Multiphase Oxidizer-Rich Environments Using Computational Fluid Dynamics and Lagrangian Particle Tracking

Particle impacts from foreign object debris (FOD) in an oxidizer-rich environment have been shown to be a credible and significant risk for ignition and failure of components and vehicles. The risk of particle ignition is dependent on many factors, including particle impact velocity/energy and fluid composition (i.e. fluid phase, species presence, and state). Prediction of particle transport and particle impacts in relevant environments, including transient, multiphase flows, is important for assessing the risk and avoiding repeated expensive system level tests. The capability to model both the fluid dynamics and particle transport in complex three-dimensional geometries and multiphase environments has been developed at Marshall Space Flight Center (MSFC) in the propulsion fluid dynamic branch (ER42) by utilizing computational fluid dynamics (CFD) simulations and Lagrangian particle tracking. Using CFD tools, including density-based Loci/CHEM with a cavitation model and pressure-based Loci/STREAM with a volume of fluids model, MSFC/ER42 has provided particle impact predictions to multiple projects for oxidizer compatibility assessments. Simulations of transient, multiphase environments such as liquid slugs flowing through tubing and cavitating flow in valves have been completed alongside one-way coupled Lagrangian particle tracking to predict particle impacts.

M.V. Fischels↗
NASA NTRS1997

Multi-Body Analysis of a Tiltrotor Configuration

The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.

Ghiringhelli, G. L.↗
DOE OSTI2024

Is the HEFT matching unique?

Physics beyond the Standard Model (BSM) can be described in a consistent and general way through the Higgs effective field theory (HEFT). Measurements of model-independent HEFT coefficients allow one to constrain the parameter space of BSM models via a matching procedure. In this work, we show that this procedure is not unique and depends on the scalings of the parameters of the Lagrangian. As examples, we consider three BSM models: the real singlet extension of the SM with a Z 2 symmetry, the complex singlet extension of the SM and the two Higgs doublet model. We discuss several physical observables, and show that different scalings of the model parameters with the UV scale in the matching to the HEFT can yield quite different results. This complicates the interpretation of HEFT measurements in terms of parameters of BSM models. Additionally, as a byproduct, we report the first matching of the complex singlet extension to the HEFT. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗
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.↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

A numerical simulation of the magnetospheric gate model for the X-ray bursters

A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.

Starrfield, S.↗
DOE OSTI2021

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Anomalous Transport in Three-Dimensional Discrete Fracture Networks: Interplay Between Aperture Heterogeneity and Injection Modes

In this work, we study how the interplay between fracture aperture heterogeneity and tracer injection mode controls fluid flow and tracer transport in three-dimensional (3D) discrete fracture networks (DFNs). The direct 3-D DFN simulations show that tracer injection mode has substantial effects on tracer spreading across all levels of aperture heterogeneity. The key controlling factor for effective transport is the initial Lagrangian velocity distribution, which is determined by the interplay between injection mode and aperture heterogeneity. The fundamental difference between initial Lagrangian velocity distribution and domain-scale Eulerian velocity distribution plays a vital role in determining anomalous transport. We effectively capture the observed anomalous transport using an upscaled transport model that incorporates initial velocity distribution, stationary velocity distribution, velocity correlation length, and average advective tortuosity. With the upscaled transport model, we accurately capture the evolution of Lagrangian velocity distribution and predict longitudinal spreading in 3-D DFN.

58 GEOSCIENCES↗