Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion couples”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hydrodynamic spin-orbit coupling in asynchronous optically driven micro-rotors

Abstract Vortical flows of rotating particles describe interactions ranging from molecular machines to atmospheric dynamics. Yet to date, direct observation of the hydrodynamic coupling between artificial micro-rotors has been restricted by the details of the chosen drive, either through synchronization (using external magnetic fields) or confinement (using optical tweezers). Here we present a new active system that illuminates the interplay of rotation and translation in free rotors. We develop a non-tweezing circularly polarized beam that simultaneously rotates hundreds of silica-coated birefringent colloids. The particles rotate asynchronously in the optical torque field while freely diffusing in the plane. We observe that neighboring particles orbit each other with an angular velocity that depends on their spins. We derive an analytical model in the Stokes limit for pairs of spheres that quantitatively explains the observed dynamics. We then find that the geometrical nature of the low Reynolds fluid flow results in a universal hydrodynamic spin-orbit coupling. Our findings are of significance for the understanding and development of far-from-equilibrium materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Electrochemical properties of poly(ethylene oxide) electrolytes above the entanglement threshold

The ion transport in electrolytes depends on three transport coefficients, conductivity (κ), salt diffusion coefficient (D), and the cation transference number with respect to the solvent velocity ($t_+^0$), and the thermodynamic factor ($T_f$). Current methods for determining these parameters involve four separate experiments, and the coupled nature of the equations used to determine them generally results in large experimental uncertainty. We present data obtained from 64 independent polymer electrolytes comprising poly(ethylene oxide) (PEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt. The molecular weights of PEO ranged from 5 to 275 kg mol -1 ; these samples are all above the entanglement threshold. We minimize the experimental uncertainty in transport and thermodynamic measurements by exploiting the fact that ion transport in entangled polymer electrolytes should be independent of molecular weight. The dependence of $κ, D, t_+^0$, and $T_f$ as a function of salt concentration in the range 0.035 ≤ r ≤ 0.30 are presented with a 95% confidence interval, where r is the molar ratio of lithium ions to ethylene oxide monomer units. While κ, D, and $T_f$ are all positive as required by thermodynamic constraints, there is no constraint on the sign of $t_+^0$. We find that $t_+^0$ is negative in the salt concentration range of 0.093 ≤ r ≤ 0.189.

25 ENERGY STORAGE↗
DOE OSTI2020

Stability of a tokamak plasma with diffuse toroidal rotation

The present work considers the stability of a high- $\beta$ , large aspect ratio, circular plasma with diffuse profiles for the safety factor and the angular toroidal frequency (López & Guazzotto, Phys. Plasmas , vol. 24, 032501). An application of the Frieman–Rotenberg formalism results in a system of scalar eigenmode equations whose coupling is retained at the plasma–vacuum transition but is disregarded across the plasma column, which is a standard practice. The solution technique consists of a multidimensional shooting method for the poloidal harmonics; robust initial guesses are constructed by solving the dispersion relation in the static scenario with vanishing magnetic shear. Flow shear appears as a high- $\beta$ toroidal contribution, and we illustrate its destabilizing influence on $n=1$ external kink modes in the presence of ideal and resistive walls. Internal resonances are avoided by means of the selection of appropriate equilibrium parameters. The stabilizing influence of a finite positive average magnetic shear is also exemplified.

Physics↗
DOE OSTI2026

Decomposition of methane diluted with inert gas in an RF discharge cell

Decomposition of methane using non-thermal plasmas is an attractive route for producing hydrogen-rich gases and valuable carbon nanomaterials. Understanding how plasma discharge modes influence methane decomposition in optimizing plasma-assisted chemical conversion remains unexplored. This study explores the coupling between the discharge structure and product selectivity in RF capacitively coupled discharges operating in methane/inert gas mixtures in the pressure range of 2–3 torr. Under our experimental conditions the discharge exhibits mode transitions from uniform to striated in Ar and Kr and from diffuse to contracted in Ar and Kr with <5% CH 4 . The discharges in He and Ne remained uniform under our operating conditions, and their mixtures with CH 4 remained diffuse. A 0-d model for Ar/CH 4 discharge established a threshold for contraction while also asserting the importance of Ar$^{*}_{\textrm{m}}$ in the dissociation and ionization processes. The highest degree of methane decomposition, >99.7% with the main products of acetylene and graphitized solid carbon was achieved in the contracted discharge mode for both Kr or Ar with ⩽5% CH 4 . We demonstrate that contraction can play a crucial role in the effective decomposition of methane with value-added products and that both the electronic and thermal properties of plasma gas are responsible for this effect.

RF plasma↗
DOE OSTI2022

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗
DOE OSTI2024

Mechanochemical topological defects in an active nematic

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±$\frac{1}{2}$ defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. Here, we posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

42 ENGINEERING↗
DOE OSTI2021

Fabrication of Ionic Covalent Triazine Framework-Linked Membranes via a Facile Sol–Gel Approach

Covalent triazine framework (CTF)-based membranes have shown unique properties and wide applications in energy-related fields, but there is still no efficient strategy capable of affording CTF membranes functionalized with ionic moieties, which will bring extra merits and application possibilities. In this work, a robust CTF membrane with ionic functionalities was fabricated via a sol–gel approach promoted by a superacid (FSO 3 H). The CTF skeleton was constructed via the trimerization of cyano groups in the monomer, and piperazine moieties were introduced as reactive sites for the formation of ammonium cations coupled with FSO 3 – anions within the backbone. The obtained transparent and red protonated polymeric CTF-based membrane (PP-CTF) exhibited significant absorption at 646 nm in the solid-state ultraviolet–visible diffuse reflectance spectrum. Comparatively, the neutralized deprotonated counterparts turned to yellow color with a significant blueshift absorption to 492 nm. The theoretical calculation demonstrated that PP-CTF and deprotonated polymeric CTF membrane (DP-CTF) were stable with AB-stacking mode, and PP-CTF had a smaller band gap (0.90 eV) than that of DP-CTF (2.17 eV). The significant optical absorption and emission behavior change of PP-CTF and its robust chemical stability endowed it with the capability to act as pH indicators. The synthetic pathway developed in this work and the performance of the resultant ionic membrane opens new opportunities in the aspects of membrane design, fabrication, and application.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Lightning Induced Interior Fields And Voltage Bounds For Coaxial Topologies

We assemble bounding formulas for the interior fields and pin voltages inside a cylindrical coaxial Faraday cage which has been struck by lightning. Approximate formulas for penetrations through a circumferential door slot with subsequent coupling to the interior center conductor structure. Fields at the opposite open end are estimated and used to drive a capped connector and estimate interior pin voltages. Finally, penetrations through small circular holes and direct diffusion through the barrier are also addressed.

42 ENGINEERING↗
DOE OSTI2020

Modeling Extreme Deformations in Lithium Ion Batteries

A simultaneously coupled modeling approach to study the electrochemical and thermal behavior of lithium-ion batteries under large mechanical deformation has been developed. The thermo-electrochemical pseudo-2D (P2D) battery model is coupled with a mechanical material model. Mechanical, thermal, and electrochemical models are implemented as user-defined sub-routines in the commercial multi-physics code LS-DYNA. The mechanical strain experienced by anode, cathode and separator results in thickness and porosity changes in each layer which in turn influences electrochemical behavior. The evolution of concentration profiles and cell potential are studied under different mechanical loading conditions. Internal short-circuits caused by mechanical deformation and corresponding physical behaviors are also elucidated. We discuss the competing effects of improved transport at higher temperature due to the internal short-circuit versus a drop in the effective ionic conductivity and electrolyte diffusivity due to mechanical deformation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Electrochemical-mechanical coupling failure mechanism of composite cathode in all-solid-state batteries

Composite cathode composed of active particles and solid electrolytes (SEs) can considerably enlarge the particle-SE contact areas and achieve high areal loadings in all-solid-state batteries (ASSBs). However, the challenging interfacial instability and particle damage problems remain unsolved. Herein, we establish a 3D electrochemical-mechanical coupled model to investigate the underlying failure mechanism by considering the governing electrochemical and physics processes. Micro-scale heterogeneous primary particles with random crystallographic orientation and size inside the LiNi 1/3 Co 1/3 Mn 1/3 O 2 (NCM111) secondary particle of the model result in the anisotropic Li diffusion and volume variation within the secondary particle, leading to significant nonuniformity of the Li concentration, and GPa-level stress distributions at primary particle boundaries, and finally causing the particle internal cracks. The particle volume shrinkage under the constraint of stiff Li 7 La 3 Zr 2 O 12 (LLZO) SE triggers the interface debonding (gap>50 nm) with increased interfacial impedance to degrade cell capacity. Higher C-rates result in larger residual stress (~100 MPa)/strain/debonding gap at dis-charging end, more likely to deteriorate the cell performance. Increasing the interfacial strength between the particle and SE can suppress the interface debonding but induces high stress (up to 10 GPa). In conclusion, results reveal the underlying mechanism of the electrochemical-mechanical coupling failure mechanism for composite cathode and provide promising guidance on the further improvement of a more robust composite cathode for ASSBs.

25 ENERGY STORAGE↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗
DOE OSTI2026

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗