Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid code”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Demand Response in Residential Energy Code: Technical Brief

As buildings account for over 75% of U.S. electricity use, effectively managing their loads can greatly facilitate the transition towards a clean, reliable grid. Grid-interactive efficient buildings (GEBs) combine efficiency and demand flexibility with smart technologies and communication to provide occupant comfort and productivity while serving the grid as a distributed energy resource (DER). In turn, GEBs can play a key role in ensuring access to an affordable, reliable, sustainable, and modern U.S. electric power system. Their national adoption could provide $\$$100-200 billion in U.S. electric power system cost savings over the next two decades. The associated reduction in CO 2 emissions is estimated at 6% per year by 2030 (DOE 2021). Building codes represent standard design practice in the construction industry and continually evolve to include advanced technologies and innovative practices. Historically, national model energy codes establish minimum efficiency requirements for new construction (ICC 2020). Expanding codes to support GEB capabilities is a pivotal step towards realizing demand flexibility in support of a clean grid by addressing capabilities to improve interoperability between smart building systems, the grid, and renewable energy resources. Realizing GEBs requires buildings with automated demand response (DR) capabilities that enable standardized communication with or control of, subject to explicit consumer consent, energy smart appliances or home energy management systems. This is achieved through direct or indirect (i.e., via an aggregator) communication between appliances and the electric grid. Energy codes can also support DR communication standardization and advance the deployment of building-integrated DERs such as energy storage, generation, and electric vehicles (EVs). Incorporating automated DR capabilities in energy codes provides many benefits to the consumers. Specifically, it aligns building electric load demand with intermittent renewable energy source availability, decreases peak load on the electric grid, allows buildings to respond to utility price signals, supports electrical network reliability and market growth of products and processes aligned with clean economic growth. The incorporation of DR into the model residential energy codes was considered for both the 2021 and 2024 International Energy Conservation Code (IECC) code development cycles. The approved DR measures in the 2021 cycle were removed in response to appeals (ICC 2020). Updated language was presented for consideration again for the 2024 IECC, where it was negotiated and again approved, and again removed in response to appeals (ICC 2024). This resulted in many sections, including sections on demand responsive controls, being moved to the credits options or an appendix as a voluntary application. This technical brief updates the proposed DR components such that they can be considered by states and local governments for direct incorporation into their codes, as well as for future IECC energy code development. The proposal refinements are intended to support consistency in approach and provide a degree of certainty for building owners, designers, contractors, manufacturers, and building and fire safety professionals. The scope of this technical brief includes three strategies for DR in residential buildings: 1) smart thermostats with demand-responsive control, 2) electric water heating incorporating demand-responsive controls and communication and 3) grid Integrated solar and energy storage systems.

2021 IECC↗
NASA NTRSDate not supplied

Cross-Validation of Computational and Experimental Distributed Surface Pressures on the Space Launch System

This paper presents a new workflow for comparing experimental pressure-sensitive paint (PSP) data to computational fluid dynamic (CFD) simulations by way of mapping data from corresponding grids utilizing interpolation methods. In addition to generating quantitative and qualitative point-to-point comparisons between PSP and CFD data, this workflow extracts sectional loading data from both grids and generates lineload comparison charts for corresponding PSP and CFD runs. Experimental PSP data presented in this paper were taken from a 2016 NASA Ames Research Center Unitary Plan Wind Tunnel 11- by 11-Foot Transonic WindTunnel Facility test of the NASA Space Launch System. CFD simulation data for comparison purposes were generated using the FUN3D code. Overall, interpolation onto PSP grids versus CFD grids yields comparable surface pressure fields. However, lineload comparisons are easier to make on the CFD grid-mapped data due to the grid topology and the current capabilities of the lineload analysis tools at NASA Langley Research Center. This workflow is written using contemporary software (Python, Tecplot, PyTecplot), is compatible with existing tools at NASA Langley, and is developed to be adaptable depending on the situation.

SLS↗
DOE OSTI2025

Logical error rates for the surface code under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions

With fault-tolerant quantum computing (FTQC) on the horizon, it is critical to understand sources of logical errors in plausible hardware implementations of quantum error-correcting codes. Detailed error modeling of computational instructions on particular FTQC architectures will enable the better prediction of error propagation in FT-encoded quantum circuits while revealing where greater attention is needed in hardware design. In this work, we consider logical error rates for the surface code implemented on a hypothetical grid-based trapped-ion quantum charge-coupled device architecture. Specifically, we construct logical channels for the idling surface code and examine its diamond error under a mixed coherent and stochastic circuit-level noise model inspired by trapped ions. We include the coherent dephasing noise that is known to accumulate during physical qubit idling and transport in these systems, determining idling and transport durations using the time-resolved output of an open-source trapped-ion surface code compiler. To estimate expectation values of logical Pauli observables following hardware circuits containing non-Clifford sources of noise, we utilize a Monte Carlo technique to sample from an underlying quasiprobability distribution of Clifford circuits that we independently simulate in a phase-sensitive fashion. We verify error suppression up to code distance 𝑑 = 11 at coherent dephasing rates near and below those of current-generation trapped-ion quantum computers and find that logical error rates align with those of analogous fully stochastic simulations in this regime. Exploring higher dephasing rates at 𝑑 = 3−5, we find evidence for growing coherent rotations about all three logical Pauli axes, increased diagonal logical error process matrix elements relative to those of stochastic simulations, and a reduced dephasing rate threshold. Overall, our work paves a way toward realistic hardware emulation of small fault-tolerant quantum processes, e.g., members of an FTQC instruction set.

Quantum benchmarking↗
NASA NTRS2004

Further Improvement in 3DGRAPE

3DGRAPE/AL:V2 denotes version 2 of the Three-Dimensional Grids About Anything by Poisson's Equation with Upgrades from Ames and Langley computer program. The preceding version, 3DGRAPE/AL, was described in Improved 3DGRAPE (ARC-14069) NASA Tech Briefs, Vol. 21, No. 5 (May 1997), page 66. These programs are so named because they generate volume grids by iteratively solving Poisson's Equation in three dimensions. The grids generated by the various versions of 3DGRAPE have been used in computational fluid dynamics (CFD). The main novel feature of 3DGRAPE/AL:V2 is the incorporation of an optional scheme in which anisotropic Lagrange-based trans-finite interpolation (ALBTFI) is coupled with exponential decay functions to compute and blend interior source terms. In the input to 3DGRAPE/AL:V2 the user can specify whether or not to invoke ALBTFI in combination with exponential-decay controls, angles, and cell size for controlling the character of grid lines. Of the known programs that solve elliptic partial differential equations for generating grids, 3DGRAPE/AL:V2 is the only code that offers a combination of speed and versatility with most options for controlling the densities and other characteristics of grids for CFD.

Alter, Stephen↗
DOE OSTI2023

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗
DOE OSTI2023

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TREERAY/RADPRESSURE, is an extension to the novel radiative transfer method TREERAY implemented in the grid-based Magneto-Hydrodynamics code FLASH. In TREERAY/RADPRESSURE, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M ⊙ . We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

(ISM:) HII regions↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Study of Geometric Porosity on Static Stability and Drag Using Computational Fluid Dynamics for Rigid Parachute Shapes

This paper explores use of computational fluid dynamics to study the e↵ect of geometric porosity on static stability and drag for NASA's Multi-Purpose Crew Vehicle main parachute. Both of these aerodynamic characteristics are of interest to in parachute design, and computational methods promise designers the ability to perform detailed parametric studies and other design iterations with a level of control previously unobtainable using ground or flight testing. The approach presented here uses a canopy structural analysis code to define the inflated parachute shapes on which structured computational grids are generated. These grids are used by the computational fluid dynamics code OVERFLOW and are modeled as rigid, impermeable bodies for this analysis. Comparisons to Apollo drop test data is shown as preliminary validation of the technique. Results include several parametric sweeps through design variables in order to better understand the trade between static stability and drag. Finally, designs that maximize static stability with a minimal loss in drag are suggested for further study in subscale ground and flight testing.

Greathouse, James S.↗
NASA NTRS1986

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗
NASA NTRS1991

Grid adaptation to multiple functions for applied aerodynamic analysis

A fast and robust method of grid adaptation to multiple functions has been developed for flow analysis in 3D space. It is based on the equidistribution principle and the alternate direction adaptation method. The grid adaptation algorithm has been fully integrated with a space marching flow solver as contained within a 3D Navier-Stokes code and PAB3D-v2. The grid adaptation strategy, details of numerical implementation, examples of application, and potential extensions of the current method are presented in this paper.

Pao, S. P.↗
NASA NTRS1993

Unstructured multigrid through agglomeration

In this work the compressible Euler equations are solved using finite volume techniques on unstructured grids. The spatial discretization employs a central difference approximation augmented by dissipative terms. Temporal discretization is done using a multistage Runge-Kutta scheme. A multigrid technique is used to accelerate convergence to steady state. The coarse grids are derived directly from the given fine grid through agglomeration of the control volumes. This agglomeration is accomplished by using a greedy-type algorithm and is done in such a way that the load, which is proportional to the number of edges, goes down by nearly a factor of 4 when moving from a fine to a coarse grid. The agglomeration algorithm has been implemented and the grids have been tested in a multigrid code. An area-weighted restriction is applied when moving from fine to coarse grids while a trivial injection is used for prolongation. Across a range of geometries and flows, it is shown that the agglomeration multigrid scheme compares very favorably with an unstructured multigrid algorithm that makes use of independent coarse meshes, both in terms of convergence and elapsed times.

Venkatakrishnan, V.↗
DOE OSTI2022

Parallel transport sweeps on two-dimensional cartesian and hexagonal grids

This paper aims to provide a proof of concept for parallel transport sweeps on two-dimensional hexagonal grids for the discrete ordinates transport equation. While the method is an extension of the popular and well-established Koch-Baker-Alcoulffe (KBA) algorithm, there are significant differences between the cartesian and hexagonal grid and thereafter sweep. The most important is the three-way connectivity of hexagons within the grid which creates greater dependencies between the elements. The KBA method in structured orthogonal grids was first implemented in the DRAGON5 code and the method is first described here. The differences in implementation for the hexagonal grid are also described. Benchmark results are also presented, showing roughly 10 times speedup in computational times with roughly 100 processors, in both cases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1991

A general multiblock Euler code for propulsion integration. Volume 2: User guide for BCON, pre-processor for grid generation and GMBE

The BCON is a menu-driven graphics interface program whose input consists of strings or arrays of points generated from a computer aided design (CAD) tool or any other surface geometry source. The user needs to design the block topology and prepare the surface geometry definition and surface grids separately. The BCON generates input files that contain the block definitions and the block relationships required for generating a multiblock volume grid with the EAGLE grid generation package. The BCON also generates the block boundary conditions file which is used along with the block relationship file as input for the general multiblock Euler (GMBE) code (GMBE, volumes 1 and 3).

Su, T. Y.↗
NASA NTRS2015

Grid Convergence for Turbulent Flows(Invited)

A detailed grid convergence study has been conducted to establish accurate reference solutions corresponding to the one-equation linear eddy-viscosity Spalart-Allmaras turbulence model for two dimensional turbulent flows around the NACA 0012 airfoil and a flat plate. The study involved three widely used codes, CFL3D (NASA), FUN3D (NASA), and TAU (DLR), and families of uniformly refined structured grids that differ in the grid density patterns. Solutions computed by different codes on different grid families appear to converge to the same continuous limit, but exhibit different convergence characteristics. The grid resolution in the vicinity of geometric singularities, such as a sharp trailing edge, is found to be the major factor affecting accuracy and convergence of discrete solutions, more prominent than differences in discretization schemes and/or grid elements. The results reported for these relatively simple turbulent flows demonstrate that CFL3D, FUN3D, and TAU solutions are very accurate on the finest grids used in the study, but even those grids are not sufficient to conclusively establish an asymptotic convergence order.

Diskin, Boris↗
NASA NTRS2019

Aerodynamic and Aerothermal Simulations of Mars Concept Vehicles Using Overset DPLR

NASA's Descent System Studies (DSS) Program is studying various concept vehicles to enable landing of heavy payloads on the surface of Mars. While it is desirable to run high-fidelity CFD simulations to accurately assess the aerodynamic and aerothermal effects of various design changes during EDL, it is usually difficult to quickly generate high-quality grids suitable for such analyses. One approach to address this bottleneck in mesh generation is through the use oversetting grids. Although the overset approach is efficient and powerful in solving partial differential equations on complex geometries, new users often find it challenging to apply overset concepts for their simulations. For example, generating hyperbolic grids with sufficient overlap; priority in hole-cutting on multiple overlapping grids; and fixes to assemble overlapping viscous grids at the body surface. The objective of this presentation is to introduce a simple process that combines the advantages of near-body, point-matched, structured grids with oversetting background grids suitable for grid alignment. This approach allows for grids that can be sequenced, reclustering of mesh spacing at the wall, and grid alignment with the bow shock. The current methodology is tested on a Mid-L/D configuration using the overset DPLR code.

Overset Grids↗
NASA NTRS1989

Aerodynamic analysis of three advanced configurations using the TranAir full-potential code

Computational results are presented for three advanced configurations: the F-16A with wing tip missiles and under wing fuel tanks, the Oblique Wing Research Aircraft, and an Advanced Turboprop research model. These results were generated by the latest version of the TranAir full potential code, which solves for transonic flow over complex configurations. TranAir embeds a surface paneled geometry definition in a uniform rectangular flow field grid, thus avoiding the use of surface conforming grids, and decoupling the grid generation process from the definition of the configuration. The new version of the code locally refines the uniform grid near the surface of the geometry, based on local panel size and/or user input. This method distributes the flow field grid points much more efficiently than the previous version of the code, which solved for a grid that was uniform everywhere in the flow field. TranAir results are presented for the three configurations and are compared with wind tunnel data.

Madson, M. D.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analyis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions.

Huff, Dennis L.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the perioidic boundary conditions. Approximate nonreflecting boundary conditons were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distribution for certain flow conditions.

Huff, Dennis L.↗