Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “material flows”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Development, construction and qualification tests of the mechanical structures of the electromagnetic calorimeter of the Mu2e experiment at Fermilab

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the Charged Lepton Flavour Violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of physics beyond the Standard Model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. The calorimeter provides excellent electron identification, complementary information to aid pattern recognition and track reconstruction, and a fast calorimetric online trigger. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1340 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements.The hostile Mu2e operational conditions, in terms of radiation levels (total ionizing dose of 12 krad and a neutron fluence of 5x1010 n/cm2 @ 1 MeVeq (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on the design of the detector mechanical structures and materials choice. The support structure of the two 670 crystal matrices employs two aluminum hollow rings and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and service front-end electronics for each crystal are assembled in a unique mechanical unit inserted in a machined copper holder. The 670 units are supported by a machined plate made of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics is hosted in aluminum custom crates positioned on the external lateral surface of the two disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system.The constraints on the calorimeter mechanical structures design, the development from the conceptual design to the specifications of all the structural components, the status of components production and the components quality assurance tests are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Cooling charge determination for a local thermal management system

A system and method for estimating the state of charge, i.e., the amount of cooling remaining, in a latent heat thermal storage module of a local environmental control system. The local environmental control system stores cooling energy at night by freezing a latent heat phase change material (PCM), and then releasing the stored cooling energy by removing heat from room air flowing through the PCM during the day. Differential temperature sensors are used measure the air temperature drop across the latent thermals storage unit and the fan voltage, which is directly related to fan air flow rate, is sampled. A controller is programmed to determine the charging state from these measurements.

Khalifa, H. Ezzat↗
DOE OSTI2021

Modelling Nuclear Thermal Propulsion Reactor Startup Transients

The National Aeronautics and Space Administration (NASA) has set the goal of a manned mission to Mars by the year 2030 [1] and charged the national academy of sci- ences "to identify primary technical and programmatic chal- lenges, merits, and risks for maturing space nuclear propulsion technologies of interest to a future human Mars exploration mission" [2]. One relevant technology, nuclear thermal propul- sion (NTP), has notable advantages over traditional chemical rockets; most important among them is the ability to produce larger specific impulse on the order of 900s. The reduction of mission time is crucial for a manned mission to Mars to reduce the risk for the crew. Due to its higher specific impulse, NTP satisfies this need and is pursued as one technology to get humans to Mars [3, 4]. The construction of an NTP sys- tem has to negotiate several challenges laid out in Ref. [2]; one of these challenges is the need to startup the NTP sys- tem from essentially cold conditions to full power within one minute. This paper focuses on studying the neutronics and thermal-hydraulics behavior of a simplified NTP model dur- ing prescribed reactivity insertions and mass flow rate (MFR) ramps. It is the goal of this paper to investigate startup, peak material temperatures, and average specific impulse for a low enriched Uranium (LEU), ceramic and metal material (CER- MET) NTP system when varying reactivity insertion and MFR ramps.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Controlling product selectivity in hybrid gas/liquid reactors using gas conditions, voltage, and temperature

For the conversion of CO 2 into fuels and chemical feedstocks, hybrid gas/liquid-fed electrochemical flow reactors provide advantages in selectivity and production rates over traditional liquid phase reactors. However, fundamental questions remain about how to optimize conditions to produce desired products. Here, using an alkaline electrolyte to suppress hydrogen formation and a gas diffusion electrode catalyst composed of copper nanoparticles on carbon nanospikes, we investigate how hydrocarbon product selectivity in the CO 2 reduction reaction in hybrid reactors depends on three experimentally controllable parameters: (1) supply of dry or humidified CO 2 gas, (2) applied potential, and (3) electrolyte temperature. Changing from dry to humidified CO 2 dramatically alters product selectivity from C 2 products ethanol and acetic acid to ethylene and C 1 products formic acid and methane. Water vapor evidently influences product selectivity of reactions that occur on the gas-facing side of the catalyst by adding a source of protons that alters reaction pathways and intermediates.

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

Performance of a Hybrid HVAC-Integrated Thermal Storage Device

Thermal equipment in buildings is a primary contributor to peak loads on the electrical grid. Thermal energy storage is a cost-effective strategy to decouple electric use from thermal loads, thus reducing grid peak costs for building owners. One method for storing thermal energy in a building is to integrate a phase change material (PCM) directly into the heating, ventilation, and air conditioning system. These systems often require additional glycol loops, pumps, valves, and heat exchangers to couple the storage to the cooling system and building space, which increases the complexity and cost. This work will discuss an alternate approach where the storage is added directly into the heat pump evaporator. A detailed two-dimensional finite difference heat transfer model of a PCM-refrigerant-glycol heat exchanger was developed to simulate the performance of this component. The fluid stream was discretized along the flow direction to capture changes in the fluid properties and local heat transfer rates, and the phase change material was discretized in both the x and y directions to capture the movement of the melt front. The model was used to understand the impact of different material and geometric properties on the charge and discharge characteristics of the device. Finally, a Ragone framework analogous to that used for electrochemical batteries was used to maximize the energy density and round-trip efficiency of the device while supplying loads appropriate for space cooling in buildings.

buildings↗
DOE OSTI2023

Final Report DE-FE0031785 Mohsen Ahmadian, Ph.D. Demonstration of Proof of Concept of a Multiphysics Approach for Real-Time Remote Monitoring of Dynamic Changes in Pressure and Salinity in Hydraulically Fractured Networks

Hydraulic fracturing has evolved into a multistep process with varying flow rates, carrier fluids (e.g., gel or slickwater), proppant loadings, and proppant grain sizes. As a result, primary recovery from a hydraulically fractured tight-oil reservoir is often a tiny fraction of the original oil in place, ranging between 5 and 10%. As stated in the FOA1990, “part of this problem is due to the inability of current well completion processes to effectively stimulate the entire reservoir volume in contact with the wellbore. Innovative technologies are needed that can help improve the effectiveness of reservoir completion methods, maximize stimulated reservoir volumes, and optimize recovery over the entire producing life span of a well”. We first need to enhance the current fracture diagnostic techniques to improve a well-completion design. However, detecting and delineating a subsurface hydraulic fracture is extremely difficult because the induced fracture network is only fractionally propped, and these propped fractures are generally very thin. Microseismic and tiltmeter monitoring techniques can provide information on the fracture extent but provide little or no information on the movement and final distribution of proppant or production fluids. On the other hand, electromagnetic (EM) imaging has shown the capability to monitor proppant distribution throughout the fracture area, especially in the presence of Electrically Active Proppants (EAPs). A previous EM survey of hydraulic fracturing at the Devine Fracture Pilot Site (DFPS) and subsequent EM code developments demonstrated this survey as a robust technique to remotely interrogate the extent of the EAP-filled hydraulic fracture during its propagation. The objectives of the project were threefold: (1) to capitalize on the material properties of an EAP to demonstrate remote monitoring of relative changes in pressure, pressure, and flow that are commonly encountered during production from a hydraulically fractured reservoir; (2) to evaluate EM imaging tools, to achieve Objective 1 in near real-time; and (3) to develop a multi-physics joint inversion approach to precisely predict flow patterns and physiochemical changes within an EAP-filled fracture network. This research project was built upon our previous work at the Devine Test Site managed by the Bureau of Economic Geology (BEG) at The University of Texas at Austin (UT-Austin). It also leveraged a significant investment from the Advanced Energy Consortium (AEC) to address the DOE's interest in subsurface flow, containment, and characterization by multiple signals. This three-year and three-month project succeeded in demonstrating the feasibility of a real-time dynamic fluid flow mapping technique at Technology Readiness Level 5 (TRL5) by utilizing a commercially available surface-based Controlled-Source Electromagnetic (CSEM) method (Objectives 1, 2). We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, the observed electric field in our study is influenced by bottomhole pressure, flow rate, and salinity, which is demonstrated by comparing these parameters with the electrical field potentials. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce the whole measured electric field magnitudes. Preliminary estimates showed that including Streaming Potential (SP) in our geophysical model is likely needed to reduce the simulation misfit.

02 PETROLEUM↗
DOE OSTI2022

CO 2 methanation reaction pathways over unpromoted and NaNO 3 -promoted Ru/Al 2 O 3 catalysts

Catalytic CO 2 sorbents, materials that adsorb and pre-concentrate CO 2 on the catalyst surface prior to subsequent conversion, are becoming important materials in CO 2 capture and utilization. In this work, a prototypical CO 2 methanation catalyst – Ru/Al 2 O 3 – and a related catalytic sorbent – NaNO 3 /Ru/Al 2 O 3 – are used for CO 2 methanation in flowing hydrogen in a fixed bed reactor at temperatures ranging from 220 to 280 °C. Activation energies for the NaNO 3 /Ru/Al 2 O 3 material are slightly higher than unpromoted Ru/Al 2 O 3 catalysts, and the reaction orders vary more significantly. In situ IR spectroscopy and steady-state isotopic kinetic analysis (SSITKA) using in situ IR/MS spectroscopy show that bicarbonate and linear carbonyl species are the likely reaction intermediates over unpromoted Ru/Al 2 O 3 , while bidentate carbonate, formate and linear carbonyl species are among likely reaction intermediates over NaNO 3 /Ru/Al 2 O 3 . Rate laws consistent with the obtained experimental data are proposed after kinetic modeling of multiple plausible reaction pathways. In conclusion, results suggest that the pathway over the NaNO 3 /Ru/Al 2 O 3 catalyst likely has an additional kinetically relevant irreversible step in the CO 2 methanation reaction pathway.

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

Active Learning Guided Computational Discovery of Plant-Based Redoxmers for Organic Nonaqueous Redox Flow Batteries

Organic nonaqueous redox flow batteries (ONRFBs) are promising energy storage devices due to their scalability and reliance on sourceable materials. However, finding suitable redox-active organic molecules (redoxmers) for these batteries remains a challenge. Using plant-based compounds as precursors for these redoxmers can decrease their costs and environmental toxicity. In this computational study, flavonoid molecules have been examined as potential redoxmers for ONRFBs. Flavone and isoflavone derivatives were selected as catholyte (positive charge carrier) and anolyte (negative charge carrier) molecules, respectively. To drive their redox potentials to the opposite extremes, in silico derivatization was performed using a novel algorithm to generate a library of > 40000 candidate molecules that penalizes overly complex structures. A multiobjective Bayesian optimization based active learning algorithm was then used to identify best redoxmer candidates in these search spaces. Furthermore, our study provides methodologies for molecular design and optimization of natural scaffolds and highlights the need of incorporating expert chemistry awareness of the natural products and the basic rules of synthetic chemistry in machine learning.

25 ENERGY STORAGE↗
DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
DOE OSTI2023

Comparative analysis of sublimation and thermal decomposition of TATB

This experimental study investigated the effects of confinement, starting mass, and heating rate on TATB thermal decomposition and sublimation using a combined Thermo-Gravimetric Analyzer and Differential Scanning Calorimetry (TGA/DSC) instrument. The confinement of volatile products was varied using different pinhole sizes with TGA/DSC pans. The measurements showed the open pan experiments without lids/pinholes resulted in complete sublimation of TATB between 320 °C and 360 °C. The heat of sublimation was determined to be 176 kJ/mol (42 kcal /mol), consistent with literature data obtained from other experimental techniques. The use of pinholes suppressed the sublimation of TATB such that the decrease in pinhole size resulted in 1) an increase in the enthalpy of reaction and an increase in the amount of carbonaceous material remaining at the end of decomposition, and 2) convergence of the two peak temperatures corresponding to maximum heat flow and maximum weight loss. Also, a transition from a two-exotherm thermal decomposition behavior towards a single-exotherm occurred as the pinhole size was decreased for a given starting mass or as the starting mass was increased for a given pinhole size. These results indicate the kinetics of TATB sublimation, TATB thermal decomposition, and gas diffusion out of a TGA/DSC pan can all compete and result in significantly different enthalpies, amounts of remaining materials, and peak temperatures depending on the pinhole size and starting mass used in the measurements. Furthermore, the results also indicate precise control of process variables (pinhole size, starting mass, and heating rate) in TGA/DSC measurements is required for thermal safety assessment of explosives.

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

Deep level defects in low-pressure chemical vapor deposition grown (010) β-Ga 2 O 3

This study provides the full-bandgap evaluation of defect state distributions in beta phase gallium oxide (β-Ga 2 O 3 ) grown by low-pressure chemical vapor deposition (LPCVD) on (010) β-Ga 2 O 3 substrates at high growth of up to 20 µm/h. Deep-level optical spectroscopy and deep-level transient spectroscopy measurements applied to Ni/β-Ga 2 O 3 Schottky diodes revealed the presence of a previously unreported defect state at E C -3.6 eV, which dominated the overall trap distribution in LPCVD grown material. However, states at E C -0.8 eV, E C -2.0 eV, and E C -4.4. eV were also detected, similar to prior studies on β-Ga 2 O 3 grown by other methods, with similar or lower concentrations for the LPCVD samples. The E C -0.8 eV and E C -2.0 eV states were previously connected to residual Fe impurities and gallium vacancies, respectively. The total concentration of traps in the LPCVD material was on par with or lower than the state-of-the-art metal–organic chemical vapor deposition-grown materials despite the much higher growth rate, and the distribution of states showed negligible dependence on SiCl 4 flow rate and doping concentration. These results demonstrate that the high growth rate of LPCVD-grown β-Ga 2 O 3 is very promising for achieving thick, low defect density, and high-quality layers needed for multi-kV device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Resin Testing and Modeling for Optimal Composite Processing

Polymer composites have properties such as high strength and stiffness, low weight, good thermal and chemical stability, as well as impact and abrasion resistance that make them ideal for high-performance applications. The chemistries of these materials are continuously improving, so determining their properties is vital for successfully producing them and achieving the desired results. Multiple methods can be employed to monitor characteristics such as heat flow, weight, dimension, and modulus as a function of time and temperature. By analyzing this information, models can be developed to predict outcomes of parameters not tested for. In one application, materials proposed for wet filament winding and the production of high pressure vessels can be analyzed to verify they will have the necessary low viscosity for good fiber wetting and long pot life for the extended handling inherent to this process. Such data about a prospective system provides valuable information on how that material could ultimately be processed to yield the desired part.

36 MATERIALS SCIENCE↗
DOE OSTI2020

X-ray and gamma-ray tomographic imaging of fuel relocation inside sodium fast reactor test assemblies during severe accidents

The present work reports on x-ray and gamma-ray high-spatial resolution computerized tomography measurements of the Pin Bundle Metallic Fuel Relocation (PBR) assemblies tested in the Metallic Uranium Safety Experiment (MUSE) facility at Argonne National Lab (ANL). The aim of the study was to characterize fuel relocation structures that develop during severe core accidents pertaining to SFR assemblies; these include but are not limited to advanced core disruption recreated in the PBR-1 assembly, and cladding breach recreated in the PBR-2 assembly. We report the x-ray tomography measurements were able to resolve small quantities of relocation fuel; with increased presence of relocation fuel, the x-ray 1measurements spatially mapped the material but could not resolve the inner regions of these. The gamma-tomography measurements showed improved results, resolving the relocation structures in great detail. The upper plenum of the PBR-1 assembly where the molten uranium was initially inserted presented high structural damage, reflected by the partial and complete disintegration of the central rods. Relocation fuel filled the subchannels, adhering to surviving cladding walls and the assembly casing. In the lower portion of the measured section, the tomogram degrades due to photon starvation effects hinting at the increased amount of relocation fuel potentially plugging the assembly; flow blockage in this section was difficult to determine due to the tomogram’s degradation from photon starvation. Small fragments were observed further down the assembly, dislodged from the initial insertion of the molten material. This section was used as an unperturbed assembly reference, with a calculated blockage of less than 1% from the present fragments. The PBR-2 assembly was characterized by columnar relocation structures propagating through the subchannels. Three of the relocation structures were captured in the measured section, with evidence of cross migration on to adjacent subchannels. The measured section in this assembly captures the leading edge of two structures. The calculated flow blockage was 16% in the planes where the three relocation structures are present, but this quickly decreases to approximately 5% past the leading edge of two of the structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Strategically Modified Ligand Incorporating Mixed Phosphonate and Carboxylate Groups to Enhance Performance in All-Iron Redox Flow Batteries

Iron redox flow batteries (Fe-RFBs) hold significant promise for achieving cost-effectiveness and utilizing abundant materials for stationary energy storage applications. Here, a design of a novel Fe complex utilizing a nitrogenous phosphonate/carboxylate mixed ligand, N,N-Bis(phosphonomethyl)glycine (BMPG), is presented to achieve high performance Fe anolyte. Compared to its all-phosphonate form, nitrilotri(methylphosphonic acid) (NTMPA), the new complex Fe(BPMG) 2 demonstrates a negatively shifted redox potential, resulting in ≈0.07 V (≈10%) increase in battery output voltage. Full battery testing paired with ferrocyanide catholyte demonstrates stable cycling (capacity degradation <0.0001%/cycle) over 730 consecutive charge/discharge cycles with Coulombic Efficiency of 100% at a current density of 20 mA cm -2 under near neutral pH (≈8). Of particular interest, density functional theory (DFT) studies and operando Raman measurements provide strong evidence supporting a molecular structure in BPMG, which reveals the mixed phosphonate/carboxylate groups in BPMG maintain the octahedral coordination of the Fe ion center with phosphonates exclusively, while leaving the carboxylate unbound for both Fe(II) and Fe(III) complexes. This structural similarity between BPMG-based Fe(II) and Fe(III) complexes effectively mitigates the slow redox reaction kinetics observed in Fe(NTMPA) 2 anolyte, where significant ligand reorientation occurs between Fe(II) and Fe(III) complexes.

25 ENERGY STORAGE↗
DOE OSTI2023

Development of a versatile electrochemical cell for in situ grazing-incidence X-ray diffraction during non-aqueous electrochemical nitrogen reduction

In situ techniques are essential to understanding the behavior of electrocatalysts under operating conditions. When employed, in situ synchrotron grazing-incidence X-ray diffraction (GI-XRD) can provide time-resolved structural information of materials formed at the electrode surface. In situ cells, however, often require epoxy resins to secure electrodes, do not enable electrolyte flow, or exhibit limited chemical compatibility, hindering the study of non-aqueous electrochemical systems. Here, a versatile electrochemical cell for air-free in situ synchrotron GI-XRD during non-aqueous Li-mediated electrochemical N 2 reduction (Li-N 2 R) has been designed. This cell not only fulfills the stringent material requirements necessary to study this system but is also readily extendable to other electrochemical systems. Under conditions relevant to non-aqueous Li-N 2 R, the formation of Li metal, LiOH and Li 2 O as well as a peak consistent with the α-phase of Li 3 N was observed, thus demonstrating the functionality of this cell toward developing a mechanistic understanding of complicated electrochemical systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2021

Bounds on the Rate-Dependent Plastic Flow of Tantalum up to 75 GPa

In this study we report on the continued development of thermodynamics-based analysis of shock waves propagation with the objective of extracting information related to materials strength at high strain rates and pressures. Building on previous results reported for peak stresses of 10 GPa and 25 GPa, we present a series of three-step gas-gun shock experiments designed to explore the pressure and strain rate dependence of plastic flow in polycrystalline tantalum. These experiments at nominal peak stresses of 50 GPa and 75 GPa show the irreversible deformation before pullback to be almost entirely confined to the shock loading, with negligible plastic relaxation on the post-shock plateau. We also add a reverse-ballistics shot at 25 GPa, which was designed to reveal the pullback response with negligible interference from free-surface effects. General thermodynamic considerations allow us to place bounds on the plastic behavior even for parts of the curve that change far too rapidly for the velocimetric time resolution of (conservatively) ~ 5 ns. To analyze the data, we found it necessary to substantially improve the interpolation/extrapolation scheme in order to improve its robustness, flexibility and range of applicability. We describe the new scheme based on splines, as well as an extension of free-surface corrections to the post-shock rarefaction waves. Reanalysis with the new scheme produces results essentially within the error bars previously reported, showing that the known systematic errors associated with free-surface effects are relatively inconsequential for determining thermodynamic paths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elucidating lipid nanoparticle properties and structure through biophysical analyses

Designing lipid nanoparticle (LNP) delivery systems with specific targeting, potency and minimal side effects is crucial for their clinical use. However, traditional characterization methods, such as dynamic light scattering, cannot accurately quantify physicochemical properties of LNPs and how these are influenced by the lipid composition and mixing method. Here, we structurally characterize polydisperse LNP formulations by applying emerging solution-based biophysical methods that have higher resolution and provide biophysical data beyond size and polydispersity. These techniques include sedimentation velocity analytical ultracentrifugation, field-flow fractionation followed by multiangle light scattering and size-exclusion chromatography in line with synchrotron small-angle X-ray scattering. Here, we show that LNPs have intrinsic polydispersity in size, RNA loading and shape, which depend on both the formulation technique and the lipid composition. Lastly, we predict LNP transfection in vitro and in vivo by examining the relationship between mRNA translation and physicochemical characteristics. Solution-based biophysical methods will be essential for determining LNP structure–function relationships, facilitating the creation of new design rules for LNPs.

36 MATERIALS SCIENCE↗