Convergence of mass transfer particle tracking schemes for the simulation of advection-diffusion-reaction equations
Not Available
97 MATHEMATICS AND COMPUTING↗
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.
Accelerating the Lagrangian Particle Tracking in Hydrologic Modeling to Continental-Scale
Unprecedented climate change and anthropogenic activities have induced increasing ecohydrological problems, which have motivated the development of large-scale modeling for solutions. Water age/quality is as important as water quantity for understanding the water cycle. However, current scientific progress in tracking water parcels at large-scale with high spatiotemporal resolutions is far behind that in water balance/quantity owing to the lack of powerful tools. EcoSLIM is a particle tracking model that works with the hydrologic model ParFlow-CLM, which couples surface-subsurface hydrology with land surface processes. Here, we demonstrate a parallel framework to accelerate EcoSLIM to continental-scale on a distributed, multi-GPU platform with CUDA-Aware MPI. In tests from catchment-, to regional-, and then to continental-scale using 25-million to 1.6-billion particles, EcoSLIM shows significant speedup and excellent parallel performance. The parallel framework is portable to atmospheric and oceanic particle tracking models, where parallelization is inadequate and a standard parallel framework is absent. Parallelized EcoSLIM is a promising tool to accelerate our understanding of the terrestrial water cycle and the upscaling of subsurface hydrology to Earth system models.
54 ENVIRONMENTAL SCIENCES↗
Linking Small-Scale Flight Manoeuvers and Density Profiles to the Vertical Movement of Insects in the Nocturnal Stable Boundary Layer
Huge numbers of insects migrate over considerable distances in the stably-stratified night-time atmosphere with great consequences for ecological processes, biodiversity, ecosystem services and pest management. We used a combination of meteorological radar and lidar instrumentation at a site in Oklahoma, USA, to take a new look at the general assistance migrants receive from both vertical and horizontal airstreams during their long-distance flights. Movement in the nocturnal boundary layer (NBL) presents very different challenges for migrants compared to those prevailing in the daytime convective boundary layer, but we found that Lagrangian stochastic modelling is effective at predicting flight manoeuvers in both cases. A key feature for insect transport in the NBL is the frequent formation of a thin layer of fast-moving air – the low-level jet. Modelling suggests that insects can react rapidly to counteract vertical air movements and this mechanism explains how migrants are retained in the jet for long periods (e.g. overnight, and perhaps for several hours early in the morning). This results in movements over much longer distances than are likely in convective conditions, and is particularly significant for the reintroduction of pests to northern regions where they are seasonally absent due to low winter temperatures.
54 ENVIRONMENTAL SCIENCES↗
Using deep machine learning to interpret proton radiography data from a pulsed power experiment
Deep machine learning is used to analyze a proton radiograph from a tin pulsed power experiment and determine density values for each pixel in the image. Two promising convolutional neural network architectures that have proven to be effective for image analysis in other applications are applied to analyze a proton radiograph and find density values. The process of creating a suitable training dataset is described, involving the Lagrangian hydrodynamic model used for simulations of the experiment, the proton radiography forward model to make synthetic images for training, and the manner in which data augmentation is used to expand the resulting image dataset. It is shown that machine learning not only produces a reasonable density field but is also able to predict features in the density field that are suggested by the proton radiograph but not captured by simulations.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Viscous Gubser flow with conserved charges to benchmark fluid simulations
We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.
Hydrodynamic models↗
Neutrino Properties and Interactions
In this chapter, we present basic concepts of neutrino physics. We start with a brief introduction to the standard model electroweak sector, followed by calculations of some relevant neutrino interaction cross sections. Here, we obtain the oscillation formalism from the Dirac equation, alongside with a self-contained derivation of matter effects from the standard model electroweak Lagrangian. Some key features of neutrino oscillation phenomenology are discussed as well. We then review the most precise oscillation measurements and finalize with some broad discussion of the current open questions in neutrino physics.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Ocean & Geohazard Analysis Tool
The Ocean & Geohazard Analysis (OGA) software tool is designed to summarize insights into key offshore hazards drawing from a diverse set of approaches, including artificial intelligence, machine learning, probabilistic and statistical, and offshore data sources. The offshore hazards that can be analyzed include submarine landslides, extreme wind/wave/current event probabilities, earthquakes, and metocean pathways (CIAM Climatological Isolation and Attraction Model–Climatological Lagrangian Coherent Structures - Submissions - EDX (doe.gov)). Currently, the tool is developed for use in the Gulf of Mexico. The data underlying the offshore hazard analyses can be found here: https://edx.netl.doe.gov/dataset/gulf-of-mexico-risk-analysis-database-gomrad This work was conducted under the Advanced Offshore Research Portfolio, FWP Number 1022409 at National Energy Technology Laboratory, U.S. Dept. of Energy. Disclaimer This project was funded by the United States Department of Energy, National Energy Technology Laboratory, in part, through a site support contract. Neither the United States Government nor any agency thereof, nor any of their employees, nor the support contractor, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.
AIML↗
Methods for Coupling a Lagrangian Discrete Element Sea Ice Model to an Eulerian Ocean Model .
Abstract not provided.
Peterson, Kara↗
Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design
Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.
36 MATERIALS SCIENCE↗
Bulk Spray and Individual Plume Characterization of LPG and Iso-Octane Sprays at Engine-Like Conditions
This study presents experimental and numerical examination of directly injected (DI) propane and iso-octane, surrogates for liquified petroleum gas (LPG) and gasoline, respectively, at various engine like conditions with the overall objective to establish the baseline with regards to fuel delivery required for future high efficiency DI-LPG fueled heavy-duty engines. Sprays for both iso-octane and propane were characterized and the results from the optical diagnostic techniques including high-speed Schlieren and planar Mie scattering imaging were applied to differentiate the liquid-phase regions and the bulk spray phenomenon from single plume behaviors. The experimental results, coupled with high-fidelity internal nozzle-flow simulations were then used to define best practices in CFD Lagrangian spray models. Optical imaging revealed that unlike iso-octane, propane’s spray propagation was fed by its flash boiling, spray collapse, and high degree of vaporization, resulting in a direct proportionality of propane’s penetration length to temperature. These unique features of propane and its variation from iso-octane’s spray pattern, contributed to its classification as an unconventional spray. Appropriate corrections to the injection and breakup models were developed to reproduce the under-expanded jet dynamics and to mimic the flash boiling-driven spray development observed with propane sprays. The simulation results were found to be sensitive to cone and inclusion angles of the blob injector. The current work represents a first assessment of the capability of the commonly available models for engine-spray simulations and highlights the fact that, despite the reasonable agreement obtained in the fuel vapor morphology, the representation of the liquid phase lacks accuracy and requires further model development.
02 PETROLEUM↗
Computer experiments on low density crossed- field electron beams
Lagrangian computer model analysis of long wave diocotron instability in crossed field electron beam
ELECTRON BEAM↗
Calculation of a hollow-cone liquid spray in a uniform air stream
Fluid dynamic computer codes for the simulation of flows in gas turbine engine combustion systems are being developed. NASA is currently sponsoring a two-phase program for the evaluation of the performance of current codes, taking into account also an improvement of accuracy, if needed. The present investigation forms a part of this program. The numerical technique used includes a Lagrangian spray model for liquid fuels. The spray model, in conjunction with the turbulence model, determines the distribution of fuel in the burning zone of the combustor. The numerical technique was applied to a hollow-cone pressure atomizer spraying water into a coflowing confined airstream.
Sturgess, G. J.↗
Spatially random models, estimation theory, and robot arm dynamics
Spatially random models provide an alternative to the more traditional deterministic models used to describe robot arm dynamics. These alternative models can be used to establish a relationship between the methodologies of estimation theory and robot dynamics. A new class of algorithms for many of the fundamental robotics problems of inverse and forward dynamics, inverse kinematics, etc. can be developed that use computations typical in estimation theory. The algorithms make extensive use of the difference equations of Kalman filtering and Bryson-Frazier smoothing to conduct spatial recursions. The spatially random models are very easy to describe and are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. The models can also be used to generate numerically the composite multibody system inertia matrix. This is done without resorting to the more common methods of deterministic modeling involving Lagrangian dynamics, Newton-Euler equations, etc. These methods make substantial use of human knowledge in derivation and minipulation of equations of motion for complex mechanical systems.
Rodriguez, G.↗
Isentropic mass exchange between the Tropics and extratropics in the stratosphere
The isentropic mass exchange between the Tropics and extratropics in the stratosphere is investigated with a semi-Lagrangian transport model for the periods from 1 June to 31 October 1992 and from 1 December 1992 to 30 April 1993 using winds from the U.K. Meteorological Office data assimilation system. Calculations with an idealized, initially zonally symmetric tracer show that in the middle and upper stratosphere the bulk of tropical air is transported into the midlatitudes of the winter hemisphere although there exist quasi-permeable barriers in the subtropics. The transport takes place in the form of planetary-scale 'tongues' of material that are drawn poleward in association with the episodic amplification of planetary-scale waves in high latitudes of the winter hemisphere. Once air of tropical origin is transported to the midlatitudes it is irreversibly mixed with the midlatitude air in the 'surf zone.' Air of tropical origin can, however, hardly penetrate into the interior of the winter polar vortex until the breakdown of the vortex. Transport of tropical air into the midlatitudes of the summer hemisphere is strongly inhibited. In the lower stratosphere, tropical air is transported into the northern and southern midlatitudes. During the period from 1 June to 31 October 1992, the amount of tropical air transported into the Northern Hemisphere is, however, much smaller than that transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the northern midlatitudes and into the southern midlatitudes, and there exist no significant transport barriers in the equatorial area.
Chen, Ping↗
Quasi-horizontal transport and mixing in the Antarctic stratosphre
The quasi-horizontal transport and mixing properties of the Antarctic stratosphere are investigated with a simi-Lagrangian transport model and a 'contour advection' technique for the winter and spring of 1992 using analyzed winds from the United Kingdom Meteorological Office data assimiliation system. Transport calculations show that passive tracers are well mixed inside the polar vortex as well as in the midlatitude 'surf zone.' A the vortex edge, strong radial gradients in the tracer fields are well preserved, and their evolutions follow that of the potential vorticity until some time after the breakdown of the polar vortex. In the middle stratosphere there is little tracer exchange across the vortex edge in August and September. Some vortex air is eroded into the surf zone in filamentary form in October, and very strong exchange of air occurs between high and middle latitudes in November. In the lower stratosphere the vortex is not so isolated from the midlatitudes as in the middle stratosphere, and there is more mass exchange across the vortex edge. Calculations of the lengthening of material contours using the contour advection technique show that in the middle stratosphere, strong stirring (i.e., stretching and folding of material elements) occurs in the inner vortex, while the strongest stirring occurs in the midlatitude surf zone and the weakest occurs at the vortex edge. In the lower strtosphere, strong stirring occurs in the inner vortex. Stirring is moderate at the vortex edge and in the midlatitudes.
Chen, Ping↗
Ignition-and-Growth Modeling of NASA Standard Detonator and a Linear Shaped Charge
The main objective of this study is to quantitatively investigate the ignition and shock sensitivity of NASA Standard Detonator (NSD) and the shock wave propagation of a linear shaped charge (LSC) after being shocked by NSD flyer plate. This combined explosive train was modeled as a coupled Arbitrary Lagrangian-Eulerian (ALE) model with LS-DYNA hydro code. An ignition-and-growth (I&G) reactive model based on unreacted and reacted Jones-Wilkins-Lee (JWL) equations of state was used to simulate the shock initiation. Various NSD-to-LSC stand-off distances were analyzed to calculate the shock initiation (or failure to initiate) and detonation wave propagation along the shaped charge. Simulation results were verified by experimental data which included VISAR tests for NSD flyer plate velocity measurement and an aluminum target severance test for LSC performance verification. Parameters used for the analysis were obtained from various published data or by using CHEETAH thermo-chemical code.
Oguz, Sirri↗
CFD Predictions of Soot & CO Emissions Generated by a Partially-Fueled 9-Element Lean-Direct Injection Combustor
A study was undertaken to investigate the CO & soot emissions generated by a partially-fueled 9- element LDI (Lean-Direct Injection) combustor configuration operating in the idle range of jet engine conditions. In order to perform the CFD analysis, several existing soot/chemistry models were implemented into the OpenNCC (Open National Combustion Code). The calculations were based on a Reynolds-Averaged Navier Stokes (RANS) simulation with standard k-epsilon turbulence model, a 62- species jet-a/air chemistry, a 2-equation soot model, & a Lagrangian spray solver. A separate transport equation was solved for all individual species involved in jet-a/air combustion. In the test LDI configuration we examined, only five of the nine injectors were fueled with the major pilot injector operating at an equivalence ratio of near one and the other four main injectors operating at an equivalence ratio near 0.55. The calculations helped to identify several reasons behind the soot & CO formation in different regions of the combustor. The predicted results were compared with the reported experimental data on soot mass concentration (SMC) & emissions index of CO (EICO). The experimental results showed that an increase in either T3 and/or F/A ratio lead to a reduction in both EICO & SMC. The predicted results were found to be in reasonable agreement. However, the predicted EICO differed substantially in one test condition associated with higher F/A ratio.
Raju, M. S.↗
Defining the upper boundary of the Asian Tropopause Aerosol Layer (ATAL) using the static stability
The Asian Tropopause Aerosol Layer (ATAL) is located in the Upper Troposphere and Lower Stratosphere (UTLS) during the Asian Summer Monsoon. However, what dynamical feature separates the ATAL from the well-known stratospheric ‘Junge layer’ is not yet clear. In this study, using the in-situ (Radiosonde, Ozonesonde, backscatter sonde and cryogenic frost-point hygrometer) observations from multiple locations in India (Gadanki (13.45° N, 79.18° E), Hyderabad (17.47° N, 78.58° E) and Varanasi (25.27° N, 82.99° E)) and multi-satellite observations ((Cloud-Aerosol Lidar and Infrared Pathfinder Observation, (CALIPSO), Atmospheric Chemistry Experiment (ACE) Fourier Transform Spectrometer (FTS) and Constellation Observation System for Meteorology, Ionosphere and Climate (COSMIC) Global Position System (GPS) Radio Occultation (RO) (COSMIC GPS-RO)) we show that the ATAL can exist up to the layer of maximum stability (LmaxS), located a few kilometers above the tropopause, determined using the square of Brunt Väisäla frequency. These in-situ observations over Indian stations collected during the ISRO-NASA Balloon Measurement Campaigns of the Asian Tropopause Aerosol Layer (BATAL) show that the ATAL top can reach up to ∼442 K potential temperature level over the Indian region. The LmaxS delineated from COSMIC GPSRO observations over the Asian Summer Monsoon Anticyclone (ASMA) region indicates that the top of ATAL can reach up to 454 K potential temperature level, which is lower than the earlier Lagrangian transport model predicted 460 K. The temperature inversion at LmaxS acts as a lid and constrains the direct transport of aerosols to higher altitudes.
S.T. Akhil Raj↗