SEARCH · Search NASA

Results for “Grid Simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Structure Damage Simulations Accounting for Inertial Effects and Impact and Optimization of Grid-Stiffened Non-Circular Shells

The goal of this research project is to develop modelling and analysis strategy for the penetration of aluminium plates impacted by titanium impactors. Finite element analysis is used to study the penetration of aluminium plates impacted by titanium impactors in order to study the effect of such uncontained engine debris impacts on aircraft-like skin panels. LS-DYNA3D) is used in the simulations to model the impactor, test fixture frame and target barrier plate. The effects of mesh refinement, contact modeling, and impactor initial velocity and orientation were studied. The research project also includes development of a design tool for optimum design of grid-stiffened non-circular shells or panels subjected to buckling.

Mei, Chuh↗

NASA NTRS2018

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on a Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability were used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids. Results showed poor agreement in predicted inlet performance on the refined grids compared to experimental data. This was true regardless of whether the refinement was adjoint-based or manual, the cell type in/near the boundary-layer regions, or the grid adaptation metric used. In addition, the 40-probe total pressure recovery was shown to decrease asymptotically as the number of adaptation cycles is increased. Solutions on the unadapted grids generally had better agreement with experimental data than their refined grid counterparts.

Inlet Performance↗

NASA NTRS2018

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Reynolds Averaged Navier-Stokes (RANS) simulations were performed on the Lockheed Martin Quiet SuperSonic Technology (QueSST) aircraft preliminary design to assess inlet performance. The FUN3D flow solver and its adjoint-based grid refinement capability was used for the simulations in hopes of determining internal "best practices" for predicting inlet performance on top-aft-mounted inlets. Several parameters were explored including tetrahedral vs. pentahedral cells in/around the boundary-layer regions, an engine axis-aligned linear pressure sensor vs. a pressure box objective as the grid adaptation metric, and the number of grid adaptation cycles performed. Additional simulations were performed on manually refined grids for comparison with the adjoint-based adapted grids

Adjoint-based Grid Adaptation↗

NASA NTRS2018

Numerical Simulations of a Quiet SuperSonic Technology (QueSST) Aircraft Preliminary Design

Adjoint-based Grid Adaptation↗

DOE OSTI2013

MAGNETOHYDRODYNAMIC MODELING OF SOLAR SYSTEM PROCESSES ON GEODESIC GRIDS

This report describes a new magnetohydrodynamic numerical model based on a hexagonal spherical geodesic grid. The model is designed to simulate astrophysical flows of partially ionized plasmas around a central compact object, such as a star or a planet with a magnetic field. The geodesic grid, produced by a recursive subdivision of a base platonic solid (an icosahedron), is free from control volume singularities inherent in spherical polar grids. Multiple populations of plasma and neutral particles, coupled via charge-exchange interactions, can be simulated simultaneously with this model. Our numerical scheme uses piecewise linear reconstruction on a surface of a sphere in a local two-dimensional 'Cartesian' frame. The code employs Haarten-Lax-van-Leer-type approximate Riemann solvers and includes facilities to control the divergence of the magnetic field and maintain pressure positivity. Several test solutions are discussed, including a problem of an interaction between the solar wind and the local interstellar medium, and a simulation of Earth's magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Quantum Simulation of the First-Quantized Pauli-Fierz Hamiltonian

We provide an explicit recursive divide-and-conquer approach for simulating quantum dynamics and derive a discrete first-quantized nonrelativistic QED Hamiltonian based on the many-particle Pauli-Fierz Hamiltonian. We apply this recursive divide-and-conquer algorithm to this Hamiltonian and compare it to a concrete simulation algorithm that uses qubitization. Our divide-and-conquer algorithm, using lowest-order Trotterization, scales for fixed grid spacing as O ~ ( Λ N 2 η 2 t 2 / ϵ ) for grid size N , η particles, simulation time t , field cutoff Λ , and error ϵ . Our qubitization algorithm scales as O ~ ( N ( η + N ) ( η + Λ 2 ) t log ( 1 / ϵ ) ) . This shows that even a naive partitioning and low-order splitting formula can yield, through our divide-and-conquer formalism, superior scaling to qubitization for large Λ . We compare the relative costs of these two algorithms on systems that are relevant for applications such as the spontaneous emission of photons and the photoionization of electrons. We observe that for different parameter regimes, one method can be favored over the other. Finally, we give new algorithmic and circuit-level techniques for gate optimization, including a new way of implementing a group of multicontrolled- X gates that can be used for better analysis of circuit cost. Published by the American Physical Society 2024

Mukhopadhyay, Priyanka (ORCID:0000000164639100)↗

DOE OSTI2017

Galactic Angular Momentum in Cosmological Zoom-in Simulations. I. Disk and Bulge Components and the Galaxy–Halo Connection

We investigate the angular momentum evolution of four disk galaxies residing in Milky-Way–sized halos formed in cosmological zoom-in simulations with various sub-grid physics and merging histories. We decompose these galaxies, kinematically and photometrically, into their disk and bulge components. The simulated galaxies and their components lie on the observed sequences in the j {sub *}–M {sub *} diagram, relating the specific angular momentum and mass of the stellar component. We find that galaxies in low-density environments follow the relation j{sub ∗} ∝ M{sub ∗}{sup α} past major mergers, with α∼0.6 in the case of strong feedback, when bulge-to-disk ratios are relatively constant, and α∼1.4 in the other cases, when secular processes operate on shorter timescales. We compute the retention factors (i.e., the ratio of the specific angular momenta of stars and dark matter) for both disks and bulges and show that they vary relatively slowly after averaging over numerous but brief fluctuations. For disks, the retention factors are usually close to unity, while for bulges, they are a few times smaller. Our simulations therefore indicate that galaxies and their halos grow in a quasi-homologous way.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2009

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗

DOE OSTI2015

Icing modelling in NSMB with chimera overset grids

In aerospace Engineering, the accurate simulation of ice accretion is a key element to increase flight safety and avoid accidents related to icing effects. The icing code developed in the NSMB solver is based on an Eulerian formulation for droplets tracking, an iterative Messinger model using a modified water runback scheme for ice thickness calculation and mesh deformation to track the ice/air interface through time. The whole process is parallelized with MPI and applied with chimera grids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Stochastic superparameterization in quasigeostrophic turbulence

In this article we expand and develop the authors' recent proposed methodology for efficient stochastic superparameterization algorithms for geophysical turbulence. Geophysical turbulence is characterized by significant intermittent cascades of energy from the unresolved to the resolved scales resulting in complex patterns of waves, jets, and vortices. Conventional superparameterization simulates large scale dynamics on a coarse grid in a physical domain, and couples these dynamics to high-resolution simulations on periodic domains embedded in the coarse grid. Stochastic superparameterization replaces the nonlinear, deterministic eddy equations on periodic embedded domains by quasilinear stochastic approximations on formally infinite embedded domains. The result is a seamless algorithm which never uses a small scale grid and is far cheaper than conventional SP, but with significant success in difficult test problems. Various design choices in the algorithm are investigated in detail here, including decoupling the timescale of evolution on the embedded domains from the length of the time step used on the coarse grid, and sensitivity to certain assumed properties of the eddies (e.g. the shape of the assumed eddy energy spectrum). We present four closures based on stochastic superparameterization which elucidate the properties of the underlying framework: a ‘null hypothesis’ stochastic closure that uncouples the eddies from the mean, a stochastic closure with nonlinearly coupled eddies and mean, a nonlinear deterministic closure, and a stochastic closure based on energy conservation. The different algorithms are compared and contrasted on a stringent test suite for quasigeostrophic turbulence involving two-layer dynamics on a β-plane forced by an imposed background shear. The success of the algorithms developed here suggests that they may be fruitfully applied to more realistic situations. They are expected to be particularly useful in providing accurate and efficient stochastic parameterizations for use in ensemble-based state estimation and prediction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

NASA NTRS2002

Numerical Simulation of the Combustion of Fuel Droplets: Finite Rate Kinetics and Flame Zone Grid Adaptation (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved hndamental understanding droplet combustion may contribute to the clean and safe utilization of fossil hels (Williams, Dryer, Haggard & Nayagam, 1997, 72). The present state of knowledge on convective extinction of he1 droplets derives fiom experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA s Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗

DOE OSTI2015

Geographic Visualization of Power-Grid Dynamics

The visualization enables the simulation analyst to see changes in the frequency through time and space. With this technology, the analyst has a bird's eye view of the frequency at loads and generators as the simulated power system responds to the loss of a generator, spikes in load, and other contingencies. The significance of a contingency to the operation of an electrical power system depends critically on how the resulting tansients evolve in time and space. Consequently, these dynamic events can only be understood when seen in their proper geographic context. this understanding is indispensable to engineers working on the next generation of distributed sensing and control systems for the smart grid. By making possible a natural and intuitive presentation of dynamic behavior, our new visualization technology is a situational-awareness tool for power-system engineers.

Sukumar, SreenivasR.↗

DOE OSTI2015

Geographic Visualization of Power-Grid Dynamics

Sukumar, Sreenivas R.↗

DOE OSTI2026

Digital Grid Twin–Direct Communication Scheme Test Bed for Assessing Relay-to-Relay Radio Antenna and Optical Fiber Performance and Misoperations

This study introduces a novel “Digital Grid Twin–Direct Communication Scheme” test bed. This advanced platform evaluates point-to-point communication between transmitter and receiver relays with optical fiber and radio omnidirectional antenna systems, implemented at the Advanced Protection lab in the Grid Research Innovation and Development Center at Oak Ridge National Laboratory. The increased diversity of energy sources has led to more protective relay misoperations. In North America, microgrid protection schemes now use point-to-point communication along distribution lines between relays to implement advanced logic in nonradial grids that include both high- and low-inertia generators. This trend challenges utilities to minimize misoperations while ensuring rapid fault clearance and accurate selectivity coordination between primary and backup relays. This study assesses relay-to-relay communication schemes by introducing an advanced testing platform based on a digital grid twin protection test bed using a synchronized time source system. The platform evaluates the communication system using radio antennas or optical fiber links by integrating protective relays that operate breakers within the digital twin and record relay events and communication signals. In the experiments, transmitter and receiver relays were configured with inverse time overcurrent and breaker trip detection logic to assess the total time of the communication protection schemes based on the sum of the relay protection element operating time, radio latency, propagation delay, baud rate delay, and relay processing time. These delays were derived from recorded relay events and communication signals from the interface of a real-time simulator set as a digital grid twin. The test bed successfully simulated various electrical faults along a distribution line while ensuring effective and reliable point-to-point communication between transmitter and receiver relays. The radio antenna communication system exhibited latency because of the radio. This latency depends on the baud rate setting and type of radio application; in general, the higher the baud rate, the lower the radio latency. The measured radio latency (for Mirrored Bits with an encryption card at 9,600 bps) was about 9–10 ms. Additionally, calculated propagation delay per mile for radio antennas and optical fiber was 5.36 µs/mi and 8.04 µs/mi, respectively. Optical fiber communication did not demonstrate radio latency. Instead, the protection element operating time depends mainly on the protection logic function set in the relay, and the relay processing time depends on the processing rate of the relay in samples per power system cycle.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗

DOE OSTI2026

Regional Earthquake Ground Motion Simulations for Southern California With EQSIM: Insights From the 2008 Chino Hills, 2024 Highland Park, and 2021 Carson Earthquakes

This study presents physics-based, 3D simulations using the EQSIM framework for several earthquakes in the Los Angeles region. The primary objective was to assess the ability of deterministic physics-based ground motion simulations to reproduce the observed motions from historical events. The selected events included the mathematical equation M w 5.4 2008 Chino Hills, the mathematical equation M w 4.4 2024 Highland Park, and the mathematical equation M w 4.3 2021 Carson events. The simulated motions were evaluated by comparing the recorded and simulated seismograms, as well as the Fourier amplitude spectra, across multiple seismic stations. The SCEC 3D velocity model, CVM-S4.26.M01, was used to represent the regional geology, and ground motion simulations were carried out with a resolution of up to 5 Hz. The results indicate that the simulated motions captured the recorded motions up to approximately 4 Hz. While careful iterations regarding source parameters and corner frequencies were required, and, for the case of the Highland Park event, some of the near-source stations had relatively low accuracy, the present study established a positive step toward the utilization of physics-based simulations in practical applications. The computational efficiencies exhibited by EQSIM, especially on GPU clusters, further supported this assertion, as wall-clock times of simulations involving more than 10 billion grid points were as low as mathematical equation minutes. This permits ensemble simulations for a considered scenario event so that modeling uncertainties (e.g., source and geology) can be bracketed.

EQSIM↗

DOE OSTI2018

CWRF performance at downscaling China climate characteristics

The performance of the regional Climate-Weather Research and Forecasting model (CWRF) for downscaling China climate characteristics is evaluated using a 1980–2015 simulation at 30 km grid spacing driven by the ECMWF Interim reanalysis (ERI). It is shown that CWRF outperforms the popular Regional Climate Modeling system (RegCM4.6) in key features including monsoon rain bands, diurnal temperature ranges, surface winds, interannual precipitation and temperature anomalies, humidity couplings, and 95th percentile daily precipitation. Even compared with ERI, which assimilates surface observations, CWRF better represents the geographic distributions of seasonal mean climate and extreme precipitation. These results indicate that CWRF may significantly enhance China climate modeling capabilities.

Meteorology & Atmospheric Sciences↗