Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Lagrangian”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

REFAME: Rain Estimation Using Forward Adjusted-Advection of Microwave Estimates

Sensors flying on satellites provide the only practical means of estimating the precipitation that falls over the entire globe, particularly across the vast unpopulated expanses of Earth s oceans. The sensors that observe the Earth using microwave frequencies provide the best data, but currently these are mounted only on satellites flying in "low Earth orbit". Such satellites constantly move across the Earth s surface, providing snapshots of any given location every 12-36 hours. The entire constellation of low-orbit satellites numbers less than a dozen, and their orbits are not coordinated, so a location will frequently go two or more hours between snapshots. "Geosynchronous Earth orbit" (GEO) satellites continuously observe the same region of the globe, allowing them to provide very frequent pictures. For example, the "satellite movies" shown on television come from GEO satellites. However, the sensors available on GEO satellites cannot match the skill of the low-orbit microwave sensors in estimating precipitation. It is perhaps obvious that scientists should try to combine these very different kinds of data, taking advantage of the strengths of each, but this simple concept has proved to be a huge challenge. The scheme in this paper is "Lagrangian", meaning we follow the storm systems, rather than being tied to a fixed grid of boxes on the Earth s surface. Whenever a microwave snapshot occurs, we gladly use the resulting precipitation estimate. Then at all the times between the microwave snapshots we force the storm system to make a smooth transition from one snapshot s values to the next. We know that a lot more changes occur between the snapshots, but this smooth transition the best we can do with the microwave data alone. The key new contribution in this paper is that we also look at the relative variations in the GEO estimates during these in-between times and force the estimated changes in the precipitation to have similar variations. Preliminary testing shows that this approach has enough promise that it should be developed and studied.

Behrangi, Ali↗
NASA NTRS1995

An Extended Lagrangian Method

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Liou, Meng-Sing↗
NASA NTRS1990

Integration of the shallow water equations on the sphere using a vector semi-Lagrangian scheme with a multigrid solver

A vector semi-Lagrangian semi-implicit two-time-level finite-difference integration scheme for the shallow water equations on the sphere is presented. A C-grid is used for the spatial differencing. The trajectory-centered discretization of the momentum equation in vector form eliminates pole problems and, at comparable cost, gives greater accuracy than a previous semi-Lagrangian finite-difference scheme which used a rotated spherical coordinate system. In terms of the insensitivity of the results to increasing timestep, the new scheme is as successful as recent spectral semi-Lagrangian schemes. In addition, the use of a multigrid method for solving the elliptic equation for the geopotential allows efficient integration with an operation count which, at high resolution, is of lower order than in the case of the spectral models. The properties of the new scheme should allow finite-difference models to compete with spectral models more effectively than has previously been possible.

Bates, J. R.↗
DOE OSTI2023

General relativistic simulations of collapsing binary neutron star mergers with Monte Carlo neutrino transport

Recent gravitational wave observations of neutron-star-neutron-star and neutron-star-black-hole binaries appear to indicate that massive neutron stars may not be too uncommon in merging systems. These discoveries have led to an increased interest in the simulation of merging compact binaries involving massive stars. In this paper, we present a first set of evolution of massive neutron star binaries using Monte Carlo radiation transport for the evolution of neutrinos. Here, we study a range of systems, from nearly symmetric binaries that collapse to a black hole before forming a disk or ejecting material, to more asymmetric binaries in which tidal disruption of the lower mass star leads to the production of more interesting postmerger remnants. For the latter type of systems, we additionally study the impact of viscosity on the properties of the outflows, and compare our results to two recent simulations of identical binaries performed with the whiskyTHC code. We find agreement on the black hole properties, disk mass, and mass and velocity of the outflows within expected numerical uncertainties, and some minor but noticeable differences in the evolution of the electron fraction when using a subgrid viscosity model, with viscosity playing a more minor role in our simulations. The method used to account for r-process heating in the determination of the outflow properties appears to have a larger impact on our result than those differences between numerical codes. We also use the simulation with the most ejected material to verify that our newly implemented Lagrangian tracers provide a reasonable sampling of the matter outflows as they leave the computational grid. We note that, given the lack of production of hot outflows in these mergers, the main role of neutrinos in these systems is to set the composition of the postmerger remnant. One of the main potential uses of our simulations is, thus, as improved initial conditions for longer evolutions of such remnants.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Toward Polydisperse Flows With MFIX-EXA

In the presence of large size disparities, single-grid neighbor search algorithms lead to inflated neighbor lists that significantly degrade the performance of Lagrangian particle solvers. If Eulerian–Lagrangian (EL) frameworks are to remain performant when simulating realistic systems, improved neighbor detection approaches must be adopted. To this end, we consider the application of a multigrid neighbor search (MGNS) algorithm in the mfix-exa software package, an exascale EL solver built upon the AMReX library. Here, details regarding the implementation and verification of MGNS are provided along with speedup curves for a bidisperse mixing layer. MGNS is shown to yield up to 15$\times$ speedup on CPU and 6$\times$ speedup on GPU for the problems considered here. The mfix-exa software is then validated for a variety of polydisperse flows. Finally, a brief discussion is given for how dynamic MGNS may be completed, with application to spatially varying particle size distributions.

42 ENGINEERING↗
NASA NTRS2014

Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

Flexible inhibitors are generally used in solid rocket motors (SRMs) as a means to control the burning of propellant. Vortices generated by the flow of propellant around the flexible inhibitors have been identified as a driving source of instabilities that can lead to thrust oscillations in launch vehicles. Potential coupling between the SRM thrust oscillations and structural vibration modes is an important risk factor in launch vehicle design. As a means to predict and better understand these phenomena, a multidisciplinary simulation capability that couples the NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This capability is crucial to the development of NASA's new space launch system (SLS). This paper summarizes the efforts in applying the coupled software to demonstrate and investigate fluid-structure interaction (FSI) phenomena between pressure waves and flexible inhibitors inside reusable solid rocket motors (RSRMs). The features of the fluid and structural solvers are described in detail, and the coupling methodology and interfacial continuity requirements are then presented in a general Eulerian-Lagrangian framework. The simulations presented herein utilize production level CFD with hybrid RANS/LES turbulence modeling and grid resolution in excess of 80 million cells. The fluid domain in the SRM is discretized using a general mixed polyhedral unstructured mesh, while full 3D shell elements are utilized in the structural domain for the flexible inhibitors. Verifications against analytical solutions for a structural model under a steady uniform pressure condition and under dynamic modal analysis show excellent agreement in terms of displacement distribution and eigenmode frequencies. The preliminary coupled results indicate that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor. This insight could have profound implications for SRM and flexible inhibitor designs for current and future launch vehicles including SLS.

Yang, H. Q.↗
NASA NTRS1989

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies was developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite-core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The viscous flow field of the two-dimensional body is calculated on an Eulerian grid. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a rotor wake with the flow about a 4 to 1 elliptic cylinder at 45 degree incidence was calculated for a Reynolds number of 3000. The results demonstrate the significant variations in the lift and drag on the elliptic cylinder in the presence of the interacting rotor wake.

Stremel, Paul M.↗
DOE OSTI2020

A massively parallel and scalable multi-CPU material point method

Harnessing the power of modern multi-GPU architectures, we present a massively parallel simulation system based on the Material Point Method (MPM) for simulating physical behaviors of materials undergoing complex topological changes, self-collision, and large deformations. Our system makes three critical contributions. First, we introduce a new particle data structure that promotes coalesced memory access patterns on the GPU and eliminates the need for complex atomic operations on the memory hierarchy when writing particle data to the grid. Second, we propose a kernel fusion approach using a new Grid-to-Particles-to-Grid (G2P2G) scheme, which efficiently reduces GPU kernel launches, improves latency, and significantly reduces the amount of global memory needed to store particle data. Finally, we introduce optimized algorithmic designs that allow for efficient sparse grids in a shared memory context, enabling us to best utilize modern multi-GPU computational platforms for hybrid Lagrangian-Eulerian computational patterns. We demonstrate the effectiveness of our method with extensive benchmarks, evaluations, and dynamic simulations with elastoplasticity, granular media, and fluid dynamics. In comparisons against an open-source and heavily optimized CPU-based MPM codebase [Fang et al. 2019] on an elastic sphere colliding scene with particle counts ranging from 5 to 40 million, our GPU MPM achieves over 100x per-time-step speedup on a workstation with an Intel 8086K CPU and a single Quadro P6000 GPU, exposing exciting possibilities for future MPM simulations in computer graphics and computational science. Moreover, compared to the state-of-the-art GPU MPM method [Hu et al. 2019a], we not only achieve 2x acceleration on a single GPU but our kernel fusion strategy and Array-of-Structs-of-Array (AoSoA) data structure design also generalizes to multi-GPU systems. Our multi-GPU MPM exhibits near-perfect weak and strong scaling with 4 GPUs, enabling performant and large-scale simulations on a 10243 grid with close to 100 million particles with less than 4 minutes per frame on a single 4-GPU workstation and 134 million particles with less than 1 minute per frame on an 8-GPU workstation.

Wang, Xinlei↗
NASA NTRS1976

Application of a novel finite difference method to dynamic crack problems

A versatile finite difference method (HEMP and HEMP 3D computer programs) was developed originally for solving dynamic problems in continuum mechanics. It was extended to analyze the stress field around cracks in a solid with finite geometry subjected to dynamic loads and to simulate numerically the dynamic fracture phenomena with success. This method is an explicit finite difference method applied to the Lagrangian formulation of the equations of continuum mechanics in two and three space dimensions and time. The calculational grid moves with the material and in this way it gives a more detailed description of the physics of the problem than the Eulerian formulation.

Chen, Y. M.↗
NASA NTRS1992

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗
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↗
DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗
NASA NTRS1953

Diffusion of Heat from a Line Source in Isotropic Turbulence

An experimental and analytical study has been made of some features of the turbulent heat diffusion behind a line heated wire stretched perpendicular to a flowing isotropic turbulence. The mean temperature distributions have been measured with systematic variations in wind speed, size of turbulence-producing grid, and downstream location of heat source. The nature of the temperature fluctuation field has been studied. A comparison of Lagrangian and Eulerian analyses for diffusion in a nondecaying turbulence yields an expression for turbulent-heat-transfer coefficient in terms of turbulence velocity and a Lagrangian "scale." the ratio of Eulerian to Lagrangian microscale has been determined theoretically by generalization of a result of Heisenberg and with arbitrary constants taken from independent sources, shows rough agreement with experimental results. A convenient form has been deduced for the criterion of interchangeability of instantaneous space and time derivatives in a flowing turbulence.

Uberoi, Mahinder S↗
NASA NTRS1993

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗
NASA NTRS1991

A new Lagrangian random choice method for steady two-dimensional supersonic/hypersonic flow

Glimm's (1965) random choice method has been successfully applied to compute steady two-dimensional supersonic/hypersonic flow using a new Lagrangian formulation. The method is easy to program, fast to execute, yet it is very accurate and robust. It requires no grid generation, resolves slipline and shock discontinuities crisply, can handle boundary conditions most easily, and is applicable to hypersonic as well as supersonic flow. It represents an accurate and fast alternative to the existing Eulerian methods. Many computed examples are given.

Loh, C. Y.↗
DOE OSTI2024

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↗
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.↗
DOE OSTI2021

Development of a Filtered CFD-DEM Drag Model with Multiscale Markers Using an Artificial Neural Network and Nonlinear Regression

Here, the accuracy of coarse-grained Euler-Lagrangian simulations of fluidized beds heavily depends on the mesoscale drag models to account for the influences of the unresolved sub-grid structures. Traditional filtered drag models are regressed with mesoscale markers such as voidage and slip velocities. In this research, a filtered drag was regressed with both mesoscale and macro-scale markers using fine grid Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) simulations. The traditional non-linear regression method was compared with machine learning regression using an Artificial Neural Network (ANN) implemented in PyTorch and coupled with MFiX. The new drag showed higher accuracy than the Wen-Yu drag and another filtered drag derived from the two-fluid model. The nonlinear regression shows slightly better results than ANN regression in cases with similar R 2 values. The utilization of the gas inlet velocity as an additional macro-scale marker reduced the errors by up to 55.3% in the tested cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗