Search NASA⌕ Search

SEARCH · Search NASA

Results for “confined diffusion”

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

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of pore size on water dynamics in mesoporous silica

Water confined in mesoporous silica plays a central role in its many uses ranging from gas sorption to nanoconfined chemical reactions. Here in this paper, the influence of pore diameter (2.5–5.4 nm) on water hydrogen bond (H-bond) dynamics in MCM41 and SBA15 mesoporous silicas is investigated using femtosecond infrared vibrational spectroscopy and molecular dynamics simulations on selenocyanate (SeCN-) anions dissolved in the pores. As shown recently, SeCN- spectral diffusion is a reliable probe of surrounding water H-bond structural motions. Additionally, the long CN stretch vibrational lifetime facilitates measurement of the full range of confined dynamics, which are much slower than in bulk water. The simulations shed light on quantitative details that are inaccessible from the spatially averaged observables. The dependence of SeCN- orientational relaxation and that of spectral diffusion on the distance from the silica interface are quantitatively described with an exponential decay and a smoothed step-function, respectively. The distance-dependence of both quantities is found to be independent of the diameter of the pores, and the spatial distribution of SeCN- is markedly non-uniform, reaching a maximum between the interface and the pore center. The results indicate that the commonly invoked two-state, or core–shell, model is a more appropriate description of spectral diffusion. Using these insights, we model the full time-dependence of the measured dynamics for all pore sizes and extract the “core” and “shell” dynamical correlation functions and SeCN- spatial probability distributions. The results are critically compared to those for water confined in reverse micelles.

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

Turbulence-reduced high-performance scenarios in Wendelstein 7-X

In the Wendelstein 7-X (W7-X) stellarator, turbulence is the dominant transport mechanism in most discharges. This leads to a 'clamping' of ion temperature over a wide range of heating power, predominantly flat density profiles where hollow profiles driven by neoclassical thermo-diffusion would be expected and by rapid impurity transport in injection experiments. Significantly reduced turbulent transport is observed in the presence of strong core density gradients found transiently after core pellet injection and irregularly after boronisation or boron pellet injection. Density peaking is also achieved in a controlled manner in purely neutral beam heated discharges where particle transport analysis reveals an abrupt reduction in the main-ion particle flux leading to significant density profile peaking not explained by the NBI particle source alone. The plasmas exhibit a heat diffusivity of around $\chi = 0.25 \pm 0.1\,\mathrm{m}^2\ \mathrm{s}^{-1}$ at mid radius, a factor of around 4 lower than ECRH dominated discharges. Despite the improved confinement, the achieved ion temperature is limited by broader heat deposition and the lower power-per-particle given the higher density. This is overcome with limited reintroduction of ECRH power, where the low heat diffusivity diffusivity is maintained, the density rise supressed and ion temperatures above the clamping limit are achieved. The applicability of these plasmas for a high performance scenario on transport relevant time scales is assessed, including initial predictions for planned heating upgrades of W7-X, based on a range of assumptions about particle transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

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

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

High-Rate, Selective Electrosynthesis of Cyclohexanone Oxime via In Situ Generation and Release of Hydroxylamine on Bismuth

Oxime compounds are key industrial intermediates for nylon precursors and commodity chemicals. However, conventional routes rely on multistep reactions and hydroxylamine (NH 2 OH) salts, raising significant safety and sustainability concerns. Although electrosynthesis offers an alternative, oxime formation on d-block transition metals suffers from poor selectivity, as nitrogen oxyanion intermediates bind strongly to the surface and are readily over-reduced to ammonia. Here, we report morphology-controlled p-block bismuth rhombic dodecahedra (Bi RDs) that promote in situ NH 2 OH generation and its desorption into the electrolyte, enabling an electrochemical-chemical decoupled route for cyclohexanone oxime (CHO) synthesis. Bi RDs deliver nearly 100% Faradaic efficiency (FE) at −0.5 V vs. RHE and a yield of 1.4 mmol h –1 cm –2 at −0.9 V vs. RHE in an H-cell, while maintaining a CHO selectivity of nearly 100% at 100 mA cm –2 in a flow cell. Under identical conditions, d-block electrodes (Cu, Pd, Ag) show FE below 30%. Density functional theory calculations reveal that Bi 6p orbital-derived surface states weaken intermediate binding and facilitate NH 2 OH desorption, suppressing over-reduction. Kinetic analysis, post-addition trapping experiments, and in situ ATR-FTIR and Raman spectroscopy suggest the following reaction mechanism: NH 2 OH is selectively generated at the electrode surface, released as a freely diffusing intermediate, and undergoes homogeneous condensation with cyclohexanone in the bulk electrolyte, bypassing the surface-confined Langmuir–Hinshelwood pathway. These findings demonstrate that regulating intermediate desorption through p-block orbital chemistry provides a general strategy for achieving high selectivity in electro-organic nitrogen synthesis.

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

Electrical control of spin relaxation anisotropy during drift transport in a two-dimensional electron gas

Spin relaxation was studied in a two-dimensional electron gas confined in a wide GaAs quantum well. Recently, the control of the spin relaxation anisotropy by diffusive motion was first shown in D. Iizasa et al., arXiv:2006.08253. Here, we demonstrate electrical control by drift transport in a system with two subbands occupied. The combined effect of in-plane and gate voltages was investigated using time-resolved Kerr rotation. Furthermore, the measured relaxation time presents strong anisotropy with respect to the transport direction. For an in-plane accelerating electric field along [110], the lifetime was strongly suppressed irrespective of the applied gate voltage. Remarkably, for transport along [$1\bara{1}0$], the data shows spin lifetime that was gate dependent and longer than in the [110] direction regardless of the in-plane voltage. In agreement, independent results of anisotropic spin precession frequencies are also presented. Nevertheless, the long spin lifetime, strong anisotropy and drift response seen in the data are beyond the existing models for spin drift and diffusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

On self-preserving, variable-density, turbulent free jets

Published experimental data on incompressible, compresible, free binary, and confined binary turbulent axisymmetric jet flows are compiled and characterized, and the effect of varying turbulent diffusivity across the mixing region of a free jet is investigated analytically, applying the similarity-solution approach of So and Hwang (1986) to the self-preserving region. It is shown that closed-form solutions, represented by Gaussian error functions and having the turbulent Reynolds number and a profile-shape factor as free parameters, can be obtained if the turbulent diffusivities of momentum, mass, or heat are assumed to be different and to vary in both the streamwise and radial directions. An entrainment function uniquely related to the turbulent Reynolds number is derived, and good agreement between theoretical predictions and experimental measurements is demonstrated in graphs.

So, R. M. C.↗
NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗
NASA NTRS2004

An Evolving Trio of Hybrid Stars: C 111

Our goal is to understand the behavior of the outer atmosphere in this intermediate stage to create a comprehensive picture of atmospheric evolution. In the hybrid phase, the large-scale magnetic dynamo activity decays and hydrodynamic processes assume importance. Some hot plasma is still confined close to the star by magnetic loops, yet the confining field is breaking open, the atmosphere can escape through these open field lines, and the diffuse corona may be warm. There may well be a more extended and variable transition process. It remains for FUSE to identify the controlling parameters of the hybrid stars. It shows the positions of our 3 targets in the color-magnitude diagram where it is seen that they are at the extreme end of the hybrid region. Originally we had been awarded the hybrid star Iota Aur, but due to newly imposed pointing constraints of FUSE, that target was not accessible. And so we substituted Iota Dra, a giant of mass similar to our other targets but less evolved. In addition, Iota Dra was recently found to harbor a sub-stellar objects, possibly a planet, and so it could reveal the stellar environment of the planet. This substitution was accepted.

Oliversen, Ronald↗
DOE OSTI2026

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗
DOE OSTI2023

Disturbed, diffuse, or just missing? A global study of the H I content of Hickson compact groups

Hickson compact groups (HCGs) are dense configurations of four to ten galaxies, whose H I morphology appears to follow an evolutionary sequence of three phases, with gas initially confined to galaxies, then significant amounts spread throughout the intra-group medium, and finally with almost no gas remaining in the galaxies themselves. It has also been suggested that several groups may harbour a diffuse H I component that is resolved out by interferometric observations. The H I deficiency of HCGs is expected to increase as the H I morphological phase progresses along the evolutionary sequence. If this is the case, H I deficiency would be a rough proxy for the age and evolutionary state of a HCG. We aim to test this hypothesis for the first time using a large sample of HCGs and to investigate the evidence for diffuse H I in HCGs. We performed a uniform reduction of all publicly available VLA H I observations (38 HCGs) with a purpose-built pipeline that also maximises the reproducibility of this study. The resulting H I data cubes were then analysed with the latest software tools to perform a manual separation of emission features into those belonging to galaxies and those extending into the intra-group medium. We thereby classified the H I morphological phase of each group as well as quantified their H I deficiency compared to galaxies in isolation. We find little evidence that H I deficiency can be used as a proxy for the evolutionary phase of a compact group in either of the first two phases, with the distribution of H I deficiency being consistent in both. However, for the final phase, the distribution clearly shifts to high H I deficiencies, with more than 90% of the expected H I content typically missing. Across all HCGs studied, we identify a few cases where there is strong evidence for a diffuse gas component in the intra-group medium, which might be detectable with improved observations. We also classify a new sub-phase where groups contain a lone H I -bearing galaxy, but are otherwise devoid of gas. The new morphological phase we have identified is likely the result of an evolved, gas-poor group acquiring a new, gas-rich member. The large spread of H I deficiencies in the first two morphological phases suggests that there is a broad range of initial H I content in HCGs, which is perhaps influenced by large-scale environment, and that the timescale for morphological changes is, in general, considerably shorter than the timescale for the destruction or consumption of neutral gas in these systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Confined Ni-In intermetallic alloy nanocatalyst with excellent coking resistance for methane dry reforming

Carbon dioxide and methane are two main greenhouse gases which are contributed to serious global warming. Fortunately, dry reforming of methane (DRM), a very important reaction developed decades ago, can convert these two major greenhouse gases into value-added syngas or hydrogen. The main problem retarding its industrialization is the seriously coking formation upon the nickel-based catalysts. Herein, a series of confined indium-nickel (In-Ni) intermetallic alloy nanocatalysts (In x Ni@SiO 2 ) have been prepared and displayed superior coking resistance for DRM reaction. The sample containing 0.5 wt.% of In loading (In 0.5 Ni@SiO 2 ) shows the best balance of carbon deposition resistance and DRM reactivity even after 430 h long term stability test. The boosted carbon resistance can be ascribed to the confinement of core–shell structure and to the transfer of electrons from Indium to Nickel in In-Ni intermetallic alloys due to the smaller electronegativity of In. Additionally, both the silica shell and the increase of electron cloud density on metallic Ni can weaken the ability of Ni to activate C–H bond and decrease the deep cracking process of methane. The reaction over the confined InNi intermetallic alloy nanocatalyst was conformed to the Langmuir-Hinshelwood (L-H) mechanism revealed by in situ diffuse reflectance infrared Fourier transform spectroscopy (in-situ DRIFTS). This work provides a guidance to design high performance coking resistance catalysts for methane dry reforming to efficiently utilize these two main greenhouse gases.

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

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Morphology and Dynamics in Hydroxide-Conducting Polysulfones

In alkaline zinc–manganese dioxide batteries, there is a need for selective polymeric separators that have good hydroxide ion conductivity but that prevent the transport of zincate (Zn(OH) 4 ) 2- . Here we investigate the nanoscale structure and hydroxide transport in two cationic polysulfones that are promising for these separators. We present the synthesis and characterization for a tetraethylammonium-functionalized polysulfone (TEA-PSU) and compare it to our previous work on an N-butylimidazolium-functionalized polysulfone (NBI-PSU). We perform atomistic molecular dynamics (MD) simulations of both polymers at experimentally relevant water contents. The MD simulations show that both polymers develop well phase separated nanoscale water domains that percolate through the polymer. Calculation of the total scattering intensity from the MD simulations reveal weak or nonexistent ionomer peaks at low wave vectors. The lack of an ionomer peak is due to a loss of contrast in the scattering. The small water domains in both polymers, with median diameters on the order of 0.5–0.7 nm, lead to hydroxide and water diffusion constants that are 1–2 orders of magnitude smaller than their values in bulk water. This confinement lowers the conductivity but also may explain the strong exclusion of zincate from the PSU membranes seen experimentally.

25 ENERGY STORAGE↗
DOE OSTI2022

Scalable Simulation of Pressure Gradient-Driven Transport of Rarefied Gases in Complex Permeable Media Using Lattice Boltzmann Method

Accurate representations of slip and transitional flow regimes present a challenge in the simulation of rarefied gas flow in confined systems with complex geometries. In these regimes, continuum-based formulations may not capture the physics correctly. This work considers a regularized multi-relaxation time lattice Boltzmann (LB) method with mixed Maxwellian diffusive and halfway bounce-back wall boundary treatments to capture flow at high Kn. The simulation results are validated against atomistic simulation results from the literature. We examine the convergence behavior of LB for confined systems as a function of inlet and outlet treatments, complexity of the geometry, and magnitude of pressure gradient and show that convergence is sensitive to all three. The inlet and outlet boundary treatments considered in this work include periodic, pressure, and a generalized periodic boundary condition. Compared to periodic and pressure treatments, simulations of complex domains using a generalized boundary treatment conserve mass but require more iterations to converge. Convergence behavior in complex domains improves at higher magnitudes of pressure gradient across the computational domain, and lowering the porosity deteriorates the convergence behavior for complex domains.

42 ENGINEERING↗
NASA NTRS1979

Effect of pressure on structure and NO sub X formation in CO-air diffusion flames

A study was made of nitric oxide formation in a laminar CO-air diffusion flame over a pressure range from 1 to 50 atm. The carbon monoxide (CO) issued from a 3.06 mm diameter port coaxially into a coflowing stream of air confined within a 20.5 mm diameter chimney. Nitric oxide concentrations from the flame were measured at two carbon monoxide (fuel) flow rates: 73 standard cubic/min and 146 sccm. Comparison of the present data with data in the literature for a methane-air diffusion flame shows that for flames of comparable flame height (8 to 10 mm) and pseudoequivalence ratio (0.162), the molar emission index of a CO-air flame is significantly greater than that of a methane-air flame.

Maahs, H. G.↗
NASA NTRS2001

The sugar model: catalytic flow reactor dynamics of pyruvaldehyde synthesis from triose catalyzed by poly-l-lysine contained in a dialyzer

The formation of pyruvaldehyde from triose sugars was catalyzed by poly-l-lysine contained in a small dialyzer with a 100 molecular weight cut off (100 MWCO) suspended in a much larger triose substrate reservoir at pH 5.5 and 40 degrees C. The polylysine confined in the dialyzer functioned as a catalytic flow reactor that constantly brought in triose from the substrate reservoir by diffusion to offset the drop in triose concentration within the reactor caused by its conversion to pyruvaldehyde. The catalytic polylysine solution (400 mM, 0.35 mL) within the dialyzer generated pyruvaldehyde with a synthetic intensity (rate/volume) that was 3400 times greater than that of the triose substrate solution (12 mM, 120 mL) outside the dialyzer. Under the given conditions the final yield of pyruvaldehyde was greater than twice the weight of the polylysine catalyst. During the reaction the polylysine catalyst was poisoned presumably by reaction of its amino groups with aldehyde reactants and products. Similar results were obtained using a dialyzer with a 500 MWCO. The dialyzer method of catalyst containment was selected because it provides a simple and easily manipulated experimental system for studying the dynamics and evolutionary development of confined autocatalytic processes related to the origin of life under anaerobic conditions.

Biogenesis↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗