SEARCH · Engineering Papers

Results for “large grids”

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 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Nonlinear soil-structure interaction calculations simulating the SIMQUAKE experiment using STEALTH 2D

Transient, nonlinear soil-structure interaction simulations of an Electric Power Research Institute, SIMQUAKE experiment were performed using the large strain, time domain STEALTH 2D code and a cyclic, kinematically hardening cap soil model. Results from the STEALTH simulations were compared to identical simulations performed with the TRANAL code and indicate relatively good agreement between all the STEALTH and TRANAL calculations. The differences that are seen can probably be attributed to: (1) large (STEALTH) vs. small (TRANAL) strain formulation and/or (2) grid discretization differences.

Tang, H. T.↗

NASA NTRS1984

Theoretical analysis of aircraft afterbody flow

The strongly interactive flow field about aircraft afterbodies is investigated using computational techniques by which the thin-shear-layer formulation of the compressible, Reynolds-averaged Navier-Stokes equations is solved. A time-dependent implicit numerical algorithm is used to obtain solutions for a variety of afterbody and nozzle geometries, within the class of bodies of revolution, for both subsonic and supersonic external flow, and for sonic and supersonic underexpanded jets. Only centered nozzles with either a sharp lip or a blunt base are considered. In all cases, computed results are compared with experimental data taken at flight Reynolds numbers for like-flow conditions. Turbulence closure is realized using algebraic eddy-viscosity concepts. A new and unique adaptive-grid technique is used to resolve flow regimes with large gradients and to improve the accuracy and efficiency of the computational scheme. Special singular point boundary conditions are used for similar purposes, and are especially effective for highly under-expanded jets. For all cases considered, except one with a very large base-to-nozzle-exit-diameter ratio, the agreement with experimental measurements is excellent. For geometries with large base regions, enhancements in the turbulence transport model are necessary to support improvements in the flow-field simulation.

Deiwert, G. S.↗

NASA NTRS1987

Theoretical analysis of aircraft afterbody flow

Computational methods solving the thin shear layer formulation of the compressible, Reynolds-averaged Navier-Stokes equations are presently used to investigate the strongly interactive flow field about aircraft afterbodies. Solutions for a variety of axisymmetric afterbody and nozzle geometries are solved by means of a time-dependent implicit numerical algorithm for both subsonic and supersonic external flows, and the results obtained are compared with experimental data. A novel adaptive-grid technique is used to resolve flow regimes having large gradients, as well as to improve the accuracy and efficiency of the computational scheme.

Deiwert, George S.↗

NASA NTRS1989

NASA Workshop on Computational Structural Mechanics 1987, part 3

Computational Structural Mechanics (CSM) topics are explored. Algorithms and software for nonlinear structural dynamics, concurrent algorithms for transient finite element analysis, computational methods and software systems for dynamics and control of large space structures, and the use of multi-grid for structural analysis are discussed.

Sykes, Nancy P.↗

NASA NTRS1989

Numerical simulation of the effect of spatial disturbances on vortex asymmetry

The steady asymmetric vortex pattern observed on slender bodies of revolution at large angle of attack was investigated using fine-grid Navier-Stokes computations. The computed results demonstrate the marked asymmetry which has been observed in experiments. To obtain asymmetry, it was found essential to introduce a space-fixed time-invariant perturbation into the computation. The computational results suggest that vortex asymmetry is forced by amplification of small disturbances, such as those due to surface roughness, occurring within the body viscous boundary layer.

Degani, David↗

NASA NTRS1990

Satellite radar altimetry over ice. Volume 4: Users' guide for Antarctica elevation data from Seasat

A gridded surface-elevation data set and a geo-referenced data base for the Seasat radar altimeter data over Greenland are described. This is a user guide to accompany the data provided to data centers and other users. The grid points are on a polar stereographic projection with a nominal spacing of 20 km. The gridded elevations are derived from the elevation data in the geo-referenced data base by a weighted fitting of a surface in the neighborhood of each grid point. The gridded elevations are useful for the creating of large-scale contour maps, and the geo-referenced data base is useful for regridding, creating smaller-scale contour maps, and examinating individual elevation measurements in specific geographic areas. Tape formats are described, and a FORTRAN program for reading the data tape is listed and provided on the tape.

Zwally, H. Jay↗

NASA NTRS1994

A linear method for analyzing lightning field changes

A constrained, least-squares method for analyzing multiple-station measurements of lightning field changes (delta Es) is introduced. Previous methods have attempted to fit the spatial pattern of lightning delta Es using nonlinear models, such as a point charge (Q) or a point dipole (P) model. With the linear method, the delta Es are described not by models but by a general volume charge distribution that is deposited on a large (40 x 40 x 20 cu km) Cartesian grid above the measuring network. A linear system of equations is used to relate the measured delta Es to the charges that are deposited at each grid point. With this approach, the information content of the measurements can be quantified by an eigenanalysis of the covariance matrix of the linear system. Constraints can be used to reduce the infinity of possible solutions to the linear system and also to reduce systematic biases that can be introduced by the method of solution. It is shown that a Landweber iterative method, derived from the general method of steepest descent, can be used to solve the linear system and that the resulting volume charge distributions are generally consistent with computer-simulated charge sources, when these sources are over the measuring network. The Landweber iteration has also provided solutions for natural lightning events that are consistent with Q- and P-model results.

Koshak, William J.↗

NASA NTRS1995

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗

NASA NTRS2014

Parallel Adaptive High-Order CFD Simulations Characterizing Cavity Acoustics for the Complete SOFIA Aircraft

This paper presents one-of-a-kind MPI-parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting of the telescope. A temporally fourth-order Runge-Kutta, and spatially fifth-order WENO-5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32,000 cores and 4 billion cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregularities caused by the highly complex geometry. Limits to scaling beyond 32K cores are identified, and targeted code optimizations are discussed.

Acoustics↗

NASA NTRS2016

Parallel Adjective High-Order CFD Simulations Characterizing SOFIA Cavity Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge-Kutta, and spatially fth-order accurate WENO- 5Z scheme was used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

CFD↗

NASA NTRS2016

Towards a Viscous Wall Model for Immersed Boundary Methods

Immersed boundary methods are frequently employed for simulating flows at low Reynolds numbers or for applications where viscous boundary layer effects can be neglected. The primary shortcoming of Cartesian mesh immersed boundary methods is the inability of efficiently resolving thin turbulent boundary layers in high-Reynolds number flow application. The inefficiency of resolving the thin boundary is associated with the use of constant aspect ratio Cartesian grid cells. Conventional CFD approaches can efficiently resolve the large wall normal gradients by utilizing large aspect ratio cells near the wall. This paper presents different approaches for immersed boundary methods to account for the viscous boundary layer interaction with the flow-field away from the walls. Different wall modeling approaches proposed in previous research studies are addressed and compared to a new integral boundary layer based approach. In contrast to common wall-modeling approaches that usually only utilize local flow information, the integral boundary layer based approach keeps the streamwise history of the boundary layer. This allows the method to remain effective at much larger y+ values than local wall modeling approaches. After a theoretical discussion of the different approaches, the method is applied to increasingly more challenging flow fields including fully attached, separated, and shock-induced separated (laminar and turbulent) flows.

Immersed↗

NASA NTRS2015

Parallel Adaptive High-Order CFD Simulations Characterizing SOFIA Cavitiy Acoustics

This paper presents large-scale MPI-parallel computational uid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-meter infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady ow eld inside and over the cavity interferes with the optical path and mounting structure of the telescope. A tempo- rally fourth-order accurate Runge-Kutta, and a spatially fth-order accurate WENO-5Z scheme were used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh re nement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32k CPU cores and 4 billion compu- tational cells shows excellent scaling. Dynamic load balancing based on execution time on individual AMR blocks addresses irregular numerical cost associated with blocks con- taining boundaries. Limits to scaling beyond 32k cores are identi ed, and targeted code optimizations are discussed.

Parallell↗

NASA NTRS2018

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid computational model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, along with the associated values of net charge, electric field, and electric potential were calculated throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the accelerating grids) to examine the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown. The results from this simulation are compared to existing literature on ion-ion plasma thrusters.

Sankaran, K.↗

NASA NTRS2023

Solar Flare Ribbon Fronts. I. Constraining Flare Energy Deposition with IRIS Spectroscopy

Spectral lines formed at lower atmospheric layers show peculiar profiles at the "leading edge" of ribbons during solar flares. In particular, increased absorption of the BBSO/GST He I λ10830 line, as well as broad and centrally reversed profiles in the spectra of the Mg II and C II lines observed by the IRIS satellite, has been reported. In this work, we aim to understand the physical origin of such peculiar IRIS profiles, which seem to be common of many, if not all, flares. To achieve this, we quantify the spectral properties of the IRIS Mg II profiles at the ribbon leading edge during four large flares and perform a detailed comparison with a grid of radiative hydrodynamic models using the RADYN+FP code. We also studied their transition region (TR) counterparts, finding that these ribbon front locations are regions where TR emission and chromospheric evaporation are considerably weaker compared to other parts of the ribbons. Based on our comparison between the IRIS observations and modeling, our interpretation is that there are different heating regimes at play in the leading edge and the main bright part of the ribbons. More specifically, we suggest that bombardment of the chromosphere by more gradual and modest nonthermal electron energy fluxes can qualitatively explain the IRIS observations at the ribbon leading front, while stronger and more impulsive energy fluxes are required to drive chromospheric evaporation and more intense TR emission in the bright ribbon. Our results provide a possible physical origin for the peculiar behavior of the IRIS chromospheric lines in the ribbon leading edge and new constraints for the flare models.

Solar flares↗

NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulations of Transonic Buffet Over a Supercritical Airfoil

A series of scale-resolving simulations of flow over the ONERA OAT15A airfoil have been performed at an angle of attack of 3.5◦, just past the onset of buffet. The focus of this study is to document the sensitivity of the wall-modeled large eddy simulation (WMLES) methodology for curvilinear structured overset grids within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework to mesh spacing, mesh distribution, and domain size. A secondary purpose of the study is to compare the results from WMLES to unsteady Reynolds-averaged Navier Stokes (URANS) simulations and hybrid RANS-LES (HRLES) within the same LAVA solver framework. The study provides a unique perspective regarding comparisons between different turbulence modeling approaches, time-integration methods, and computational performance since many of the same numerical routines are used for all three types of simulations. The results are compared with experiments and previous numerical studies of the same geometry and flow conditions.

TTT↗

DOE OSTI2024

Utility-Scale Shared Energy Storage Deployment: Challenges, Research Gaps, and Opportunities

Although community energy storage (CES) and behind-the-meter (BTM) energy storage systems have been widely used to offer homeowners and communities a variety of localized benefits, their scalability and grid support functionalities are limited. On the other hand, utility-scale shared energy storage (USES) systems may offer a number of benefits for grid integration, scalability, and economic viability. When compared to BTM and CES alternatives, these large-scale systems provide more storage capacity, more efficient operations, and more economically viable options. The deployment of USES presents opportunities for optimizing grid performance, integrating renewable energy resources, and improving energy security at the community level. However, significant research gaps exist in optimizing the integration and operation of these systems, especially to allocate energy for consumer use, grid services, and enhancing energy resilience. This paper reviews the literature in this regard, focusing on the opportunities, research gaps, and challenges associated with USES deployment. Firstly, the paper provides an overview of USES systems and emphasizes their benefits. Secondly, the key challenges are identified, and research gaps associated with the operation and integration of these systems are highlighted. Lastly, some potential solutions and opportunities that can be adopted to facilitate the rapid deployment and management of USES are presented. Technological, economic, regulatory, and environmental aspects are also discussed in this paper, providing an overview of the current state and future prospects of this technology.

Gautam, Mukesh [BATTELLE (PACIFIC NW LAB)] (ORCID:↗

NASA NTRS1984

Identification of the dynamics of a two-dimensional grid structure using least square lattice filters

It is doubtful whether the dynamics of large space structures (LSS) can be predicted well enough for control system design applications. Hence, dynamic modeling from on-orbit measurements followed by a modification of the control system is of interest, taking into account the utilization of adaptive control concepts. The present paper is concerned with the model determination phase of the adaptive control problem. Using spectral decoupling to determine mode shapes, mode frequency and damping data can be obtained with the aid of an equation error parameter identification method. This method employs a second-order auto-regressive moving average (ARMA) model to represent the natural mode amplitudes. The discussed procedure involves an extension of the application of the least square lattice filter in system identification to a nonintegral, two-dimensional grid structure made of overlapping bars.

Montgomery, R. C.↗

DOE OSTI2025

A novel approach for large-scale wind energy potential assessment

Increasing wind energy generation is central to grid decarbonization, yet methods to estimate wind energy potential are not standardized, leading to inconsistencies and even skewed results. This study aims to improve the fidelity of wind energy potential estimates through an approach that integrates geospatial analysis and machine learning (i.e., Gaussian process regression). We demonstrate this approach to assess the spatial distribution of wind energy capacity potential in the Contiguous United States (CONUS). We find that the capacity-based power density ranges from 1.70 MW/km2 (25th percentile) to 3.88 MW/km2 (75th percentile) for existing wind farms in the CONUS. The value is lower in agricultural areas (2.73 ± 0.02 MW/km2, mean ± 95 % confidence interval) and higher in other land cover types (3.30 ± 0.03 MW/km2). Notably, advancements in turbine manufacturing could reduce power density in areas with lower wind speeds by adopting low specific-power turbines, but improve power density in areas with higher wind speeds (>8.35 m/s at 120m above the ground), highlighting opportunities for repowering existing wind farms. Wind energy potential is shaped by wind resource quality and is regionally characterized by land cover and physical conditions, revealing significant capacity potential in the Great Plains and Upper Texas. The results indicate that areas previously identified as hot spots using existing approaches (e.g., the west of the Rocky Mountains) may have a limited capacity potential due to low wind resource quality. Improvements in methodology and capacity potential estimates in this study could serve as a new basis for future energy systems analysis and planning.

Dai, Tao↗