Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SURROGATE MODELS”

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 1,315 records · Page 73

Enhanced Data Efficiency Using Deep Neural Networks and Gaussian Processes for Aerodynamic Design Optimization

Adjoint-based optimization methods are attractive for aerodynamic shape design primarily due to their computational costs being independent of the dimensionality of the input space and their ability to generate high-fidelity gradients that can then be used in a gradient-based optimizer. This makes them very well suited for high-fidelity simulation based aerodynamic shape optimization of highly parametrized geometries such as aircraft wings. However, the development of adjoint-based solvers involve careful mathematical treatment and their implementation require detailed software development. Furthermore, they can become prohibitively expensive when multiple optimization problems are being solved, each requiring multiple restarts to circumvent local optima. In this work, we propose a machine learning enabled, surrogate-based framework that replaces the expensive adjoint solver, without compromising on predicting predictive accuracy. Specifically, we first train a deep neural network (DNN) from training data generated from evaluating the high-fidelity simulation model on a model-agnostic design of experiments on the geometry shape parameters. The optimum shape may then be computed by using a gradient-based optimizer coupled with the trained DNN. Subsequently, we also perform a gradient-free Bayesian optimization, where the trained DNN is used as the prior mean. We observe that the latter framework (DNN-BO) improves upon the DNN-only based optimization strategy for the same computational cost. Overall, this framework predicts the true optimum with very high accuracy, while requiring far fewer high-fidelity function calls compared to the adjoint-based method. Furthermore, we show that multiple optimization problems can be solved with the same machine learning model with high accuracy, to amortize the offline costs associated with constructing our models. Our methodology finds applications in the early stages of aerospace design. (C) 2021 Published by Elsevier Masson SAS.

Renganathan, S. Ashwin↗

A scalable planning framework of energy storage systems under frequency dynamics constraints

As the penetration of renewables increases in power systems, the declining system inertia can cause frequency stability issues. Battery energy storage systems (BESSs) respond fast and therefore can relieve the low inertia difficulty but need to be appropriately sized considering the associated cost. This paper presents a novel stochastic optimization model for economically planning BESS capacity while considering the spatial–temporal correlation of wind generation and generator outages under frequency stability constraints, which include the rate-of-change of frequency (RoCoF), frequency nadir (FN), and quasi-steady-state (QSS) frequency. A set of new FN constraints that can be easily linearized is developed. To account for renewable uncertainties, a realistic uncertainty modeling approach, Random Field, is adopted to generate wind generation scenarios by considering both spatial and temporal evolutions of wind speed profiles. The ESS sizing is formulated as a mixed-integer linear programming problem and solved by using a scalable decomposition-and-coordination approach, Surrogate Absolute Value Lagrangian Relaxation (SAVLR). To further improve the scalability and reduce computational burdens, a rolling-horizon-based update is developed and incorporated into SAVLR for providing a practical solution to the long-term planning of very large-scale power systems. Finally, a modified IEEE 118-bus system and the Polish system are used to validate the effectiveness and scalability of the model and solution methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluation and improvement of the parameterization of aerosol hygroscopicity in global climate models using in-situ surface measurements (Final Report)

Aerosols are tiny particles suspended on the atmosphere that can interact with incoming solar radiation and affect the Earth radiative budget. They do so by scattering and absorbing solar radiation, and these properties vary depending on the aerosol size and chemical composition. Moreover, by taking up water from the surrounding air, hygroscopic aerosol particles will grow in size and change their chemical composition, thus modifying their optical properties (scattering and absorption) and their final impact on radiative forcing calculations. An accurate knowledge of aerosol hygroscopicity is crucial for estimating the net radiative impact of aerosols. We took a three-pronged approach to improve our understanding of aerosol hygroscopicity and how it is implemented in Earth system models. In the first part of our project, we developed a benchmark dataset from existing aerosol hygroscopic growth measurements made by tandem nephelometer humidogram systems. We analyzed, using a standardized methodology, observations from 26 in-situ stations around the globe. Measurement data was collected from multiple data providers, reviewed and harmonized to create a consistent dataset of the scattering enhancement factor due to aerosol water uptake. This dataset is archived in several publicly available databases for use by interested researchers. In the second part of the project, we used the benchmark hygroscopicity dataset to perform a global study on aerosol hygroscopicity and aerosol optical properties. Measurements show a global picture of scattering enhancement with larger values for Arctic and marine sites and lower for urban and desert sites. We assessed the RH dependence of aerosol radiative forcing and showed that the overall effect of aerosol hygroscopicity on DARF is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 compared to dry conditions (RH<40%). Finally, we explored using aerosol single scattering albedo (SSA) and scattering Angstrom exponent (SAE) as possible proxies for aerosol hygroscopicity. SSA showed more promise as a surrogate for the scattering enhancement factor than SAE, but neither was ideal. In the third part of the project we evaluated the output of ten Earth system models (ESMs) against the benchmark hygroscopicity dataset. ESMs utilize various schemes for aerosol hygroscopicity which had not been previously tested against observations on a global scale. We found that ESMs currently overestimate scattering enhancement due to hygroscopic growth. Model parameterizations of hygroscopicity and model chemistry are two main factors driving the observed diversity in hygroscopicity simulations among the models. In addition, our study makes several suggestions for modelers, including improving the parameterizations of organic and sea salt aerosol hygroscopicity. Future hygroscopicity model evaluation experiments should include the model data related to particles size which was not available for our study.

54 ENVIRONMENTAL SCIENCES↗

Comparative analysis of SARS-CoV-2 neutralization titers reveals consistency between human and animal model serum and across assays

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires ongoing monitoring to judge the ability of newly arising variants to escape the immune response. A surveillance system necessitates an understanding of differences in neutralization titers measured in different assays and using human and animal serum samples. We compared 18 datasets generated using human, hamster, and mouse serum and six different neutralization assays. Datasets using animal model serum samples showed higher titer magnitudes than datasets using human serum samples in this comparison. Fold change in neutralization of variants compared to ancestral SARS-CoV-2, immunodominance patterns, and antigenic maps were similar among serum samples and assays. Most assays yielded consistent results, except for differences in fold change in cytopathic effect assays. Hamster serum samples were a consistent surrogate for human first-infection serum samples. These results inform the transition of surveillance of SARS-CoV-2 antigenic variation from dependence on human first-infection serum samples to the utilization of serum samples from animal models.

Cell Biology↗

Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems

The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.

CFD/DEM↗

Estimating and Forecasting Time-Varying Groundwater Recharge in Fractured Rock: A State-Space Formulation with Preferential and Diffuse Flow to the Water Table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

groundwater, Recharge, infiltration, recursive est↗

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗

Blended fuel property analysis of butyl-exchanged polyoxymethylene ethers as renewable diesel blendstocks

Methyl-terminated polyoxymethylene ethers (MM-POMEs), having the formula CH 3 O-(CH 2 O)n-CH 3 (n = 3-6), are a class of oxygenates with desirable diesel-like fuel properties including high cetane number and low soot formation. However, their low energy density and high water-solubility present barriers to their adoption. Both concerns were recently addressed by our research group by synthesizing a mixture of POME structures having butyl end-groups and n = 1-6, termed B*POME1-6. B*POME1-6 maintained the advantageous properties of the parent MM-POMEs, and exhibited improved energy density and most notably, dramatically decreased water solubility. For evaluation against a set of criteria for a blended diesel blendstock, a 20 vol% blend of B*POME1-6 with a base diesel fuel was investigated here. Oxidation stability, cetane number, sooting tendency, lubricity and conductivity were improved in the B*POME1-6 blend compared with the base diesel, while also maintaining the flash point, cloud point, energy density, viscosity, and boiling point requirements. The B*POME1-6 product demonstrated a synergistic blending behavior at 10 vol% and a linear blending behavior at 20-30 vol% blends, in agreement with similar POME blends at comparable blend levels. Finally, common environmental and toxicity models performed on B*POME1-6 component molecules suggested they have a greater propensity to partition into the water compartment compared to a common diesel surrogate, but with a lower tendency to bioaccumulate.

33 ADVANCED PROPULSION SYSTEMS↗

Validation of Pronghorn for Natural-Circulation Molten Salt Loops

This paper presents the development and validation of a high-fidelity thermal-hydraulic model of a molten salt natural circulation flow loop, designed for integration within a digital twin framework. The study evaluates the performance of Idaho National Laboratory’s Pronghorn against experimental data from Texas A&M University Molten Salt Flow Loop (MSFL) four Hitec-salt test benchmark data. Natural circulation of high-Prandtl-number fluids exhibits complex, counter-intuitive flow patterns that make pointwise thermocouple readings unreliable. Experimental work at TAMU’s MSFL provides benchmark data, including flow visualization at a test-section and centerline steady-state temperature measurements along the loop. Validation includes four single-phase natural circulation test cases with Hitec salt. Key metrics include flow profile agreement and steady-state temperature accuracy. Pronghorn results for two-dimensional single-phase agree qualitatively with the experimental flow profile. This paper illustrates the importance of Computational Fluid Dynamics (CFD) in elucidating the behavior of high-Prandtl-number thermal-hydraulics, along with how misleading centerline temperature measurements can be. Pronghorn reproduces the axial and radial stratification that makes single thermocouple readings unreliable. Future research will focus on reduced-order modeling techniques to enable rapid simulation suitable for real-time digital twin applications. The validated cases provide a basis for developing reduced-order surrogates aimed at real-time digital-twin applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Binding of radionuclides and surrogate to 18-crown-6 ether by density functional theory

For this work, we use density functional theory to investigate the interactions of cerium, americium, and curium cations with crown ethers. Our calculations reveal that the modeled structure of cerium integrated within the crown ether is in good agreement with experimental data, with the negative binding energy indicating that capturing the cerium nitrates is thermodynamically favorable. Our results demonstrate that crown ethers can also bind americium and curium, providing insights into the potential applications of crown ether in radionuclide sequestration. Finally, we explore the impact of the skeleton modification of different crown ethers through by substitution of nitrogen atoms in the core of the crown ether for oxygen atoms and find that this structural modification significantly increases the radionuclide binding energies. These findings provide insights on the potential for the use of organic linkers such as crown ethers to address the urgent needs in radionuclide sequestration, separation and sensing.

36 MATERIALS SCIENCE↗

Improved Compaction Experiments and Modeling of Waste Isolation Pilot Plant Standard, Non-degraded, Waste Containers

A credible simulation of disposal room porosity at the Waste Isolation Pilot Plant (WIPP) requires a tenable compaction model for the 55-gallon waste containers within the room. A review of the legacy waste material model, however, revealed several out-of-date and untested assumptions that could affect the model’s compaction behavior. For example, the legacy model predicted non-physical tensile out-of-plane stresses under plane strain compression. (Plane strain compression is similar to waste compaction in the middle of a long drift.) Consequently, a suite of new compaction experiments were performed on containers filled with surrogate, non-degraded, waste. The new experiments involved uniaxial, triaxial, and hydrostatic compaction tests on quarter-scale and full-scale containers. Special effort was made to measure the volume strain during uniaxial and triaxial tests, so that the lateral strain could be inferred from the axial and volume strain. These experimental measurements were then used to calibrate a pressure dependent, viscoplastic, constitutive model for the homogenized compaction behavior of the waste containers. This new waste material model’s predictions agreed far better with the experimental measurements than the legacy model’s predictions, especially under triaxial and hydrostatic conditions. Under plane strain compression, the new model predicted reasonable compressive out-of-plane stresses, instead of tensile stresses. Moreover, the new model’s plane strain behavior was substantially weaker for the same strain, yet substantially stronger for the same porosity, than the legacy model’s behavior. Although room for improvement exists, the new model appears ready for prudent engineering use.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Experimental Measurements of Soot and Soot Precursors of Sustainable Aviation Fuels

Sustainable Aviation Fuels (SAFs) are promising for decarbonizing the challenging-to-electrify aerospace sector and reducing yearly greenhouse gas emissions by thousands of tons. Developing predictive detailed chemical kinetic models for SAF (and traditional) fuel mixtures necessitates a robust experimental database to elucidate the intricate pyrolysis and oxidation chemistry of SAF fuel blend components. Two critical components of SAF surrogate mixtures are iso-dodecane isomers and trans-decalin. This report characterizes the structure of three slightly sooting, laminar, non-premixed, nitrogen-diluted Planar Mixing Layer Flames (PMLFs) fueled by either pure ethylene, ethylene doped with 1500ppm of either 2,2,4,6,6-pentamethyl-heptane or trans-decalin. The flames have the same stoichiometric mixture fraction and total mole fraction of fuels in the fuel stream (X F,F = X C2H4,F + X Dopant,F = 0.26), resulting in approximately the same maximum temperature (T max ≈ 1800 K) to highlight the effect of doping on soot and its gaseous precursors. The PMLF is unconfined and established between adjacent planar jets of the fuel and oxidizer streams surrounded by annular shielding nitrogen. The flow is stabilized by exhausting the hot buoyant flame products through a slot in a downstream plate onto which the cold flow impinges. Importantly, any horizontal cross-section of the PMLF has a self-similar thermochemical structure of a reactive boundary layer that grows at increasing Height Above the Burner (HAB). Therefore, the horizontal cross-sections of a PMLF can be modeled as a One-Dimensional Counterflow Flame (1D-CF) with very low strain rates which are unachievable in traditional CFs and yield a structure that is several millimeters thick. In brief, the PMLF is ideal for performing spatially resolved measurements with minimal effects of the sampling-induced perturbations, to develop and validate detailed chemical kinetic models on time scales of a few tens of milliseconds that are relevant in aviation applications. The horizontal profiles of C0-C18 gaseous species at HAB = 50 mm and soot volume fraction at 25 and 50 mm are measured using capillary sampling followed by GC-MS analyses and Laser Induced Emission Spectroscopy (LIES), respectively. Additionally, the Elastic Laser Light Scattering (E-LLS) coefficient is measured at HAB = 50 mm and 25 mm to determine the strain rates of the 1D-CFs equivalent to the characterized PMLF horizontal cross sections and the profile of the E-LLS equivalent diameter of soot. The comparison of the results in the base PMLF fueled only by ethylene with those in either doped PMLF in which 1500 ppm of ethylene are replaced with 2,2,4,6,6-pentamethyl-heptane or trans-decalin, reveals that the doping increases the mole fraction of several Polycyclic Aromatic Hydrocarbons (PAHs) at both HABs by a factor of approximately 2.5, and the soot volume fraction at both HABs by factors of ~1.5 for 2,2,4,6,6-pentamethyl-heptane and ~1.25 for trans-decalin doping, respectively. Conversely, doping causes a near doubling of the soot equivalent diameter only at HAB = 50 mm, without significant effects at 25 mm. The experimental results partially validate the tested state-of-the-art detailed chemical kinetic model and also point toward the further improvement of its predictive capabilities.

10 SYNTHETIC FUELS↗

sMF-BO-2CoGP: A sequential multi-fidelity constrained Bayesian optimization framework for design applications

Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.

42 ENGINEERING↗

Localization of Nonblinking Point Sources Using Higher-Order-Mode Detection and Optical Heterodyning: Developing a Strategy for Extending the Scope of Molecular, Super-resolution Imaging

While the stochastic, “blinking” nature of fluorescent systems has enabled the super-resolution of their localization by the fitting of their point-spread functions (PSFs), this strategy cannot be exploited for similar resolution of “nonblinking” systems, such as those that might be encountered in a coherent Raman experiment. An alternative method for subdiffraction-limited imaging lies in the exploitation of optical heterodyning. For example, if a Gaussian PSF (a TEM00 mode) of a point emitter is displaced with respect to the origin of the optical system, photons in the higher-order TEM modes carry information about that displacement. Information concerning the displacement can be extracted from photons in these higher-order modes. These photons can be collected by optical heterodyning, which exploits the large gain in a detector’s response to an optical signal from an emitter coupled to a local oscillator, which is prepared in the TEM of interest, e.g., TEM10. We have generalized and developed the heterodyning technique to localize point emitters via the detection of higher-order spatial modes. We have developed a theoretical approach to find a practical estimation limit of the localization parameters using a realistic model that accounts for shot noise, background noise, and Gaussian noise. To demonstrate the applicability of the method, we designed experiments in which a laser is a surrogate for one and two point emitters. Using the Fisher information and its accompanying Cramér-Rao lower bound, we demonstrate super-resolution localization in these cases: we show that objects can be localized to roughly 2–3 orders of magnitude of their point-spread function’s size for a given optical system. Finally and most importantly, it is suggested that the results will ultimately be generalizable to multiple emitters and, most importantly, to “nonblinking” molecular systems, which will be essential for broadening the scope of super-resolution measurements beyond the limits of fluorescence-based techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Use of Satellite Nightlights for Power Outages Prediction

Hurricanes are a dominant disaster in the Caribbean, always causing serious power outages throughout the islands. Hurricane Maria was a prime example, causing unimaginable destruction of the power infrastructure of Puerto Rico (PR). Consequently, one month after the hurricane landfall, approximately 80% of the population was still without power. After an event of such massive destruction, the electric power restoration process progresses very slowly. This timeline can be improved using power outage (PO) forecast models that help identify the vulnerable places before the hurricane landfall. Generally, these models are trained with historical power outages records, associated data on weather conditions, and additional information about the natural and built environments. However, PO records are often difficult to acquire, and, in many instances, the power utility companies may not record them. This study utilizes a satellite-based Visible Infrared Imaging Radiometer Suite (VIIRS) night light data product as a surrogate for the power delivery to predict hurricane-induced PO in areas having limited to nonexistent historical data records. The processed satellite data is then used along with geographic variables, and simulated weather data to formulate machine learning-based algorithms to predict PO for future hurricane events. These models are applied and validated in the context of the PR catastrophic storm, Hurricane Maria.

54 ENVIRONMENTAL SCIENCES↗

Dakota A Multilevel Parallel Object-Oriented Framework for Design Optimization Parameter Estimation Uncertainty Quantification and Sensitivity Analysis: Version 6.12 Theory Manual

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗