Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D transport”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Application of laminar flow control to the High Speed Civil Transport - The NASA Supersonic Laminar Flow Control Program

The NASA Supersonic Laminar Flow Control (SLFC) program encompasses the development of refined CFD methods and boundary layer stability codes for the highly 3D supersonic flow conditions encountered by the F-16XL technology demonstration aircraft and the prospective High Speed Civil Transport (HSCT). While the F-16XL-1 aircraft continues to gather SLFC data, work is under way on the F-16XL-2 aircraft: which will furnish attach-line design criteria, code-calibration data, and an improved understanding of the flowfield over a wing that will add confidence to the design of HSCTs' boundary layer-controlling air-suction panels.

Fischer, Michael C.↗
DOE OSTI2023

Diffusion–convection model of runaway electrons due to large magnetohydrodynamic perturbations in post-thermal quench plasmas

Systematic test particle tracing simulations for runaway electrons (REs) are performed for six post-thermal quench equilibria from DIII-D and ITER, where large scale, kink-like n = 1 (n is the toroidal mode number) magnetohydrodynamic (MHD) instabilities are found. The modeled particle guiding center orbits allow extraction of the effective diffusion–convection coefficients of REs in the presence of large three-dimensional (3D) perturbations up to 10% of the equilibrium toroidal field. With a fixed spatial distribution of the field perturbation, the RE transport coefficients along the plasma radial coordinate track reasonably well with the surface-averaged perturbation level. A substantial variation in the value of the transport coefficients—by three orders of magnitude in most cases, however, occurs with varying launching location of REs along the plasma radius. Large 3D perturbations almost always lead to comparable diffusion and convection processes, meaning that diffusion alone is insufficient to describe the particle motion. At lower (but still high) level of perturbation, the RE convection is found to be dominant over diffusion. A similar observation is made when the perturbation is too strong. In the presence of large perturbation, the dependence of the RE transport on the particle energy is sensitive to the spatial distribution of the perturbation. Based on numerically obtained RE transport coefficients, an analytic fitting model is proposed to quantify the particle diffusion and convection processes due to large MHD events in post-thermal quench plasmas. The model is shown to reasonably well reproduce the direct test particle tracing results for the RE loss fraction and can, thus, be useful for incorporating into other kinetic RE codes in order to simulate the RE beam evolution in the presence of large 3D perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Collisionless plasma transport mechanisms in stochastic open magnetic field lines in tokamaks

The collisionless plasma transport in stochastic magnetic fields connecting to wall boundary has been studied to understand the underlying mechanisms of thermal quench physics in tokamak disruption. In this study, we present a comprehensive understanding of the plasma dynamics in three-dimensional stochastic open magnetic field lines, taking into account the consistent coupling of electron and ion dynamics through the ambipolar electric field. The open stochastic field lines act as a 3D magnetic mirror consisting of the magnetic well and hill regions. It was found that the magnetic hill regions play a critical role in enhancing electron thermal transport. When the plasma collapses to the wall, the 3D ambipolar potential arises in the stochastic layer to maintain the quasi-neutrality of the plasma. The E × B vortices induced by the 3D ambipolar potential mix the plasma across stochastic field lines and enhance the radial transport. Particularly, the E × B mixing between the magnetic well and hill regions provides a collisionless detrapping mechanism that plays a major role in the loss of high-energy trapped electrons. As a result, the electron temperature steadily decreases at a rate of about -0.5 keV ms -1 , comparable to the rate observed in thermal quench experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗
DOE OSTI2022

Entropy Stabilization Effects and Ion Migration in 3D “Hollow” Halide Perovskites

A recently discovered new family of 3D halide perovskites with the general formula (A) 1–x (en) x (Pb) 1–0.7x (X) 3–0.4x (A = MA, FA; X = Br, I; MA = methylammonium, FA = formamidinium, en = ethylenediammonium) is referred to as “hollow” perovskites owing to extensive Pb and X vacancies created on incorporation of en cations in the 3D network. The “hollow” motif allows fine tuning of optical, electronic, and transport properties and bestowing good environmental stability proportional to en loading. To shed light on the origin of the apparent stability of these materials, we performed detailed thermochemical studies, using room temperature solution calorimetry combined with density functional theory simulations on three different families of “hollow” perovskites namely en/FAPbI 3 , en/MAPbI 3 , and en/FAPbBr 3 . We found that the bromide perovskites are more energetically stable compared to iodide perovskites in the FA-based hollow compounds, as shown by the measured enthalpies of formation and the calculated formation energies. The least stable FAPbI 3 gains stability on incorporation of the en cation, whereas FAPbBr 3 becomes less stable with en loading. This behavior is attributed to the difference in the 3D cage size in the bromide and iodide perovskites. Configurational entropy, which arises from randomly distributed cation and anion vacancies, plays a significant role in stabilizing these “hollow” perovskite structures despite small differences in their formation enthalpies. With the increased vacancy defect population, we have also examined halide ion migration in the FA-based “hollow” perovskites and found that the migration energy barriers become smaller with the increasing en content.

36 MATERIALS SCIENCE↗
NASA NTRS2017

A Dose of Reality: Radiation Analysis for Realistic Human Spacecraft

INTRODUCTION As with most computational analyses, a tradeoff exists between problem complexity, resource availability and response accuracy when modeling radiation transport from the source to a detector. The largest amount of analyst time for setting up an analysis is often spent ensuring that any simplifications made have minimal impact on the results. The vehicle shield geometry of interest is typically simplified from the original CAD design in order to reduce computation time, but this simplification requires the analyst to "re-draw" the geometry with a limited set of volumes in order to accommodate a specific radiation transport software package. The resulting low-fidelity geometry model cannot be shared with or compared to other radiation transport software packages, and the process can be error prone with increased model complexity. The work presented here demonstrates the use of the DAGMC (Direct Accelerated Geometry for Monte Carlo) Toolkit from the University of Wisconsin, to model the impacts of several space radiation sources on a CAD drawing of the US Lab module. METHODS The DAGMC toolkit workflow begins with the export of an existing CAD geometry from the native CAD to the ACIS format. The ACIS format file is then cleaned using SpaceClaim to remove small holes and component overlaps. Metadata is then assigned to the cleaned geometry file using CUBIT/Trelis from csimsoft (Registered Trademark). The DAGMC plugin script removes duplicate shared surfaces, facets the geometry to a specified tolerance, and ensures that the faceted geometry is water tight. This step also writes the material and scoring information to a standard input file format that the analyst can alter as desired prior to running the radiation transport program. The scoring results can be transformed, via python script, into a 3D format that is viewable in a standard graphics program. RESULTS The CAD model of the US Lab module of the International Space Station, inclusive of all the racks and components, was simplified to remove holes and volume overlaps. Problematic features within the drawing were also removed or repaired to prevent runtime issues. The cleaned drawing was then run through the DAGMC workflow to prepare for analysis. Pilot tests modeling transport of 1GeV proton and 800MeV/A oxygen sources show that reasonable results are converged upon in an acceptable amount of overall computation time from drawing preparation to data analysis. The FLUKA radiation transport code will next be used to model both a GCR and a trapped radiation source. These results will then be compared with measurements that have been made by the radiation instrumentation deployed inside the US Lab module. DISCUSSION Early analyses have indicated that the DAGMC workflow is a promising toolkit for running vehicle geometries of interest to NASA through multiple radiation transport codes. In addition, recent work has shown that a realistic human phantom, provided via a subcontract with the University of Florida, can be placed inside any vehicle geometry for a combinatorial analysis. This added functionality gives the user the ability to score various parameters at the organ level, and the results can then be used as input for cancer risk models.

Barzilla, J. E.↗
DOE OSTI2024

3D Printing of Functional Hydrogel Devices for Screenings of Membrane Permeability and Selectivity

Developing a fundamental understanding of the effects of varying ligand chemistries on mass transport rates is key to designing membranes with solute-specific selectivity. While permeation cells offer a robust method to characterize membrane performance, they are limited to assessing a single membrane chemistry or salt solution per test. As a result, investigating the effects of varying ligand chemistries on membrane performance can be a tedious process, involving both the preparation of multiple samples and numerous, time-consuming permeation tests. This study uses digital light processing (DLP) 3D printing to fabricate a millifluidic flow-based permeation device made from a hydrogel active ester network that can be easily functionalized with ion-selective ligands. Without the need for bonding or assembly steps, ligands can be introduced and tested in the permeation device by simply injecting a small volume of a ligand solution. Various salt concentrations and molecular species can be cycled through a single device by switching the solution feeding into the salt reservoir, thereby reducing the number of samples needed for permeability and selectivity screenings. This research sets the groundwork for formulation development and postprocessing methods to 3D-print functional millifluidic devices capable of assessing solute selectivity in membranes and polymer adsorbents for aqueous separations. In this work, comparable salt permeability trends were observed with both 3D-printed devices and traditional assays. Devices were functionalized with an imidazole ligand to investigate salt permeability and selectivity of monovalent and divalent salts. Measurements showed increasing permeability for monovalent salts (NaCl) relative to divalent salts (MgCl 2 , CuCl 2 ) in functionalized membranes, with higher monovalent/divalent selectivity at increasing imidazole grafting densities. Here, the methods and findings described here represent a step toward developing higher-throughput methods with 3D-printed devices for screening the effects of ligand chemistry on mass transport rates in membrane materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Increasing the Pressure‐Free Stripping Capacity of the Lithium Metal Anode in Solid‐State‐Batteries by Carbon Nanotubes

Abstract Lithium metal is the key anode material for solid‐state‐batteries as its successful implementation will drastically increase their energy and power densities. However, anode contact loss during stripping leads to dendrites upon plating and subsequent cell failure. Design strategies to mitigate these issues are crucial to enable the use of lithium metal anodes. This paper reports the dissolution kinetics of composite anodes made of lithium metal and carbon nanotubes (CNTs) with a garnet‐type solid electrolyte (SE). In addition to an enhancement of the effective diffusion within the anode, its dissolution is fundamentally changed from being 2D to 3D. By maintaining contact with the SE, the CNTs facilitate lithium transport to the interface, which yields more than 20 mAh cm −2 discharge capacity at 100 µ A cm −2 without the application of external stack pressure ( > 1 MPa). Conclusions drawn from electrochemical data on the anode microstructure are validated using cryo‐focused‐ion‐beam scanning electron microscopy and correlated with the mechanical properties. Micro‐indentation, acoustic analysis, and stress–strain testing show that mechanical properties of the anode, like yield strength and hardness, are adjustable. Overall, it is shown that the mechanical and electrochemical properties of Li–CNT composite electrodes can be tailored to suit the requirements of a practical cell.

Fuchs, Till↗
DOE OSTI2024

Enhancement of Conjugated Polymer Microstructure and Mixed-Conducting Properties via Chalcogenophene Heteroatom Substitution

Heteroatom substitution is a powerful tool to tune the intra- and intermolecular structure of conjugated polymers as well as their resulting optoelectronic and electrochemical properties. Here, a series of oligoethylene glycol bithiophene chalcogenophene polymers (p(g3T2-X)) with systematically varied furan, thiophene, selenophene, and tellurophene comonomers have been synthesized for mixed ionic-electronic conducting applications. Their microstructures have been thoroughly characterized ex situ and in situ with X-ray scattering, and their mixed conducting properties have been probed in electrochemical transistor testbeds. Chalcogenophene heteroatom choice was found to clearly dictate the polymer microstructure (crystallite dimensionality and orientation) and tune mixed conducting properties. Proceeding down Group 16, from O to Se systematically directed the molecular ordering of 2D polymer crystallites from face-on (O) to mixed (S) to edge-on (Se) orientations, with Te driving the polymer to form well-oriented edge-on 3D crystallites. Heteroatom dictated crystallite quality, and orientation tuned relative ionic transport by 2 orders of magnitude. Hole mobility (μ hole ) and mixed conducting figure of merit (μC*) were each tuned over an order of magnitude depending on heteroatom choice, with the Te-containing polymer reaching μ hole = 3.60 cm 2 V –1 s –1 and μC* = 483 F cm –1 V –1 s –1 , due to improved molecular ordering. Insights from this polymer series highlight target microstructures for enhanced mixed conduction in future conjugated polymers.

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

Xylem network connectivity and embolism spread in grapevine( Vitis vinifera L.)

Abstract Xylem networks are vulnerable to the formation and spread of gas embolisms that reduce water transport. Embolisms spread through interconduit pits, but the three-dimensional (3D) complexity and scale of xylem networks means that the functional implications of intervessel connections are not well understood. Here, xylem networks of grapevine (Vitis vinifera L.) were reconstructed from 3D high-resolution X-ray micro-computed tomography (microCT) images. Xylem network performance was then modeled to simulate loss of hydraulic conductivity under increasingly negative xylem sap pressure simulating drought stress conditions. We also considered the sensitivity of xylem network performance to changes in key network parameters. We found that the mean pit area per intervessel connection was constant across 10 networks from three, 1.5-m stem segments, but short (0.5 cm) segments fail to capture complete network connectivity. Simulations showed that network organization imparted additional resistance to embolism spread beyond the air-seeding threshold of pit membranes. Xylem network vulnerability to embolism spread was most sensitive to variation in the number and location of vessels that were initially embolized and pit membrane vulnerability. Our results show that xylem network organization can increase stem resistance to embolism spread by 40% (0.66 MPa) and challenge the notion that a single embolism can spread rapidly throughout an entire xylem network.

Wason, Jay↗
DOE OSTI2024

Anomalous thermal effect in Zr Te 5 observed via photothermal measurements

In this study, we explore the magneto-thermoelectric power (MTP) of $ZrTe$ 5 , a canonical Dirac semimetal, through a novel photothermal technique. Unlike conventional thermoelectric studies that rely on on-chip heaters and are limited by fabrication processes, especially for stress-sensitive materials, our approach utilizes photothermal effects to induce temperature gradients. Our experiments, applying a magnetic field approximately parallel and transverse to the photocurrent detection direction, reveal that the photothermal method efficiently and reliably extracts both diagonal and off-diagonal components of the thermoelectric coefficient of $ZrTe$ 5 . Here, we observe that the longitudinal MTP reproduces features previously reported in thermal transport studies, while the photoinduced transverse MTP confirms the anomalous Nernst effect. This photothermal measurement technique opens new avenues for investigating transport properties in a wide range of quantum materials, both in 3D and 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analysis of Approximations in Modeling of BWR Bundle Void Distributions

In boiling water reactors, complex heterogeneous bundle designs, control blades adjacent to the corner of bundles, and the presence of boiling can lead to complex internal void distributions. A few approximations exist to model these void distributions. They could be modeled using a 1D axial solver in which each axial node is assumed to be at an average void, or each pin cell could be modeled with its own void concentration. In the latter case, the void could be discretized in pin-centered or coolant-centered channels. The goal of this project was to quantify the effect of using the different approximations for modeling internal void distributions on neutronics calculations. Using 3D void distributions calculated with CTF, Monte Carlo Neutral Particle (MCNP) transport code models were created for GE-9 and GE-14 lattices. For each model, the internal void distribution from CTF at a given axial node was selected, and a lattice calculation was carried out with MCNP. Comparisons between models using a lattice-averaged void, or using a void distribution in coolant-centered channels, showed large differences in reactivity which in some cases were well above 1,000 pcm, and it also showed differences in normalized fission rates greater than 20%. It was also found that using a lattice average void can lead to a significant difference in the estimation of the worth of a control blade. The differences found when comparing results from models using pin-centered and coolant-centered channels were up to 200 pcm in reactivity and up to 1.4% in the normalized fission rates. In addition to these two sets of comparisons, MCNP models were set up so that each subchannel had a saturated liquid component around the fuel pins and a saturated vapor component in the center to approximate annular flow. In comparison to the models using coolant-centered subchannels, up to 1–3% differences in normalized fission rates could be found.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2001

Large-Scale Stratospheric Transport Processes

The paper discusses the following: 1. The Brewer-Dobson circulation: tropical upwelling. 2. Mixing into polar vortices. 3. The latitudinal structure of "age" in the stratosphere. 4. The subtropical "tracer edges". 5. Transport in the lower troposphere. 6. Tracer modeling during SOLVE. 7. 3D modeling of "mean age". 8. Models and measurements II.

Plumb, R. Alan↗
NASA NTRS1979

Tank Insulation

For NASA's Apollo program, McDonnell Douglas Astronautics Company, Huntington Beach, California, developed and built the S-IVB, uppermost stage of the three-stage Saturn V moonbooster. An important part of the development task was fabrication of a tank to contain liquid hydrogen fuel for the stage's rocket engine. The liquid hydrogen had to be contained at the supercold temperature of 423 degrees below zero Fahrenheit. The tank had to be perfectly insulated to keep engine or solar heat from reaching the fuel; if the hydrogen were permitted to warm up, it would have boiled off, or converted to gaseous form, reducing the amount of fuel available to the engine. McDonnell Douglas' answer was a supereffective insulation called 3D, which consisted of a one-inch thickness of polyurethane foam reinforced in three dimensions with fiberglass threads. Over a 13-year development and construction period, the company built 30 tanks and never experienced a failure. Now, after years of additional development, an advanced version of 3D is finding application as part of a containment system for transporting Liquefied Natural Gas (LNG) by ship.

Source record↗
DOE OSTI2025

Microstructure-based modeling of inner oxygen pressure in solid oxide electrolysis cells: Analysis of electrode delamination and mitigation

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here a microstructure-based electrochemical model, which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

25 ENERGY STORAGE↗