Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Greening the 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 127 records · Page 7

New large-strain FFT-based formulation and its application to model strain localization in nano-metallic laminates and other strongly anisotropic crystalline materials

This paper presents a new robust large-strain (LS) elasto-viscoplastic (EVP) formulation based on Fast Fourier Transforms (FFTs) for the prediction of the micro-mechanical response and microstructure evolution of polycrystalline and multiphase materials, with emphasis on the effect of strong crystallographic and/or morphologic anisotropy on localization of plastic deformation. In this work, the novel LS-EVPFFT formulation allows treatment of complex initial geometries and large deformations considering three grids of material points: a regular grid in the reference configuration, where FFTs can be performed; an irregular grid in the initial configuration, created by applying a stress-free displacement field to the reference regular grid; and an irregular grid in the current configuration, undergoing large strains and rotations as the material is loaded. Further numerical stability of the new formulation also required the use of a novel expression for the discrete modified Green’s operator, which reduces spurious field oscillations. After presenting and validating the new formulation by comparison with preexisting implementations and analytical solutions, LS-EVPFFT is applied to the prediction of slip and kink bands formation in polycrystalline columnar ice, and kink bands in single crystal zinc wires, showing good agreement with classic experiments. Finally, the model is used to study kink band formation during compression of Cu–Nb nano-metallic laminates (NMLs), in which accurate treatment of the complex geometry associated with the tortuosity of interfaces and large deformations become critical, showing consistency with corresponding micropillar experiments.

36 MATERIALS SCIENCE↗

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗

Parallel Computation of Sensitivity Derivatives with Application to Aerodynamic Optimization of a Wing

This paper focuses on the parallel computation of aerodynamic derivatives via automatic differentiation of the Euler/Navier-Stokes solver CFL3D. The comparison with derivatives obtained by finite differences is presented and the scaling of the time required to obtain the derivatives relative to the number of processors employed for the computation is shown. Finally, the derivative computations are coupled with an optimizer and surface/volume grid deformation tools to perform an optimization to reduce the drag of a three-dimensional wing.

Biedron, Robert T.↗

Assessment of UMUSCL Scheme for DNS of Turbulent Flows

Direct numerical simulation (DNS) is performed using an unstructured node-based finite-volume methodology for canonical turbulent flow configurations to assess the effect of the upwinding parameter (κ) in the unstructured monotonic upstream scheme for conservation laws (UMUSCL) reconstruction scheme. The NASA FUN3D solver is used, and four turbulent validation cases are considered: periodic Taylor-Green vortex, (2) flow through a plane channel, (3) flow through a channel with a constriction, and (4) flow over a flat plate. The results are presented for κ = 0.5, 0.9 and 1 on hexahedral grids, with the last value corresponding to a zero-numerical dissipation scheme. Detailed comparisons with available reference data indicate that for the smallest length-scales simulated, κ = 0.5 is too dissipative while κ = 1 could lead to an energy build-up. Overall, κ =0.9 yields satisfactory results in terms of accuracy and robustness for all the cases considered in this study.

Direct Numerical Simulation↗

Spatial and temporal variations of the fluxes of carbon dioxide and sensible and latent heat over the FIFE site

Airborne measurements of flux densities of carbon dioxide CO2, sensible heat, and latent heat (H2O) obtained over the First ISLSCP Field Experiment (FIFE) site during three intensive field campaigns in 1987 and one in 1989 are examined to characterize the spatial and temporal variations of CO2 and energy transfer processes. These data were collected by the National Research Council Twin Otter using low-level flight patterns, all flown at constant pressure altitude during relatively clear days. The spatial variations are larger in 1989 than in 1987 and a higher correlation is observed between the fluxes and the surface features. The temporal patterns are easier to characterize with the relatively homogeneous situation of 1987. Functional relationships obtained between fluxes of CO2 and latent heat, CO2 fluxes and greenness index, latent heat fluxes and greenness index, and between sensible heat fluxes and surface air temperature differences are presented for one day in 1987 and one in 1989 as an example of the kind of information that can be obtained from grid flights at constant pressure altitude.

Desjardins, R. L.↗

Deconvolution of wide field-of-view radiometer measurements of earth-emitted radiation. II - Analysis of first year of Nimbus 6 ERB data

The theory of deconvolution considered by Smith and Green (1981) is applied to Nimbus 6 data in order to interpret the data with the maximum obtainable resolution. Attention is given to the data sample, sun-contaminated measurements, the measurement model, the deconvolution, the degree variance, the spherical harmonic coefficients, the geographical distribution of longwave radiation, time histories of zonal coefficients, and the effect of a grid system. Degree variance plots for 12 months of longwave radiation data show that the limit for a spherical harmonic representation of the Nimbus 6 wide field-of-view longwave data is 12th degree. The degree variance plots also reveal that most of the power is in the lower degree terms. Contour maps of the radiation field show the geographical distribution of earth emitted radiant exitance for each month.

Bess, T. D.↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

Direct Numerical Simulations↗

The Effects of Atmospheric Models on the Estimation of Infrasonic Source Functions at the Source Physics Experiment

We invert infrasound signals for an equivalent seismoacoustic source function using different atmospheric models to produce the necessary Green’s functions. The infrasound signals were produced by a series of underground chemical explosions as part of the Source Physics Experiment (SPE). In a previous study, we inverted the infrasound data using so-called predictive atmospheric models, which were based on historic, regional-scaled, publicly available weather observations interpolated onto a 3D grid. For the work presented here, we invert the same infrasound data, but using atmospheric models based on weather data collected in a time window that includes the approximate time of the explosion experiments, which we term postdictive models. We build two versions of the postdictive models for each SPE event: one that is based solely on the regional scaled observations, and one that is based on both regional scaled observations combined with on-site observations obtained by a weather sonde released at the time of the SPE. We then invert the observed data set three times, once for each atmospheric model type. We find that the estimated seismoacoustic source functions are relatively similar in waveform shape regardless of which atmospheric model that we used to construct the Green’s functions. However, we find that the amplitude of the estimated source functions is systematically dependent on the atmospheric model type: using the predictive atmospheric models to invert the data generally yields estimated source functions that are larger in amplitude than those estimated using the postdictive models.

54 ENVIRONMENTAL SCIENCES↗

Grid Optimization Competition on Synthetic and Industrial Power Systems

This paper summarizes a grid optimization (GO) competition effort in the United States to find the best solution strategies for up to interconnect-scale power system networks with around 32,000 buses. The optimization problem is a mixedinteger, non-convex non-linear problem, (MINLP) and includes discrete variables such as unit commitment and line switching, control settings (transformer taps and phase shifters with impedance correction tables), and bus shunts. The case study includes six actual industry grids as well as 16 realistic synthetic grids created by three different dataset teams. The winners are selected and ranked based on scoring criteria, which consider the solution quality (such as objective functions) within time limits. Nine winner teams are selected from 26 competitor teams. The results achieved by different teams are described and the performance of different algorithms on synthetic grids and actual industry grids are compared and analyzed.

mixed-integer non-linear programming↗

Modeling and control of nuclear–renewable integrated energy systems: Dynamic system model for green electricity and hydrogen production

The need for decarbonization and diversification of energy resources has led to the development of integrated energy systems (IESs), where multiple resources supply more than one energy sector. Here, one such IES with small modular nuclear reactors and renewables (wind and solar) as generating resources, catering to the demand of the electric grid while producing hydrogen for industries, is modeled in this paper. The physics-based component models are represented using the Modelica language and interconnected to form the IES. The control and coordination of the overall system are ensured by designing a suitable control architecture composed of individual subsystem-level controls and supervisory control. The dynamic performance and the load-following capability of the IES are evaluated, while satisfying the safe operational limits of the components. Different configurations and modes of IES operation are considered, where the adaptability of the control system in the presence of varying demands and renewable generations is validated. The simulation results indicate that hydrogen as a flexible load facilitates the supply of varying grid demand. Additionally, the renewables are also accommodated into the IES owing to the flexibility of the balance of plant associated with the nuclear reactors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Technique for Solving Electrically Small to Large Structures for Broadband Applications

Fast iterative algorithms are often used for solving Method of Moments (MoM) systems, having a large number of unknowns, to determine current distribution and other parameters. The most commonly used fast methods include the fast multipole method (FMM), the precorrected fast Fourier transform (PFFT), and low-rank QR compression methods. These methods reduce the O(N) memory and time requirements to O(N log N) by compressing the dense MoM system so as to exploit the physics of Green s Function interactions. FFT-based techniques for solving such problems are efficient for spacefilling and uniform structures, but their performance substantially degrades for non-uniformly distributed structures due to the inherent need to employ a uniform global grid. FMM or QR techniques are better suited than FFT techniques; however, neither the FMM nor the QR technique can be used at all frequencies. This method has been developed to efficiently solve for a desired parameter of a system or device that can include both electrically large FMM elements, and electrically small QR elements. The system or device is set up as an oct-tree structure that can include regions of both the FMM type and the QR type. The system is enclosed with a cube at a 0- th level, splitting the cube at the 0-th level into eight child cubes. This forms cubes at a 1st level, recursively repeating the splitting process for cubes at successive levels until a desired number of levels is created. For each cube that is thus formed, neighbor lists and interaction lists are maintained. An iterative solver is then used to determine a first matrix vector product for any electrically large elements as well as a second matrix vector product for any electrically small elements that are included in the structure. These matrix vector products for the electrically large and small elements are combined, and a net delta for a combination of the matrix vector products is determined. The iteration continues until a net delta is obtained that is within the predefined limits. The matrix vector products that were last obtained are used to solve for the desired parameter. The solution for the desired parameter is then presented to a user in a tangible form; for example, on a display.

Jandhyala, Vikram↗

Techno-Economic Analysis of a Cryogenic Flux Capacitor for Grid Storage

Abstract The Cryogenic Flux Capacitor (CFC) is a cold, dense fluid storage core with integrated design features that afford the designer flexibility and provide new possibilities for the storage and discharge of energy. The stored energy, in this case, is represented by hydrogen physically bonded within the nanoscale pores within the aerogel composite blanket material, and the process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (∼1,000 m2/g) allows for storage densities close to, or in some cases exceeding, that of normal boiling point liquids. Its performance easily exceeds what can be achieved via ambient temperature, high-pressure gas storage for an equivalent volume. CFC storage is predicted to be easily scalable, is constructed from readily available commercial materials, lends itself to a range of pressure applications, and is geometry insensitive. The high energy density of CFC-stored hydrogen allows long durations of storage, such as daily to monthly cycling, which corresponds to approximately 10 to 100 hours of duration, respectively. The techno-economic analysis examines different types of hydrogen storage for commercial viability. The assessment includes the use of hydrogen blended with natural gas in a combined cycle gas turbine. The study provides results for a range of blends from 100% natural gas to 100% hydrogen. The required fuel usages are estimated for the intended application in a grid scenario that features large amounts of renewable penetration, 10 hours, and as a baseload, 100 hours. The study presents preliminary estimates of the capital cost of storage and operating cost of storage. In addition, the cost of firing natural gas and capturing the carbon dioxide (CO2) in a carbon capture and sequestration (CCS) system is assessed. Based on the study, it is shown that CFC can competitively meet performance and cost needs for grid-based energy storage.

08 HYDROGEN↗

Continuum models for static and dynamic analysis of repetitive lattices

A simple, rational approach is presented for developing continuum models for large repetitive lattice structures subjected to static and dynamic loadings. The procedure involves introducing kinematic assumptions to reduce the dimensionality of the lattice, and equating the strain and kinetic energies of the continuum model to those of the original lattice structure. The proposed procedure is applied to obtain effective elastic and dynamic characteristics of continuum plate models for double-layered grids. Numerical results are presented of stress and free vibration problems for two double-layered grids. These problems demonstrate the high accuracy of the continuum plate models developed. Also, an assessment is made of different approximations to these models.

Noor, A. K.↗

Orbiter Entry Aeroheating Working Group Viscous CFD Boundary Layer Transition Trailblazer Solutions

Boundary layer transition correlations for the Shuttle Orbiter have been previously developed utilizing a two-layer boundary layer prediction technique. The particular two-layer technique that was used is limited to Mach numbers less than 20. To allow assessments at Mach numbers greater than 20, it is proposed to use viscous CFD to the predict boundary layer properties. This report addresses if the existing Orbiter entry aeroheating viscous CFD solutions, which were originally intended to be used for heat transfer rate predictions, adequately resolve boundary layer edge properties and if the existing two-layer results could be leveraged to reduce the number of needed CFD solutions. The boundary layer edge parameters from viscous CFD solutions are extracted along the wind side centerline of the Space Shuttle Orbiter at reentry conditions, and are compared with results from the two-layer boundary layer prediction technique. The differences between the viscous CFD and two-layer prediction techniques vary between Mach 6 and 18 flight conditions and Mach 6 wind tunnel conditions, and there is not a straightforward scaling between the viscous CFD and two-layer values. Therefore: it is not possible to leverage the existing two-layer Orbiter flight boundary layer data set as a substitute for a viscous CFD data set; but viscous CFD solutions at the current grid resolution are sufficient to produce a boundary layer data set suitable for applying edge-based boundary layer transition correlations.

Wood, William A.↗

A Mariner 9 Atlas of the moons of Mars

During its active lifetime Mariner 9 obtained 32 high-resolution pictures of Phobos, and 9 of Deimos. The spacecraft orbit was adequate for about 80% complete coverage of Phobos. However, the coverage of the outer satellite was restricted to the Mars-facing side. The approaches used in processing the pictures are discussed. For each satellite picture a coordinate grid of the satellite was obtained with the aid of a computer.

Veverka, J.↗

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗

A global climatology of albedo, roughness length and stomatal resistance for atmospheric general circulation models as represented by the Simple Biosphere Model (SiB)

Components of the Simple Biosphere Model (SiB) of Sellers et al. (1986) were used to generate global monthly fields of surface albedo (0.4-4.0 microns), roughness length and minimum surface (stomatal) resistance. SiB consists of three submodels which describe the roles of radiative transfer, turbulent transfer and surface resistance in determining the energy balance of the vegetated land surface. These three submodels were detached from SiB and used on the SiB parameter set (total and green leaf area index, leaf angle orientation, canopy dimensions, etc.) to calculate global monthly fields of albedo, roughness length and minimum stomatal resistance at 1 x 1 deg resolution. Time series of various parameters are also displayed for each vegetation type for specified grid points. The SiB results compare reasonably well with appropriate measurements obtained from the literature and have the additional merit of being mutually consistent; the three submodels use many common parameters, which ensures that, for each grid area, the calculated surface properties are closely interrelated as is the case in nature. The derived fields provide a check on the operation of the submodels and the correctness of the parameter set. They can also be used as prescribed fields for GCMs that do not have biophysically based land surface parameterizations.

Dorman, J. L.↗