SEARCH · Engineering Papers

Results for “dependability modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Cloud response to co-condensation of water and organic vapors over the boreal forest

Abstract. Accounting for the condensation of organic vapors along with water vapor (co-condensation) has been shown in adiabatic cloud parcel model (CPM) simulations to enhance the number of aerosol particles that activate to form cloud droplets. The boreal forest is an important source of biogenic organic vapors, but the role of these vapors in co-condensation has not been systematically investigated. In this work, the environmental conditions under which strong co-condensation-driven cloud droplet number enhancements would be expected over the boreal biome are identified. Recent measurement technology, specifically the Filter Inlet for Gases and AEROsols (FIGAERO) coupled to an iodide-adduct chemical ionization mass spectrometer (I-CIMS), is utilized to construct volatility distributions of the boreal atmospheric organics. Then, a suite of CPM simulations initialized with a comprehensive set of concurrent aerosol observations collected in the boreal forest of Finland during spring 2014 is performed. The degree to which co-condensation impacts droplet formation in the model is shown to be dependent on the initialization of temperature, relative humidity, updraft velocity, aerosol size distribution, organic vapor concentration, and the volatility distribution. The predicted median enhancements in cloud droplet number concentration (CDNC) due to accounting for the co-condensation of water and organics fall on average between 16 % and 22 %. This corresponds to activating particles 10–16 nm smaller in dry diameter that would otherwise remain as interstitial aerosol. The highest CDNC enhancements (ΔCDNC) are predicted in the presence of a nascent ultrafine aerosol mode with a geometric mean diameter of ∼ 40 nm and no clear Hoppel minimum, indicative of pristine environments with a source of ultrafine particles (e.g., via new particle formation processes). Such aerosol size distributions are observed 30 %–40 % of the time in the studied boreal forest environment in spring and fall when new particle formation frequency is the highest. To evaluate the frequencies with which such distributions are experienced by an Earth system model over the whole boreal biome, 5 years of UK Earth System Model (UKESM1) simulations are further used. The frequencies are substantially lower than those observed at the boreal forest measurement site (< 6 % of the time), and the positive values, peaking in spring, are modeled only over Fennoscandia and the western parts of Siberia. Overall, the similarities in the size distributions between observed and modeled (UKESM1) are limited, which would limit the ability of this model, or any model with a similar aerosol representation, to project the climate relevance of co-condensation over the boreal forest. For the critical aerosol size distribution regime, ΔCDNC is shown to be sensitive to the concentrations of semi-volatile and some intermediate-volatility organic compounds (SVOCs and IVOCs), especially when the overall particle surface area is low. The magnitudes of ΔCDNC remain less affected by the more volatile vapors such as formic acid and extremely low- and low-volatility organic compounds (ELVOCs and LVOCs). The reasons for this are that most volatile organic vapors condense inefficiently due to their high volatility below the cloud base, and the concentrations of LVOCs and ELVOCs are too low to gain significant concentrations of soluble mass to reduce the critical supersaturations enough for droplet activation to occur. A reduction in the critical supersaturation caused by organic condensation emerges as the main driver of the modeled ΔCDNC. The results highlight the potential significance of co-condensation in pristine boreal environments close to sources of fresh ultrafine particles. For accurate predictions of co-condensation effects on CDNC, also in larger-scale models, an accurate representation of the aerosol size distribution is critical. Further studies targeted at finding observational evidence and constraints for co-condensation in the field are encouraged.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗

DOE OSTI2021

Study of passively stable, fully detached divertor plasma regimes attained in innovative long-legged divertor configurations

Numerical modeling of divertor configurations with radially or vertically extended, tightly baffled, outer divertor legs has demonstrated the existence of a passively-stable fully detached divertor regime. In the simulations, long-legged divertors provide up to an order-of- magnitude increase in peak power handling capability compared to conventional divertors. The key physics for attaining the passively stable, fully detached regime in these simulations involves the interplay of strong convective plasma transport to the divertor leg outer sidewall, confinement of neutral gas in the divertor volume, geometric effects including a secondary X-point, and atomic radiation. New analysis shows that in this regime the detachment front location is set by the balance between the power entering the divertor leg and the losses to the walls of the divertor channel. Correspondingly, the maximum power that can be accommodated by the divertor, while still staying detached, increases with the poloidal length of the leg. The detached regime access window in terms of input power, density and impurity seeding concentration varies quantitatively depending on divertor geometry and modeling assumptions most specifically, cross-field transport to the side walls.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Light sterile neutrinos, left-right symmetry, and 0νββ decay

We investigate neutrinoless double beta (0νββ) decay rates in minimal left-right symmetric models in presence of relatively light right-handed neutrinos. By use of an effective field theory approach, we systematically include all contributions in the model as well as the dependence of the decay amplitude on the masses of right-handed neutrinos. In type-I and type-II seesaw scenarios, we analyze the impact of right-handed neutrinos heavier than about 10 MeV, showing that this effect can lead to a detection of 0νββ decay in the next-generation experiments even for the normal hierarchy and a relatively large right-handed scale set by the mass of hypothetical right-handed gauge bosons. Finally, we comment on a possible connection between light right-handed neutrinos and the strong CP problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

High-precision mass measurements of the isomeric and ground states of 44 V: Improving constraints on the isobaric multiplet mass equation parameters of the A=44, 0 + quintet

Background: The quadratic isobaric multiplet mass equation (IMME) has been very successful at predicting the masses of isobaric analog states in the same multiplet, while its coefficients are known to follow specific trends as functions of mass number. The Atomic Mass Evaluation 2016 [Chin. Phys. C 41, 030003 (2017)] 44 V mass value results in an anomalous negative c coefficient for the IMME quadratic term; a consequence of large uncertainty and an unresolved isomeric state. The b and c coefficients can provide useful constraints for construction of the isospin-nonconserving Hamiltonians for the pf shell. In addition, the excitation energy of the 0 + , T = 2 level in 44 V is currently unknown. This state can be used to constrain the mass of the more exotic 44 Cr. Purpose: The aim of the experimental campaign was to perform high-precision mass measurements to resolve the difference between 44 V isomeric and ground states, to test the IMME using the new ground state mass value and to provide necessary ingredients for the future identification of the 0 + , T = 2 state in 44 V. Method: High-precision Penning trap mass spectrometry was performed at LEBIT, located at the National Superconducting Cyclotron Laboratory, to measure the cyclotron frequency ratios of [ 44g,m VO] + versus [ 32 SCO] + , a well-known reference mass, to extract both the isomeric and ground state masses of 44 V. Results: The mass excess of the ground and isomeric states in 44 V were measured to be -23 804.9(80) keV/c 2 and -23 537.0(55) keV/c 2 , respectively. This yielded a new proton separation energy of S p = 1 773(10) keV. Conclusion: The new values of the ground state and isomeric state masses of 44 V have been used to deduce the IMME b and c coefficients of the lowest 2 + and 6 + triplets in A = 44. The 2 + c coefficient is now verified with the IMME trend for lowest multiplets and is in good agreement with the shell-model predictions using charge-dependent Hamiltonians. The mirror energy differences were determined between 44 V and 44 Sc, in line with isospin-symmetry for this multiplet. Finally, the new value of the proton separation energy determined, to an uncertainty of 10 keV, will be important for the determination of the 0 + , T = 2 state in 44 V and, consequently, for prediction of the mass excess of 44 Cr.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

PFLOTRAN-SIP: A PFLOTRAN Module for Simulating Spectral-Induced Polarization of Electrical Impedance Data

Spectral induced polarization (SIP) is a non-intrusive geophysical method that collects chargeability information (the ability of a material to retain charge) in the time domain or its phase shift in the frequency domain. Although SIP is a temporal method, it cannot measure the dynamics of flow and solute/species transport in the subsurface over long times (i.e., 10–100 s of years). Data collected with the SIP technique need to be coupled with fluid flow and reactive-transport models in order to capture long-term dynamics. To address this challenge, PFLOTRAN-SIP was built to couple SIP data to fluid flow and solute transport processes. Specifically, this framework couples the subsurface flow and transport simulator PFLOTRAN and geoelectrical simulator E4D without sacrificing computational performance. PFLOTRAN solves the coupled flow and solute-transport process models in order to estimate solute concentrations, which were used in Archie’s model to compute bulk electrical conductivities at near-zero frequency. These bulk electrical conductivities were modified while using the Cole–Cole model to account for frequency dependence. Using the estimated frequency-dependent bulk conductivities, E4D simulated the real and complex electrical potential signals for selected frequencies for SIP. These frequency-dependent bulk conductivities contain information that is relevant to geochemical changes in the system. This study demonstrated that the PFLOTRAN-SIP framework is able to detect the presence of a tracer in the subsurface. SIP offers a significant benefit over ERT in the form of greater information content. It provided multiple datasets at different frequencies that better constrained the tracer distribution in the subsurface. Consequently, this framework allows for practitioners of environmental hydrogeophysics and biogeophysics to monitor the subsurface with improved resolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

ADAM: A web platform for graph-based modeling and optimization of supply chains

Modeling and optimization are essential tasks that arise in the analysis and design of supply chains (SCs). SC models are essential for understanding emergent behavior such as transactions between participants, inherent value of products exchanged, as well as impact of externalities (e.g., policy and climate) and of constraints. Unfortunately, most users of SC models have limited expertise in mathematical optimization, and this hinders the adoption of advanced decision-making tools. Here, in this work, we present ADAM, a web platform that enables the modeling and optimization of SCs. ADAM facilitates modeling by leveraging intuitive and compact graph-based abstractions that allow the user to express dependencies between locations, products, and participants. ADAM model objects serve as repositories of experimental, technology, and socio-economic data; moreover, the graph abstractions facilitate the organization and exchange of models and provides a natural framework for education and outreach. Here, we discuss the graph abstractions and software design principles behind ADAM, its key functional features and workflows, and application examples.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Denoising Autoencoder for Reconstructing Sensor Observation Data and Predicting Evapotranspiration: Noisy and Missing Values Repair and Uncertainty Quantification

Abstract Machine learning (ML) methods applied in scientific research often deal with interrelated features in high‐dimensional data. Reducing data noise and redundancy is needed to increase prediction accuracy and efficiency especially when dealing with data from field sensors. We explored an unsupervised learning method, the denoising autoencoder (DAE), to extract the underlying data structure from noisy raw data in the context of predicting hydrologic quantities from multiple field sensors. These sensors have intrinsic instrumental noise and occasional malfunctions that cause missing values. Our DAE neural network reconstructed meteorological sensor data containing noise and missing values to predict evapotranspiration in a mountainous watershed. The DAE reconstructed the sensor variables with a mean coefficient of determination value of 0.77 across 15 dimensions representing individual sensors. It reduced variance and bias uncertainties compared to a classical autoencoder model. The reconstruction quality varied across dimensions depending on their cross‐correlation and alignment with the underlying data structure. Uncertainties arising from the model structure were overall higher than those resulting from data corruption. We attached the DAE structure to a downstream ET‐prediction neural network in three formats and achieved reasonably accurate ET predictions . The use of the DAE notably reduced variance uncertainty in ET prediction. However, excessive variance reduction may be accompanied by an increase in bias due to the intrinsic bias‐variance tradeoff. Our method of evaluating and reducing uncertainties in aggregated data from different sources can be used to improve predictive models, process understanding, and uncertainty quantification for better water resource management. Plain Language Summary We present a machine learning method, namely the denoising autoencoder, which reduces the effects of data noise and missing values typically present in scientific data sets collected through sensor measurements. This method selects the most relevant information from noisy raw data collected by the instruments and fills in missing values. To demonstrate the effectiveness of our method, we applied it to predict evapotranspiration, a hydrologic variable that represents the water moved from the land surface to the atmosphere through a combination of evaporation and plant water use (transpiration). We also used a random sampling technique (the Monte Carlo method) to compare the uncertainty in the predictions when using the raw and noisy data versus the reconstructed data. The denoising process produced more accurate predictions of evapotranspiration with less uncertainty. Improved predictions of evapotranspiration can lead to a better understanding and accounting of water budgets. This ML approach is broadly suitable for a wide variety of applications that involve noisy sensor data with missing values. Key Points We used a denoising autoencoder (DAE) neural network to reduce noise in meteorological and soil sensor observations by on average We used Monte Carlo sampling to estimate the bias and variance of all model outputs, including uncertainty sources from data and the model We attached the DAE component to a downstream neural network to predict ET with the variance reduced by , compared to that without the DAE

denoising autoencoder↗

DOE OSTI2021

Substantial hysteresis in emergent temperature sensitivity of global wetland CH4 emissions

Abstract Wetland methane (CH 4 ) emissions ( $${F}_{{{CH}}_{4}}$$ F C H 4 ) are important in global carbon budgets and climate change assessments. Currently, $${F}_{{{CH}}_{4}}$$ F C H 4 projections rely on prescribed static temperature sensitivity that varies among biogeochemical models. Meta-analyses have proposed a consistent $${F}_{{{CH}}_{4}}$$ F C H 4 temperature dependence across spatial scales for use in models; however, site-level studies demonstrate that $${F}_{{{CH}}_{4}}$$ F C H 4 are often controlled by factors beyond temperature. Here, we evaluate the relationship between $${F}_{{{CH}}_{4}}$$ F C H 4 and temperature using observations from the FLUXNET-CH 4 database. Measurements collected across the globe show substantial seasonal hysteresis between $${F}_{{{CH}}_{4}}$$ F C H 4 and temperature, suggesting larger $${F}_{{{CH}}_{4}}$$ F C H 4 sensitivity to temperature later in the frost-free season (about 77% of site-years). Results derived from a machine-learning model and several regression models highlight the importance of representing the large spatial and temporal variability within site-years and ecosystem types. Mechanistic advancements in biogeochemical model parameterization and detailed measurements in factors modulating CH 4 production are thus needed to improve global CH 4 budget assessments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Varied phenomenology of models displaying dynamical large-deviation singularities

Singularities of dynamical large-deviation functions are often interpreted as the signal of a dynamical phase transition and the coexistence of distinct dynamical phases, by analogy with the correspondence between singularities of free energies and equilibrium phase behavior. Here we study models of driven random walkers on a lattice. These models display large-deviation singularities in the limit of large lattice size, but the extent to which each model's phenomenology resembles a phase transition depends on the details of the driving. We also compare the behavior of ergodic and nonergodic models that present large-deviation singularities. Finally, we argue that dynamical large-deviation singularities indicate the divergence of a model timescale, but not necessarily one associated with cooperative behavior or the existence of distinct phases.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

A Review of Bayesian Networks for Spatial Data

We report Bayesian networks are a popular class of multivariate probabilistic models as they allow for the translation of prior beliefs about conditional dependencies between variables to be easily encoded into their model structure. Due to their widespread usage, they are often applied to spatial data for inferring properties of the systems under study and also generating predictions for how these systems may behave in the future. We review published research on methodologies for representing spatial data with Bayesian networks and also summarize the application areas for which Bayesian networks are employed in the modeling of spatial data. We find that a wide variety of perspectives are taken, including a GIS-centric focus on efficiently generating geospatial predictions, a statistical focus on rigorously constructing graphical models controlling for spatial correlation, as well as a range of problem-specific heuristics for mitigating the effects of spatial correlation and dependency arising in spatial data analysis. Special attention is also paid to potential future directions for integration of Bayesian networks with spatial processes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Flux effects in precipitation under irradiation simulation of Fe-Cr alloys

Radiation-enhanced precipitation of Cr-rich α’ in irradiated Fe-Cr alloys, which results in hardening and embrittlement, depends on the irradiating particle and the displacement per atom (dpa) rate. Here, we utilize a Cahn-Hilliard phase-field based approach, that includes simple models for nucleation, irradiating particle and rate dependent radiation-enhanced diffusion and cascade mixing to simulate α’ evolution under neutrons, heavy ions, and electron irradiations. Different irradiating particles manifest very different cascade mixing efficiencies. The model was calibrated using neutron data. For cascade inducing neutron/heavy-ion dpa rates at 300 °C between 10-8 and 10-6 dpa/s the model predicts approximately constant number density, decreasing radius, decreasing α’ Cr composition, and lower α’ volume fraction. The model then predicts a dramatic transition to no α’ formation above approximately 10-5 dpa/s, while electron irradiation, with weak mixing, had little effect at dpa rates up to 10-3 dpa/s. These model predictions are consistent with experiments. We explain the results in terms of the flux dependence of the radiation enhanced diffusion, cascade mixing, and their ratio, which all vary significantly in relevant flux ranges for neutron and cascade inducing ion irradiations. These results show that both cascade mixing and radiation enhanced diffusion must be accounted for when attempting to emulate neutron-irradiation effects using accelerated ion irradiations.

Ke, Jia-Hong↗

DOE OSTI2023

Observation of Weak Counterion Size Dependence of Thermoelectric Transport in Ion Exchange Doped Conducting Polymers Across a Wide Range of Conductivities

Conducting polymers are of interest for a broad range of applications from bioelectronics to thermoelectrics. The factors that govern their complex charge transport physics include the structural disorder present in these highly doped polymer films and the Coulombic interactions between the electronic charge carriers and the dopant counterions. Previous studies have shown that at low doping levels carriers are strongly trapped in the vicinity of the counterions, while at high doping levels charge transport is not limited by Coulombic trapping, which manifests itself in the conductivity being independent of the size of the dopant counterion. Here a recently developed ion exchange doping method is used to investigate the ion size dependence of a semicrystalline polythiophene–based model system across a wide range of conductivities. It is found that the regime in which the charge and thermoelectric transport is not or only weakly dependent on ion size, extends to surprisingly low conductivities. This surprising observation is explained by a heterogeneous doping that involves doping of the amorphous domains to high doping levels first before doping of the ordered, crystalline domains occurs. The study provides new insights into how the thermoelectric physics of conducting polymers evolves as a function of doping level.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Process prediction and detection of faults using probabilistic bidirectional recurrent neural networks on real plant data

Attaining Industry 4.0 for manufacturing operations requires advanced monitoring systems and real-time data analytics of plant data, among other topics. We propose a Probabilistic Bidirectional Recurrent Network (PBRN) for industrial process monitoring for the early detection of faults. The model is based on a Gated Recurrent Unit (GRU) neural network that allows the model to retain long-term dependencies between sensor data along a time horizon, hence learning the dynamic behavior of the process. To reduce the false-positive detection rate of the model, we compel the model to learn from a highly noisy sensor reading while outputting noise-free sensor outputs. The performance of the proposed model is compared to other data-driven statistical process monitoring schemes using real plant data from an industrial Air Separations Unit (ASU) containing noisy sensor readings. We show that the model can learn from noisy data without reducing its performance. Using two different fault cases, we demonstrate the model’s ability to carry out early fault detection with average false-positive rates of 2.9% and 4.9% for both fault cases. The missed detection rates are 0.1% and 0.2%, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Selective atomic sieving across metal/oxide interface for super-oxidation resistance

Abstract Surface passivation, a desirable natural consequence during initial oxidation of alloys, is the foundation for functioning of corrosion and oxidation resistant alloys ranging from industrial stainless steel to kitchen utensils. This initial oxidation has been long perceived to vary with crystal facet, however, the underlying mechanism remains elusive. Here, using in situ environmental transmission electron microscopy, we gain atomic details on crystal facet dependent initial oxidation behavior in a model Ni-5Cr alloy. We find the (001) surface shows higher initial oxidation resistance as compared to the (111) surface. We reveal the crystal facet dependent oxidation is related to an interfacial atomic sieving effect, wherein the oxide/metal interface selectively promotes diffusion of certain atomic species. Density functional theory calculations rationalize the oxygen diffusion across Ni(111)/NiO(111) interface, as contrasted with Ni(001)/NiO(111), is enhanced. We unveil that crystal facet with initial fast oxidation rate could conversely switch to a slow steady state oxidation.

(S)TEM, oxidizing metals↗

DOE OSTI2021

Understanding Hepatitis B Virus Dynamics and the Antiviral Effect of Interferon Alpha Treatment in Humanized Chimeric Mice

Whereas the mode of action of lamivudine (LAM) against hepatitis B virus (HBV) is well established, the inhibition mechanism(s) of interferon alpha (IFN-α) is less completely defined. To advance our understanding, we, for this study, mathematically modeled HBV kinetics during 14-day pegylated IFN-α-2a (pegIFN), LAM, or pegIFN-plus-LAM (pegIFN+LAM) treatment of 39 chronically HBV-infected humanized uPA/SCID chimeric mice. Serum HBV DNA and intracellular HBV DNA were measured frequently. We developed a multicompartmental mathematical model and simultaneously fit it to the serum and intracellular HBV DNA data. Unexpectedly, even in the absence of an adaptive immune response, a biphasic decline in serum HBV DNA and intracellular HBV DNA was observed in response to all treatments. Kinetic analysis and modeling indicate that the first phase represents inhibition of intracellular HBV DNA synthesis and secretion, which was similar under all treatments with an overall mean efficacy of 98%. In contrast, there were distinct differences in HBV decline during the second phase, which was accounted for in the model by a time-dependent inhibition of intracellular HBV DNA synthesis, with the steepest decline observed during pegIFN+LAM treatment (1.28/day) and the slowest (0.1/day) during pegIFN monotherapy. Reminiscent of observations in patients treated with pegIFN and/or LAM, a biphasic HBV decline was observed in treated humanized mice in the absence of an adaptive immune response. Interestingly, combination treatment did not increase the initial inhibition of HBV production but rather enhanced second-phase decline, providing insight into the dynamics of HBV treatment response and the mode of action of IFN-α against HBV.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

Stochastic Unit Commitment: Model Reduction via Learning

As weather-dependent renewable generation increases its share in the generation mix of most electric energy systems, a stochastic unit commitment becomes the natural day-ahead scheduling tool. However, such a tool is generally computationally intractable if a detailed uncertainty description is considered. Taking this into account, we proposed a learning method to make the stochastic unit commitment problem tractable. Here, recent advances in statistical learning and machine learning to address optimization problems can be advantageously applied to the rather intractable stochastic unit commitment problem. Considering these advances, we explore simple learning techniques to drastically reduce the size of a stochastic unit commitment problem without significantly altering its optimal solution. The considered stochastic unit commitment problem is formulated as a two-stage stochastic programming problem. The first stage represents commitment decisions, while the second one represents the operation conditions under different scenarios. Taking into account historical solved instances (or proxies for them), we reduce the size (measured by numbers of constraints and variables) of the stochastic unit commitment problem by (i) fixing unchanged binary variables and by (ii) eliminating inactive inequality constraints. Our numerical results show that the reduced problem generally requires significantly less time to solve while obtaining high-quality solutions, which are very close to or indistinguishable from the one obtained by solving the original problem. We use an Illinois 200-bus system to illustrate and characterize the performance of the proposed problem-reduction method.

42 ENGINEERING↗