Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed model predictive control”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Forest reorganisation effects on fuel moisture content can exceed changes due to climate warming in wet temperate forests

Abstract The distributions of vegetation and fire activity are changing rapidly in response to climate warming. In many regions, climate effects on dead fuel moisture content (FMC) are expected to increase future wildfire activity. However, forest FMC is largely driven by microclimate conditions, which are moderated from open weather by vegetation canopies. As shifts in vegetation increase under climate warming, the extent to which future fire activity will be driven by climate directly or associated vegetation shifts remains unresolved. Here, we present a study aimed at quantifying the relative magnitudes of (i) direct climate warming, and (ii) vegetation change, on FMC. Field sites to evaluate these effects were established in a natural laboratory of altered forest states to mature wet temperate forest in south‐eastern Australia. FMC was estimated using a process‐based model and 48 years of reconstructed climate data. Canopy effects on microclimate were captured by transferring inputs from climate to microclimate using models parameterised with field observations. To evaluate the relative magnitude of climate and vegetation effects, we calculated the maximum difference in mean annual FMC across annual climate replicates and compared this to FMC differences across reorganising forest sites. Our results show vegetation effects on FMC can exceed those related to expected climate change. Changes to forest structure and composition increased (+15.7%) and decreased (−12.3%) mean annual FMC, with a larger negative effect when forest cover was completely removed (−18.5%). In contrast, the largest climate effect on FMC was −6.6% across 48‐years of data. Our study demonstrates that the magnitude of vegetation effects on FMC can exceed expected climate change effects. Models of future fire activity that do not account for changing vegetation effects on microclimate are omitting a key biophysical control on FMC and therefore may not be accurately predicting future fire activity.

Brown, Tegan P.↗
NASA NTRS2002

Study of Fluid Flow Control in Protein Crystallization using Strong Magnetic Fields

An important component in biotechnology, particularly in the area of protein engineering and rational drug design is the knowledge of the precise three-dimensional molecular structure of proteins. The quality of structural information obtained from X-ray diffraction methods is directly dependent on the degree of perfection of the protein crystals. As a consequence, the growth of high quality macromolecular crystals for diffraction analyses has been the central focus for biochemists, biologists, and bioengineers. Macromolecular crystals are obtained from solutions that contain the crystallizing species in equilibrium with higher aggregates, ions, precipitants, other possible phases of the protein, foreign particles, the walls of the container, and a likely host of other impurities. By changing transport modes in general, i.e., reduction of convection and sedimentation, as is achieved in "microgravity", researchers have been able to dramatically affect the movement and distribution of macromolecules in the fluid, and thus their transport, formation of crystal nuclei, and adsorption to the crystal surface. While a limited number of high quality crystals from space flights have been obtained, as the recent National Research Council (NRC) review of the NASA microgravity crystallization program pointed out, the scientific approach and research in crystallization of proteins has been mainly empirical yielding inconclusive results. We postulate that we can reduce convection in ground-based experiments and we can understand the different aspects of convection control through the use of strong magnetic fields and field gradients. Whether this limited convection in a magnetic field will provide the environment for the growth of high quality crystals is still a matter of conjecture that our research will address. The approach exploits the variation of fluid magnetic susceptibility with concentration for this purpose and the convective damping is realized by appropriately positioning the crystal growth cell so that the magnetic susceptibility force counteracts terrestrial gravity. The general objective is to test the hypothesis of convective control using a strong magnetic field and magnetic field gradient and to understand the nature of the various forces that come into play. Specifically we aim to delineate causative factors and to quantify them through experiments, analysis and numerical modeling. Once the basic understanding is obtained, the study will focus on testing the hypothesis on proteins of pyruvate dehydrogenase complex (PDC), proteins E1 and E3. Obtaining high crystal quality of these proteins is of great importance to structural biologists since their structures need to be determined. Specific goals for the investigation are: 1. To develop an understanding of convection control in diamagnetic fluids with concentration gradients through experimentation and numerical modeling. Specifically solutal buoyancy driven convection due to crystal growth will be considered. 2. To develop predictive measures for successful crystallization in a magnetic field using analyses and numerical modeling for use in future protein crystal growth experiments. This will establish criteria that can be used to estimate the efficacy of magnetic field flow damping on crystallization of candidate proteins. 3. To demonstrate the understanding of convection damping by high magnetic fields to a class of proteins that is of interest and whose structure is as yet not determined. 4. To compare quantitatively, the quality of the grown crystals with and without a magnetic field. X-ray diffraction techniques will be used for the comparative studies. In a preliminary set of experiments, we studied crystal dissolution effects in a 5 Tesla magnet available at NASA Marshall Space Flight Center (MSFC). Using a Schlieren setup, a 1mm crystal of Alum (Aluminum-Potassium Sulfate) was introduced in a 75% saturated solution and the resulting dissolution plume was observed. The experiment was conducted both in the presence and absence of a magnetic field gradient. The magnet produces a gradient field of approx. 1 Tesla2/cm. Image analysis of the recorded images indicated an enhanced plume velocity that was of the order of the measurement limit. For this experiment, both the gradient and gravity fields are in the same direction resulting in an enhanced effective gravity that tends to accelerate the observed plume velocity. While the results are not conclusive, pending further tests, it clearly points out the inadequacy of the MSFC magnet for conducting protein crystallization experiments and the need for a stronger magnet. In spacebased experiments, however, where the gravitational effects are small, only a weak magnetic field will be required to control or mitigate the effects of convective contamination.

Ramachandran, Narayanan↗
NASA NTRSDate not supplied

Implementing Atmospheric Infrared Sounder (AIRS) and Cross-Track Infrared Sounder (CrIS) Cloud-Clearing Algorithm into the NASA GEOS: Focus on the 2017 Atlantic Tropical Cyclone Season

Numerical Weather Prediction (NWP) centers assimilate cloud-free infrared (IR) radiances because the assimilation of all-sky IR radiances is not yet operationally achievable. The cloud-clearing procedure offers a simpler, but effective strategy that produces cloud-affected radiances suitable for assimilation in partially cloudy regions. Several studies conducted by this team have demonstrated that IR Cloud-Cleared Radiances (CCRs), if thinned more aggressively than clear-sky radiances, can improve analysis and forecasts, particularly in meteorologically active areas. However, CCRs are not used by operational centers due partly to the thought that the process of cloud-clearing may affect latency and introduce difficult-to-control external dependencies. This study presents the results of implementing an Atmospheric Infrared Sounder (AIRS) and Cross-Track Infrared Sounder (CrIS) cloud-clearing procedure into the NASA Goddard Earth Observing System (GEOS) to demonstrate the portability of the procedure. The AIRS and CrIS cloud-clearing algorithms have been deprived of external dependencies, made customizable to any specific model, and the computational efficiency has been improved via parallelization. The revised AIRS and CrIS cloud-clearing algorithms allow a customized choice of channel selection, the use of a user-specified model's fields as first guess, and can perform in real time. Data assimilation experiments with the hybrid 4DEnVar GEOS system were successfully performed for the 2017 tropical cyclones (TC) season with a focus on three major hurricanes (Harvey, Irma, and Maria). This study shows that assimilation of locally-generated CCRs have a positive impact on both global skill and TC representation, compared to the assimilation of AIRS and CrIS clear-sky radiances, and a comparable or slightly improved impact compared to assimilation of CCRs produced by external sources, such as NASA's Distributed Active Archive Centers and NOAA’s Comprehensive Large Array-data Stewardship System. The customization and computational efficiency of the revised procedure would enable its usability in a real-time forecast context.

Niama Boukachaba↗
NASA NTRS2000

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. 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 cm/s, respectively. Five models participated. 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 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. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results 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, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) 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, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗
NASA NTRSDate not supplied

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗
NASA NTRS2015

Aircraft Measurements of BrO, IO, Glyoxal, NO2, H2O, O2-O2 and Aerosol Extinction Profiles in the Tropics: Comparison with Aircraft-/Ship-Based in Situ and Lidar Measurements

Tropospheric chemistry of halogens and organic carbon over tropical oceans modifies ozone and atmospheric aerosols, yet atmospheric models remain largely untested for lack of vertically resolved measurements of bromine monoxide (BrO), iodine monoxide (IO) and small oxygenated hydrocarbons like glyoxal (CHOCHO) in the tropical troposphere. BrO, IO, glyoxal, nitrogen dioxide (NO2), water vapor (H2O) and O2-O2 collision complexes (O4/ were measured by the University of Colorado Airborne Multi-AXis Differential Optical Absorption Spectroscopy (CU AMAXDOAS) instrument, aerosol extinction by high spectral resolution lidar (HSRL), in situ aerosol size distributions by an ultra high sensitivity aerosol spectrometer (UHSAS) and in situ H2O by vertical-cavity surface-emitting laser (VCSEL) hygrometer. Data are presented from two research flights (RF12, RF17) aboard the National Science Foundation/ National Center for Atmospheric Research Gulfstream V aircraft over the tropical Eastern Pacific Ocean (tEPO) as part of the "Tropical Ocean tRoposphere Exchange of Reactive halogens and Oxygenated hydrocarbons" (TORERO) project (January/February 2012). We assess the accuracy of O4 slant column density (SCD) measurements in the presence and absence of aerosols. Our O4-inferred aerosol extinction profiles at 477 nm agree within 6% with HSRL in the boundary layer and closely resemble the renormalized profile shape of Mie calculations constrained by UHSAS at low (sub-Rayleigh) aerosol extinction in the free troposphere. CU AMAX-DOAS provides a flexible choice of geometry, which we exploit to minimize the SCD in the reference spectrum (SCDREF, maximize signal-to-noise ratio) and to test the robustness of BrO, IO and glyoxal differential SCDs. The RF12 case study was conducted in pristine marine and free tropospheric air. The RF17 case study was conducted above the NOAA RV Ka'imimoana (TORERO cruise, KA-12-01) and provides independent validation data from ship-based in situ cavity-enhanced DOAS and MAX-DOAS. Inside the marine boundary layer (MBL) no BrO was detected (smaller than 0.5 pptv), and 0.2-0.55 pptv IO and 32-36 pptv glyoxal were observed. The near-surface concentrations agree within 30% (IO) and 10% (glyoxal) between ship and air-craft. The BrO concentration strongly increased with altitude to 3.0 pptv at 14.5 km (RF12, 9.1 to 8.6 deg N; 101.2 to 97.4 deg W). At 14.5 km, 5-10 pptv NO2 agree with model predictions and demonstrate good control over separating tropospheric from stratospheric absorbers (NO2 and BrO). Our profile retrievals have 12-20 degrees of freedom (DoF) and up to 500m vertical resolution. The tropospheric BrO vertical column density (VCD) was 1.5 x 10(exp 13) molec cm(exp -2) (RF12) and at least 0.5 x 10(exp 13) molec cm(exp -2) (RF17, 0- 10 km, lower limit). Tropospheric IO VCDs correspond to 2.1 x 10(exp 12) molec cm(exp -2) (RF12) and 2.5 x 10(exp 12) molec cm(exp -2) (RF17) and glyoxal VCDs of 2.6 x 10(exp 14) molec cm(exp -2) (RF12) and 2.7 x 10(exp 14) molec cm(exp -2) (RF17). Surprisingly, essentially all BrO as well as the dominant IO and glyoxal VCD fraction was located above 2 km (IO: 58 plus or minus 5 %, 0.1-0.2 pptv; glyoxal: 52 plus or minus 5 %, 3-20 pptv). To our knowledge there are no previous vertically resolved measurements of BrO and glyoxal from aircraft in the tropical free troposphere. The atmospheric implications are briefly discussed. Future studies are necessary to better understand the sources and impacts of free tropospheric halogens and oxygenated hydrocarbons on tropospheric ozone, aerosols, mercury oxidation and the oxidation capacity of the atmosphere.

Volkamer, R.↗
DOE OSTI2023

Quantification of hydraulic trait control on plant hydrodynamics and risk of hydraulic failure within a demographic structured vegetation model in a tropical forest (FATES–HYDRO V1.0)

Abstract. Vegetation plays a key role in the global carbon cycle and thus is an important component within Earth system models (ESMs) that project future climate. Many ESMs are adopting methods to resolve plant size and ecosystem disturbance history, using vegetation demographic models. These models make it feasible to conduct more realistic simulation of processes that control vegetation dynamics. Meanwhile, increasing understanding of the processes governing plant water use, and ecosystem responses to drought in particular, has led to the adoption of dynamic plant water transport (i.e., hydrodynamic) schemes within ESMs. However, the extent to which variations in plant hydraulic traits affect both plant water stress and the risk of mortality in trait-diverse tropical forests is understudied. In this study, we report on a sensitivity analysis of an existing hydrodynamic scheme (HYDRO) model that is updated and incorporated into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (FATES–HYDRO V1.0). The size- and canopy-structured representation within FATES is able to simulate how plant size and hydraulic traits affect vegetation dynamics and carbon–water fluxes. To better understand this new model system, and its functionality in tropical forest systems in particular, we conducted a global parameter sensitivity analysis at Barro Colorado Island, Panama. We assembled 942 observations of plant hydraulic traits on 306 tropical plant species for stomata, leaves, stems, and roots and determined the best-fit statistical distribution for each trait, which was used in model parameter sampling to assess the parametric sensitivity. We showed that, for simulated leaf water potential and loss of hydraulic conductivity across different plant organs, the four most important traits were associated with xylem conduit taper (buffers increasing hydraulic resistance with tree height), stomatal sensitivity to leaf water potential, maximum stem hydraulic conductivity, and the partitioning of total hydraulic resistance above vs. belowground. Our analysis of individual ensemble members revealed that trees at a high risk of hydraulic failure and potential tree mortality generally have a lower conduit taper, lower maximum xylem conductivity, lower stomatal sensitivity to leaf water potential, and lower resistance to xylem embolism for stem and transporting roots. We expect that our results will provide guidance on future modeling studies using plant hydrodynamic models to predict the forest responses to droughts and future field campaigns that aim to better parameterize plant hydrodynamic models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2021

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗
NASA NTRS1994

FIFE data analysis: Testing BIOME-BGC predictions for grasslands

The First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) was conducted in a 15 km by 15 km research area located 8 km south of Manhattan, Kansas. The site consists primarily of native tallgrass prairie mixed with gallery oak forests and croplands. The objectives of FIFE are to better understand the role of biology in controlling the interactions between the land and the atmosphere, and to determine the value of remotely sensed data for estimating climatological parameters. The goals of FIFE are twofold: the upscale integration of models, and algorithm development for satellite remote sensing. The specific objectives of the field campaigns carried out in 1987 and 1989 were the simultaneous acquisition of satellite, atmospheric, and surface data; and the understanding of the processes controlling surface energy and mass exchange. Collected data were used to study the dynamics of various ecosystem processes (photosynthesis, evaporation and transpiration, autotrophic and heterotrophic respiration, etc.). Modelling terrestrial ecosystems at scales larger than that of a homogeneous plot led to the development of simple, generalized models of biogeochemical cycles that can be accurately applied to different biomes through the use of remotely sensed data. A model was developed called BIOME-BGC (for BioGeochemical Cycles) from a coniferous forest ecosystem model, FOREST-BGC, where a biome is considered a combination of a life forms in a specified climate. A predominately C4-photosynthetic grassland is probably the most different from a coniferous forest possible, hence the FIFE site was an excellent study area for testing BIOME-BGC. The transition from an essentially one-dimensional calculation to three-dimensional, landscape scale simulations requires the introduction of such factors as meteorology, climatology, and geomorphology. By using remotely sensed geographic information data for important model inputs, process-based ecosystem simulations at a variety of scales are possible. The second objective of this study is concerned with determining the accuracy of the estimated fluxes from BIOME-BGC, when extrapolated spatially over the entire 15-km by 15-km FIFE site. To accomplish this objective, a topographically distributed map of soil depth at the FIFE site was developed. These spatially-distributed fluxes were then tested with data from aircraft by eddy-flux correlation obtained during the FIFE experiment.

Hunt, E. Raymond, Jr.↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗
DOE OSTI2025

Synchronous Machine Governor Upgrade

Conventional generation sources play a critical role in the stability and reliability of the electrical grid, particularly as we transition towards more renewable energy sources. To understand and accurately emulate their behavior for optimizing grid operations and ensuring seamless integration with renewable technologies, it is essential to better emulate the grid- and plant-level impacts of conventional generation sources, such as natural gas (NG) driven heat recovery steam generators (HRSGs) and combustion turbines (CTs). Therefore, a governor model is developed in a programmable logic controller (PLC) to investigate the performance of the conventional generator under various dynamic operating conditions and to identify the impact on grid stability in a controlled environment. The governor model aims to enable the hardware-in-the-loop (HIL) based emulation of these conventional generation sources using the existing 2 MVA synchronous machine/generator that is driven by a flexible 2.5 MW variable speed drive. This setup will allow us to replicate the dynamic characteristics and response behaviors of NG-driven HRSGs and CTs. The controls for the emulated conventional plants follow the industry standard and are adjustable, ensuring they accurately reflect the operational capabilities and limitations of real-world systems. These controls include load-following capabilities, ramp rates, startup and shutdown sequences, and emissions characteristics. By incorporating these adjustable controls, we aim to capture the nuanced impacts of conventional generation, such as their ability to provide ancillary services like frequency regulation, voltage support, and spinning reserve. In this report, we simulate two types of dynamic operations: grid-connected and islanding. For each dynamic operation, representative starting sequences are tested, including turbine purge, ignition, speed ramping up, generator excitation and synchronizing, and breaker close. The HIL based tests provides insights for field deployment, specifically the high-fidelity governor model provides results to predict the potential stability and reliability risk and suggest possible integration measures (e.g., generation and load balancing, tuning of governor control parameters). Ultimately, this enhanced emulation capability will be integrated into our Advanced Research on Integrated Energy Systems (ARIES), enabling us to conduct comprehensive studies on the interactions between conventional and renewable energy sources. By better understanding these interactions, we can develop strategies to optimize the overall performance and reliability of the grid. This will support the deployment of advanced grid management techniques, such as demand response, grid-forming inverters, and energy storage systems. The main contributions are summarized as follows: (1) This report introduces a PLC-based governor model for gas turbines. This model accurately simulates the dynamic behavior of conventional generation sources under various operational scenarios; (2) The model is integrated with an HIL testbed that includes a 2.5 MW variable speed drive and a 2 MVA synchronous machine. This setup enables realistic, real-time emulation of conventional power plants, particularly NG driven HRSGs and CTs; (3) The developed model is adaptable to various gas turbine configurations and allows for precise control over parameters such as MW ramp rates. This flexibility makes it a valuable tool for future research and industry collaboration; and (4) By incorporating the model into the National Renewable Energy Laboratory's Advanced Research on Integrated Energy Systems, the report lays the groundwork for future studies on interactions between conventional and renewable energy sources, enhancing the ability to develop advanced grid management strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Peak Reduction Using Mode Adjustment of Heat Pump Water Heaters in a Residential Neighborhood

Building electrification is putting pressure on distribution grid worldwide. Peak reduction is an important concern that can help reduce the growing stress and allow to defer investments in new capacity. Water heaters represent a convenient way of reducing peak because they are less dependent on weather, and their storage volume allows for asynchronous water heating and hot water use. Previous empirical studies investigated the ability of water heaters to reduce peak through the adjustment of the temperature setpoint. However, not all equipment vendors offer this option. This study aims at understanding the feasibility of peak reduction with an alternative configuration available in the market - by adjusting the device mode rather than temperature setpoint. The peak reduction methodology is tested in an occupied 46-townhome neighborhood located in Atlanta, GA. We find that peak shifting is possible with the adjustable mode approach, with the change in the peak load by 30-60%.

demand response↗
DOE OSTI2019

Roll-to-Roll Advanced Materials Manufacturing DOE Laboratory Collaboration - FY2018 Final Report

High-value R2R processing is used to support a wide range of products in applications which span many industrial business sectors. The overall R2R process technology can be considered as “mature” as the process methodology and has been in use for decades. This continuous processing technique traditionally involves deposition of material(s) onto moving webs or carriers or other continuous R2R, belt-fed, or conveyor-based processes that enable successive steps to build a final version which serves to support the deposited materials. Current process technologies, which typify “roll to roll”, include tape casting, silk-screen printing, reel-to-reel vacuum deposition/coating and R2R lithography. Products supported by R2R manufacturing include micro-electronics, electro-chromic window films, PVs, fuel cells for energy conversion, battery electrodes for energy storage, and barrier materials. Due to innovation in materials and process equipment, high-quality yet very low-cost multi-layer technologies can be manufactured on a very cost-competitive basis. To move energy-related products from high-cost niche applications to the commercial sector, a means must be available to enable manufacture of these products in a cost-competitive manner that is affordable by the general consumer. Fortunately, products such as fuel cells, thin- and mid-film PVs, batteries, electrochromic and piezoelectric films, water separation membranes, and other energy saving technologies readily lend themselves to manufacture using R2R approaches.Within the DOE EERE AMO, it was recognized that establishment of a program supported at the DOE National Laboratories, along with their immense design of materials and equipment modelling capabilities enabled with use of high-performance computing, could take advantage of available R2R infrastructure to manufacture new technologies. In FY 2016, a R2R Consortium was established and provided with initial “seed” funding to take an approach that was envisioned to be supportable by advanced manufacturing R2R processes. This collaborative approach was designed to foster identification and development of materials and processes related to R2R for clean-energy product development. Using computational and experimental capabilities by acknowledged subject matter experts within the supported National Laboratory system, this collaborative project would leverage the capabilities and expertise at each of four laboratories to further the development of an enabling high-volume cost-competitive platform technology. The collaboration team that is comprised of ORNL, ANL, NREL and LBNL, coordinating with EPB and other selected industry partners, was formed in April 2016 to initially address enhancing battery electrode performance and R2R manufacturing challenges. The research efforts were to predict and measure changes and results in electrode morphology and performance based on process condition changes; to evaluate mixed, active, particle size deposition and drying for novel electrode materials; to model various process condition changes and the resulting morphology and electrode performance; and to develop and validate NDE techniques for in-line measurement of battery electrode material properties. These efforts carried through FY 2017 and completed at the end of FY 2018.The approach was to look at compositions of materials with different particle sizes to make electrode samples using a R2R manufacturing process. The shape, size, and morphology of the materials, the chemistry of the formulation, the nature of slurries, their coating rate, the rate of drying all play a role in determining the final coating architecture, quality, and performance. A commercial cathode material was selected to make a series of cathodes and anodes by single pass, dual pass and slot die methods. Analysis of all the compiled results for this battery electrode development were that the best performing cathodes included: dual-pass electrode with large particles near the foil, mixed small and large particles, and small particles only. Whereas, the best performing anodes were with a mix of small and large graphite particles. An additional core project was added in FY 2017 to conduct studies of fuel cell materials that can be produced using R2R processes. The goal of this project is to explore, understand and optimize material and process parameters to support increased throughput, increased quality, and reduced cost for high volume production of gas-diffusion electrodes (GDEs) for PEMFCs. Project work in FY 2018 were to use a R2R process to fabricate electrodes without ionomer overlayer that can produce the equivalent mass activity as spray-coated electrodes with an ionomer overlayer. Oxygen-limiting current measurements were utilized to optimize oxygen mass transport. Alcohol-rich solvents were investigated, and results were that the alcohol and water ratio can be tuned to control ionomer distribution. Water-rich solvents produce a more dispersed ink which results in better high-current density performance.

36 MATERIALS SCIENCE↗
NASA NTRS2002

The NASA Low-Pressure Turbine Flow Physics Program: A Review

An overview of the NASA Glenn Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The flow in the LPT is unique for the gas turbine. It is characterized by low Reynolds number and high freestream turbulence intensity and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The flow of most interest is on the suction surface, where large losses are generated due to separation. The LPT is a large, multistage, heavy, jet engine component that suffers efficiency degradation between takeoff and cruise conditions due to decrease in Reynolds number with altitude. The performance penalty is around 2 points for large commercial bypass engines and as much as 7 points for small, high cruise altitude, military engines. The gas-turbine industry is very interested in improving the performance of the LPT and in reducing its weight, part count and cost. Many improvements can be accomplished by improved airfoil design, mainly by increasing the airfoil loading that can yield reduction of airfoils and improved performance. In addition, there is a strong interest in reducing the design cycle time and cost. Key enablers of the needed improvements are computational tools that can accurately predict LPT flows. Current CFD tools in use cannot yet satisfactorily predict the unsteady, transitional and separated flow in the LPT. The main reasons are inadequate transition & turbulence models and incomplete understanding of the LPT flow physics. NASA Glenn has established its LPT program to answer these needs. The main goal of the program is to develop and assess models for unsteady CFD of LPT flows. An approach that consists of complementing and augmenting experimental and computational work elements has been adopted. The work is performed in-house and by several academic institutions, in cooperation and interaction with industry. The program was reviewed at the Minnowbrook II meeting in 1997. This review will summarize the progress that was made since and will introduce newly started projects. The LPT program is focused on three areas: acquisition of experimental and numerical databases and on modeling and computation. Priority was initially given to experiments. There are three classes of experiments: simulated LPT passages, linear cascade, both with and without wakes, and low-speed rotating rig. They are being conducted as follows: At NASA GRC on a flat surface with blade pressure distribution, at the US Naval Academy on a curved surface. The addition of wakes is studied at the University of Minnesota in a curved passage with a retractable wake generator, and at Texas A&M University in a linear cascade with continuously running wake generator. The pressure distribution of the Pratt & Whitney blade 'Pak B' is used in all these experiments. Experiments have been performed also in the GEAE Low-Speed Rotating Turbine (LSRT) rig with GE-designed airfoils. Work on numerically generated database is in progress at the University of Kentucky, using the DNS/LES code LESTool developed there. Turbulence/transition model assessment and development is performed also at the University of Kentucky, where a new intermittency transport model was developed and many experimental test cases have been numerically computed. Assessments of models using simulations of multistage LPT experiments were performed at Virginia Commonwealth University using the Corsair code. Work on suction surface separation delay, using passive and active flow-control, has also been initiated. Following the overview, Principal Investigators attending the workshop will present in detail several of the projects supported by NASA.

Ashpis, David E.↗