Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 199 records · Page 11

Three-dimensional cooperative inversion of airborne magnetic and gravity gradient data using deep-learning techniques

Using multiple geophysical methods has become a prevailing approach in numerous geophysical applications to investigate subsurface structures and parameters. These multimethod-based exploration strategies have the potential to greatly diminish uncertainties and ambiguities encountered during geophysical data analysis and interpretation. One of the applications is the cooperative inversion of airborne magnetic and gravity gradient data for the interpretation of data obtained in mineral, oil and gas, and geothermal explorations. In this paper, a unified cooperative inversion framework is designed by combining the standard separate inversions with a deep neural network (DNN), which serves as the link between different types of data. A well-trained DNN takes the separately inverted susceptibility and density models as the inputs and provides improved models that will be used as the initial models of deterministic inversions. A two-round iteration strategy is adopted to guarantee the reasonability of the recovered models and overall efficiency of the inversion. In addition, this deep-learning (DL)-based framework demonstrates excellent generalization abilities when tested on models that are entirely distinct from the training data sets. The framework can easily incorporate multiphysics without necessitating any structural changes to the network. Synthetic experiments validate that our DL-based method outperforms conventional separate inversions and cross-gradient-based joint inversion in view of the accuracy of the recovered models and inversion efficiency. Successful application to field data further verifies the effectiveness of our DL-based method.

Geochemistry & Geophysics↗

Chemistry Simulations Using MERRA-2 Reanalysis with the GMI CTM and Replay in Support of the Atmospheric Composition Community

Simulations using reanalyzed meteorological conditions have been long used to understand causes of atmospheric composition change over the recent past. Using the new Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) meteorology, chemistry simulations are being conducted to create products covering 1980-2016 for the atmospheric composition community. These simulations use the Global Modeling Initiative (GMI) chemical mechanism in two different models: the GMI Chemical Transport Model (CTM) and the GEOS-5 model developed Replay mode. Replay mode means an integration of the GEOS-5 general circulation model that is incrementally adjusted each time step toward the MERRA-2 analysis. The GMI CTM is a 1 x 1.25 simulation and the MERRA-2 GMI Replay simulation uses the native MERRA-2 approximately horizontal resolution on the cubed sphere. The Replay simulations is driven by the online use of key MERRA-2 meteorological variables (i.e. U, V, T, and surface pressure) with all other variables calculated in response to those variables. A specialized set of transport diagnostics is included in both runs to better understand trace gas transport and changes over the recent past.

Chemistry-Climate Model↗

Chemistry Simulations using the MERRA-2 Reanalysis with the GMI CTM and Replay in Support of the Atmospheric Composition Community

Simulations using reanalysis meteorological fields have long been used to understand the causes of atmospheric composition change in the recent past. Using the new MERRA-2 reanalysis, we are conducting chemistry simulations to create products covering 1980-2016 for the atmospheric composition community. These simulations use the Global Modeling Initiative (GMI) chemical mechanism in two different models: the GMI Chemical Transport Model (CTM) and the GEOS-5 model in Replay mode. Replay mode means an integration of the GEOS-5 general circulation model that is incrementally adjusted each time step toward the MERRA-2 reanalysis. The GMI CTM is a 1 deg x 1.25 deg simulation and the MERRA-2 GMI Replay simulation uses the native MERRA-2 grid of approximately 1/2 deg horizontal resolution on the cubed sphere. A specialized set of transport diagnostics is included in both runs to better understand trace gas transport and its variability in the recent past.

MERRA-2↗

AeroCom Phase III Multi-Model Evaluation of the Aerosol Life Cycle and Optical Properties Using Ground and Space-Based Remote Sensing as Well as Surface In Situ Observations

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

space-based remote sensing↗

Assessment of NO2 Observations During DISCOVER-AQ and KORUS-AQ Field Campaigns

NASA’s Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ, conducted in 2011–2014) campaign in the United States and the joint NASA and National Institute of Environmental Research (NIER) Korea–United States Air Quality Study (KORUS-AQ, conducted in 2016) in South Korea were two field study programs that provided comprehensive, integrated data sets of airborne and surface observations of atmospheric constituents, including nitrogen dioxide (NO2), with the goal of improving the interpretation of spaceborne remote sensing data. Various types of NO2 measurements were made, including in situ concentrations and column amounts of NO2 using ground- and aircraft-based instruments, while NO2 column amounts were being derived from the Ozone Monitoring Instrument (OMI) on the Aura satellite. This study takes advantage of these unique datasets by first evaluating in situ data taken from two different instruments on the same aircraft platform, comparing coincidently sampled profile-integrated columns from aircraft spirals with remotely sensed column observations from ground-based Pandora spectrometers, intercomparing column observations from the ground (Pandora), aircraft (in situ vertical spirals), and space (OMI), and evaluating NO2 simulations from coarse Global Modeling Initiative (GMI) and high-resolution regional models. We then use these data to interpret observed discrepancies due to differences in sampling and deficiencies in the data reduction process. Finally, we assess satellite retrieval sensitivity to observed and modeled a priori NO2 profiles. Contemporaneous measurements from two aircraft instruments that likely sample similar air masses generally agree very well but are also found to differ in integrated columns by up to 31.9 %. These show even larger differences with Pandora, reaching up to 53.9 %, potentially due to a combination of strong gradients in NO2 fields that could be missed by aircraft spirals and errors in the Pandora retrievals. OMI NO2 values are about a actor of 2 lower in these highly polluted environments due in part to inaccurate retrieval assumptions (e.g., a priori pro-files) but mostly to OMI’s large footprint (>312 km2).

Nitrogen dioxide↗

Diabatic Initialization of Mesoscale Models in the Southeastern United States: Can 0 to 12h Warm Season QPF be Improved?

It is well known that numerical warm season quantitative precipitation forecasts lack significant skill for numerous reasons. Some are related to the model--it may lack physical processes required to realistically simulate convection or the numerical algorithms and dynamics employed may not be adequate. Others are related to initialization-mesoscale features play an important role in convective initialization and atmospheric observation systems are incapable of properly depicting the three-dimensional stability structure at the mesoscale. The purpose of this study is to determine if a mesoscale model initialized with a diabatic initialization scheme can improve short-term (0 to 12h) warm season quantitative precipitation forecasts in the Southeastern United States. The Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory is used to diabatically initialize the Pennsylvania State University/National center for Atmospheric Research (PSUNCAR) Mesoscale Model version 5 (MM5). The SPORT Center runs LAPS operationally on an hourly cycle to produce analyses on a 15 km covering the eastern 2/3 of the United States. The 20 km National Centers for Environmental Prediction (NCEP) Rapid Update Cycle analyses are used for the background fields. Standard observational data are acquired from MADIS with GOES/CRAFT Nexrad data acquired from in-house feeds. The MM5 is configured on a 140 x 140 12 km grid centered on Huntsville Alabama. Preliminary results indicate that MM5 runs initialized with LAPS produce improved 6 and 12h QPF threat scores compared with those initialized with the NCEP RUC.

Lapenta, William M.↗

Region-based convolutional neural network for wind turbine wake characterization from scanning lidars

A convolutional neural network is applied to lidar scan images from three experimental campaigns to identify and characterize wind turbine wakes. Initially developed as a proof-of-concept model and applied to a single data set in complex terrain, the model is now improved and generalized and applied to two other unique lidar data sets, one located near an escarpment and one located offshore. The model, initially developed using lidar scans collected in predominantly westerly flow, exhibits sensitivity to wind flow direction. The model is thus successfully generalized through implementing a standard rotation process to scan images before input into the convolutional neural network to ensure the flow is westerly. The sample size of lidar scans used to train the model is increased, and along with the generalization process, these changes to the model are shown to enhance accuracy and robustness when characterizing dissipating and asymmetric wakes. Applied to the offshore data set in which nearly 20 wind turbine wakes are included per scan, the improved model exhibits a 95% success rate in characterizing wakes and a 74% success rate in characterizing dissipating wake fragments. The improved model is shown to generalize well to the two new data sets, although an increase in wake characterization accuracy is offset by an increase in model sensitivity and false positive wake identifications.

17 WIND ENERGY↗

Development of the INL Thermal Energy Distribution System (TEDS) in the Modelica Eco-System for Validation and Verification

This report provides an overview on the development of the Thermal Energy Distribution System (TEDS) in the Modelica process model eco-system at Idaho National Laboratory as part of the integrated energy system initiative. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Engineering Science Laboratory (ESL) at INL. The model was then used during the pre-construction phase of the experimental effort to inform experimental design (insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit including: a 200kW Chromalox heater, a single-tank packed bed thermal energy storage system filled with 0.125 inch Alumina (Al 2 O 3 ), an ethylene-glycol to Therminol-66 heat exchanger, system piping, six control valves, and all associated temperature, pressure, and mass flow sensors. The model does not include nitrogen fill gas tanks, or associated overfill tanks, as these are not part of standard system control. Two simulation sets were run. The first is a five-hour test that operates as a shakedown test for the facility. This simulation puts the facility through all five potential operating modes and showcases the ability of valving, control sensors, and component controllers to meet the system demands. The second case imposes a typical summer day demand on the system from a region with mixed commercial and residential electrical needs where the generator alone cannot meet peak demand but instead requires the thermal storage unit to act as a peaking unit. Over the course of the fifteen-hour discharge cycle heat losses account for a heat loss of approximately 15°C in outlet temperature. Such heat losses would lead to a decrease in overall system efficiency and may lead to additional design changes for such designs in the future that may include a topping heater. Through commencement of this work, a systems level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model has been utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Once TEDS begins operation, this model will be refined, tuned, and used for validation and verification purposes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Observation of a single protein by ultrafast X-ray diffraction

* runs/ - The raw data from the pnCCD * histograms.h5 - Histogram of the values recorded by each pixel throughout the data acquisition. Used for pedestal and gain correction. * mask.h5 - Detector bad pixels mask. * 1SS8.pdb - PDB model of the assembled GroEL 14-mer, from the 1SS8 PDB entry. * 1SS8_H2O.pdb - Hydrated version of the assembled GroEL 14-mer, derived from the 1SS8 PDB entry. The various hydrated models are created by removing water from this initial model. * rotations.h5 - Quaternions of the rotations used for template matching

European XFEL↗

Loss terms in free-piston Stirling engine models

Various models for free piston Stirling engines are reviewed. Initial models were developed primarily for design purposes and to predict operating parameters, especially efficiency. More recently, however, such models have been used to predict engine stability. Free piston Stirling engines have no kinematic constraints and stability may not only be sensitive to the load, but also to various nonlinear loss and spring constraints. The present understanding is reviewed of various loss mechanisms for free piston Stirling engines and how they have been incorporated into engine models is discussed.

Gordon, Lloyd B.↗

Photolysis Rate Coefficient Calculations in Support of SOLVE Campaign

The objectives for this SOLVE project were 3-fold. First, we sought to calculate a complete set of photolysis rate coefficients (j-values) for the campaign along the ER-2 and DC-8 flight tracks. En route to this goal, it would be necessary to develop a comprehensive set of input geophysical conditions (e.g., ozone profiles), derived from various climatological, aircraft, and remotely sensed datasets, in order to model the radiative transfer of the atmosphere accurately. These j-values would then need validation by comparison with flux-derived j-value measurements. The second objective was to analyze chemistry along back trajectories using the NASA/Goddard chemistry trajectory model initialized with measurements of trace atmospheric constituents. This modeling effort would provide insight into the completeness of current measurements and the chemistry of Arctic wintertime ozone loss. Finally, we sought to coordinate stellar occultation measurements of ozone (and thus ozone loss) during SOLVE using the Midcourse Space Experiment(MSX)/Ultraviolet and Visible Imagers and Spectrographic Imagers (UVISI) satellite instrument. Such measurements would determine ozone loss during the Arctic polar night and represent the first significant science application of space-based stellar occultation in the Earth's atmosphere.

Lloyd, Steven A.↗

Photolysis Rate Coefficient Calculations in Support of SOLVE Campaign

The objectives for this SOLVE project were 3-fold. First, we sought to calculate a complete set of photolysis rate coefficients (j-values) for the campaign along the ER-2 and DC-8 flight tracks. En route to this goal, it would be necessary to develop a comprehensive set of input geophysical conditions (e.g., ozone profiles), derived from various climatological, aircraft, and remotely sensed datasets, in order to model the radiative transfer of the atmosphere accurately. These j-values would then need validation by comparison with flux-derived j-value measurements. The second objective was to analyze chemistry along back trajectories using the NASA/Goddard chemistry trajectory model initialized with measurements of trace atmospheric constituents. This modeling effort would provide insight into the completeness of current measurements and the chemistry of Arctic wintertime ozone loss. Finally, we sought to coordinate stellar occultation measurements of ozone (and thus ozone loss) during SOLVE using the MSX/UVISI satellite instrument. Such measurements would determine ozone loss during the Arctic polar night and represent the first significant science application of space-based stellar occultation in the Earth's atmosphere.

Lloyd, Steven A.↗

Acoustic-Modal Testing of the Ares I Launch Abort System Attitude Control Motor Valve

The Attitude Control Motor (ACM) is being developed for use in the Launch Abort System (LAS) of NASA's Ares I launch vehicle. The ACM consists of a small solid rocket motor and eight actuated pintle valves that directionally allocate.thrust_- 1t.has-been- predicted-that significant unsteady. pressure.fluctuations.will.exist. inside the-valves during operation. The dominant frequencies of these oscillations correspond to the lowest several acoustic natural frequencies of the individual valves. An acoustic finite element model of the fluid volume inside the valve has been critical to the prediction of these frequencies and their associated mode shapes. This work describes an effort to experimentally validate the acoustic finite model of the valve with an acoustic modal test. The modal test involved instrumenting a flight-like valve with six microphones and then exciting the enclosed air with a loudspeaker. The loudspeaker was configured to deliver broadband noise at relatively high sound pressure levels. The aquired microphone signals were post-processed and compared to results generated from the acoustic finite element model. Initial comparisons between the test data and the model results revealed that additional model refinement was necessary. Specifically, the model was updated to implement a complex impedance boundary condition at the entrance to the valve supply tube. This boundary condition models the frequency-dependent impedance that an acoustic wave will encounter as it reaches the end of the supply tube. Upon invoking this boundary condition, significantly improved agreement between the test data and the model was realized.

Davis, R. Benjamin↗

Multidecadal Changes in the UTLS Ozone from the MERRA-2 Reanalysis and the GMI Chemistry Model

Long-term changes of ozone in the UTLS (Upper Troposphere / Lower Stratosphere) reflect the response to decreases in the stratospheric concentrations of ozone-depleting substances as well as changes in the stratospheric circulation induced by climate change. To date, studies of UTLS ozone changes and variability have relied mainly on satellite and in-situ observations as well as chemistry-climate model simulations. By comparison, the potential of reanalysis ozone data remains relatively untapped. This is despite evidence from recent studies, including detailed analyses conducted under SPARC (Scalable Processor Architecture) Reanalysis Intercomparison Project (S-RIP), that demonstrate that stratospheric ozone fields from modern atmospheric reanalyses exhibit good agreement with independent data while delineating issues related to inhomogeneities in the assimilated observations. In this presentation, we will explore the possibility of inferring long-term geographically and vertically resolved behavior of the lower stratospheric (LS) ozone from NASA's MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications -2) reanalysis after accounting for the few known discontinuities and gaps in its assimilated input data. This work builds upon previous studies that have documented excellent agreement between MERRA-2 ozone and ozonesonde observations in the LS. Of particular importance is a relatively good vertical resolution of MERRA-2 allowing precise separation of tropospheric and stratospheric ozone contents. We also compare the MERRA-2 LS ozone results with the recently completed 37-year simulation produced using Goddard Earth Observing System in "replay" mode coupled with the GMI (Global Modeling Initiative) chemistry mechanism. Replay mode dynamically constrains the model with the MERRA-2 reanalysis winds, temperature, and pressure. We will emphasize the areas of agreement of the reanalysis and replay and interpret differences between them in the context of our increasing understanding of model transport driven by assimilated winds.

Ozone↗

A Multi-Satellite Assimilation and Modeling Platform to Construct a Global and Complete view of the Hydrologic Cycle

What captivates me about the research field of satellite remote sensing of the earth system science is the opportunity to help human beings by addressing challenging questions related to the current and future availability of natural resources for domestic, agricultural, and industrial needs. For example, a continuously varying climate poses a threat to the water resource, not only in developing countries, but also in economically well-developed regions, creating the urgent need to characterize and predict its spatial and temporal availability. Specifically, my curiosity is driven by the challenge of merging cutting-edge space technology and earth observations (i.e., remote sensing) with state-of-the-art models for the purpose of improving our scientific knowledge about the variability and the change of the hydrologic cycle.The only practical way to observe the land surface processes (e.g., hydrologic cycle) on continental to global scales is via satellite remote sensing. Though remote sensing can make spatially comprehensive measurements of various components of the land surface system, it cannot provide direct information on the entire system, and the measurements represent only a snapshot in time. Physical based models may be used to continuously predict the temporal and spatial earth processes, but these predictions are often poor, due to model initialization, parameter and meteorological forcing errors, and inadequate model physics and/or resolution. Thus, an attractive prospect is to combine the strengths of land surface models and observations (and minimize the weaknesses of both) to provide a superior land surface state estimate. This is the goal of land surface data assimilation. Data assimilation provides a better estimate of the environmental states than either models or observations could individually do. The broad concept of data assimilation can be applied to various disciplines such as hydrology, ecology, environmental hazards, agriculture and economy.My research is targeted to integrating multiple satellite information having multi-sensor and multi-resolution assimilation within today's state-of-the-art hydrologic models. Multi-sensor and multi-resolution assimilation techniques represent a necessary milestone in future earth science applications because they offer the capability of comprehensively integrating disparate types of earth observations via deep-learning techniques. A multi-satellite assimilation and modeling platform will ultimately provide a robust and complete dataset to more fully understand earth's dynamics. In fact, the multi-satellite assimilation platforms build a comprehensive description of the earth processes that is geared toward accurately representing the complexity of natural and anthropogenic interactions in land surface processes. Most importantly, these platforms are suitable as decision support tools for applications across different aspects of earth system science.My recent work has focused on recent multi-sensor data assimilation techniques targeted at improving snow, soil moisture, groundwater, and terrestrial water storage hydrological states.

Girotto, Manuela↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – An Uncertainty Analysis

This report summarizes spreading and heat transfer results for molten salt pouring onto a flat stainless steel substrate that were calculated using MELTSPREAD to gain insight into bulk salt behavior during a molten salt spill accident. MELTSPREAD was developed at Argonne to model the one-dimensional flow and freezing behavior of molten corium and is being applied to model molten salts for the first time as part of this work. An uncertainty analysis of thermophysical properties and spill conditions was performed using MELTSPREAD to determine which properties and conditions have the greatest impact on model outcome. The model was run with and without the inclusion of decay heat. Eutectic FLiNaK was used as the salt composition because it is well-characterized and appropriate for model development. Small salt volumes were used for the initial model runs described in this report so that the model results can be compared to the results from ongoing experiments on molten salt spreading and heat transfer currently being conducted at Argonne at a benchtop scale. The spreading of a small volume of FLiNaK (50 mL) was found to be limited by the balance between the molten salt surface tension and gravity and not by freezing for the pour conditions and assumed mechanisms of heat transfer in the model. Heat transfer from the salt was highly inefficient to due to its low thermal conductivity, high heat capacity, and high heat of fusion and varying these properties had no effect on the spreading behavior. Changing the viscosity and density of the salt slightly affected the spreading behavior. Setting the salt surface tension to approximately zero to simulate the salt wetting the substrate had the greatest effect on the spreading behavior and this increased both the spreading velocity and overall spreading area. Including decay heat in the model had a significant effect on the long-term cooling behavior of the salt but had no effect on the overall spreading behavior. Future effort should involve modeling a larger volume of molten salt to represent a reactor case spill scenario because the spreading behavior of large volumes of salt may differ from the small volumes that were the focus of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Effects of diffusion and mass flows on C IV and Si IV lines formed in the solar atmosphere

A model for the transition region is derived from an initial interpretation of EUV observations, assuming ionization equilibrium and constant elemental abundance with height. The effects of diffusion and mass flows are then included in the initial model and the emergent profiles of several C IV and Si IV lines are computed. It is found that diffusion and mass flows have a strong effect on both the emergent intensity and spectral shape of these lines. Diffusion acts to deplete the transition region of heavy ions to an extent which depends on the detailed temperature and density structure. The net effect is a weighting of the coronal emission relative to that in the transition. On the basis of the results, it is suggested that the downflows observed in the network in lines of C IV and Si IV could be due to gravitational settling of the ions following their injection, via diffusion, into the corona from spicules. When flows are superposed on the basic diffusion model, the rapid change in elemental abundance characteristic of the model is virtually eliminated and the coronal contribution to the emission in C IV and Si IV lines becomes negligible relative to that in the transition region. Flows have a strong effect on the computed line intensities and introduce large asymmetries into the line shapes.

Roussel-Dupre, R.↗