Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 397 records · Page 22

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Super-swelling behavior of stacked lipid bilayer systems

Bilayer systems comprising lipid mixtures are the most well-studied model of biological membranes. While the plasma membrane of the cell is a single bilayer, many intra- and extra-cellular biomembranes comprise stacks of bilayers. Most bilayer stacks in nature are periodic, maintaining a precise water layer separation between bilayers. That equilibrium water separation is governed by multiple inter-bilayer forces and is highly responsive. Biomembranes re-configure inter-bilayer spacing in response to temperature, composition, or mass transport cues. In synthetic bilayer systems for applications in cosmetics or topical treatments, control of the hydration level is a critical design handle. Herein we investigate a binary lipid system that leverages key inter-bilayer forces leading to unprecedented levels of aqueous swelling while maintaining a coherent multilamellar form. We found that combining cationic lipids with bicontinuous cubic phase-forming lipids (lipids with positive Gaussian modulus), results in the stabilization of multilamellar phases against repulsive steric forces that typically lead to bilayer delamination at high degrees of swelling. Using ultra-small-angle X-ray scattering alongside confocal laser scanning microscopy, we characterized various super-swelled states of 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP) and glycerol monooleate (GMO) lipids, as well as other analogous systems, at varied concentration and molar ratios. Through these experiments we established swelling profiles of various binary lipid systems that were near-linear with decreasing lipid volume fraction, showing maximum swelling with periodicity well above 200 nanometers.

42 ENGINEERING↗

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗

Not all types of secondary organic aerosol mix: two phases observed when mixing different secondary organic aerosol types

Abstract. Secondary organic aerosol (SOA) constitutes a large fraction of atmospheric aerosol. To assess its impacts on climate and air pollution, knowledge of the number of phases in internal mixtures of different SOA types is required. Atmospheric models often assume that different SOA types form a single phase when mixed. Here, we present visual observations of the number of phases formed after mixing different anthropogenic and biogenic SOA types. Mixing SOA types generated in environmental chambers with oxygen-to-carbon (O/C) ratios between 0.34 and 1.05, we found 6 out of 15 mixtures of two SOA types to result in two phase particles. We demonstrate that the number of phases depends on the difference in the average O/C ratio between the two SOA types (Δ(O/C)). Using a threshold Δ(O/C) of 0.47, we can predict the phase behavior of over 90 % of our mixtures, with one- and two-phase particles predicted for Δ(O/C)<0.47 and Δ(O/C)≥0.47, respectively. This threshold ΔO/C value provides a simple parameter to predict whether mixtures of fresh and aged SOA form one- or two-phase particles in the atmosphere. In addition, we show that phase-separated SOA particles form when mixtures of volatile organic compounds emitted from real trees are oxidized.

Mahrt, Fabian (ORCID:0000000270596765)↗

Milestone 1.2.15: Feasibility of In Situ Accelerated Aluminum Coupon Radiolysis by Synchrotron X-Rays

Extended dry storage of aluminum-clad spent nuclear fuel (ASNF) requires an assessment of the radiolytic generation of molecular hydrogen gas (H2) from the ASNF’s corrosion layers. H2 accumulation can potentially lead to storage canister embrittlement/rupture and the accumulation of flammable gas mixtures in accident scenarios. To date, computational models have been developed for the prediction of radiolytic H2 generation from samples exposed to lower absorbed dose regimes (up to ~3 MGy). However, greater accuracy is needed for higher absorbed doses (> 25 MGy) where the concentration of H2 ultimately reaches a steady-state. Data in this area is limited due to the amount of time taken (months) to accumulate such high gamma doses. Furthermore, a recent study observed radiation-induced damage to the surface of pre-corroded aluminum alloy coupons at high absorbed gamma doses. Despite this observation and its implications, there is currently no computational connection between the radiolytic formation of H2 and changes in the composition and morphology of the cladding’s corrosion layers with absorbed dose. To develop a better understanding of the processes and effects described above, the present study aimed to evaluate the feasibility of leveraging synchrotron x-ray capabilities for coupled accelerated irradiation and in situ surface characterization of ASNF alloys. To that end, National Synchrotron Light Source II (NSLS-II) beamtime was secured and a series of aluminum alloy 1100 (AA1100) wires, prepared under a variety of conditions, were interrogated using ex situ x-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques at Idaho National Laboratory (INL) and ex situ and in situ XRD capabilities at the NSLS-II x-ray powder diffraction (XPD) beamline. Results indicated that performing synchrotron XRD using the available XPD beamline configuration can provide sufficient radiation dose but cannot discern changes in the composition and morphology of sample corrosion layers. Additionally, ex situ XRD and SEM data from Idaho National Laboratory (INL) indicated that the pre-corrosion of AA1100 wires did not generate an appreciably thick corrosion layer, as compared with previous aluminum coupon studies. These thinner surface corrosion layers were not sufficient for detection by both ex situ and in situ synchrotron XRD, specifically under consideration of the NSLS-II XPD beamline’s configuration primarily capturing information from the sample’s bulk material, i.e., aluminum metal. In summary, the use of the NSLS-II XPD beamline for promoting accelerated radiolysis and in situ surface characterization is not appropriate for this program’s goals.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

High-Fidelity Arc-Discharge Model for Hydrogen-Plasma-Smelting-Reduction of Iron Ore

Electrification and use of renewable hydrogen is currently a necessity for decarbonizing the iron-and-steel industry. In this regard, hydrogen plasma smelting reduction (HPSR) is a novel pathway that is being explored for reduction of iron ore. HPSR provides several decarbonization merits compared to conventional blast furnaces. Firstly, the use of renewable hydrogen drastically reduces the CO2 emissions compared to the use of coke. Secondly, renewable electricity in the form of a thermal plasma for making reactive hydrogen species (radicals, ions) are more efficient at reducing iron ore compared to neutral H2. Thirdly, a molten product compatible with downstream processes is obtained from the intense heat transfer from the plasma. However, the scale-up of this technology requires fundamental exploration of hydrogen plasma dynamics and its interaction with complex solid material that include phase changing iron-ore and slag. In this work, we present a first principles continuum scale model for thermal plasmas in Ar/H2 gas mixtures typically used for HPSR. The thermal plasma governing equations for mass, momentum and energy with Lorentz force and Joule heating source terms are solved along with electromagnetic equations for electrostatic and magnetic vector potential. Our solver will be based on Pele, a suite of reacting flow solvers designed for advanced scientific computing architectures (Henry De Frahan et al., Proceedings of SIAM Parallel Processing, 13-25, 2024), and will utilize adaptive mesh generation for enhanced resolutions at locations of intense physicochemical interactions. This study will present the impact of Ar to H2 ratios on excited/dissociated hydrogen species concentrations, plasma temperature and conductivity along with the impact of outgassed species (water, metal vapor, O, OH radicals) from ore surface on gas phase chemistry. Furthermore, the heat and species flux to the surface will be quantified as a function of applied voltages in a transferred arc configuration.

hydrogen plasma↗

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗

Extraction of Furfural and Furfural/5-Hydroxymethylfurfural from Mixed Lignocellulosic Biomass-Derived Feedstocks

We present a combined in silico and experimental study on the extraction of furfural and 5-hydroxymethylfurfural (HMF) in aqueous–organic biphasic systems. We predict the liquid–liquid equilibria and furfural/HMF partition coefficients of over 2200 water-organic biphasic systems using the multiscale COSMO-RS model and measure experimentally single-component (furfural) and mixture (furfural and HMF) partition coefficients at room and dehydration reaction-relevant temperatures in 28 solvents. We find the experimental data to be within a factor of 2 from the COSMO-RS predictions. Even though furfural and HMF have chemical similarity, the slight differences in molecular structure render the separation of furfural easier by the supply of more solvents of higher partition coefficient for extraction. We leverage this molecular difference and experimentally demonstrate that with an additional extraction step, using dichloromethane or toluene, we can selectively extract furfural from furfural-HMF mixtures, which can coexist in lignocellulosic biomass dehydration products, despite their partition coefficients being generally correlated. We complement solvent selection criteria for biphasic lignocellulosic biomass processes with a simple mass balance extraction model for determining volume ratios in multistage extraction. Lastly, the molecular nature of the preferential furfural extraction is rationalized using COSMO-RS σ-profile analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Engineering Escherichia coli to produce medium chain oleochemicals from C 2 substrates

Many strategies to create a circular bioeconomy have been proposed. To be successful, CO 2 must be reduced with renewable energy into chemical building blocks, from which the chemical industry can be supported. Circular strategies include leveraging photosynthesis to produce sugar and lipid intermediates or renewable electricity to produce hydrogen or other electron carriers to support CO 2 reduction. Acetogens can anaerobically reduce CO 2 with H 2 to produce mixtures of small organic molecules in gas fermentations. We previously demonstrated that acetate, a common product of gas fermentation, can be converted to the model oleochemical dodecanol in engineered Escherichia coli. Here, we explored the conversion of ethanol and mixtures of ethanol and acetate to the same model oleochemicals. Co-feeding ethanol can supply both carbon and additional reducing power relative to acetate alone. In this work, we engineered E. coli to catabolize ethanol and expressed two distinct ethanol metabolism pathways in different operons and combined them with improved engineered acetate activation. We evaluated the performance of these operons in dodecanol-producing strains when fed ethanol or acetate and found ethanol to be a better carbon source when judged by product titers. The engineered strains fed ethanol produced about two-fold more dodecanol than the strains fed acetate. This increase was in part, due to change in product distribution. Cells fed ethanol produced predominantly dodecanol, whereas cells fed acetate generated a mixture of dodecanol and dodecanoic acid. Dodecanol titers were further improved by employing feeding strategies in controlled bioreactors.

Escherichia coli↗

Reaction mechanism for atmospheric pressure plasma treatment of cysteine in solution

Mechanisms for the cold atmospheric plasma (CAP) treatment of cells in solution are needed for more optimum design of plasma devices for wound healing, cancer treatment, and bacterial inactivation. However, the complexity of organic molecules on cell membranes makes understanding mechanisms that result in modifications (i.e. oxidation) of such compounds difficult. As a surrogate to these systems, a reaction mechanism for the oxidation of cysteine in CAP activated water was developed and implemented in a 0-dimensional (plug-flow) global plasma chemistry model with the capability of addressing plasma-liquid interactions. Reaction rate coefficients for organic reactions in water were estimated based on available data in the literature or by analogy to gas-phase reactions. The mechanism was validated by comparison to experimental mass-spectrometry data for COST-jets sustained in He/O 2 , He/H 2 O and He/N 2 /O 2 mixtures treating cysteine in water. Results from the model were used to determine the consequences of changing COST-jet operating parameters, such as distance from the substrate and inlet gas composition, on cysteine oxidation product formation. Results indicate that operating parameters can be adjusted to select for desired cysteine oxidation products, including nitrosylated products.

60 APPLIED LIFE SCIENCES↗

Geochemical narrowing of cement fracture aperture during multiphase flow of supercritical CO2 and brine

For carbon capture and storage operations, wells are used as a necessary conduit for injecting CO2 at depth, but they can also act as leakage conduits if the cement seals are compromised. The specific objective of this research was to investigate the nature of self-healing of a fracture within cement under multiphase flow of CO2 and brine, and to compare this to a newly developed model. This was accomplished through a laboratory investigation, flowing reacting supercritical CO2 and brine through a cement fracture with multiphase flow. The influent end of the fracture was degraded as evidenced by the formation of cracks across the surface. At the effluent end of the fracture, the initial aperture of 137 µm was reduced to 50 µm. This reduction by 87 µm was similar to results from a recently developed model tested in this study, which predicted an aperture reduction of 80 µm. Our results demonstrated that the model accurately predicted the reactions which occurred in this multiphase flow system.

wellbore, cement, fracture, multiphase-flow, super↗

Distance Dependence of Förster Resonance Energy Transfer Rates in 2D Perovskite Quantum Wells via Control of Organic Spacer Length

Two-dimensional (2D) semiconductors are attractive candidates for a variety of optoelectronic applications owing to the unique electronic properties that arise from quantum confinement along a single dimension. Incorporating nonradiative mechanisms that enable directed migration of bound charge carriers, such as Forster resonance energy transfer (FRET), could boost device efficiencies provided that FRET rates outpace undesired relaxation pathways. However, predictive models for FRET between distinct 2D states are lacking, particularly with respect to the distance d between a donor and acceptor. We approach FRET in systems with binary mixtures of donor and acceptor 2D perovskite quantum wells (PQWs), and we synthetically tune distances between donor and acceptor by varying alkylammonium spacer cation lengths. FRET rates are monitored using transient absorption spectroscopy and ultrafast photoluminescence, revealing rapid picosecond lifetimes that scale with spacer cation length. We theoretically model these binary mixtures of PQWs, describing the emitters as classical oscillating dipoles. Finally, we find agreement with our empirical lifetimes and then determine the effects of lateral extent and layer thickness, establishing fundamental principles for FRET in 2D materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bubble Point Measurements of cis-1,1,1,4,4,4-Hexafluorobutene [R-1336mzz(Z)] + trans-1,2-Dichloroethene [R-1130(E)] mixtures

Saturation pressures of pure R-1336mzz(Z) and R-1130(E) and bubble point pressures of three R-1336mzz(Z)/1130(E) blends were measured from 265 K to 360 K. For each pure refrigerant or refrigerant blend, a total of twenty unique saturation pressures or bubble points were measured. In total 100 unique state points were obtained. Presently, no Helmholtz-energy-explicit type equation of state (EoS) is available for R-1130(E). While an extended corresponding states EoS for R-1130(E) is available to estimate the properties of the R-1336mzz(Z)/1130(E) blend, this model does not resolve the azeotropic behavior of the mixture. Therefore, the perturbed-chain statistical-associating fluid theory (PC-SAFT) EoS is used to model the vapor–liquid equilibria of the R-1336mzz(Z)/1130(E) blend. PC-SAFT model parameters are reported, and the overall performance of the model is characterized by deviations from the experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanistic modeling of copper corrosions in data center environments

Air-side economizers are increasingly used to take advantage of “free-cooling” in data centers with the intent of reducing the carbon footprint of buildings. However, they can introduce outdoor pollutants to indoor environment of data centers and cause corrosion damage to the information technology equipment. Here, to evaluate the reliability of information technology equipment under various thermal and air-pollution conditions, a mechanistic model based on multi-ion transport and chemical reactions was developed. The model was used to predict Cu corrosion caused by Cl2-containing pollutant mixtures. It also accounted for the effects of temperature (25 °C and 28 °C), relative humidity (50%, 75%, and 95%), and synergism. It also identified higher air temperature as a corrosion barrier and higher relative humidity as a corrosion accelerator, which agreed well with the experimental results. The average root mean square error of the prediction was 13.7 Å. The model can be used to evaluate the thermal guideline for data centers design and operation when Cl2 is present based on pre-established acceptable risk of corrosion in data centers’ environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dissolved gas recovery from water using a sidestream hollow-fiber membrane module: First principles model synthesis and steady-state validation

This paper presents a first-principles model for the recovery of dissolved gases from liquids using a sidestream hollow-fiber membrane module. The model avoids the use of new empirical coefficients, thus providing a parametric understanding of the process behavior for future design and optimization of membrane modules. This type of first-principles model could be particularly useful when gas recovery is beneficial to biological or chemical reactions of interest, such as the acetogenesis reactions in two-stage anaerobic digesters. The steady-state behavior of the model was validated against both new experimental data for the recovery of H 2 , CH 4 and H 2 –CH 4 mixtures from pure water, as well as existing published data. The modeled gas recovery predictions agreed with experimental data to an absolute average error of 13%, and an average R value of 0.98. Parametric analysis of mixed-gas recovery suggests possible key transition points in the composition of the recovered gases. For example, at 40 °C, increasing trans-membrane pressure while keeping hydraulic residence time (HRT) under 0.5 s will result in an increase in the ratio of H 2 to CH 4 recovered. Otherwise, increasing trans-membrane pressure will instead decrease the ratio of H 2 to CH 4 recovered. The model has potential to be extended to transient analysis, but has yet to be validated with transient experimental data. Furthermore, this model was successfully implemented in both Python and MATLAB, and provides valuable insights for future net-energy optimization for anaerobic digestion systems with in-situ gas recovery.

Anaerobic Digestion↗