SEARCH · Engineering Papers

Results for “mixed effects”

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.

SourceResearch connectionFrom yearThrough year

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

At least 973 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Computational Fluid Dynamics Combustion Modeling for Rotating Detonation Engines

This paper focuses on the development and validation of a combustion model for Computational Fluid Dynamics (CFD) modeling of Rotating Detonation Engines. A zero-dimensional Partially Stirred Reactor (PaSR) with a detailed chemical kinetic mechanism for hydrogen and air is used to model turbulent combustion. The model is computationally efficient and is based on the notion of partial mixing at the sub-grid level with turbulent exchange between mixed and unmixed regions. The ability of the PaSR model to accurately represent both detonative and deflagrative combustion is assessed by validating the results against experimental data. The effects of mesh resolution on the solution are also studied in order to determine if a mesh independent solution is obtainable with the Large Eddy Simulation (LES) approach to modeling turbulence. A comparison is made between the PaSR model and simply ignoring turbulence chemistry interactions which assumes that all species are perfectly mixed at the sub-grid level.

Strakey, Peter↗

DOE OSTI2024

Directional Superradiance in a Driven Ultracold Atomic Gas in Free Space

Ultracold atomic systems are among the most promising platforms that have the potential to shed light on the complex behavior of many-body quantum systems. One prominent example is the case of a dense ensemble illuminated by a strong coherent drive while interacting via dipole-dipole interactions. Despite being subjected to intense investigations, this system retains many open questions. A recent experiment carried out in a pencil-shaped geometry [Ferioli Nat. Phys. 19, 1345 (2023)] has reported measurements that have seemed consistent with the emergence of strong collective effects in the form of a “superradiant” phase transition in free space, when looking at the light-emission properties in the forward direction. Motivated by the experimental observations, we carry out a systematic theoretical analysis of the steady-state properties of the system as a function of the driving strength and atom number N . We observe signatures of collective effects in the weak-driving regime, which disappear with increasing drive strength as the system evolves into a single-particle-like mixed state comprised of randomly aligned dipoles. Although the steady state features some similarities to the reported superradiant-to-normal nonequilibrium transition, also known as cooperative resonance fluorescence, we observe significant qualitative and quantitative differences, including a different scaling of the critical drive parameter (from N to N ). We validate the applicability of a mean-field treatment to capture the steady-state dynamics under currently accessible conditions. Furthermore, we develop a simple theoretical model that explains the scaling properties by accounting for interaction-induced inhomogeneous effects and spontaneous emission, which are intrinsic features of interacting disordered arrays in free space. Published by the American Physical Society 2024

Agarwal, Sanaa (ORCID:0000000336863804)↗

DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

DOE OSTI2026

Studies of localized mix and asymmetry caused by the stalk in inertial fusion experiments at OMEGA

In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Intermediate Principal Stress Effects on the 3D Cracking Behavior of Flawed Rocks Under True Triaxial Compression

Crack initiation, growth, and coalescence in flawed rocks have been extensively studied for 2D (planar, penetrating) flaws under uniaxial/biaxial compression. However, little is known as to the mechanisms and processes of cracking from 3D flaws under true triaxial compression, where the intermediate principal stress (σ 2 ) is distinguished from the major and minor principal stresses. In this work, we systematically investigate the effects of σ 2 on the 3D cracking behavior of rock specimens with preexisting flaws, through the use of mechanistic simulations of mixed-mode fracture in rocks. We explore how two characteristics of σ 2 , namely, (i) its orientation with respect to the flaw and (ii) its magnitude, affect two aspects of the cracking behavior, namely, (i) the cracking pattern and (ii) the peak stress. Results show that the orientation of σ 2 exerts more control over the cracking pattern than the flaw inclination angle. The peak stress becomes highest when σ 1 is parallel to the flaw, whereas it becomes lowest when σ 2 is parallel to the flaw. Also, the effects of σ 2 magnitude are more significant when σ 2 becomes more oblique to the flaw plane. On the basis of our observations, we propose mechanisms underlying the cracking behavior of 3D flawed rocks under true triaxial compression.

58 GEOSCIENCES↗

DOE OSTI2025

Insights into determining pore size properties of ultrafiltration membranes

The selectivity of porous membranes is often characterized using solute rejection tests, where membranes are challenged with dilute aqueous solutions of neutral solutes at operating conditions that minimize concentration polarization and fouling. In single solute tests, a membrane is challenged with one molecular weight (MW) solute at a time from low to high MW. Since single solute methods are time-intensive, mixed solute tests have become more common, where a mixture of several MW solutes challenges a membrane at once. However, the presence of large solutes in a mixture increases the rejection of smaller solutes. Furthermore, there are no universally accepted operating conditions or standard methods used by membrane manufacturers or researchers for the experiments, leading to difficulty in pore size and pore characteristic comparisons. In this paper, commercial ultrafiltration membranes were challenged with single and mixed solute polyethylene glycol (PEG) and dextran aqueous solutions. First, rejection values determined using total organic carbon (TOC) and high-performance liquid chromatography (HPLC) from single solute filtration experiments are compared. Differences in rejection curves obtained by the two techniques are attributed to solute polydispersity. Mixed solute filtration experiments with binary mixtures of solutes showcased co-solute interactions, which increase with both the size and weight percent of large solute in the mixture. Mixed solute filtration experiments at varying operating conditions (i.e., stir speed and flux) were conducted to determine operating conditions that mitigate co-solute interactions. Stir speed had a minimal effect on co-solute interactions. In contrast, low flux conditions can help minimize co-solute interactions, leading to pore size distributions that closely resemble results observed in single solute filtration using narrowly dispersed solutes. Additionally, at low flux conditions, the predicted membrane pore size distributions utilizing mixed solute experiments with PEG and dextran were similar.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Compositional patterning in irradiated alloys: Effective potentials and effective interfacial energy

Compositional patterning (CP) in binary alloys during energetic particle irradiation is studied using a kinetic model that considers two competing kinetic processes, a thermally activated one promoting macroscopic phase separation (MPS) of the concentration field c(r,t) and a forced one resulting in finite-range random atomic mixing. The forced mixing is modeled by a Gaussian relocation distribution with a characteristic distance R. A series of approximate kinetic models are introduced by expanding the mixing function into a series of n terms, thus replacing the non-local evaluations of the concentration field c(r'-r,t) by local derivatives of c(r,t). This approach makes it possible to obtain exact effective potentials and build steady-state diagrams for each order-n model. Phase field (PF) simulations using these order-n models reveal that near the onset of patterning, phase evolution is accurately described using an order-3 model, which changes smoothly from an extended Cahn-Hilliard free energy in the MPS regime to a one-mode Swift-Hohenberg functional in the CP regime. Deeper into the patterning regime, higher-order models are required to achieve convergence, yielding square-like concentration profiles characteristic of a strong segregation regime. These higher-order effective free energies are analogous to multimodal Swift-Hohenberg functionals. Here, a new definition for the effective interfacial energy is proposed in the CP regime, since the interfacial area is no longer an excess quantity in that regime, precluding the use of the standard thermodynamic definition of interfacial energy.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Influence of pilot-fuel mixing on the spatio-temporal progression of two-stage autoignition of diesel-sprays in low-reactivity ambient fuel-air mixture

The spatial and temporal locations of autoignition for direct-injection compression-ignition engines depend on fuel chemistry, temperature, pressure, and mixing trajectories in the fuel jets. Dual-fuel systems can provide insight into both fuel-chemistry and physical effects by varying fuel reactivities and engine operating conditions. In this context, the spatial and temporal progression of two-stage autoignition of a diesel-fuel surrogate, n-heptane, in a lean-premixed charge of synthetic natural-gas (NG) and air is imaged in an optically accessible heavy-duty diesel engine. The lean-premixed charge of NG is prepared by fumigation upstream of the engine intake manifold. Optical diagnostics include high-speed (15kfps) cool-flame chemiluminescence-imaging as an indicator of low-temperature heat-release (LTHR) and OH* chemiluminescence-imaging as an indicator high-temperature heat-release (HTHR). NG prolongs the ignition delay of the pilot fuel and increases the combustion duration. Zero-dimensional chemical-kinetics simulations provide further understanding by replicating a Lagrangian perspective for mixtures evolving along streamlines originating either at the fuel nozzle or in the ambient gas, for which the pilot-fuel concentration is either decreasing or increasing, respectively. The zero-dimensional simulations predict that LTHR initiates most likely on the air streamlines before transitioning to HTHR, either on fuel-streamlines or on air-streamlines in regions of near-constant φ. Due to the relatively short pilot-fuel injection-durations, the transient increase in entrainment near the end of injection (entrainment wave) is important for quickly creating auto-ignitable mixtures. To achieve desired combustion characteristics, e.g., multiple ignition-kernels and favorable combustion phasing and location (e.g., for reducing wall heat-transfer or optimizing charge stratification), adjusting injection parameters could tailor mixing trajectories to offset changes in fuel ignition chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

In‐situ Analysis of Paste Properties in Resonant Acoustic Mixers for Quality Monitoring

Formulation control is key to achieving consistent target properties of energetic materials, as feedstock variations and slight deviations in the ratios of different ingredients can have major effects on final product properties, particularly in dense pastes with high particle loading >65 vol.%. In large‐scale operations, it is imperative to either correct or remove batches of material that perform outside baseline property specifications as early as possible to avoid unnecessary processing of suboptimal material. Quality monitoring is the practice of measuring material properties during processing using process analytical technologies as opposed to only testing the properties of the final product; it is a key principle in the quality‐by‐design frameworks used for designing formulations and manufacturing processes. Herein, a process analytical technology method for correlating material properties of dense pastes directly after mixing in a Resonant Acoustic Mixer to motor data is developed and used to detect differences in the particle content of dense paste formulations. This method was also capable of detecting variations in powder feedstock properties, such as particle packing efficiency, and is sensitive enough to detect changes of 2 wt.% in the total solids content of the formulation. The techniques presented herein show excellent promise for use as a process analytical technology capable of quantifying formulation effects on material movement modes during resonant acoustic mixing.

Materials science↗

DOE OSTI2021

Understanding interfacial segregation in polymer blend films with random and mixed side chain bottlebrush copolymer additives

Bottlebrush polymers are complex macromolecules with tunable physical properties dependent on the chemistry and architecture of both the side chains and the backbone. Prior work has demonstrated that bottlebrush polymer additives can be used to control the interfacial properties of blends with linear polymers but has not specifically addressed the effects of bottlebrush side chain microstructures. Here, using a combination of experiments and self-consistent field theory (SCFT) simulations, we investigated the effects of side chain microstructures by comparing the segregation of bottlebrush additives having random copolymer side chains with bottlebrush additives having a mixture of two different homopolymer side chain chemistries. Specifically, we synthesized bottlebrush polymers with either poly(styrene-ran-methyl methacrylate) side chains or with a mixture of polystyrene (PS) and poly(methyl methacrylate) (PMMA) side chains. Furthermore, the bottlebrush additives were matched in terms of PS and PMMA compositions, and they were blended with linear PS or PMMA chains that ranged in length from shorter to longer than the bottlebrush side chains. Experiments revealed similar behaviors of the two types of bottlebrushes, with a slight preference for mixed side-chain bottlebrushes at the film surface. SCFT simulations were qualitatively consistent with experimental observations, predicting only slight differences in the segregation of bottlebrush additives driven by side chain microstructures. Specifically, these slight differences were driven by the chemistries of the bottlebrush polymer joints and side chain end-groups, which were entropically repelled and attracted to interfaces, respectively. Using SCFT, we also demonstrated that the interfacial behaviors were dominated by entropic effects with high molecular weight linear polymers, leading to enrichment of bottlebrush near interfaces. Surprisingly, the SCFT simulations showed that the chemistry of the joints connecting the bottlebrush backbones and side chains played a more significant role compared with the side chain end groups in affecting differences in surface excess of bottlebrushes with random and mixed side chains. This work provides new insights into the effects of side chain microstructure on segregation of bottlebrush polymer additives.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dark magnetohydrodynamics: Black hole accretion in superradiant dark photon clouds

Black holes threaded by massive vector fields can be subject to a superradiant instability, growing a cloud of massive vector particles around it. In this work, we consider what happens if such a dark matter candidate field mimicking a dark photon interacts with an accretion flow onto the black hole. By including a kinetic mixing term with the standard model photon, we extend the commonly used equations of general-relativistic magnetohydrodynamics to a dark photon constituent. The coupling to the dark photon then appears as an effective dynamo term together with a dark Lorentz force acting on the accreting matter. We numerically study the interactions between the superradiant dark photon cloud and the inner accretion flow by solving the coupled system in full numerical relativity. By parameterically varying the mixing parameter between the dark and standard model sector, we provide a first investigation of how the accretion flow could be modified. In conclusion, depending on the coupling strength, our solutions exhibit increased wind launching, as well as oscillation modes in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind Effects

Here, the physical release behavior of hydrogen is important to understand from a safety and design perspective. The consequences of unignited pressurized gaseous hydrogen plumes exiting vent stacks were considered by extending and modifying existing general hydrogen plume models. Entrainment, vent stack backpressure, and the flow regime of hydrogen exiting the vent were found to be significant factors affecting plume shape and size, but further investigation and validation with unchoked, low-Froude-number flows is recommended to improve the model’s robustness. Additionally, models for the effects of wind on unignited plume momentum and entrainment were added to explore this behavior. Wind was assumed to increase mixing of hydrogen with the ambient air, and to affect the momentum of the released jet. Introducing wind into the plume model led to a shorter plume for all wind and jet directions. A high counter-flowing wind led to non-physical results and challenges in interpreting the visualization. The proposed jet plume wind sub-models (specifically entrainment coefficients) were fit and compared to experimental data of different releases of hydrogen into a wind tunnel, but the quantity of data available and experimental conditions were limited. Thus, collection of more empirical data and for a wider range of conditions is recommended for improvement of the proposed computational models. Developing reduced-order models for these physical phenomena can improve accessibility to predicted physical behavior and the rate at which hydrogen systems can be safely designed and deployed.

entrainment↗

DOE OSTI2020

An Observational Comparison of Level of Neutral Buoyancy and Level of Maximum Detrainment in Tropical Deep Convective Clouds

Tropical deep convective clouds are important drivers of large-scale atmospheric circulation representing the main vertical transport pathway through the depth of the troposphere for heat, momentum, water, and chemical species. The strength and depth of this transport are impacted by the convective updraft size and intensity that are driven by buoyancy, dynamical forcing, and mixing of environmental air, i.e., entrainment. In this paper, we identify tropical deep convective systems with well-defined forward anvils using Atmospheric Radiation Measurement (ARM) ground-based profiling radars, at three ARM fixed-sites in the Tropical Western Pacific (TWP; i.e., Manus, Nauru, Darwin) and three ARM Mobile Facility deployments in Niamey, Niger; Gan Island, Maldives; and Manacapuru, Brazil. We use the difference between the level of neutral buoyancy (LNB) and the level of maximum detrainment (LMD) as a proxy for the effective bulk convective entrainment (εproxy). The LNB, the theoretical height that a parcel raised above the level of free convection would reach with no mixing, is calculated based on pre-convection radiosonde measurements using parcel theory. The LMD is the height of the maximum reflectivity observed in forward anvil clouds by profiling radars. Deep convective systems over the TWP show higher LNBs that extend to 16.3 km on average and larger εproxy (median LNB minus LMD up to 6.5 km) compared to their continental counterparts in the Amazon and West Africa. Oceanic conditions show larger convective available potential energy (CAPE) coupled with higher moisture at low levels which favors larger εproxy. In contrast, continental cases initiate and develop under high convective inhibition, steeper environmental lapse rate, and high wind shear conditions, which show smaller offset between LNB and LMD. Deep convective cases that promote significant cold pools at the surface experience less εproxy. Using a Random Forest regression algorithm, CAPE is associated with the highest feature importance score for predicting convective εproxy, followed by low-level relative humidity. For continental cases, the low-level wind shear also indicates higher importance.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Time-lens photon Doppler velocimetry (TL-PDV)

Here, we describe a time lens (TL) to expand the dynamic range of photon Doppler velocimetry (PDV) systems. The principle and preliminary design of a TL-PDV system are explained and shown to be feasible through simulations. In a PDV system, an interferometer is used for measuring frequency shifts due to the Doppler effect from the target motion. However, the sampling rate of the electronics could limit the velocity range of a PDV system. A four-wave-mixing (FWM) TL applies a quadratic temporal phase to an optical signal within a nonlinear FWM medium (such as an integrated photonic waveguide or a highly nonlinear optical fiber). By spectrally isolating the mixing product, termed the idler, and with appropriate lengths of dispersion prior to and after this FWM TL, a temporally magnified version of the input signal is generated. Therefore, the frequency shifts of PDV can be “slowed down” with the magnification factor M of the TL. M = 1 corresponds to a regular PDV system without a TL. M = 10 has been shown to be feasible for a TL-PDV system. The use of this effect for PDV can expand the velocity measurement range and allow for the use of lower bandwidth electronics. TL-PDV will open up new avenues for various dynamic material experiments.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Towards precise and accurate calculations of neutrinoless double-beta decay

We present the results of a National Science Foundation Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a roadmap for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The roadmap includes specific projects in theoretical and computational physics as well as the use of Bayesian methods to quantify both intra- and inter-model uncertainties. Here, the goal of this ambitious program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Nanopore Confinement of C-H-O Mixed-Volatile Fluids Relevant to Subsurface Energy Systems

The overarching goal of OU’s contribution to the project was to use novel computational approaches to quantify the effect of variable nanopore features (pore size, volume, topology, and chemistry) and the presence of water or salt solutions on the sorption and transport of carbon-bearing fluids. To achieve this goal OU collaborators tested test two related hypotheses: • The transport of water and aqueous electrolytes in nanopores is controlled by pore size, pore wall composition, and the type of salt – structure-maker (e.g., CaCl 2 ) versus structure-breaker (e.g., NaCl) that affect the hydrogen-bonding network. • The structure, solubility, and transport of carbon-bearing molecules in water or aqueous electrolyte-filled nanopores are controlled by the substrate type including degree of hydration, and pore features, and regulated by the hydration structure of the guest molecules.

58 GEOSCIENCES↗