Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Final Technical Report

The Department of Energy is interested in technologies that support the sustainable production of fuels, chemicals, and other bioproducts from plant biomass, to offset the nation’s reliance on fossil resources. The plant cell wall of energy crops provides the largest reservoir of raw materials for bioproducts. However, the widespread use of plant cell walls is hampered by their complexity and resistance to breakdown. To improve the productivity and cost-effectiveness of using energy crops to generate bioproducts, the fundamental problem of deconstructing plant cell walls must be addressed. This project developed and evaluated an innovative genetic modification technology to produce strategically designed enzymes that specifically accumulate in the plant cell wall. The resulting enzyme-engineered energy crops are expected to grow normally under natural conditions but break down more quickly and easily under high temperature during the production of biobased products. As such, this plant cell wall targeting enzyme engineering effort will reduce the cost of plant cell wall deconstruction and ultimately improve the economics of bioproducts. The overall objective of this project is to develop and evaluate the in-planta enzyme engineering technology to reduce lignocellulose deconstruction cost. The concept was first validated using tobacco plant, a model plant system that is typically used in lab testing for initial concept validation. Then the enzyme optimization was validated using switchgrass, the energy crop to be used to produce bioproducts. There are three specific objectives in this Phase I project: (1) validate the enzyme optimization concept using tobacco plant, a model plant system. (2) validate the enzyme optimization concept using switchgrass. (3) techno-economic analysis (TEA) for further scale-up application. By the end of this project, in-planta enzyme engineering was validated in both tobacco and switchgrass plants, with improved enzyme activity and saccharification efficiency. The in-planta enzyme engineering in Tabacco didn’t have a significant impact on plant growth and development. Transgenic tobacco plants with in-planta cellulose degrading enzymes showed higher biomass digestibility than wild type. Gene construction and transformation in switchgrass was much longer than expected, which delayed the research progress. Besides, in-planta engineering of lignin degrading enzyme is more challenging than cellulose degrading enzyme, in terms of expression detection. Expression of lignin degrading enzyme and cellulose degrading enzyme improved biomass yield and saccharification efficiency of switchgrass, respectively. It is promising to express both genes in switchgrass for optimized overall performance. According to the results of TEA, switchgrass biomass production cost is mainly attributed to by fertility and harvesting. Biomass production profit can increase up to 10-fold depending on biomass price. The PHA production profit is also sensitive to the biomass price. The proposed technology could potentially reduce the biomass deconstruction cost from 33% to 9% of PHA revenue, making the biomass-based PHA competitive to petroleum-based polymers even in case of relatively high biomass price of biomass. Therefore, cultivation of the genetically engineered self-deconstruction switchgrass for Polyhydroxyalkanoate (PHA) production could benefit switchgrass grower and PHA producer with attractive profits for both sectors. This new enzyme optimization approach will be beneficial for bioindustries that use energy crops as feedstocks. It will improve the economic viability of converting energy crops to renewable products that support a sustainable society and helps address the Nation’s long-term strategic needs for renewable products and reduction of reliance on fossil resources.

42 ENGINEERING↗
DOE OSTI2021

A classical, non-singular, bouncing universe

We present a model for a classical, non-singular bouncing cosmology without violation of the null energy condition (NEC). The field content is General Relativity plus a real scalar field with a canonical kinetic term and only renormalizable, polynomial-type self-interactions for the scalar field in the Jordan frame. The universe begins vacuum-energy dominated and is contracting at t=-∞. We consider a closed universe with a positive spatial curvature, which is responsible for the universe bouncing without any NEC violation. An Rφ 2 coupling between the Ricci scalar and the scalar field drives the scalar field from the initial false vacuum to the true vacuum during the bounce. The model is sub-Planckian throughout its evolution and every dimensionful parameter is below the effective-field-theory scale M P , so we expect no ghost-type or tachyonic instabilities. This model solves the horizon problem and extends co-moving particle geodesics to past infinity, resulting in a geodesically complete universe without singularities. We solve the Friedman equations and the scalar-field equation of motion numerically, and analytically under certain approximations.

Astronomy & Astrophysics↗
DOE OSTI2023

Series FACTS Devices for Increasing Resiliency in Severe Weather Conditions

Severe weather conditions are low-probability, high-impact events that affect grid operations. The majority of power outages are caused by severe weather conditions. Grid resiliency to weather events can be enhanced by decreasing the reliance on its affected sections. One way to do this is to reduce the power flow through lines vulnerable to severe weather. If a line is disconnected, its initial power flow is distributed through the neighbor lines, which may cause congestion in the grid. FACTS devices can be used to control the power flow of lines that have a higher chance of power outages. Most previous works do not consider weather events in power flow control. In this work, a linearized optimal power flow (OPF)–based algorithm is developed to minimize the real power flow of vulnerable lines considering the thermal limits of lines to prevent infeasible solutions; the simulation is fast, making it suitable for large-scale systems. The proposed optimization problem is presented as a mixed-integer linear program (MILP), making it capable of using short-term load forecasting due to its high solution speed. The proposed optimization problem considers multiple lines with different outage probabilities and the uncertainties of the weather forecast. Moreover, it estimates the power reduction in vulnerable lines due to changes in the series FACTS devices. The performance of the proposed optimization problem is tested on IEEE 14-, 30-, and 118-bus systems for several scenarios. The results are validated with the AC power flow results from MATPOWER.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

CO 2 Adsorption and Hydrogenation on Inverse InO x /Cu(111) Catalysts: Active Role of the Oxide–Metal Interface

The direct conversion of carbon dioxide (CO 2 ) into methanol via hydrogenation is essential for industrial applications. Recent studies on catalysts that contain an inverse oxide/metal configuration have shown very good catalytic performance for the CO 2 hydrogenation to methanol process. Here, in this study, we investigated the behavior of indium oxide-Cu(111) interfaces under pure CO 2 and CO 2 /H 2 mixtures using synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS). Initially, a single layer of copper oxide (Cu x O) was grown on the Cu(111) surface by controlled oxidation. On this surface, indium was deposited at room temperature. Oxygen atoms transferred from Cu x O/Cu(111) to the indium metal upon deposition, forming In-O-Cu bonds and active interfaces. Although Cu(111) is not very active for the binding and activation of CO 2 , the formed InO x -Cu(111) interfaces had no problem adsorbing and dissociating the molecule at room temperature. Reaction of CO 2 with H 2 on InO x -Cu(111) yielded surface-bound H 3 CO, CO 2 δ− , CO 3 , and CH x species that are typical intermediates in the production of methanol and other oxygenates. The InO x -Cu(111) interface underwent dynamic chemical changes under reaction conditions, forming In-Cu alloys at low indium coverages (< 0.05 monolayer), while at higher indium coverages a mixture of an In-Cu and InO x was detected in XPS. These findings indicate that InO x /In-Cu interfaces can play a key role in processes aimed at the trapping and valorization of CO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES↗
DOE OSTI2022

Field demonstration and implementation analysis of model predictive control in an office HVAC system

Model Predictive Control (MPC) is a promising technique to address growing needs for heating, ventilation, and air-conditioning (HVAC) systems to operate more efficiently and with greater flexibility. However, due to a number of factors, including the required implementation expertise, lack of high quality data, and a risk-adverse industry, MPC has yet to gain widespread adoption. While many previous studies have shown the advantages of MPC, few analyzed the implementation effort and associated practical challenges. In addition, previous work has developed an open-source, Modelica-based tool-chain that automatically generates optimal control, parameter estimation, and state estimation problems aimed at facilitating MPC implementation. Therefore, this study demonstrates usage of this tool-chain to implement MPC in a real office building, discusses practical challenges of implementing MPC, and estimates the implementation effort associated with various tasks in order to inform the development of future workflows and serve as an initial benchmark for their impact on reducing implementation effort. This study finds that the implemented MPC saves approximately 40% of HVAC energy over the existing control during a two-month trial period and that tasks related to data collection and controller deployment activities can each require as much effort as model generation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2023

Understanding and Mitigating the Dissolution and Delamination Issues Encountered with High-Voltage LiNi 0.5 Mn 1.5 O 4

In our initial study on the high-voltage 5 V cobalt-free spinel LiNi 0.5 Mn 1.5 O 4 (LNMO) cathode, we discovered a severe delamination issue in the laminates when cycled at a high upper cut-off voltage (UCV) of 4.95 V, especially when a large cell format was used. This delamination problem prompted us to investigate further by studying the transition metal (TM) dissolution mechanism of cobalt-free LNMO cathodes, and as a comparison, some cobalt-containing lithium nickel manganese cobalt oxides (NMC) cathodes, as the leachates from the soaking experiment might be the culprit for the delamination. Unlike other previous reports, we are interested in the intrinsic stability of the cathode in the presence of a baseline Gen2 electrolyte consisting of 1.2 M of LiPF 6 in ethylene carbonate/ethyl methyl carbonate (EC/EMC), similar to a storage condition. The electrode laminates (transition metal oxides, transition metal oxides, TMOs, coated on an Al current collector with a loading level of around 2.5 mAh/cm 2 ) or the TMO powders (pure commercial quality spinel LNMO, NMC, etc.) were stored in the baseline solution, and the transition metal dissolution was studied through nuclear magnetic resonance, such as 1 H NMR, 19 F NMR, scanning electron microscope (SEM), X-ray photoelectron spectroscopy (XPS) and inductively coupled plasma mass spectrometry (ICP-MS). Significant electrolyte decomposition was observed and could be the cause that leads to the TM dissolution of LNMO. To address this TM dissolution, additives were introduced into the baseline electrolyte, effectively alleviating the issue of TM dissolution. The results suggest that the observed delamination is caused by electrolyte decompositions that lead to etching, and additives such as lithium difluorooxalato borate and p-toluenesulfonyl isocyanate can alleviate this issue by forming a firm cathode electrolyte interface. This study provides a new perspective on cell degradation induced by electrode/electrolyte interactions under storage conditions.

25 ENERGY STORAGE↗
DOE OSTI2021

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗
DOE OSTI2020

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Research on Deformation Law of Deep Foundation Pit of Station in Core Region of Saturated Soft Loess Based on Monitoring

Saturated soft loess has a large pore structure, high compressibility, low strength, fluid plastic state, and poor engineering properties. It is still one of the key problems that engineering needs to solve. In order to study the influence of deep foundation pit excavation in the saturated soft loess area on the deformation of foundation pits, the deformation laws of ground settlement, enclosure structure, and supporting axial force were clarified based on the field tests on the deformation characteristics of deep foundation pits in areas with the high-water level in saturated soft loess, combined with geological conditions and on-site construction procedures. The results indicate the following: water supply and construction process were found to be the main factor in changing the surface settlement curve of deep foundation pits in saturated soft loess; increasing the construction speed of the pit bottom floor, inverted braces, floor frame beams, and sidewall frame beams to close the structure, which is conducive to restraining the deformation of the continuous underground wall and foundation pits in similar areas. In the initial stage of support layout, the axial force of steel support tends to increase too fast or even exceed the standard control value. Therefore, a reasonable preadding axial force is an effective means to control the deformation of the continuous underground wall and the axial force of the steel support. The current research results may provide a reference for constructing deep foundation pits in similar areas.

Wang, Xueyan↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗
DOE OSTI2023

Maritime Fuel Cell Generator Project: 2018 – 2023

This report summarizes activity in the Maritime Fuel Cell (MarFC) Generator Project from 2018 – 2023. FY 2018 saw the implementation of upgrades and repairs, making the unit more reliable and operator friendly. In FY2019 the team engaged the Scripps Institution of Oceanography (SIO) to use the MarFC to provide zero-emission shore power to the research vessel R/V Robert Gordon Sproul while in port at the Nimitz Marine Facility in San Diego, CA. In FY2020, the MarFC unit was shipped to San Diego, CA. A fueling contract with IGX was established to support MarFC operations at SIO, with renewable hydrogen provided by the California State University Los Angeles (CSULA) hydrogen station. The project team (Sandia, Cummins/Hydrogenics) provided training to SIO staff on the technical details of the MarFC, the safe use of hydrogen in general and the MarFC in particular. The first fueling of the unit at the SIO pier was successfully completed by IGX. The first powering of a vessel with fuel-cell shore power was conducted with the R/V Robert Gordon Sproul. While the mechanical systems (lights, AC, ventilation, hydraulics, pumps and cranes) were powered without incident, problems arose when powering the computer systems. Inspections carried out in FY2021 revealed the MarFC needed routine maintenance. Maintenance was performed and the unit was upgraded. The MarFC was turned on after the year pause, and initial test data on power levels and stability were collected. FY2022 was a year spent repairing, upgrading and testing the MarFC unit. Spikes in power and voltage were observed above 60 kW that could potentially extend below 60 kW with time. Such spikes could cause problems with the Sproul electrical systems. These age-related problems, the extended time for the Sproul spent in dry dock for scheduled upgrades, and the associated need to reschedule the vessel’s high-priority science missions made it no longer possible to deploy the unit at SIO. After due consideration, the decision was made by DOE, MARAD and the project team to cease the deployment, remove the MarFC from SIO, and not pursue further deployment activities. On December 2, 2022, the MarFC unit was removed from the Scripps Nimitz Marine Facility and shipped to Fridley, Michigan. The Cummins/Hydrogenics plan for the unit is to assess the condition of the MarFC subcomponents, and then use it as a training/learning system for technical employees new to hydrogen fuel-cell technology. After summarizing project activity from 2018 – 2023, this report provides a review of lessons learned. This report provides next steps in contemplating a follow-on project that would further advance the use of fuel-cell-based shore power in a marine setting. A comparison is made of the project results to the original objectives. This report ends with an accounting of presentations stemming from the project, and a list of references.

08 HYDROGEN↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2023

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A continuum damage approach to spallation and the role of microinertia

Spall failure is of interest due to its prevalence in high strain rate problems in which the spallation is driven by the interaction of release waves. In this article, a porosity-based damage model that includes microinertial effects is used to examine spall failure. The model is successfully calibrated to plate impact-driven experiments and then used to evaluate experimental conditions producing more extreme strain rate conditions, such as those in laser-driven experiments. The incorporation of microinertia allows us to better understand the increase in apparent macroscopic spall strength seen at high strain rates. Correspondingly, we conclude that the incorporation of microinertial effects improves the model’s predictive capabilities. Microinertial effects result in more severe local tensile stresses that affect the damage evolution, and microinertia can play a significant role in the approach to the material’s ideal strength at extreme loading rates. Additionally, a preliminary parametric study is also carried out to investigate the role of microstructural aspects such as nucleation volume fraction and initial pore radius. One counter-intuitive result from the microinertial effects is that, for a given nucleation site volume fraction, having larger initial pore nucleation sites can lead to an improved macroscopic spall strength.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of a hybrid neural network and transfer learning model for optimized ICP-MS/MS operation

Correct function and calibration of instrumentation is a crucial assumption for any scientific experiment. One such instrument, tandem inductively coupled plasma mass spectrometer (ICP-MS/MS), has in-depth calibration settings that range across 30+ different parameters, making it difficult to determine optimal conditions without expertise and some degree of trial and error. Often, these settings are hand-tuned, a time-intensive process prone to local maxima and human error. While some automation is available, the automation also may favor local optimizations over a global optimum. In addition to these difficulties, day to day instrument variability can further complicate the calibration process. We propose a solution to this problem as a machine learning (ML) algorithm that learns how each parameter helps determine the calibration sensitivity across several elements, and re-weights parameters over time as instrument variability changes (e.g., a global neural network (NN) with a time-dependent transfer learning (TL) component). This model would be able to generate a surface of predicted calibration sensitivities and their respective parameters, and a simple multivariate algorithm would be able to pull out the optimum results with the settings associated with them. Here-in, we describe our initial findings in working towards this goal, including data extraction from historical files, exploratory data analysis, and some initial model building to better describe the data and the feasibility of our goal.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗