Search NASA⌕ Search

SEARCH · Search NASA

Results for “DECOUPLING”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2025

Understanding neutron capture processes in uranium deposits using combined U-Sm-Nd isotopic compositions

Valuable insights into the history and evolution of a geologic deposit can be found by investigating neutron capture reactions. Thermal neutron capture reactions occur within both the samarium (Sm) and the uranium (U) systems, where 149 Sm and 235 U can capture neutrons to become 150 Sm and 236 U, respectively. Although largely unexplored, paired measurements of 150 Sm and 236 U could be important for understanding neutron capture effects within uranium ore bodies, and such measurements are potentially useful in nuclear forensics for assessing a material's provenance or mineral exploration. In this work, we refined measurement procedures of Sm isotope compositions utilizing MC-ICPMS. While geologic reference materials were found to have indistinguishable Sm isotope compositions, we found significant isotope variations consistent with nuclear field shift among synthetic Sm standards. Here, this observation highlights that future high-precision Sm isotope investigations need to carefully evaluate synthetic standard(s) against geologic reference materials until an unfractionated and agreed-upon standard is identified. Here, we applied this method to a set of nine uranium ores from the South Australian Beverley North uranium deposits. Although 236 U excesses had been previously reported for these U ores, we found no measurable isotopic shifts in 149 Sm- 150 Sm at the current level of precision (±5 parts per million). One possible explanation for this disparity in the observed neutron capture signatures between U and Sm is that the source(s) of the U and Sm in these ores may be decoupled. This is consistent with the finding that these ores have variable 143 Nd/ 144 Nd, thus demonstrating that diverse sources were involved in the formation of the Beverley North deposits. Alternatively, this deposit may be too young (<50Ma) to have accumulated measurable neutron capture effects in Sm to be detected with the methods employed here.

Mineral exploration↗
DOE OSTI2026

Experimental study of a thermal energy storage-integrated heat pump system for load shifting in space cooling and heating

This study experimentally investigates an all-season thermal energy storage-integrated heat pump (TES-HP) system developed to enhance building energy efficiency and support grid-interactive operation through load shifting in both cooling and heating. A 4-ton commercial rooftop air-source heat pump was modified by integrating a hydronic thermal energy storage (TES) unit containing a phase change material (PCM) with a melting temperature of 22.0 °C. The system employs two reversing valves and three electronic expansion valves to enable six operating modes, including normal, TES charging, and TES discharging in both cooling and heating seasons. A subcooling-controlled electronic expansion valve was implemented to mitigate refrigerant maldistribution caused by unequal internal volumes among the heat exchangers. Compared with a baseline heat pump, the TES-HP reduced power consumption by 30–50% during cooling at an ambient temperature of 40.6 °C, and by up to 60% during heating at −15.0 °C, while maintaining high performance under cold-climate conditions. Additionaly to experimental evaluation, a simplified annualized on-peak analysis was conducted to compare the TES-HP with the baseline system, indicating an on-peak electricity reduction of approximately 876 kWhₑ per unit per year and associated on-peak cost savings under time-of-use pricing. These results demonstrate that the TES-HP effectively decouples HP operation from adverse ambient conditions, improving flexibility, resilience, and energy efficiency. This work advances prior research by experimentally proving the integration of PCM-based TES for a year-round, load-flexible HVAC operation.

42 ENGINEERING↗
DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Neural architecture search via similarity adaptive guidance

Evolutionary neural network architecture search (ENAS) has attracted the attention of many experts due to its global optimization capabilities to automatically search for convolutional neural network architectures based on the target task. The current search space for ENAS is not to design a fully structured network, but to search for smaller cell architectures to reduce search costs. However, blind search strategies do not effectively utilize the potential experience of the population. In order to utilize the potential experience learned by the current population to guide the evolutionary search of the population, we propose a similarity guided neural network architecture search algorithm based on cell architecture, which utilizes the similarity between pairwise architectures in the population as empirical knowledge learned by the population. Our proposed algorithm provides a novel method for calculating architecture similarity, which calculates architecture similarity separately from the cell and macro-structure. Then we decouple the connections and operations in the cell and calculate connection and operation similarity separately. In addition, we propose adaptive similarity selection and binary tournament selection strategies to enhance the algorithm’s global and local search capabilities and effectively explore the search space. Finally, we design an improved single-point crossover operator to enhance the local search ability of the evolutionary operator. The experimental results show that SAGNAS is a competitive algorithm that achieves 97.44% and 81.60% in CIFAR10 and CIFAR100 with only 1.9 GPU-days spent.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Macromolecular gelatin properties affect fibrin microarchitecture and tumor spheroid behavior in fibrin-gelatin gels

The biophysical properties of extracellular matrices (ECM) are known to regulate cell behavior, however decoupling cell behavior changes due to the relative contributions of material microstructure versus biomechanics or nutrient permeability remains challenging, especially within complex, multi-material matrices. We developed four gelatin-fibrin interpenetrating network (IPN) formulations which are identical in composition but possess variable gelatin molecular weight distributions, and display differences in microstructure, biomechanics, and diffusivity. In this work we interrogate the response of multicellular tumor spheroids to these IPN formulations and found that a high stiffness, gelatin-network dominated IPNs impeded remodeling and invasion of multicellular tumor spheroids; whereas relatively lower stiffness, fibrin-network dominated IPNs permitted protease-dependent remodeling and spheroid invasion. Cell proliferation correlated to nutrient diffusivity across tested IPN formulations. These findings demonstrate the complexity of ECM IPNs, relative to single polymer matrices, and highlight that cell response does not derive from a single aspect of the ECM, but rather from the interplay of multiple biomechanical properties. The methodology developed in this study represents a framework for future studies which aim to characterize cellular phenotypic responses to biophysical cues present within complex, multi-material matrices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A novel method for quantifying irradiation damage in nuclear graphite using Raman spectroscopy

Raman spectroscopy has long been used in studying irradiation damage in carbon/graphite materials. It is however unclear if the measurements from different types of materials are directly comparable. Further, decoupling the contribution from irradiation temperature and fluence to the total damage in nuclear graphite possessing irradiation damage gradient is currently not readily feasible as it requires a complete set of test reactor irradiation experiments over the relevant temperature and fluence range. A novel methodology has therefore been proposed and developed to quantify total irradiation damage evolution at crystal level based on graphite Raman G-band position shift. Specifically, G-band positions derived from a large number of spectra collected from the fractured surface of microfine-grained POCO ZXF-5Q graphite possessing proton irradiation damage gradient were plotted with open literature Raman data on HOPG, BEPO (AXGP graphite), PCEA and IG-110 graphite as a function of dpa. The total damage level within 2σ beam radius in this POCO graphite was estimated to be equivalent to ~2–5 dpa at ~350–370 °C. Derived G-band positions were then mapped to the three-stage amorphization trajectory model indicating beam centre area has entered the second stage, i.e., transitioning from nanocrystalline graphite into amorphous carbon. G-band position relative shift (ΔG) curve with a ‘turn-around’ peak as a function of total damage can be used for in-service lifetime prediction. Importantly, the developed methodology has the potential to ‘unify’ total damage levels across different grades of nuclear graphite caused by ion, neutron and proton irradiation at different temperatures.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Optimizing complexity in the kinetic modelling of integrated flue gas purification for pressurized oxy-combustion

A 5-step mechanism and a corresponding kinetic model are presented for integrated flue gas purification, where SOx and NOx are removed during pressurized oxy-coal combustion. Three coupled subsystems (‘blocks”) of gas–liquid transformations are considered: 1) oxidation of NO to NO 2 and dissolution of NO 2 in water; 2) dissolution of SO 2 in the water; and 3) interactions between S- and N- containing substances in the liquid phase. The model has been formulated in terms of only species that can be experimentally measured (model of optimal complexity). Subsystem I is limiting regarding NO removal. Subsystem II is very fast in comparison with the other subsystems and can be treated separately. Analytical expressions for all independent subsystems have been found and compared with computational results for the whole system. The temporal domain of validity of all these subsystems is presented. Finally, the validity of the NO 2 quasi-steady-state assumption is verified, and it is reported that variables such as NO, NO 2 and HNO 3 can be considered decoupled from the whole system since they appear in only Subsystem I, while HSO 3 - and HNO 2 are coupled. The presented kinetic model was verified experimentally by determining the 1:1 ratio of HNO 2 and HNO 3 , in gas–liquid NO-O 2 interaction. Additionally, new experimental facts on kinetic orders and stoichiometry between HNO 2 and HSO 3 - have been used to optimize the model. It is concluded that all NO x - SO x interactions can be essentially grasped by a sequence of only 4 irreversible reactions considering, in addition, the 5th reaction, i.e. the fast ‘SO 2 -water’ equilibrium.

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

Microwave-driven heterogeneous catalysis for activation of dinitrogen to ammonia under atmospheric pressure

This paper presents an innovative approach to producing energy-dense, carbon–neutral liquid ammonia as a means for carrying energy. This approach synergistically integrates microwave reaction chemistry with novel heterogeneous catalysis that decouples dinitrogen activation from high-temperature and high-pressure reactions, altering reaction pathways and increasing ammonia formation rate. Results presented here demonstrate that ammonia synthesis can be conducted at 280 ? and ambient pressure to achieve ~1 mmol ammonia/g cat/h over supported ruthenium catalyst systems utilizing microwave irradiation. It is further shown that adding promoter ions such as potassium, cerium, and barium significantly improves the ammonia production rate over undoped ruthenium-based catalysts. This effect could be attributed to enhanced dielectric loss processes that lead to stronger microwave absorption by the catalyst. Measurement of the equilibrium constant under microwave conditions showed a higher ammonia yield than under thermal equilibrium conditions for both the iron- and ruthenium-based catalysts. Finally, this study also illustrates the advantages of using a variable-frequency microwave reactor for ambient-pressure ammonia synthesis. Mechanistically, investigators believed that the oscillating electric fields of the radiation can couple with adsorbed nitrogen on the surface and accelerate its dissociation. Since dinitrogen dissociation on the surface is rate limiting, this effectively accelerates the reaction.

Hu, Jianli (John)↗
DOE OSTI2022

Cross-kingdom expression of synthetic genetic elements promotes discovery of metabolites in the human microbiome

Small molecules encoded by biosynthetic pathways mediate cross-species interactions and harbor untapped potential, which has provided valuable compounds for medicine and biotechnology. Since studying biosynthetic gene clusters in their native context is often difficult, alternative efforts rely on heterologous expression, which is limited by host-specific metabolic capacity and regulation. Here, in this work, we describe a computational-experimental technology to redesign genes and their regulatory regions with hybrid elements for cross-species expression in Gram-negative and -positive bacteria and eukaryotes, decoupling biosynthetic capacity from host-range constraints to activate silenced pathways. These synthetic genetic elements enabled the discovery of a class of microbiome-derived nucleotide metabolites—tyrocitabines—from Lactobacillus iners. Tyrocitabines feature a remarkable orthoester-phosphate, inhibit translational activity, and invoke unexpected biosynthetic machinery, including a class of “Amadori synthases” and “abortive” tRNA synthetases. Our approach establishes a general strategy for the redesign, expression, mobilization, and characterization of genetic elements in diverse organisms and communities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Cation uptake by fine aggregate in hardened cement mortar and its effect on electrical properties

Electrical resistivity tests can be used to evaluate the transport properties of concrete and provide a durability assessment. However, the electrical resistivity is largely dependent on the pore solution composition and recent work suggests that some aggregates have the capacity for cation uptake. This study first aims to provide further evidence for adsorption of cations on aggregate surfaces, without formation of reaction products (e.g. alkali-silica reaction). Secondly, hardened mortar samples were prepared using a fine aggregate with a high alkali affinity and a non-reactive fine aggregate as a control. The electrical resistivity of mortars was measured, and the pore solution of these mortars was obtained through high-pressure extraction. The effect of aggregate moisture dilution on the pore solution was decoupled by using a pore partitioning model. The results indicate that aggregate minerology can influence the pore solution composition through cation uptake. Specific minerals of minor quantity, like biotite, may be responsible for cation exchange. While adsorbed cations strongly affected pore solution and formation factor measurement, the bulk resistivity measurements on hardened mortar were only marginally influenced. Research on other implications of similar aggregate interactions with pore solutions are an intriguing area for future research.

Alkali-aggregate reaction↗
DOE OSTI2021

Sectoral and growth rate control on elemental uptake by individual calcite crystals

Heterogeneous distribution of trace elements (impurities) within individual calcite crystals is a phenomenon commonly observed in natural and laboratory systems. Changes in thermodynamic intensive parameters (mostly chemical potential and temperature) cannot always explain the inhomogeneous impurity patterns in calcite crystals and it has been suggested that growth rate and crystallographic orientation may exert strong effects on elemental incorporation into calcite. In addition, there are a number of experimental studies on micro-scale element (E) distribution between non-equivalent pairs of calcite vicinal faces (known as sectoral zoning); however, the variability of partition coefficients (e.g., K E = (E/Ca) calcite /(E/Ca) fluid ) within individual crystals remains undetermined. In this study, we have extended the work on elemental distribution between crystal sectors to evaluation of partition coefficients of trace and minor elements (Li, B, Mg, and Sr) in calcite crystal faces (10–14) and (01-12), whose growth rates were assessed. Growth entrapment model (GEM) and lattice strain theory were applied to explain K E heterogeneity by varying near-surface diffusivity of Mg and Sr and by varying surface enrichment factor for Li, B, Mg, and Sr. Decoupling of sectoral and growth rate effects reveals that sectoral zoning plays a key role in elemental distribution. More specifically, K Li and K B vary by more than one order of magnitude and K Mg varies by a factor of two within individual crystal faces. Strontium sectoral distribution is different from those of Li, B, and Mg and K Sr varies by up to a factor of two. Finally, these behaviors likely reflect different mechanisms of incorporation of Li, B, Mg, and Sr into calcite.

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

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

A comparative study of two numerical approaches for solving Kim–Kim–Suzuki phase-field models

Among the standard multi-phase multi-component phase-field (PF) methods, the Kim–Kim–Suzuki (KKS) method has the advantage of decoupling interfacial energy from bulk energy and solving concentration as the conserved variable. There are two approaches to numerically solving a KKS method: the global solution approach (GSA) solves all variables in a global system simultaneously, and the local solution approach (LSA) solves phase concentrations locally using a Newton solver. This work compares the performance of LSA and GSA for solving four KKS models of increasing complexity with the finite element method using the MOOSE framework. The solution accuracy, degrees of freedom (DOFs), memory usage, and computational efficiency are compared. We find that GSA and LSA generate similar solutions, with a maximum difference of only 0.34%. For each model, LSA has a lower number of DOFs, utilizes less memory, and less wall time. Additionally, the savings of memory and wall time in LSA increase with increasing mesh density of the same model and are more pronounced in models with higher dimensionality and more nodes. However, GSA is easier to implement in existing codes and can better solve highly nonlinear systems by utilizing sophisticated solvers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗
DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗