Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed grid”

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

Hestia-SWIFL: hourly anthropogenic fossil fuel CO2 and heat on the 2km WRF grid, version 1.1

The Hestia-SWIFL version 1.1 anthropogenic heat (AH) and fossil fuel CO2 (FFCO2) emissions data product represent emissions due to the combustion of fossil fuel and cement production within the state of Arizona from 2019 to 2022. This product was developed as part of the Southwest Urban Corridor Integrated Field Laboratory (SW-IFL) project, which aims to provide new knowledge and tools that address extreme heat, air quality, climate change and related urban environmental issues by integrating high-resolution observations, modeling, and civic engagement. The emissions are generated using a bottom-up/engineering approach and are tied to results generated by the Vulcan Project version 4, an effort to quantify space/time-resolved FFCO2 & AH emissions for the entire United States landscape. A large number of data sources are combined to best estimate the emissions at fine scales such as air quality emissions data, traffic flow data, building information, sociodemographic information, and fuel statistics. The AH product provides emissions for two emissions sources (transportation and point source emissions) in units of Watts per hour per square meter (W/m2) per year (annual files) or per hour (hourly files). The FFCO2 product provides emissions from nine individual emission sectors as well as the total, and in units of tons of carbon (tC) per grid cell per year or per hour. The output made available here places the native spatial resolution of the Hestia FFCO2 & AH emissions data product (points, lines, and polygons) into a regularized 2km x 2km grid at hourly and annual temporal resolutions, and stored in netCDF files. The exact spatial extent is defined by the ASU Weather Research Forecast (WRF) simulation grid. All data are processed using R/Python pm high-performance computing system. 2-27-2026 updates: Bugs in airport hourly profile (both AH and FFCO2) and building spatial patterns (FFCO2 only) were fixed. Hourly emissions are reprocessed for all years to reflect those changes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Flight evaluation of Loran-C for general aviation area navigation

This paper reports on a flight evaluation of Loran-C which was part of a long range study of area navigation systems for general aviation. Tests involved two different Loran-C receivers, 6 different aircraft, and a variety of antennas. Uncorrected position fixes were typically accurate to one quarter mile. With measured corrections, repeatability was good to within 200 ft. Signal reliability was 99.7%. The receiver was not sensitive to atmospheric noise. The time difference grid demonstrated a long term stability of 0.3 microsecond. Vertical whip and ADF E-field antennas were found suitable for airborne use. Loran-C was found satisfactory for instrument approaches to runways at general aviation airports where published latitude-longitude coordinates were available. Accuracy was further improved by using locally measured Loran-C time difference coordinates.

Hollister, W. M.↗
NASA NTRS1998

Computational Investigation and Validation of Twin-Tail Buffet Response Including Dynamics and Control

Multidisciplinary tools for prediction of single rectangular-tail buffet are extended to single swept-back-tail buffet in transonic-speed flow, and multidisciplinary tools for prediction and control of twin-tail buffet are developed and presented. The configuration model consists of a sharp-edged delta wing with single or twin tails that are oriented normal to the wing surface. The tails are treated as cantilevered beams fixed at the root and allowed to oscillate in both bending and torsion. This complex multidisciplinary problem is solved sequentially using three sets of equations on a dynamic single or multi-block grid structure. The first set is the unsteady, compressible, Reynolds-averaged Navier-Stokes equations which are used for obtaining the flow field vector and the aerodynamic loads on the tails. The Navier-Stokes equations are solved accurately in time using the implicit, upwind, flux-difference splitting, finite volume scheme. The second set is the coupled bending and torsion aeroelastic equations of cantilevered beams which are used for obtaining the bending and torsion deflections of the tails. The aeroelastic equations'are solved accurately in time using, a fifth-order-accurate Runge-Kutta scheme. The third set is the grid-displacement equations and the rigid-body dynamics equations, which are used for updating the grid coordinates due to the tail deflections and rigid-body motions. The tail-buffet phenomenon is predicted for highly-swept, single vertical tail placed at the plane of geometric symmetry, and for highly-swept, vertical twin tails placed at three different spanwise separation distances. The investigation demonstrates the effects of structural inertial coupling and uncoupling of the bending and torsion modes of vibration, spanwise positions of the twin-tail, angle of attack, and pitching and rolling dynamic motions of the configuration model on the tail buffet loading and response. The fundamental issue of twin-tail buffet alleviation is addressed using two active flow-control methods. These methods are the tangential leading-edge blowing and the flow suction from the leading-edge vortex cores along their paths. Qualitative and quantitative comparisons with the available experimental data are presented. The comparisons indicate that the present multidisciplinary aeroelastic analysis tools are robust, accurate and efficient.

Kandil, Osama A.↗
NASA NTRS1978

Mounting structure

A mounting platform for heat producing instruments operated in a narrow equilibrium temperature range comprises a grid-like structure with relatively large openings therein. The instruments are secured to and thermally coupled with the grid surface facing the instruments. Excess heat from the instruments is selectively radiated to the ambient through openings in the grid, the grid surfaces at these openings exhibiting low thermal emissivity and adsorptivity. The remainder of the grid is maintained at the equilibrium temperature and is covered with a thermal insulating blanket. Thus, the entire system including the platform and instruments is maintained substantially isothermal, whereby the instruments remain in fixed physical relationship to one another.

Ganssle, Eugene Robert↗
NASA NTRS2014

Aeroelastic Analyses of the SemiSpan SuperSonic Transport (S4T) Wind Tunnel Model at Mach 0.95

Detailed aeroelastic analyses of the SemiSpan SuperSonic Transport (S4T) wind tunnel model at Mach 0.95 with a 1.75deg fixed angle of attack are presented. First, a numerical procedure using the Computational Fluids Laboratory 3-Dimensional (CFL3D) Version 6.4 flow solver is investigated. The mesh update method for structured multi-block grids was successfully applied to the Navier-Stokes simulations. Second, the steady aerodynamic analyses with a rigid structure of the S4T wind tunnel model are reviewed in transonic flow. Third, the static analyses were performed for both the Euler and Navier-Stokes equations. Both the Euler and Navier-Stokes equations predicted a significant increase of lift forces, compared to the results from the rigid structure of the S4T wind-tunnel model, over various dynamic pressures. Finally, dynamic aeroelastic analyses were performed to investigate the flutter condition of the S4T wind tunnel model at the transonic Mach number. The condition of flutter was observed at a dynamic pressure of approximately 75.0-psf for the Navier-Stokes simulations. However, it was observed that the flutter condition occurred a dynamic pressure of approximately 47.27-psf for the Euler simulations. Also, the computational efficiency of the aeroelastic analyses for the S4T wind tunnel model has been assessed.

Hur, Jiyoung↗
DOE OSTI2019

Hybrid approach for estimating dynamic states of synchronous generators

Accurate and robust estimation of dynamic states is essential for monitoring and controlling a power grid. To achieve higher accuracy and robustness, this paper proposes a hybrid estimation approach to estimate the dynamic states of synchronous generators using multiple models, classified into a set of local models and a wide-area model. The proposed estimation approach, first, estimates dynamic states using the extended Kalman filter (EKF) based on a local model and a wide-area model. Then, at each time step, probability indexes, which quantify the likelihood of both local and wide-area models, are determined using hypothesis testing based on measurements innovations. Finally, using the weighted fix approach, the states estimated from the local and wide-area models are combined based on their probability indexes. Simulation studies using the two-area four-machine system shows that the proposed approach can improve estimation accuracy and increase robustness against model errors.

Akhlaghi, Shahrokh↗
DOE OSTI2024

Updated sustainability status of cadmium telluride thin‐film photovoltaic systems and projections

This paper provides a comprehensive assessment of the up-to-date life-cycle sustainability status of cadmium-telluride based photovoltaic (PV) systems. Current production modules (Series 6 and Series 7) are analyzed in terms of their energy performance and environmental footprint and compared with the older series 4 module production and current single-crystalline Silicon (sc-Si) module production. For fixed-tilt systems with Series 6 modules operating under average US irradiation of 1800 kWh/m 2 /year, the global warming potential (GWP) is reduced from 16 g CO 2eq /kWh in Series 4 systems to 10 CO 2 eq /kWh in Series 6 systems. For operation in US-SW irradiation of 2300 kWh/m 2 /year, the GWP is reduced from 11 to 8 CO 2eq /kWh and for 1-axis tracking systems operating in Phoenix, Arizona, with point-of array irradiation of 3051 kWh/m 2 /year the GWP is reduced to 6.5 CO 2eq /kWh. Similar reductions have happened in all environmental indicators. Energy payback times (EPBT) of currently installed systems range from 0.6 years for fixed–tilt ground–mounted installations at average US irradiation at latitude tilt installations to 0.3 years for one-axis trackers at high US-SW irradiation, considering average fossil-fuel dominated electricity grids with fuel to electricity conversion efficiency of 0.3. The resulting energy return on energy investment (EROI) also depends on the conversion efficiency of the electricity grid and on the operation life expectance. For a 30-year operational life and grid conversion efficiency of 0.3, EROI ranges from 50 (at US average irradiation) to 70 for US-SW irradiation. The EROI declines with increased grid conversion efficiency; for CdTe PV operating in south California with grid conversion efficiency of 49%, the EROI is about 50 and is projected to fall to 30 when the state's 2030 target of 80% renewable energy penetration materializes. Material alternatives that show a potential of further reductions in degradation rates and materials for enhanced encapsulation that would enable longer operation lives have also been investigated. A degradation rate of 0.3%/year, which has been verified by accelerated testing, is assumed in 30-year scenarios; this is projected to be reduced to 0.2%/year in the near-term and potentially to 0.1%/year in the longer term. With such low degradation rates and enhanced edge-sealing, modules can last 40- to 50-years. Consequently, all impact indicators will be proportionally reduced while EROI will increase. This detailed LCA was conducted according to ISO standards and IEA PVPS Task 12 guidelines. Furthermore, the study revealed that the choices of system models, methods and temporal system boundaries can significantly impact the results and points out to the need to include assumptions regarding these choices in the “transparency in reporting” requirements listed in the IEA PVPS Task 12 Guidelines.

14 SOLAR ENERGY↗
DOE OSTI2017

Gravitational wave production from preheating: parameter dependence

Parametric resonance is among the most efficient phenomena generating gravitational waves (GWs) in the early Universe. The dynamics of parametric resonance, and hence of the GWs, depend exclusively on the resonance parameter q . The latter is determined by the properties of each scenario: the initial amplitude and potential curvature of the oscillating field, and its coupling to other species. Previous works have only studied the GW production for fixed value(s) of q . We present an analytical derivation of the GW amplitude dependence on q , valid for any scenario, which we confront against numerical results. By running lattice simulations in an expanding grid, we study for a wide range of q values, the production of GWs in post-inflationary preheating scenarios driven by parametric resonance. We present simple fits for the final amplitude and position of the local maxima in the GW spectrum. Our parametrization allows to predict the location and amplitude of the GW background today, for an arbitrary q . The GW signal can be rather large, as h {sup 2Ω}{sub GW}( f {sub p} ) ∼< 10{sup −11}, but it is always peaked at high frequencies f {sub p} ∼> 10{sup 7} Hz. We also discuss the case of spectator-field scenarios, where the oscillatory field can be e.g. a curvaton, or the Standard Model Higgs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Integrated Model, Batch and Domain Parallelism in Training Neural Networks

We propose a new integrated method of exploiting model, batch and domain parallelism for the training of deep neural networks (DNNs) on large distributed-memory computers using minibatch stochastic gradient descent (SGD). Our goal is to find an efficient parallelization strategy for a fixed batch size using $P$ processes. Our method is inspired by the communication-avoiding algorithms in numerical linear algebra. We see $P$ processes as logically divided into a $P_r \times P_c$ grid where the $P_r$ dimension is implicitly responsible for model/domain parallelism and the $P_c$ dimension is implicitly responsible for batch parallelism. In practice, the integrated matrix-based parallel algorithm encapsulates these types of parallelism automatically. We analyze the communication complexity and analytically demonstrate that the lowest communication costs are often achieved neither with pure model nor with pure data parallelism. We also show how the domain parallel approach can help in extending the theoretical scaling limit of the typical batch parallel method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Design optimization of transonic airfoils

Numerical optimization procedures were considered for the design of airfoils in transonic flow based on the transonic small disturbance (TSD) and Euler equations. A sequential approximation optimization technique was implemented with an accurate approximation of the wave drag based on the Nixon's coordinate straining approach. A modification of the Euler surface boundary conditions was implemented in order to efficiently compute design sensitivities without remeshing the grid. Two effective design procedures producing converged designs in approximately 10 global iterations were developed: interchanging the role of the objective function and constraint and the direct lift maximization with move limits which were fixed absolute values of the design variables.

Joh, C.-Y.↗
DOE OSTI2010

Grid-based methods for diatomic quantum scattering problems. III. Double photoionization of molecular hydrogen in prolate spheroidal coordinates

Our previously developed finite-element discrete-variable representation in prolate spheroidal coordinates is extended to two-electron systems with a study of double ionization of H{sub 2} with fixed nuclei. Particular attention is paid to the development of fast and accurate methods for treating the electron-electron interaction. The use of exterior complex scaling in the implementation offers a simple way of enforcing Coulomb boundary conditions for the electronic double continuum. While the angular distributions calculated in this study are found to be completely consistent with our earlier treatments that employed single-center expansions in spherical coordinates, we find that the magnitude of the integrated cross sections are sensitive to small changes in the initial-state wave function. The present formulation offers significant advantages with respect to convergence and efficiency and opens the way to calculations on more complicated diatomic targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
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.↗
DOE OSTI2018

ICON-A, the Atmosphere Component of the ICON Earth System Model: I. Model Description

ICON–A is the new icosahedral nonhydrostatic (ICON) atmospheric general circulation model in a configuration using the Max Planck Institute physics package, which originates from the ECHAM6 general circulation model, and has been adapted to account for the changed dynamical core framework. The coupling scheme between dynamics and physics employs a sequential updating by dynamics and physics, and a fixed sequence of the physical processes similar to ECHAM6. To allow a meaningful initial comparison between ICON–A and the established ECHAM6–LR model, a setup with similar, low resolution in terms of number of grid points and levels is chosen. The ICON–A model is tuned on the base of the Atmospheric Model Intercomparison Project (AMIP) experiment aiming primarily at a well balanced top–of atmosphere energy budget to make the model suitable for coupled climate and Earth system modeling. The tuning addresses first the moisture and cloud distribution to achieve the top–of–atmosphere energy balance, followed by the tuning of the parameterized dynamic drag aiming at reduced wind errors in the troposphere. The resulting version of ICON–A has overall biases, which are comparable to those of ECHAM6. Problematic specific biases remain in the vertical distribution of clouds and in the stratospheric circulation, where the winter vortices are too weak. Biases in precipitable water and tropospheric temperature are, however, reduced compared to the ECHAM6. As a result, ICON–A will serve as the basis of further development and as the atmosphere component to the coupled model, ICON–Earth system model (ESM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Comparison of Linear Solvers for Resolving Flow in Three-Dimensional Discrete Fracture Networks

We compare various methods for resolving steady flow within three-dimensional discrete fracture networks, including direct methods, Krylov subspace methods with and without preconditioning, and multi-grid methods. We compared the performance of the methods based on compute times and scaling of the solution as a function of the number of grid nodes and log-variance of the hydraulic aperture. The methods are applied to three test cases: (a) variable density of networks with a truncated power-law distribution of fracture lengths, (b) a fixed network composed of monodisperse fracture sizes but varied permeability/aperture heterogeneity, (c) and a network based on field site in Nevada, US. We chose these cases to allow us to study the impact of the mesh size and flow properties, as well as to demonstrate our conclusions on a large-scale, realistic problem (more than 40 million mesh nodes). A direct solution using Cholesky factorization outperformed other methods for every example but was closely followed in performance by some algebraic multigrid (AMG) preconditioned Krylov subspace methods. Among the Krylov methods, conjugate gradients (CG) with an AMG preconditioner performs the best. Generally, Cholesky factorization is recommended, but CG with an AMG preconditioner may be suitable for very large problems beyond 40 million nodes where the entire linear system cannot reside in memory.

58 GEOSCIENCES↗
NASA NTRS1991

Reynolds number effects on supersonic asymmetrical flows over a cone at high angle of attack

The supersonic viscous flow over a 5-degree half-angle cone at an angle of attack of four times the cone half-angle is studied computationally using both the conical and the three-dimensional Navier-Stokes equations. The numerical solutions were obtained with an implicit, upwind-biased algorithm. Asymmetrical flowfields of the absolute-instability type are found using the conical-flow equations which agree with published results. However, the absolute instabilities of the originally symmetric flow found with the conical equations do not occur in the three-dimensional simulations, although spurious asymmetric three-dimensional flows for symmetric bodies arise if the grid resolution is insufficient in the nose region. The asymmetric flows computed with the three-dimensional equations are convective instabilities and are possible if the local Reynolds number exceeds a critical value and a fixed geometric asymmetry is imposed. A continuous range of asymmetries can be developed, depending on the size of the disturbance and the Reynolds number. As the Reynolds number is increased, the asymmetries demonstrate a bistable behavior at levels of side force consistent with those predicted using the conical equations. Below a certain critical Reynolds number, any flow asymmetries arising from a geometrical asymmetry are damped with increasing distance downstream from the geometrical asymmetry.

Thomas, J. L.↗
NASA NTRS1992

Analysis of an advanced ducted propeller subsonic inlet

A time marching Navier-Stokes code called PARC (PARC2D for 2-D/axisymmetric and PARC3D for 3-D flow simulations) was validated for an advanced ducted propeller (ADP) subsonic inlet. The code validation for an advanced ducted propeller (ADP) subsonic inlet. The code validation was implemented for a non-separated flow condition associated with the inlet operating at angles-of-attack of 0 and 25 degrees. The inlet test data were obtained in the 9 x 15 ft Low Speed Wind Tunnel at NASA Lewis Research Center as part of a cooperative study with Pratt and Whitney. The experimental study focused on the ADP inlet performance for take-off and approach conditions. The inlet was tested at a free stream Mach number of 0.2, at angles-of-attack between O and 35 degrees, and at a maximum propeller speed of 12,000 RPM which induced a corrected air flow rate of about 46 lb/sec based on standard day conditions. The computational grid and flow boundary conditions (BC) were based on the actual inlet geometry and the funnel flow conditions. At the propeller face, two types of BC's were applied: a mass flow BC and a fixed flow properties BC. The fixed flow properties BC was based on a combination of data obtained from the experiment and calculations using a potential flow code. Comparison of the computational results with the test data indicates that the PARC code with the propeller face fixed flow properties BC provided a better prediction of the inlet surface static pressures than the predictions when the mass flow BC was used. For an angle-of-attack of 0 degrees, the PARC2D code with the propeller face mass flow BC provided a good prediction of inlet static pressures except in the region of high pressure gradient. With the propeller face fixed flow properties BC, the PARC2D code provided a good prediction of the inlet static pressures. For an angle-of-attack of 25 degrees with the mass flow BC, the PARC3D code predicted statis pressures which deviated significantly from the test data; however, with the fixed flow properties BC, a good comparison with the test data was obtained.

Iek, Chanthy↗
DOE OSTI2013

Hurricanes in an Aquaplanet World: Implications of the Impacts of External Forcing and Model Horizontal Resolution

High-resolution climate models have been shown to improve the statistics of tropical storms and hurricanes compared to low-resolution models. The impact of increasing horizontal resolution in the tropical storm simulation is investigated exclusively using a series of Atmospheric Global Climate Model (AGCM) runs with idealized aquaplanet steady-state boundary conditions and a fixed operational storm-tracking algorithm. The results show that increasing horizontal resolution helps to detect more hurricanes, simulate stronger extreme rainfall, and emulate better storm structures in the models. However, increasing model resolution does not necessarily produce stronger hurricanes in terms of maximum wind speed, minimum sea level pressure, and mean precipitation, as the increased number of storms simulated by high-resolution models is mainly associated with weaker storms. The spatial scale at which the analyses are conducted appears to have more important control on these meteorological statistics compared to horizontal resolution of the model grid. When the simulations are analyzed on common low-resolution grids, the statistics of the hurricanes, particularly the hurricane counts, show reduced sensitivity to the horizontal grid resolution and signs of scale invariant.

aquaplanet↗
DOE OSTI2014

Motion as perturbation. II. Development of the method for dosimetric analysis of motion effects with fixed-gantry IMRT

Purpose: In this work, the feasibility of implementing a motion-perturbation approach to accurately estimate volumetric dose in the presence of organ motion—previously demonstrated for VMAT-–is studied for static gantry IMRT. The method's accuracy is improved for the voxels that have very low planned dose but acquire appreciable dose due to motion. The study describes the modified algorithm and its experimental validation and provides an example of a clinical application. Methods: A contoured region-of-interest is propagated according to the predefined motion kernel throughout time-resolved 4D phantom dose grids. This timed series of 3D dose grids is produced by the measurement-guided dose reconstruction algorithm, based on an irradiation of a staticARCCHECK (AC) helical dosimeter array (Sun Nuclear Corp., Melbourne, FL). Each moving voxel collects dose over the dynamic simulation. The difference in dose-to-moving voxel vs dose-to-static voxel in-phantom forms the basis of a motion perturbation correction that is applied to the corresponding voxel in the patient dataset. A new method to synchronize the accelerator and dosimeter clocks, applicable to fixed-gantry IMRT, was developed. Refinements to the algorithm account for the excursion of low dose voxels into high dose regions, causing appreciable dose increase due to motion (LDVE correction). For experimental validation, four plans using TG-119 structure sets and objectives were produced using segmented IMRT direct machine parameters optimization in Pinnacle treatment planning system (v. 9.6, Philips Radiation Oncology Systems, Fitchburg, WI). All beams were delivered with the gantry angle of 0°. Each beam was delivered three times: (1) to the static AC centered on the room lasers; (2) to a static phantom containing a MAPCHECK2 (MC2) planar diode array dosimeter (Sun Nuclear); and (3) to the moving MC2 phantom. The motion trajectory was an ellipse in the IEC XY plane, with 3 and 1.5 cm axes. The period was 5 s, with the resulting average motion speed of 1.45 cm/s. The motion-perturbed high resolution (2 mm voxel) volumetric dose grids on the MC2 phantom were generated for each beam. From each grid, a coronal dose plane at the detector level was extracted and compared to the corresponding moving MC2 measurement, using gamma analysis with both global (G) and local (L) dose-error normalization. Results: Using the TG-119 criteria of (3%G/3 mm), per beam average gamma analysis passing rates exceeded 95% in all cases. No individual beam had a passing rate below 91%. LDVE correction eliminated systematic disagreement patterns at the beams’ aperture edges. In a representative example, application of LDVE correction improved (2%L/2 mm) gamma analysis passing rate for an IMRT beam from 74% to 98%. Conclusions: The effect of motion on the moving region-of-interest IMRT dose can be estimated with a standard, static phantom QA measurement, provided the motion characteristics are independently known from 4D CT or otherwise. The motion-perturbed absolute dose estimates were validated by the direct planar diode array measurements, and were found to reliably agree with them in a homogeneous phantom.

62 RADIOLOGY AND NUCLEAR MEDICINE↗