SEARCH · Engineering Papers

Results for “consistent scheme”

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,477 records · Page 82

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Low-Cost, User-Friendly, Rapid Analysis of Dynamic Data System Established

An issue of primary importance to the development of new jet and certain other airbreathing combined-cycle powered aircraft is the advancement of airframe-integrated propulsion technologies. Namely, engine inlets and their systems and subsystems are required to capture, convert, and deliver the atmospheric airflow demanded by such engines across their operating envelope in a form that can be used to provide efficient, stable thrust. This must be done while also minimizing aircraft drag and weight. Revolutionary inlet designs aided by new technologies are needed to enable new missions. An unwanted byproduct of pursuing these inlet technologies is increased time-variant airflow distortion. Such distortions reduce propulsion system stability, performance, operability, and life. To countermand these limitations and fully evaluate the resulting configurations, best practices dictate that this distortion be experimentally measured at large scale and analyzed. The required measurements consist of those made by an array of high-response pressure transducers located in the flow field at the aerodynamic interface plane (AIP) between the inlet and engine. Although the acquisition of the necessary pitot-pressure time histories is relatively straight-forward, until recent years, the analysis has proved to be very time-consuming, tedious, and expensive. To transform the analysis of these data into a tractable and timely proposition, researchers at the NASA Glenn Research Center created and established the Rapid Analysis of Dynamic Data (RADD) system. The system provides complete, near real-time analysis of time-varying inlet airflow distortion datasets with report quality output. This fully digital approach employs Institute of Electrical and Electronics Engineers (IEEE) binary data file format standardization to establish data-acquisition-system-independent processing on low cost personal computers. Features include invalid instrumentation code-out, logging, and multiple replacement schemes as needed for each channel of instrumentation. The AIP pressure distribution can be interpolated to simulate measurements by alternate AIP probe arrays, if desired. In addition, the RADD system provides for the application of filters that can be used to focus the analysis on the frequency range of interest.

Arend, David J.↗

DOE OSTI2025

Interactive Gas Chemistry for Enhanced Science Capabilities of the Energy Exascale Earth System Model Version 3

Atmospheric chemistry plays a crucial role in Earth system models (ESMs), controlling atmospheric composition and radiative balance; it is highly interactive with the physical climate, biogeochemical cycles, and human systems. However, it often imposes computational challenges in an ESM. Here we develop a full troposphere‐stratosphere interactive chemistry module for the US Department of Energy's Energy Exascale Earth System Model (E3SM). We intentionally build a streamlined module based on E3SM version 2 that interacts with other components and maintains all of major chemical and chemistry‐climate feedbacks. The module incorporates a new, highly efficient tracer advection scheme; linearization of stratospheric chemistry; and abridged tropospheric chemical mechanism with 28 reactive tracers. This new model, E3SM‐chem, can readily perform century‐long climate simulations of ozone, methane, and nitrous oxide based on emission scenarios as well as provide hourly budgets for the gas‐phase radicals that drive aerosol chemistry. We evaluate E3SM‐chem with an atmosphere‐only simulation as in the recent climate model intercomparison project (CMIP6) finding results similar to the other CMIP6 models. For the present‐day, E3SM‐chem matches the standard measurement metrics for stratospheric and tropospheric ozone, surface air quality, other key reactive gases like carbon monoxide, and the methane lifetime. Overall, E3SM‐chem maintains the climate fidelity of the baseline model while adding at most 20% to the computational cost of the atmosphere model. Hence, interactive chemistry can be a default configuration for long climate simulations at resolutions of 1° or finer, which is crucial for producing self‐consistent chemistry‐climate feedbacks that alter the climate system.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1997

Current Collection in a Magnetic Field

It is found that the upper-bound limit for current collection in the case of strong magnetic field from the current is close to that given by the Parker-Murphy formula. This conclusion is consistent with the results obtained in laboratory experiments. This limit weakly depends on the shape of the wire. The adiabatic limit in this case will be easily surpassed due to strong magnetic field gradients near the separatrix. The calculations can be done using the kinetic equation in the drift approximation. Analytical results are obtained for the region where the Earth's magnetic field is dominant. The current collection can be calculated (neglecting scattering) using a particle simulation code. Dr. Singh has agreed to collaborate, allowing the use of his particle code. The code can be adapted for the case when the current magnetic field is strong. The needed dm for these modifications is 3-4 months. The analytical description and essential part of the program is prepared for the calculation of the current in the region where the adiabatic description can be used. This was completed with the collaboration of Drs. Khazanov and Liemohn. A scheme of measuring the end body position is also proposed. The scheme was discussed in the laboratory (with Dr. Stone) and it was concluded that it can be proposed for engineering analysis.

Krivorutsky, E. N.↗

NASA NTRS2000

The Cirrus Parcel Model Comparison Project

The cirrus Parcel Model Comparison Project involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 centimeters per second, respectively. Five models are participating in the project. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found amongst the models for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Systematic bias exists between results of a model based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies. This information is necessary, but not sufficient, to construct consistent formulae for the two approaches. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period.

Lin, Ruei-Fong↗

NASA NTRS2006

Wire-Mesh-Based Sorber for Removing Contaminants from Air

A paper discusses an experimental regenerable sorber for removing CO2 and trace components principally, volatile organic compounds, halocarbons, and NH3 from spacecraft cabin air. This regenerable sorber is a prototype of what is intended to be a lightweight alternative to activated-carbon and zeolite-pellet sorbent beds now in use. The regenerable sorber consists mainly of an assembly of commercially available meshes that have been coated with a specially-formulated washcoat containing zeolites. The zeolites act as the sorbents while the meshes support the zeolite-containing washcoat in a configuration that affords highly effective surface area for exposing the sorbents to flowing air. The meshes also define flow paths characterized by short channel lengths to prevent excessive buildup of flow boundary layers. Flow boundary layer resistance is undesired because it can impede mass and heat transfer. The total weight and volume comparison versus the atmosphere revitalization equipment used onboard the International Space Station for CO2 and trace-component removal will depend upon the design details of the final embodiment. However, the integrated mesh-based CO2 and trace-contaminant removal system is expected to provide overall weight and volume savings by eliminating most of the trace-contaminant control equipment presently used in parallel processing schemes traditionally used for spacecraft. The mesh-based sorbent media enables integrating the two processes within a compact package. For the purpose of regeneration, the sorber can be heated by passing electric currents through the metallic meshes combined with exposure to space vacuum. The minimal thermal mass of the meshes offers the potential for reduced regeneration-power requirements and cycle time required for regeneration compared to regenerable sorption processes now in use.

Perry, Jay↗

DOE OSTI2025

Efficient screening of rare large pit anomalies on polished surfaces using a minimalist sampling scheme

Lawrence Livermore National Laboratory (LLNL) has made significant strides in generating clean energy through its inertial confinement fusion (ICF) experiments. These experiments rely on high-density carbon (HDC) coated shells to encapsulate the fusion fuel. The success of these experiments is heavily dependent on the surface quality of these shells, as even minor imperfections, such as deep pits, can negatively impact fusion yield. Ensuring the required smoothness involves an extensive surface-finishing process that spans approximately 20 stages, making it both time-intensive and resource-demanding. A critical challenge in this process is the need for high-resolution scans to detect rare deep pits, which can be costly and impractical if performed on every shell. This highlights the necessity of developing more efficient scanning methods to optimize time and cost without compromising accuracy. To address these challenges, we introduce a novel approach that employs the multivariate Dvoretzky–Kiefer–Wolfowitz (DKW) inequality to provide a probabilistic upper bound on the error in estimating pit distribution characteristics via a Kernel Density Estimator (KDE). This error bound enables efficient and reliable estimation of pit distribution characteristics at a specified statistical confidence level using a minimal number of surface scans. The integrated DKW-KDE approach was validated through surface-finishing experiments across two batches of HDC-coated shells, demonstrating consistent and robust performance across multiple stages of the surface-finishing experiments. The validation studies suggest that the integrated DKW-KDE approach achieves comparable accuracy in estimating the risk of deleterious large pits with six scans, thus conserving time and resources. Further evaluations show that performance remains consistent across batches and over multiple polishing stages. In conclusion, based on these findings, one can leverage the minimal-scan insights to strategically improve the bottleneck inspection process, thus enhancing the productivity and quality of shell polishing and similar challenging manufacturing processes.

Inertial confinement fusion↗

DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗

NASA NTRS1993

An analytical/numerical correlation study of the multiple concentric cylinder model for the thermoplastic response of metal matrix composites

The utility of a recently developed analytical micromechanics model for the response of metal matrix composites under thermal loading is illustrated by comparison with the results generated using the finite-element approach. The model is based on the concentric cylinder assemblage consisting of an arbitrary number of elastic or elastoplastic sublayers with isotropic or orthotropic, temperature-dependent properties. The elastoplastic boundary-value problem of an arbitrarily layered concentric cylinder is solved using the local/global stiffness matrix formulation (originally developed for elastic layered media) and Mendelson's iterative technique of successive elastic solutions. These features of the model facilitate efficient investigation of the effects of various microstructural details, such as functionally graded architectures of interfacial layers, on the evolution of residual stresses during cool down. The available closed-form expressions for the field variables can readily be incorporated into an optimization algorithm in order to efficiently identify optimal configurations of graded interfaces for given applications. Comparison of residual stress distributions after cool down generated using finite-element analysis and the present micromechanics model for four composite systems with substantially different temperature-dependent elastic, plastic, and thermal properties illustrates the efficacy of the developed analytical scheme.

Pindera, Marek-Jerzy↗

NASA NTRS1987

Electromagnetic fields backscattered from an s-shaped inlet cavity with an absorber coating on its inner walls

The EM backscatter from a two-dimensional S-shaped inlet cavity is analyzed using three different techniques, namely a hybrid combination of asymptotic high frequency and modal methods, an integral equation method, and the geometrical optics ray method, respectively. This inlet has a thin absorber coating on its perfectly conducting inner walls and the planar interior termination is made perfectly conducting. The effect of the absorber on the inner wall is treated via a perturbation scheme in the hybrid approach where it is assumed that the loss is sufficiently small for the method to be valid. The results are compared with the backscatter from a straight inlet cavity to evaluate the effect of offsetting the termination in the S-bend configuration such that it is not visible from the open end of the inlet. The envelope of the backscatter pattern for the straight inlet is always seen to peak around the forward axis due to the large return from the directly visible termination, and the pattern envelope tapers off away from the forward axis. Offsetting the termination causes the envelope of the backscatter pattern to flatten out, thereby reducing the return near the forward axis by several dB. The absorber coating reduces the pattern level of the straight inlet in directions away from the forward axis but has little effect on the peak near the axis; furthermore, the absorber coating is seen to consistently reduce the backscatter from the S-bend inlet for almost all incidence angles. The hybrid method gives excellent agreement with experimental data and with the integral equation solution, whereas, the geometrical optics ray tracing method is able to generally predict the average of the bachscatter pattern but not the pattern details.

Burkholder, R. J.↗

NASA NTRS1996

Simulations of Active Vortex Generators

We are interested in the study, via numerical simulations, of active vortex generators. Vortex generators may be used to modify the inner part of the boundary layer or to control separation thus enhancing the performance and maneuverability of aerodynamic configurations. we consider generators that consist of a surface cavity elongated in the streamwise direction and partially covered with a moving lid that at rest lies flush with the boundary. Streamwise vorticity is generated and ejected due to the oscillatory motion of the lid. The present simulations complement relevant experimental investigations of active vortex generators that have been conducted at NASA Ames Research Center and Stanford University used a piezoelectric device in water, allowing for small amplitude high frequency oscillations. They placed the lid asymmetrically on the cavity and observed a strong outward velocity at the small gap of the cavity. Saddoughi used a larger mechanically driven device in air to investigate this flow and observed a jet emerging from the wide gap of the configuration, contrary to the findings of Jacobson and Reynolds We are simulating the flows generated by these devices and we are conducting a parametric study that would help us elucidate the physical mechanisms present in the flow. Conventional computational schemes encounter difficulties when simulating flows around complex configurations undergoing arbitrary motions. Here we present a formulation that achieves this task on a purely Lagrangian frame by extending the formulation presented by Koumoutsakos, Leonard and Pepin (1994). The viscous effects are taken into account by modifying the strength Of the particles, whereas fast multipole schemes employing hundreds of thousands of particles allow for high resolution simulations. We shall present simulation results of an oscillating plate at various Reynolds numbers and Strouhal frequencies.

Mansour, N. N.↗

NASA NTRS2021

Dynamical and Trace Gas Responses of the Quasi-Biennial Oscillation to Increased CO2

Long-term projections of the Quasi-Biennial Oscillation (QBO) remain uncertain. Using the high-top \Middle Atmosphere" (MA) version of the NASA Goddard Institute for Space Studies (GISS) ModelE (E2.2-AP), we show that increased CO2 reduces the QBO period and weakens its amplitude, consistent with results from recent multimodel studies. The amplitude response is asymmetric, with the easterly phase weakening more than the westerly phase. Results from both coupled atmosphere{ocean and climatologically fixed sea surface temperature simulations link the period response to changes in parameterized convective gravity wave fluxes, whereas the amplitude changes are more closely tied to changes in tropical upwelling. Comparisons of simulations with and without ozone-dynamical feedbacks reveal that ozone feedbacks increase the QBO amplitude and period in pre-industrial control simulations, and reduce the response of the amplitude to increased CO2. The response of stratospheric tracers to projected changes in the QBO is also investigated. Similarity between results using fully interactive (trace gas and aerosol) versus linearized ozone chemical mechanisms suggests that the QBO response to CO2 does not depend on additional chemical feedbacks beyond those associated with ozone, which implies that projected QBO changes may be credibly represented using simpli fied stratospheric ozone schemes.

GISS Model E2.2-AP↗

NASA NTRS2019

POLARRIS: A POLArimetric Radar Retrieval and Instrument Simulator

This paper introduces a synthetic polarimetric radar simulator and retrieval package, POLArimetric Radar Retrieval and Instrument Simulator (POLARRIS), for evaluating cloud‐resolving models (CRMs). POLARRIS is composed of forward (POLARRIS‐f) and inverse (retrieval and diagnostic) components (iPOLARRIS) to generate not only polarimetric radar observables (Zh, Zdr, Kdp, ρhv) but also radar‐consistent geophysical parameters such as hydrometeor identification, vertical velocity, and rainfall rates retrieved from CRM data. To demonstrate its application and uncertainties, POLARRIS is applied to simulations of a mesoscale convective system over the Southern Great Plains on 23 May 2011, using the Weather Research and Forecasting model with both spectral bin microphysics (SBM) and the Goddard single‐moment bulk 4ICE microphysics. Statistical composites reveal a significant dependence of simulated polarimetric observables (Zdr, Kdp) on the assumptions of the particle axis ratio (oblateness) and orientation angle distributions. The simulated polarimetric variables differ considerably between the SBM and 4ICE microphysics in part due to the differences in their ice particle size distributions as revealed by comparisons with aircraft measurements. Regardless of these uncertainties, simulated hydrometeor identification distributions overestimate graupel and hail fractions, especially from the simulation with SBM. To minimize uncertainties in forward model, the particle shape and orientation angle distributions of frozen particles should be predicted in a microphysics scheme in addition to the size distributions and particle densities.

Toshi Matsui↗

DOE OSTI2022

From Salt to Electronics: Heteroepitaxy and GaAs Solar Cells

This paper presents work on the heteroepitaxy of salts, specifically fluorides, on semiconductor and heteroepitaxy of semiconductor on salts. Fluorides layers are deposited on commercial Gallium Arsenide (GaAs) wafers followed by the heteroepitaxial growth of GaAs using metal organic chemical vapor deposition (MOCVD). The fluoride layers consist of 2 lattice-engineered layers of alkaline-earth compounds to match with GaAs, and are used to sandwich another alkaline-earth compound with higher water-solubility as sacrificial layer. The triple fluoride layers enable liftoff of free-standing semiconductor films which can be further transferred to desirable substrates. 2D-X-ray Diffraction (2D-XRD) measurements confirm epitaxial growth of both the fluorides and the subsequently grown GaAs films. Single junction (SJ) solar cell devices based on thus prepared films show a power conversion efficiency (PCE) of 10.3% under 1 sun illumination. After the completion of device fabrications, the GaAs film is lifted off from the substrate by a novel water-assisted epitaxial liftoff (H 2 O-ELO) technique and transferred to a cheaper substrate. The original GaAs wafer is recycled and reused twice. Devices based on reused substrates show no significant degradation in performance. Lastly, the semiconductor-salt-semiconductor scheme has great implications in high performance, flexible, and large area electronics.

14 SOLAR ENERGY↗

NASA NTRS2001

An Assessment of the Ability of Potential Space-Borne Instruments to Resolve Spatial and Temporal Variability of Atmospheric Carbon Dioxide

Mounting concern regarding the possibility that increasing carbon dioxide concentrations will initiate climate change has stimulated interest in the feasibility of measuring CO2 mixing ratios from satellites. Currently, the most comprehensive set of atmospheric CO2 data is from the NOAA CMDL cooperative air sampling network, consisting of more than 40 sites where flasks of air are collected approximately weekly. Sporadic observations in the troposphere and stratosphere from airborne in situ and flask samplers are also available. Although the surface network is extensive, there is a dearth of data in the Southern Hemisphere and most of the stations were intentionally placed in remote areas, far from major sources. Sufficiently precise satellite observations with adequate spatial and temporal resolution would substantially increase our knowledge of the atmospheric CO2 distribution and would undoubtedly lead to improved understanding of the global carbon budget. We use a 3-D chemical transport model to investigate the ability of potential satellite instruments with a variety of orbits, horizontal resolution and vertical weighting functions to capture the variation in the modeled CO2 fields. The model is driven by analyzed winds from the Goddard Data Assimilation Office. Simulated CO2 fields are compared with existing surface and aircraft data, and the effects of the model convection scheme and representation of the planetary boundary layer are considered.

Andrews, Arlyn E.↗

NASA NTRS2001

An Assessment of the Ability of Potential Spaceborne Instruments to Resolve Spatial and Temporal Variability of Atmospheric Carbon Dioxide

Andrews, Arlyn E.↗

DOE OSTI2020

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A taxonomy of systems that combine utility-scale renewable energy and energy storage technologies

Interest in renewable energy-based hybrid systems—which combine multiple renewable energy and/or energy storage technologies—is apparent and growing. Such systems could play an important role in the grid of the future, but our ability to quantify their potential impacts is hindered by (a) the inconsistent treatment of hybrid systems and (b) an incomplete understanding of which aspects of hybridization will have the greatest influence. In this article, we review the literature and propose an organization scheme for utility-scale hybrid systems in which electricity is generated from commercially available, renewable energy technologies, and electricity is the only output. Based on this literature review, we propose three categories for delineating systems that link multiple renewable energy and/or energy storage technologies: co-located resources (locational linkage only), virtual power plants (operational linkage only), and full hybrids (location and operational linkages). Proposing a categorization based on whether the systems involve locational and/or operational linkages reflects the fact that technology combinations do not represent a meaningful delineation. In addition, the categories are ideal for informing model development efforts, by identifying the unique characteristics that must be captured in adequate representations of systems that fall under each category. Our categories can also help inform policy considerations, as they delineate system characteristics that would introduce distinct permitting, siting, interconnection, and policy challenges. Therefore, we believe that our proposed taxonomy will help organize thinking and research efforts in a way that will enable consistent evaluation of the costs, values, and potential of renewable energy-based hybrids.

25 ENERGY STORAGE↗

NASA NTRS1998

Multigrid Strategies for Viscous Flow Solvers on Anisotropic Unstructured Meshes

Unstructured multigrid techniques for relieving the stiffness associated with high-Reynolds number viscous flow simulations on extremely stretched grids are investigated. One approach consists of employing a semi-coarsening or directional-coarsening technique, based on the directions of strong coupling within the mesh, in order to construct more optimal coarse grid levels. An alternate approach is developed which employs directional implicit smoothing with regular fully coarsened multigrid levels. The directional implicit smoothing is obtained by constructing implicit lines in the unstructured mesh based on the directions of strong coupling. Both approaches yield large increases in convergence rates over the traditional explicit full-coarsening multigrid algorithm. However, maximum benefits are achieved by combining the two approaches in a coupled manner into a single algorithm. An order of magnitude increase in convergence rate over the traditional explicit full-coarsening algorithm is demonstrated, and convergence rates for high-Reynolds number viscous flows which are independent of the grid aspect ratio are obtained. Further acceleration is provided by incorporating low-Mach-number preconditioning techniques, and a Newton-GMRES strategy which employs the multigrid scheme as a preconditioner. The compounding effects of these various techniques on speed of convergence is documented through several example test cases.

Movriplis, Dimitri J.↗