Impact Modeling with Explicit Eulerian to Explicit Lagrangian Coupling Using the Modified Immersed Finite Element Method
Explore the source record for details and available documents.
Thomas, Jesse David [Sandia National Laboratories ↗
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.
Modeling of Cavity Residence Time Via Lagrangian Particle Tracking
Explore the source record for details and available documents.
Langevin↗
Realistic fuel spray modeling for gasoline direct injection engine applications
Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.
33 ADVANCED PROPULSION SYSTEMS↗
Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method
Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.
33 ADVANCED PROPULSION SYSTEMS↗
LES on unstructured deforming meshes: Towards reciprocating IC engines
A variable explicit/implicit characteristics-based advection scheme that is second-order accurate in space and time has been developed recently for unstructured deforming meshes (O'Rourke & Sahota 1996a). To explore the suitability of this methodology for Large-Eddy Simulation (LES), three subgrid-scale turbulence models have been implemented in the CHAD CFD code (O'Rourke & Sahota 1996b): a constant-coefficient Smagorinsky model, a dynamic Smagorinsky model for flows having one or more directions of statistical homogeneity, and a Lagrangian dynamic Smagorinsky model for flows having no spatial or temporal homogeneity (Meneveau et al. 1996). Computations have been made for three canonical flows, progressing towards the intended application of in-cylinder flow in a reciprocating engine. Grid sizes were selected to be comparable to the coarsest meshes used in earlier spectral LES studies. Quantitative results are reported for decaying homogeneous isotropic turbulence, and for a planar channel flow. Computations are compared to experimental measurements, to Direct-Numerical Simulation (DNS) data, and to Rapid-Distortion Theory (RDT) where appropriate. Generally satisfactory evolution of first and second moments is found on these coarse meshes; deviations are attributed to insufficient mesh resolution. Issues include mesh resolution and computational requirements for a specified level of accuracy, analytic characterization of the filtering implied by the numerical method, wall treatment, and inflow boundary conditions. To resolve these issues, finer-mesh simulations and computations of a simplified axisymmetric reciprocating piston-cylinder assembly are in progress.
Haworth, D. C.↗
Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint
An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.
Li, Y.↗
Space-based Ornithology-Studying Bird Migration and Environmental Change in North America
Natural fluctuations in the availability of critical stopover sites coupled with anthropogenic destruction of wetlands, land-use change, and anticipated losses due to climate change present migratory birds with a formidable challenge. We have developed an individual-based, spatially explicit bird migration model that simulates the migration routes, timing and energy budgets of individual birds under dynamic weather and land surface conditions. Our model incorporates biophysical constraints, individual bird energy status, bird behavior, and flight aerodynamics. We model the speed, direction, and timing of individual birds moving through a user specified Lagrangian grid. The model incorporates environmental properties including wind speed and direction, topography, dynamic hydrologic properties of the landscape, and environmental suitability. The model is driven by important variables estimated from satellite observations of the land surface, by data assimilation products from weather and climate models, and biological field data. We illustrate the use of the model to study the impact of both short- and long-term environmental variatios, e.g. climate, drought, anthropogenic, on migration timing (phenology), spatial pattern, and fitness (survival and reproductive success). We present several theoretical simulations of the spring migration of Pectoral Sandpiper (Calidris melanotos) in North America with emphasis on the Central flyway from the Gulf of Mexico to Alaska.
Smith, James↗
Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids
The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.
Raju, M. S.↗
Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows
The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.
Source record↗
A relation between the Lagrangian and Eulerian turbulent velocity autocorrelations
Direct computation of the Lagrangian autocorrelation is not feasible since generally one cannot measure the turbulent velocity of each fluid particle. A model whereby the Lagrangian autocorrelation is determined in terms of a domain integral of a set of regular Eulerian autocorrelations is advanced. The Eulerian autocorrelations are to be acquired concurrently at all positions in the flow field. Three novel averaging procedures are utilized for obtaining the relationship between the Lagrangian and Eulerian autocorrelations. This relationship is not constrained to either homogeneous or isotropic turbulence.
Koper, C. A., Jr.↗
Evaporating Spray in Supersonic Streams Including Turbulence Effects
Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.
Balasubramanyam, M. S.↗
A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes
Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.
42 ENGINEERING↗
Learning effective physical laws for generating cosmological hydrodynamics with Lagrangian deep learning
The goal of generative models is to learn the intricate relations between the data to create new simulated data, but current approaches fail in very high dimensions. When the true data-generating process is based on physical processes, these impose symmetries and constraints, and the generative model can be created by learning an effective description of the underlying physics, which enables scaling of the generative model to very high dimensions. In this work, we propose Lagrangian deep learning (LDL) for this purpose, applying it to learn outputs of cosmological hydrodynamical simulations. The model uses layers of Lagrangian displacements of particles describing the observables to learn the effective physical laws. The displacements are modeled as the gradient of an effective potential, which explicitly satisfies the translational and rotational invariance. The total number of learned parameters is only of order 10, and they can be viewed as effective theory parameters. Further, we combine N-body solver fast particle mesh (FastPM) with LDL and apply it to a wide range of cosmological outputs, from the dark matter to the stellar maps, gas density, and temperature. The computational cost of LDL is nearly four orders of magnitude lower than that of the full hydrodynamical simulations, yet it outperforms them at the same resolution. We achieve this with only of order 10 layers from the initial conditions to the final output, in contrast to typical cosmological simulations with thousands of time steps. This opens up the possibility of analyzing cosmological observations entirely within this framework, without the need for large dark-matter simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
Cross-correlation of Planck CMB lensing with DESI-like LRGs
ABSTRACT Cross-correlations between the lensing of the cosmic microwave background (CMB) and other tracers of large-scale structure provide a unique way to reconstruct the growth of dark matter, break degeneracies between cosmology and galaxy physics, and test theories of modified gravity. We detect a cross-correlation between Dark Energy Spectroscopic Instrument (DESI)-like luminous red galaxies (LRGs) selected from DECam Legacy Survey imaging and CMB lensing maps reconstructed with the Planck satellite at a significance of S/N = 27.2 over scales ℓmin = 30, ℓmax = 1000. To correct for magnification bias, we determine the slope of the LRG cumulative magnitude function at the faint limit as s = 0.999 ± 0.015, and find corresponding corrections of the order of a few per cent for $C^{\kappa g}_{\ell }, C^{gg}_{\ell }$ across the scales of interest. We fit the large-scale galaxy bias at the effective redshift of the cross-correlation zeff ≈ 0.68 using two different bias evolution agnostic models: a HaloFit times linear bias model where the bias evolution is folded into the clustering-based estimation of the redshift kernel, and a Lagrangian perturbation theory model of the clustering evaluated at zeff. We also determine the error on the bias from uncertainty in the redshift distribution; within this error, the two methods show excellent agreement with each other and with DESI survey expectations.
79 ASTRONOMY AND ASTROPHYSICS↗
Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models
In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.
fusion plasma↗
Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport
In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.
97 MATHEMATICS AND COMPUTING↗
Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models
We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the onset of hierarchical clustering. This success is found at a considerable higher non-linearity than is usual for perturbation theory. Whether a truncation of the initial power-spectrum in hierarchical models retains this improvement will be analyzed in a forthcoming work.
Buchert, T.↗
Constraining the Sulfur Dioxide Degassing Flux from Turrialba Volcano, Costa Rica Using Unmanned Aerial System Measurements
Observed sulfur dioxide (SO2)mixing ratios onboard unmanned aerial systems (UAS) duringMarch 11-13, 2013 are used to constrain the three-day averaged SO2 degassing flux fromTurrialba volcanowithin a Bayesian inverse modeling framework. A mesoscale model coupled with Lagrangian stochastic particle backward trajectories is used to quantify the source-receptor relationships at very high spatial resolutions (i.e., b1 km). The model shows better performance in reproducing the near-surface meteorological properties and observed SO2 variations when using a first-order closure non-local planetary boundary layer (PBL) scheme. The optimized SO2 degassing fluxes vary from 0.59 +/- 0.37 to 0.83 +/- 0.33 kt d−1 depending on the PBL scheme used. These fluxes are in good agreement with ground-based gas flux measurements, and correspond to corrective scale factors of 8-12 to the posteruptive SO2 degassing rate in the AeroCom emission inventory. The maximum a posteriori solution for the SO2 flux is highly sensitive to the specification of prior and observational errors, and relatively insensitive to the SO2 loss term and temporal averaging of observations. Our results indicate relatively low degassing activity but sustained sulfur emissions from Turrialba volcano to the troposphere during March 2013. This study demonstrates the utility of low-cost small UAS platforms for volcanic gas composition and flux analysis.
Xi, Xin↗