Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid steps”

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 307 records · Page 17

A Semi-Analytical Solution Approach for Solving Constant-Coefficient First-Order Partial Differential Equations

Simulation and control of many dynamic systems involve solving partial differential equations (PDE). This letter proposes a semi-analytical solution (SAS) approach for fast and high-quality solution of first-order PDEs. The region of interest of the studied PDE is divided into a grid, and an SAS is derived for each grid cell in the form of the multivariate polynomials, of which the coefficients are identified using initial value and boundary value conditions. The solutions are solved in a “time-stepping” manner, i.e. within one time step, the coefficients of the SAS are identified and the initial value of the next time step is evaluated. This approach achieves a significantly larger grid cell than the widely used finite difference method, and thus enhances the computational efficiency significantly. Furthermore, the simulation result on the natural gas pipeline model demonstrates the advantages of SAS in accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗

Empowering Rural Electrification in Honduras: An Integrated Assessment of PV/BESS and Productive Uses of Electricity in Gracias a Dios

Honduras faces significant challenges in its energy sector, particularly in rural areas where access to reliable and affordable electricity remains limited. The flagship rural electrification initiative for Honduras Secretary of Energy (SEN) is the Politica de Acceso Universal a la Electricidad (PAUEH - Universal Electricity Access Policy), a key solution for addressing this energy poverty is the deployment of more than 1700 distributed solar and hybrid mini-grid solutions. In late 2023 as a first step towards supporting SEN's electrification efforts, the National Renewable Energy Laboratory (NREL) developed a literature review of SEN electrification policy documents and conducted a series of technical capacity-building workshops with SEN and other energy sector stakeholders in Honduras focused on using NREL's open-source REopt tool to conduct techno-economic assessments and develop least-cost optimizations for potential solar + storage mini-grid systems. Building on some initial capacity building on mini-grid modeling, the National Renewable Energy Laboratory (NREL) worked with SEN to develop a detailed techno-economic assessment for electrifying two schools, a healthcare clinic, and a hospital in a hypothetical community within the department of Gracias a Dios. The analysis also evaluates the business case for cold storage productive use of energy (PUE) applications for the fisheries value chain and how the incorporation of these PUE loads potentially impacts both the viability of the PV+BESS solutions as well as local economic development. By assessing the potential for deployment of integrated PV/BESS systems to both support critical community services like education and healthcare, as well as potential for downstream enterprise and economic development, this analysis represents a first step that can help to inform specific strategies for development of pilot PV+BESS projects aligned with national priorities and sector level planning under PAUEH.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Advanced in Visualization of 3D Time-Dependent CFD Solutions

Numerical simulations of complex 3D time-dependent (unsteady) flows are becoming increasingly feasible because of the progress in computing systems. Unfortunately, many existing flow visualization systems were developed for time-independent (steady) solutions and do not adequately depict solutions from unsteady flow simulations. Furthermore, most systems only handle one time step of the solutions individually and do not consider the time-dependent nature of the solutions. For example, instantaneous streamlines are computed by tracking the particles using one time step of the solution. However, for streaklines and timelines, particles need to be tracked through all time steps. Streaklines can reveal quite different information about the flow than those revealed by instantaneous streamlines. Comparisons of instantaneous streamlines with dynamic streaklines are shown. For a complex 3D flow simulation, it is common to generate a grid system with several millions of grid points and to have tens of thousands of time steps. The disk requirement for storing the flow data can easily be tens of gigabytes. Visualizing solutions of this magnitude is a challenging problem with today's computer hardware technology. Even interactive visualization of one time step of the flow data can be a problem for some existing flow visualization systems because of the size of the grid. Current approaches for visualizing complex 3D time-dependent CFD solutions are described. The flow visualization system developed at NASA Ames Research Center to compute time-dependent particle traces from unsteady CFD solutions is described. The system computes particle traces (streaklines) by integrating through the time steps. This system has been used by several NASA scientists to visualize their CFD time-dependent solutions. The flow visualization capabilities of this system are described, and visualization results are shown.

Lane, David A.↗

PathTrace and MPVEASI: A Path Analysis Comparative Validation Study

Developed in 2018, PathTrace is a software package built with the intention of making path analysis simple and intuitive. PathTrace is a top-down pathway analysis software where a user is able to explore vulnerable pathways into a facility. The intention of utilizing a software tool like PathTrace is to characterize an existing physical protection system (PPS) and to upgrade the system to achieve a high level of response interruption, or probability of interruption (P I ) of the adversary. There are four steps for conducting path analysis using PathTrace. The first step is to identify an image to use to build the model and scale the model within PathTrace using a section of known distance (wall or fence perimeter, for example). The scaling process will produce a grid of cells through which the user is able to build a model. The second step is to fill out the grid of cells with four categories of materials: Barriers, Detection Areas, Jumps, and Targets. These materials apply associated delay and detection values to the cells in which they are applied. The third step is to represent the adversary and response forces. The adversaries are represented by their capabilities in interacting with the materials identified in step two, and the response is represented by how quickly they will be able to respond to an adversary attack. Finally, the user is able to take all of the information from the previous three steps and perform a Most Vulnerable Path (MVP) analysis. In this stage, the user is able to visualize vulnerable adversary pathways and reason about how to upgrade these pathways to provide a high level of P I .

97 MATHEMATICS AND COMPUTING↗

Fortalecimiento de la Electrificación Rural en Honduras: Una Evaluación Integrada de PV+BESS y Usos Productivos de la Electricidad en Gracias a Dios, Honduras [Empowering Rural Electrification in Honduras: An Integrated Assessment of PV/BESS and Productive Uses of Electricity in Gracias a Dios] (Spanish Translation)

Honduras faces significant challenges in its energy sector, particularly in rural areas where access to reliable and affordable electricity remains limited. The flagship rural electrification initiative for Honduras Secretary of Energy (SEN) is the Politica de Acceso Universal a la Electricidad (PAUEH - Universal Electricity Access Policy), a key solution for addressing this energy poverty is the deployment of more than 1700 distributed solar and hybrid mini-grid solutions. In late 2023 as a first step towards supporting SEN's electrification efforts, the National Renewable Energy Laboratory (NREL) developed a literature review of SEN electrification policy documents and conducted a series of technical capacity-building workshops with SEN and other energy sector stakeholders in Honduras focused on using NREL's open-source REopt tool to conduct techno-economic assessments and develop least-cost optimizations for potential solar + storage mini-grid systems. Building on some initial capacity building on mini-grid modeling, the National Renewable Energy Laboratory (NREL) worked with SEN to develop a detailed techno-economic assessment for electrifying two schools, a healthcare clinic, and a hospital in a hypothetical community within the department of Gracias a Dios. The analysis also evaluates the business case for cold storage productive use of energy (PUE) applications for the fisheries value chain and how the incorporation of these PUE loads potentially impacts both the viability of the PV+BESS solutions as well as local economic development. By assessing the potential for deployment of integrated PV/BESS systems to both support critical community services like education and healthcare, as well as potential for downstream enterprise and economic development, this analysis represents a first step that can help to inform specific strategies for development of pilot PV+BESS projects aligned with national priorities and sector level planning under PAUEH. This is the Spanish translation of NREL/TP-7A40-90865.

14 SOLAR ENERGY↗

Infrared-Bolometer Arrays with Reflective Backshorts

Integrated circuits that incorporate square arrays of superconducting-transition- edge bolometers with optically reflective backshorts are being developed for use in image sensors in the spectral range from far infrared to millimeter wavelengths. To maximize the optical efficiency (and, thus, sensitivity) of such a sensor at a specific wavelength, resonant optical structures are created by placing the backshorts at a quarter wavelength behind the bolometer plane. The bolometer and backshort arrays are fabricated separately, then integrated to form a single unit denoted a backshort-under-grid (BUG) bolometer array. In a subsequent fabrication step, the BUG bolometer array is connected, by use of single-sided indium bump bonding, to a readout device that comprises mostly a superconducting quantum interference device (SQUID) multiplexer circuit. The resulting sensor unit comprising the BUG bolometer array and the readout device is operated at a temperature below 1 K. The concept of increasing optical efficiency by use of backshorts at a quarter wavelength behind the bolometers is not new. Instead, the novelty of the present development lies mainly in several features of the design of the BUG bolometer array and the fabrication sequence used to implement the design. Prior to joining with the backshort array, the bolometer array comprises, more specifically, a square grid of free-standing molybdenum/gold superconducting-transition-edge bolometer elements on a 1.4- m-thick top layer of silicon that is part of a silicon support frame made from a silicon-on-insulator wafer. The backshort array is fabricated separately as a frame structure that includes support beams and contains a correspond - ing grid of optically reflective patches on a single-crystal silicon substrate. The process used to fabricate the bolometer array includes standard patterning and etching steps that result in the formation of deep notches in the silicon support frame. These notches are designed to interlock with the support beams on the backshort-array structure to provide structural support and precise relative positioning. The backshort-array structure is inserted in the silicon support frame behind the bolometer array, and the notches in the frame serve to receive the support beams of the backshort-array structure and thus determine the distance between the backshort and bolometer planes. The depth of the notches and, thus, the distance between the backshort and bolometer planes, can be tailored to a value between 25 to 300 m adjusting only a few process steps. The backshort array is designed so as not to interfere with the placement of indium bumps for subsequent indium bump-bonding to the multiplexing readout circuitry

Miller, Timothy M.↗

Earthquake Phase Association Using a Bayesian Gaussian Mixture Model

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual seismic events and play an important role in earthquake monitoring and research. Dense seismic networks and improved phase picking methods produce massive seismic phase datasets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. Here, we present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with a hyperbolic moveout of arrival times and a decay of amplitude with distance. We use the multivariate Gaussian distribution to model the collection of phase picks of an event; and the mean of the multivariate Gaussian distribution is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment to each earthquake and determine earthquake source parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization algorithm. The GaMMA method does not require typical association steps of other algorithms, such as grid-search or supervised training. The results for both synthetic tests and for the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

58 GEOSCIENCES↗

Spatial coupling of gyrokinetic simulations, a generalized scheme based on first-principles

In this paper, we present a scheme that spatially couples two gyrokinetic codes using first-principles. Coupled equations are presented and a necessary and sufficient condition for ensuring accuracy is derived. This new scheme couples both the field and the particle distribution function. The coupling of the distribution function is only performed once every few time-steps, using a five-dimensional (5D) grid to communicate the distribution function between the two codes. This 5D grid interface enables the coupling of different types of codes and models, such as particle and continuum codes, or delta-f and total-f models. Transferring information from the 5D grid to the marker particle weights is achieved using a new resampling technique. Demonstration of the coupling scheme is shown using two XGC gyrokinetic simulations for both the core and edge. We also apply the coupling scheme to two continuum simulations for a one-dimensional advection–diffusion problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experiences with explicit finite-difference schemes for complex fluid dynamics problems on STAR-100 and CYBER-203 computers

Several two- and three-dimensional external and internal flow problems solved on the STAR-100 and CYBER-203 vector processing computers are described. The flow field was described by the full Navier-Stokes equations which were then solved by explicit finite-difference algorithms. Problem results and computer system requirements are presented. Program organization and data base structure for three-dimensional computer codes which will eliminate or improve on page faulting, are discussed. Storage requirements for three-dimensional codes are reduced by calculating transformation metric data in each step. As a result, in-core grid points were increased in number by 50% to 150,000, with a 10% execution time increase. An assessment of current and future machine requirements shows that even on the CYBER-205 computer only a few problems can be solved realistically. Estimates reveal that the present situation is more storage limited than compute rate limited, but advancements in both storage and speed are essential to realistically calculate three-dimensional flow.

Kumar, A.↗

A verification of unsteady Navier-Stokes solutions around oscillating airfoils

A finite difference solution code for the two dimensional Navier-Stokes equations was combined with a moving-grid system. The thin layer Navier-Stokes equations with a turbulence model are solved in a time-accurate manner in order to study the unsteady aerodynamics around airfoils undergoing small amplitude pitching or heaving motions in the transonic regime. The accuracy of the solutions obtained by the use of the present moving-grid technqiue is investigated. The effects of the minimum grid size and the integrating time-step size on the solutions are also checked. Some of the solutions obtained by the present method are compared with experimental results. It is demonstrated that the unsteady aerodynamics around oscillating airfoils can be predicted fairly well by the present code for cases in which the dynamic angle of attack or displacement is small.

Nakamichi, J.↗

Knowledge-based flow field zoning

Automation flow field zoning in two dimensions is an important step towards easing the three-dimensional grid generation bottleneck in computational fluid dynamics. A knowledge based approach works well, but certain aspects of flow field zoning make the use of such an approach challenging. A knowledge based flow field zoner, called EZGrid, was implemented and tested on representative two-dimensional aerodynamic configurations. Results are shown which illustrate the way in which EZGrid incorporates the effects of physics, shape description, position, and user bias in a flow field zoning.

Andrews, Alison E.↗

Knowledge-based flow field zoning

Automation flow field zoning in two dimensions is an important step towards easing the three-dimensional grid generation bottleneck in computational fluid dynamics. A knowledge based approach works well, but certain aspects of flow field zoning make the use of such an approach challenging. A knowledge based flow field zoner, called EZGrid, was implemented and tested on representative two-dimensional aerodynamic configurations. Results are shown which illustrate the way in which EZGrid incorporates the effects of physics, shape description, position, and user bias in a flow field zoning.

Andrews, Alison E.↗

Temporal-adaptive Euler/Navier-Stokes algorithm for unsteady aerodynamic analysis of airfoils using unstructured dynamic meshes

A temporal adaptive algorithm for the time-integration of the two-dimensional Euler or Navier-Stokes equations is presented. The flow solver involves an upwind flux-split spatial discretization for the convective terms and central differencing for the shear-stress and heat flux terms on an unstructured mesh of triangles. The temporal adaptive algorithm is a time-accurate integration procedure which allows flows with high spatial and temporal gradients to be computed efficiently by advancing each grid cell near its maximum allowable time step. Results indicate that an appreciable computational savings can be achieved for both inviscid and viscous unsteady airfoil problems using unstructured meshes without degrading spatial or temporal accuracy.

Kleb, William L.↗

Temporal-adaptive Euler/Navier-Stokes algorithm for unsteady aerodynamic analysis of airfoils using unstructured dynamic meshes

A temporal adaptive algorithm for the time-integration of the two-dimensional Euler or Navier-Stokes equations is presented. The flow solver involves an upwind flux-split spatial discretization for the convective terms and central differencing for the shear-stress and heat flux terms on an unstructured mesh of triangles. The temporal adaptive algorithm is a time-accurate integration procedure which allows flows with high spatial and temporal gradients to be computed efficiently by advancing each grid cell near its maximum allowable time step. Results indicate that an appreciable computational savings can be achieved for both inviscid and viscous unsteady airfoil problems using unstructured meshes without degrading spatial or temporal accuracy.

Kleb, William L.↗

Automated domain decomposition for computational fluid dynamics

Automation of flow-field zoning in two-dimensions is an important step towards easing the three-dimensional grid generation bottleneck in computational fluid dynamics. A knowledge-based approach works well, but several aspects of flow-field zoning make the use of such an approach challenging. A proposed model and language to describe the process of zoning a flow field are presented, followed by a discussion of the implementation of EZGrid, a knowledge-based two-dimensional (2-D) flow-field zoner. Results are shown for representative two-dimensional aerodynamic configurations. Finally, an approach to the evaluation of flow-field zonings is described and used to compare the performance of EZGrid with that of a human expert.

Vogel, Alison Andrews↗

Development of iterative techniques for the solution of unsteady compressible viscous flows

A research effort was initiated at Georgia Tech in February 1991 on the development of efficient techniques for the computation of 2-D and 3-D unsteady compressible flow problems. It was found that in 2-D unsteady viscous flow applications, the generalized minimal residual (GMRES) scheme was able to significantly improve the accuracy and stability characteristics of an existing 2-D ADI (Alternating Direction Implicit) time marching scheme. That is, the GMRES/ADI combination allowed 10 to 20 times larger time steps compared to an ADI scheme. Because the GMRES algorithm requires 5 to 10 times the CPU work compared to the ADI scheme, the combined GMRES/ADI scheme yields a net factor of 2 savings in CPU cost. During the past year, we also experimented with GMRES/multigrid/ADI combination. The purpose of this combination was to compute the low frequency components of the change in the flow properties from one time step to the next on a coarse grid. This strategy reduces the memory requirements of the GMRES method roughly by a factor of 4-8 for steady flow problems.

Sankar, Lakshmi N.↗

Mixing and combustion effects in a sliding-wedge ram accelerator with hydrogen injection

The ram accelerator concept has been envisioned as the basis of a ground-based flight-test facility concept in which the test section follows the ram-accelerated projectile, and allows hydrogen mixing/combustion experiments to be conducted when hydrogen carried by the projectile is injected into the freestream. A numerical simulation is presently conducted for such mixing and combustion; the chemical reactions in question are modeled using a seven-step, seven-species model. A grid-adaptation procedure is used to resolve flow features in areas of high fluid and species gradients.

Ramakrishnan, R.↗