SEARCH · Engineering Papers

Results for “Grid computing”

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,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Grid smoothing and orthogonalization procedures were developed and implemented in the construction of two and three dimensional grids. The procedures are based on the variational methods of grid generation. The two-dimensional examples were computed using the MSU IRIS Graphics Workstation. It was demonstrated that the elliptic grid generation equations, with arbitrary forcing functions, can be solved, in their variational formulation, using a gradient method. Since gradient methods have a global convergence property, the divergence problems often encountered when using SOR iterative methods can be avoided. It is not to be concluded, however, that SOR methods should be abandoned, since gradient methods tend to converge very slowly. In fact, slow convergence was the major problem encountered in the three-dimensional grids. Further progress was made on the continuing effort to develop conservative interpolation formulas for overlapping grids.

Mastin, C. Wayne↗

NASA NTRS1993

The three-dimensional Multi-Block Advanced Grid Generation System (3DMAGGS)

As the size and complexity of three dimensional volume grids increases, there is a growing need for fast and efficient 3D volumetric elliptic grid solvers. Present day solvers are limited by computational speed and do not have all the capabilities such as interior volume grid clustering control, viscous grid clustering at the wall of a configuration, truncation error limiters, and convergence optimization residing in one code. A new volume grid generator, 3DMAGGS (Three-Dimensional Multi-Block Advanced Grid Generation System), which is based on the 3DGRAPE code, has evolved to meet these needs. This is a manual for the usage of 3DMAGGS and contains five sections, including the motivations and usage, a GRIDGEN interface, a grid quality analysis tool, a sample case for verifying correct operation of the code, and a comparison to both 3DGRAPE and GRIDGEN3D. Since it was derived from 3DGRAPE, this technical memorandum should be used in conjunction with the 3DGRAPE manual (NASA TM-102224).

Alter, Stephen J.↗

NASA NTRS1994

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗

NASA NTRS1985

Computation of unsteady transonic aerodynamics with steady state fixed by truncation error injection

A novel technique is introduced for efficient computations of unsteady transonic aerodynamics. The steady flow corresponding to body shape is maintained by truncation error injection while the perturbed unsteady flows corresponding to unsteady body motions are being computed. This allows the use of different grids comparable to the characteristic length scales of the steady and unsteady flows and, hence, allows efficient computation of the unsteady perturbations. An example of typical unsteady computation of flow over a supercritical airfoil shows that substantial savings in computation time and storage without loss of solution accuracy can easily be achieved. This technique is easy to apply and requires very few changes to existing codes.

Fung, K.-Y.↗

NASA NTRS1996

Research in Parallel Algorithms and Software for Computational Aerosciences

Phase I is complete for the development of a Computational Fluid Dynamics parallel code with automatic grid generation and adaptation for the Euler analysis of flow over complex geometries. SPLITFLOW, an unstructured Cartesian grid code developed at Lockheed Martin Tactical Aircraft Systems, has been modified for a distributed memory/massively parallel computing environment. The parallel code is operational on an SGI network, Cray J90 and C90 vector machines, SGI Power Challenge, and Cray T3D and IBM SP2 massively parallel machines. Parallel Virtual Machine (PVM) is the message passing protocol for portability to various architectures. A domain decomposition technique was developed which enforces dynamic load balancing to improve solution speed and memory requirements. A host/node algorithm distributes the tasks. The solver parallelizes very well, and scales with the number of processors. Partially parallelized and non-parallelized tasks consume most of the wall clock time in a very fine grain environment. Timing comparisons on a Cray C90 demonstrate that Parallel SPLITFLOW runs 2.4 times faster on 8 processors than its non-parallel counterpart autotasked over 8 processors.

Domel, Neal D.↗

NASA NTRS2000

DRAGON Grid: A Three-Dimensional Hybrid Grid Generation Code Developed

Because grid generation can consume 70 percent of the total analysis time for a typical three-dimensional viscous flow simulation for a practical engineering device, payoffs from research and development could reduce costs and increase throughputs considerably. In this study, researchers at the NASA Glenn Research Center at Lewis Field developed a new hybrid grid approach with the advantages of flexibility, high-quality grids suitable for an accurate resolution of viscous regions, and a low memory requirement. These advantages will, in turn, reduce analysis time and increase accuracy. They result from an innovative combination of structured and unstructured grids to represent the geometry and the computation domain. The present approach makes use of the respective strengths of both the structured and unstructured grid methods, while minimizing their weaknesses. First, the Chimera grid generates high-quality, mostly orthogonal meshes around individual components. This process is flexible and can be done easily. Normally, these individual grids are required overlap each other so that the solution on one grid can communicate with another. However, when this communication is carried out via a nonconservative interpolation procedure, a spurious solution can result. Current research is aimed at entirely eliminating this undesired interpolation by directly replacing arbitrary grid overlapping with a nonstructured grid called a DRAGON grid, which uses the same set of conservation laws over the entire region, thus ensuring conservation everywhere. The DRAGON grid is shown for a typical film-cooled turbine vane with 33 holes and 3 plenum compartments. There are structured grids around each geometrical entity and unstructured grids connecting them. In fiscal year 1999, Glenn researchers developed and tested the three-dimensional DRAGON grid-generation tools. A flow solver suitable for the DRAGON grid has been developed, and a series of validation tests are underway.

Liou, Meng-Sing↗

NASA NTRS2001

Computation of Flow Over a Drag Prediction Workshop Wing/Body Transport Configuration Using CFL3D

A Drag Prediction Workshop was held in conjunction with the 19th AIAA Applied Aerodynamics Conference in June 2001. The purpose of the workshop was to assess the prediction of drag by computational methods for a wing/body configuration (DLR-F4) representative of subsonic transport aircraft. This report details computed results submitted to this workshop using the Reynolds-averaged Navier-Stokes code CFL3D. Two supplied grids were used: a point-matched 1-to-1 multi-block grid, and an overset multi-block grid. The 1-to-1 grid, generally of much poorer quality and with less streamwise resolution than the overset grid, is found to be too coarse to adequately resolve the surface pressures. However, the global forces and moments are nonetheless similar to those computed using the overset grid. The effect of three different turbulence models is assessed using the 1-to-1 grid. Surface pressures are very similar overall, and the drag variation due to turbulence model is 18 drag counts. Most of this drag variation is in the friction component, and is attributed in part to insufficient grid resolution of the 1-to-1 grid. The misnomer of 'fully turbulent' computations is discussed; comparisons are made using different transition locations and their effects on the global forces and moments are quantified. Finally, the effect of two different versions of a widely used one-equation turbulence model is explored.

Rumsey, Christopher L.↗

NASA NTRS1995

Parallel volume ray-casting for unstructured-grid data on distributed-memory architectures

As computing technology continues to advance, computational modeling of scientific and engineering problems produces data of increasing complexity: large in size and unstructured in shape. Volume visualization of such data is a challenging problem. This paper proposes a distributed parallel solution that makes ray-casting volume rendering of unstructured-grid data practical. Both the data and the rendering process are distributed among processors. At each processor, ray-casting of local data is performed independent of the other processors. The global image composing processes, which require inter-processor communication, are overlapped with the local ray-casting processes to achieve maximum parallel efficiency. This algorithm differs from previous ones in four ways: it is completely distributed, less view-dependent, reasonably scalable, and flexible. Without using dynamic load balancing, test results on the Intel Paragon using from two to 128 processors show, on average, about 60% parallel efficiency.

Ma, Kwan-Liu↗

NASA NTRS2010

Application of the FUN3D Unstructured-Grid Navier-Stokes Solver to the 4th AIAA Drag Prediction Workshop Cases

FUN3D Navier-Stokes solutions were computed for the 4th AIAA Drag Prediction Workshop grid convergence study, downwash study, and Reynolds number study on a set of node-based mixed-element grids. All of the baseline tetrahedral grids were generated with the VGRID (developmental) advancing-layer and advancing-front grid generation software package following the gridding guidelines developed for the workshop. With maximum grid sizes exceeding 100 million nodes, the grid convergence study was particularly challenging for the node-based unstructured grid generators and flow solvers. At the time of the workshop, the super-fine grid with 105 million nodes and 600 million elements was the largest grid known to have been generated using VGRID. FUN3D Version 11.0 has a completely new pre- and post-processing paradigm that has been incorporated directly into the solver and functions entirely in a parallel, distributed memory environment. This feature allowed for practical pre-processing and solution times on the largest unstructured-grid size requested for the workshop. For the constant-lift grid convergence case, the convergence of total drag is approximately second-order on the finest three grids. The variation in total drag between the finest two grids is only 2 counts. At the finest grid levels, only small variations in wing and tail pressure distributions are seen with grid refinement. Similarly, a small wing side-of-body separation also shows little variation at the finest grid levels. Overall, the FUN3D results compare well with the structured-grid code CFL3D. The FUN3D downwash study and Reynolds number study results compare well with the range of results shown in the workshop presentations.

Lee-Rausch, Elizabeth M.↗

DOE OSTI2026

End-To-End Decentralized Transmission Line Protection in IBR-Dominated Weak Grids Using Interpretable Data-Driven Methods

Traditional transmission line protection relies on predictable synchronous-based fault signatures, which frequently fail under the non-standard, current-limited fault characteristics of Inverter-Based Resources (IBRs). This study investigates how to achieve secure, communication-free fault isolation in IBR-dominated weak grids without relying on opaque, computationally heavy "black-box" machine learning algorithms. To address this, we propose a novel, standalone, and inherently interpretable data-driven protection framework. Unlike centralized methods requiring multi-terminal communication, this decentralized approach relies solely on local measurements using a hierarchical linear-kernel Support Vector Machine (SVM). The methodology decomposes the protection task into four sequential stages that mimic traditional protection elements: fault detection and fault direction identification, fault type classification, zone classification, and location estimation. This multi-stage architecture allows for specialized feature engineering at each stage, combining high computational efficiency with logic traceability. The framework's end-to-end performance was validated via C-code and PSCAD/EMTDC co-simulation, utilizing a real-world utility network and an OEM black-box IBR model. The proposed relay achieves 97.2% overall accuracy and provides a reliable trip decision within a 2.5-cycle window. The results confirm 100% accuracy in fundamental fault detection, reliable zone selectivity across low to moderate fault resistances, and robust security against non-fault transients, proving its immediate viability for integration into commercial numerical relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Simulations of divertor designs that spatially separate power and particle exhaust using mid-leg divertor particle pumping

Predictive design modeling of a Dissipation-Focused Divertor for future operation in DIII-D reveals that increasing the poloidal distance of the pump duct entrance from the target surface along the low-field side divertor baffle increases neutral compression and modifies the spatial distribution of power dissipation. With a divertor pump located mid-leg between the target and the X-point, SOLPS-ITER boundary plasma simulations without drifts predict the formation of a dense neutral cloud near the target with > 30x higher neutral compression in detachment, a more stable detachment front located further from the target, and ~25% lower outer midplane separatrix density required for detachment onset, compared to a pump located in the scrape-off layer at the target surface. Up to 19 MW of power flowing into the divertors is modeled using the following two numerical implementations for particle pumping: a specified fraction of particles incident on variable wall sections of the plasma grid is removed from the computational domain (so-called albedo pumping), and a pump duct is modeled which includes dynamics of kinetic neutrals in the duct. The simulations show that the detachment front is located between the divertor target and the X-point and is relatively stable near the pump entrance, without a strong dependence on gas puff rate or injected power. The mid-leg pump design spatially separates the two primary functions of a divertor (power handling and particle exhaust), with the majority of power dissipation occurring near the target plate and particle exhaust taking place further upstream. The benefit of enhanced dissipation using mid-leg pumping comes at the cost of a higher outer midplane separatrix density for a given amount of particle injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Emulating climate extreme indices

We use simple pattern scaling and time-shift to emulate changes in a set of climate extreme indices under future scenarios, and evaluate the emulators' accuracy. We propose a metric for the error in emulation in the context of initial condition ensembles, to specifically characterize the role of internal variability in the emulation performance. Our metric separates systematic emulation errors from unavoidable discrepancies between emulated and target values due to internal variability. We compute the metricis at grid-point scale, and we show geographically resolved results, or aggregate them at global scale. We demonstrate the use of our error metric in the emulation of a suite of temperature and precipitation extreme indices. We test and compare simple pattern scaling and time-shift using a range of trajectories spanning targets inspired by the Paris agreement -- warming to 1.5C and 2.0C from the pre-industrial baseline -- and two of the longer-established trajectories, RCP4.5 and RCP8.5. With this suite of scenarios we can test the effects on the performance of the size of the temperature gap between emulation origin and target. We find that for most indices emulation the dominant source of discrepancy is internal variability. For at least one index, however, counting exceedances of a high temperature threshold, significant portions of the globally aggregated discrepancy and its regional pattern originate from the systematic emulation error. This error exceeds internal variability of both the target and the emulated quantities in large coherent regions at low latitudes, and the explanation can be found in the differential behavior of temperature distributions across latitudes. The metric also highlights a fundamental difference in the two methods related to the simulation of internal variability, which is dampened significantly by simple pattern scaling. This aspect is of consequence when using these methods for specific applications, where preserving variability for uncertainty quantification is deemed important. With this study we offer our metric as a diagnostic tool, facilitating the formulation of scientific hypotheses on the reasons for the error. In the meantime, we show that for many impact relevant indices by now traditional emulation techniques can be accurate within the variations unavoidably introduced by internal variability, establishing the fundamental condition for using their emulation in impact modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Department of Energy Project DE-SC0012198: Probing Electron-Molecule Scattering with HHG: Theory and Experiment (Final Report)

The study of the motion of atoms in molecules is important to understanding many areas of physical and life sciences. Such motion occurs on many different times scales, with electronic motion occurring on a sub-femtosecond time scale, and simple vibrational motion in tens to hundreds of femtoseconds. One way to follow such processes in real time is by the use of shortpulsed lasers, and in particular by studying time-resolved photoionization and the related process of high-harmonic generation (HHG). Thus there has been much effort to develop the tools necessary to probe molecular systems using short pulse lasers and understanding the sensitivity of the different possible probes to the time dependent geometric structure as well as the electronic structure of molecules. Our research has particularly focused on the connection between high-field processes and the more traditional weak field photoionization processes. Strong field and weak field processes can be connected through models that involve the same matrix elements. We have demonstrated in our study of HHG from SF 6 that the spectrum is sensitive to the interplay between the angular dependence of the ionization step and recombination step. In our study of rescattering spectroscopy, we have shown that with a combination of experiment and theory, we can use this high-field spectroscopy to determine molecular structure in molecules such as C 2 H 4 . We have also developed new computational tools based on overset grids to enable studies on larger molecular systems which use much more robust numerical approaches so that the resulting code can be a tool that non-specialists can use to study related systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

VPP Cybersecurity: Stakeholder Roles and Responsibilities

This presentation was presented at the Aggregation and Grid Security Workshop - held on June 17-18, 2025, at NREL in Golden, Colorado. The goal of the two-day workshop was to address the critical cybersecurity challenges for the future electric grid. This presentation outlines a framework for assessing the roles and responsibilities of virtual power plant stakeholders in mitigating cyber risk to virtual power plant projects and the overall power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1973

Analytical method for predicting the pressure distribution about a nacelle at transonic speeds

The formulation and development of a computer analysis for the calculation of streamlines and pressure distributions around two-dimensional (planar and axisymmetric) isolated nacelles at transonic speeds are described. The computerized flow field analysis is designed to predict the transonic flow around long and short high-bypass-ratio fan duct nacelles with inlet flows and with exhaust flows having appropriate aerothermodynamic properties. The flow field boundaries are located as far upstream and downstream as necessary to obtain minimum disturbances at the boundary. The far-field lateral flow field boundary is analytically defined to exactly represent free-flight conditions or solid wind tunnel wall effects. The inviscid solution technique is based on a Streamtube Curvature Analysis. The computer program utilizes an automatic grid refinement procedure and solves the flow field equations with a matrix relaxation technique. The boundary layer displacement effects and the onset of turbulent separation are included, based on the compressible turbulent boundary layer solution method of Stratford and Beavers and on the turbulent separation prediction method of Stratford.

Keith, J. S.↗

NASA NTRS1976

PECHCV, PECHFV, PEFHCV and PEFHFV: A set of atmospheric, primitive equation forecast models for the Northern Hemisphere, volume 3

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric prediction models were developed. The models, which use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical, are: (1) PECHCV - five sigma layers and a 63 x 63 horizontal grid, (2) PECHFV - ten sigma layers and a 63 x 63 horizontal grid, (3) PEFHCV - five sigma layers and a 187 x 187 horizontal grid, and (4) PEFHFV - ten sigma layers and a 187 x 187 horizontal grid. The models and associated computer programs are described.

Wellck, R. E.↗

NASA NTRS1977

Implicit finite difference simulation of flow about arbitrary geometries with application to airfoils

Finite difference procedures are used to solve either the Euler equations or the 'thin layer' Navier-Stokes equations subject to arbitrary boundary conditions. An automatic grid generation program is employed, and because an implicit finite difference algorithm for the flow equations is used, time steps are not severely limited when grid points are finely distributed. Computational efficiency and compatibility to vectorized computer processors is maintained by use of approximate factorization techniques. Computed results for both inviscid and viscous flow about airfoils are described and compared to various known solutions.

Steger, J. L.↗

NASA NTRS1980

Mass loss from UW Canis Majoris

The present analysis is an application of the theory described by Lucy (1971) for the calculation of P Cygni resonance line profiles formed by isotropic and coherent scattering in spherically symmetric expanding circumstellar envelopes. Copernicus satellite measurements of resonance features in the FUV spectrum of the O 7 supergiant UW Canis Majoris (= HD 57060) are compared with theoretical P Cygni profiles. Grids of line profiles are computed using four free parameters which contain information about the velocity law, ionization equilibrium, temperature, and envelope density. Thus, with the assumption of a spherically symmetric steady flow, and of solar element abundances, the stellar mass loss rate and the electron temperature of the expanding shell can be derived. The mass flow is treated in a fully transonic way, i.e., the Sobolev approximation is applied.

Drechsel, H.↗