Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural settings”

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Machine learning and serving of discrete field theories

A method for machine learning and serving of discrete field theories in physics is developed. The learning algorithm trains a discrete field theory from a set of observational data on a spacetime lattice, and the serving algorithm uses the learned discrete field theory to predict new observations of the field for new boundary and initial conditions. The approach of learning discrete field theories overcomes the difficulties associated with learning continuous theories by artificial intelligence. The serving algorithm of discrete field theories belongs to the family of structure-preserving geometric algorithms, which have been proven to be superior to the conventional algorithms based on discretization of differential equations. The effectiveness of the method and algorithms developed is demonstrated using the examples of nonlinear oscillations and the Kepler problem. In particular, the learning algorithm learns a discrete field theory from a set of data of planetary orbits similar to what Kepler inherited from Tycho Brahe in 1601, and the serving algorithm correctly predicts other planetary orbits, including parabolic and hyperbolic escaping orbits, of the solar system without learning or knowing Newton’s laws of motion and universal gravitation. The proposed algorithms are expected to be applicable when the effects of special relativity and general relativity are important.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Polymer-based microfluidic device for on-chip counter-diffusive crystallization and in situ X-ray crystallography at room temperature

Proteins are long chains of amino acid residues that perform a myriad of functions in living organisms, including enzymatic reactions, signalling, and maintaining structural integrity. Protein function is determined directly by the protein structure. X-ray crystallography is the primary technique for determining the 3D structure of proteins, and facilitates understanding the effects of protein structure on function. The first step towards structure determination is crystallizing the protein of interest. We have developed a centrifugally-actuated microfluidic device that incorporates the fluid handling and metering necessary for protein crystallization. Liquid handling takes advantage of surface forces to control fluid flow and enable metering, without the need for any fluidic or pump connections. Our approach requires only the simple steps of pipetting the crystallization reagents into the device followed by either spinning or shaking to set up counter-diffusive protein crystallization trials. The use of thin, UV-curable polymers with a high level of X-ray transparency allows for in situ X-ray crystallography, eliminating the manual handling of fragile protein crystals and streamlining the process of protein structure analysis. We demonstrate the utility of our device using hen egg white lysozyme as a model system, followed by the crystallization and in situ, room temperature structural analysis of the hub domain of calcium–calmodulin dependent kinase II (CaMKIIβ).

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

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

5.1.6 Demo #6: NanoStitch (NAWAStitch) Multifunctional Composites Trade Validation

The overall objective of the NanoStitch Multifunctional Composites Trade Validation Demonstration project was to determine the value proposition of adding NanoStitch to rapid-cure composite systems for structures used in the expanding sector of advanced air mobility (AAM) and specifically electric vertical take-off and landing (eVTOL) aircraft. The expected value proposition was stronger composite parts, with new multifunctional attributes, that are lower cost to fabricate and maintain. NanoStitch is a new composite innovation where a contiguous vertically aligned carbon nanotube (VACNT) forest is located between each ply of a laminate composite resulting in improved mechanical properties as well as other multifunctional properties such as enhanced thermal and electrical conductivity. NanoStitch was originally distributed by N12, later acquired by NAWA Technologies, now called NAWAStitch, and henceforth referred to in this document as NAWAStitch. Two significant evaluations were framed to help meet the overall objective. First, was Task 1 to evaluate the effects that adding NAWAStitch to a rapid cure material system would have on the mechanical and processing properties of the material. Second, was Task 2 to evaluate the workability of adding NAWAStitch to practical flight demonstration articles. From previous research, it was already known that the application of NAWAStitch to a composite structure can enhance the mechanical properties of that structure. However, this project will help the engineering community understand those benefits and the basic structural strength values to apply when incorporating NAWAStitch material into composite structure that utilizes a rapid-cure resin system, which is novel. Testing showed that by adding NAWAStitch to the prepreg, most material properties were improved because the VACNTs strengthened the resin to fiber interface, and the interface itself, between plies. For example compression strength improved by more than 20%, and erosion resistance by more than 50%. The UDRI team identified two flight demonstration articles through discussions with Aurora Flight Sciences (Aurora) and reviewing literature on air mobility and the eVTOL market. The demonstration articles were a thin skin honeycomb panel and a section from an eVTOL propeller blade. The manufacturing and testing trials of the demonstration articles were as follows: • NAWAStitch added no manufacturing complexity or detriment to the article physical qualities, while at the same time added multifunctional properties such as enhanced rain erosion and electrothermal de-icing besides enhanced mechanical properties • VACNTs did not increase or hurt the mechanical properties of the interface between honeycomb core materials and the thin skin face sheet • Demonstrated successful de-icing function in an eVTOL propeller meeting real-world specifications in uniformity and endurance • Multifunctional eVTOL propeller cost analysis showed a cost savings of at least 67 % for a single propeller blade with multifunctional de-icing capability and improved mechanical properties which equates to a potential savings of $10,116 per eVTOL aircraft with six three-bladed propellers Overall, the results of this project can provide additional tools for the composite industry to optimize energy efficiency through faster and stronger composite parts that are lower cost to fabricate and maintain. Additionally, the results of this effort can help US industry like NAWA America set the stage for rapid growth and establishment of a manufacturing facility in Ohio. The project has given NAWA America the basic understanding of how the addition of VACNTs to a composite structure can affect the mechanical properties and hence point them to potential commercial applications. The project has also provided Teijin Carbon US with another potential material application for their rapid cure materials that are manufactured at the Renegade facility in Miamisburg, OH.

36 MATERIALS SCIENCE↗
DOE OSTI2024

50 MW Segmented Ultralight Morphing Rotors for Wind Energy

A multi-institutional team designed a 50 Megawatt (MW) rated wind turbine featuring downwind aeroelastic morphing to reduce blade loads and allow an ultralight segmented rotor (mass reduction of about 25% compared to a conventional upwind rotor). The team used a control co-design approach with state-of-the art simulations including for the rotor and tower design using non-linear fluid-structure interactions and control algorithms, and the team also designed, built, and field-tested an aeroelastically-scaled downwind rotor to demonstrate this novel technology and validate the design tool fidelity. In a follow-on phase, these results were used along with an updated Levelized Cost Of Energy (LCOE) methodology to co-design a more detailed set of 25 MW rated turbine designs (including individual pitch control) based on minimum LCOE with highly flexible blades for an Atlantic Ocean offshore fixed-bottom design targeted towards market technology evaluation. A set of upwind and downwind designs at 25 MW rated scales were found to provide the best LCOE, with strong improvements over all previous offshore reference turbines (NREL 5MW, DTU 10 MW and IEA 15 MW). The optimized design would represent the world’s largest offshore turbine design, which combines several state-of-the-art structural, aerodynamic, and control technologies into a new and optimized system concept.

17 WIND ENERGY↗
DOE OSTI2022

On the Rapid Calculation of Binding Affinities for Antigen and Antibody Design and Affinity Maturation Simulations

The accurate and efficient calculation of protein-protein binding affinities is an essential component in antibody and antigen design and optimization, and in computer modeling of antibody affinity maturation. Such calculations remain challenging despite advances in computer hardware and algorithms, primarily because proteins are flexible molecules, and thus, require explicit or implicit incorporation of multiple conformational states into the computational procedure. The astronomical size of the amino acid sequence space further compounds the challenge by requiring predictions to be computed within a short time so that many sequence variants can be tested. In this study, we compare three classes of methods for antibody/antigen (Ab/Ag) binding affinity calculations: (i) a method that relies on the physical separation of the Ab/Ag complex in equilibrium molecular dynamics (MD) simulations, (ii) a collection of 18 scoring functions that act on an ensemble of structures created using homology modeling software, and (iii) methods based on the molecular mechanics-generalized Born surface area (MM-GBSA) energy decomposition, in which the individual contributions of the energy terms are scaled to optimize agreement with the experiment. When applied to a set of 49 antibody mutations in two Ab/HIV gp120 complexes, all of the methods are found to have modest accuracy, with the highest Pearson correlations reaching about 0.6. In particular, the most computationally intensive method, i.e., MD simulation, did not outperform several scoring functions. The optimized energy decomposition methods provided marginally higher accuracy, but at the expense of requiring experimental data for parametrization. Within each method class, we examined the effect of the number of independent computational replicates, i.e., modeled structures or reinitialized MD simulations, on the prediction accuracy. We suggest using about ten modeled structures for scoring methods, and about five simulation replicates for MD simulations as a rule of thumb for obtaining reasonable convergence. We anticipate that our study will be a useful resource for practitioners working to incorporate binding affinity calculations within their protein design and optimization process.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

The emergence of valency in colloidal crystals through electron equivalents

Colloidal crystal engineering of complex, low-symmetry architectures is challenging when isotropic building blocks are assembled. In this study we describe an approach to generating such structures based upon programmable atom equivalents (nanoparticles functionalized with many DNA strands) and mobile electron equivalents (small particles functionalized with a low number of DNA strands complementary to the programmable atom equivalents). Under appropriate conditions, the spatial distribution of the electron equivalents breaks the symmetry of isotropic programmable atom equivalents, akin to the anisotropic distribution of valence electrons or coordination sites around a metal atom, leading to a set of well-defined coordination geometries and access to three new low-symmetry crystalline phases. All three phases represent the first examples of colloidal crystals, with two of them having elemental analogues (body-centred tetragonal and high-pressure gallium), while the third (triple double-gyroid structure) has no known natural equivalent. This approach enables the creation of complex, low-symmetry colloidal crystals that might find use in various technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Covalent Organic Framework (COF) Derived Ni–N–C Catalysts for Electrochemical CO 2 Reduction: Unraveling Fundamental Kinetic and Structural Parameters of the Active Sites

Electrochemical CO 2 reduction is a potential approach to convert CO 2 into valuable chemicals using electricity as feedstock. Abundant and affordable catalyst materials are needed to upscale this process in a sustainable manner. Nickel-nitrogen-doped carbon (Ni-N-C) is an efficient catalyst for CO 2 reduction to CO, and the single-site Ni-N x motif is believed to be the active site. However, critical metrics for its catalytic activity, such as active site density and intrinsic turnover frequency, so far lack systematic discussion. In this work, we prepared a set of covalent organic framework (COF)-derived Ni-N-C catalysts, for which the Ni-N x content could be adjusted by the pyrolysis temperature. The combination of high-angle annular dark-field scanning transmission electron microscopy and extended X-ray absorption fine structure evidenced the presence of Ni single-sites, and quantitative X-ray photoemission addressed the relation between active site density and turnover frequency.

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

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗
DOE OSTI2026

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Unconventional Pathway in the Gas‐Phase Synthesis of 9 H ‐Fluorene (C 13 H 10 ) via the Radical–Radical Reaction of Benzyl (C 7 H 7 ) with Phenyl (C 6 H 5 )

Abstract The simplest polycyclic aromatic hydrocarbon (PAH) carrying a five‐membered ring—9 H ‐fluorene (C 13 H 10 )—is produced isomer‐specifically in the gas phase by reacting benzyl (C 7 H 7 ⋅) with phenyl (C 6 H 5 ⋅) radicals in a pyrolytic reactor coupled with single photon ionization mass spectrometry. The unconventional mechanism of reaction is supported by theoretical calculations, which first produces diphenylmethane and unexpected 1‐(6‐methylenecyclohexa‐2,4‐dienyl)benzene intermediates (C 13 H 12 ) accessed via addition of the phenyl radical to the ortho position of the benzyl radical. These findings offer convincing evidence for molecular mass growth processes defying conventional wisdom that radical‐radical reactions are initiated through recombination at their radical centers. The structure of 9 H ‐fluorene acts as a molecular building block for complex curved nanostructures like fullerenes and nanobowls providing fundamental insights into the hydrocarbon evolution in high temperature settings.

He, Chao↗
DOE OSTI2022

Unconventional Pathway in the Gas-Phase Synthesis of 9 H -Fluorene (C 13 H 10 ) via the Radical–Radical Reaction of Benzyl (C 7 H 7 ) with Phenyl (C 6 H 5 )

The simplest polycyclic aromatic hydrocarbon (PAH) carrying a five-membered ring - 9H-fluorene (C 13 H 10 ) - is produced isomer-specifically in the gas phase by reacting benzyl (C 7 H 7 •) with phenyl (C 6 H 5 •) radicals in a pyrolytic reactor coupled with single photon ionization mass spectrometry. The unconventional mechanism of reaction is supported by theoretical calculations, which first produces diphenylmethane and unexpectedly 1-(6-methylenecyclohexa-2,4-dienyl)benzene intermediates (C 13 H 12 ) accessed via addition of the phenyl radical to the ortho position of the benzyl radical. These findings offer convincing evidence for molecular mass growth processes defying conventional wisdom that radical–radical reactions are initiated through recombination at their radical centers. Furthermore, the structure of 9H-fluorene acts as a molecular building block for complex curved nanostructures like fullerenes and nanobowls providing fundamental insights into the hydrocarbon evolution in high temperature settings.

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

An S=1 Iron(IV) Intermediate Revealed in a Non-Heme Iron Enzyme-Catalyzed Oxidative C—S Bond Formation

Ergothioneine (ESH) and ovothiol A (OSHA) are two natural thiol-histidine derivatives. ESH has been implicated as a longevity vitamin and OSHA inhibits the proliferation of hepatocarcinoma. The key biosynthetic step of ESH and OSHA in the aerobic pathways is the O 2 -dependent C—S bond formation catalyzed by non-heme iron enzymes (e.g., OvoA in ovothiol biosynthesis), but due to the lack of identification of key reactive intermediate the mechanism of this novel reaction is unresolved. In this study, we report the identification and characterization of a kinetically competent S=1 iron(IV) intermediate supported by a four-histidine ligand environment (three from the protein residues and one from the substrate) in enabling C—S bond formation in OvoA from Methyloversatilis thermotoleran, which represents the first experimentally observed intermediate spin iron(IV) species in non-heme iron enzymes. Results reported in this study thus set the stage to further dissect the mechanism of enzymatic oxidative C—S bond formation in the OSHA biosynthesis pathway. They also afford new opportunities to study the structure-function relationship of high-valent iron intermediates supported by a histidine rich ligand environment.

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

Microstructure Clones

Background: A material’s microstructure drives its material performance. Contemporary crystal plasticity experiments compare full-field strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which features of the observed deformation are deterministic vs statistical; thus, differences between model and experiment may or may not be significant. Objective: This paper introduces the invention of microstructure clones. Microstructure clones are 2D oligocrystal specimens that have nearly identical microstructures to remedy the aforementioned experimental limitations. Having specimens with nearly identical microstructures will allow for multiple destructive tests of a microstructure (either as repeats or intentionally different experiments), an ability to “see the future” by providing insight into how a specimen will deform, variability quantification, and experimental investigations of response to small microstructural changes. Methods: This work introduces microstructure clones. Repeatability of these clones is demonstrated in tensile bars of pure nickel. Local strain measurements from digital image correlation are compared between clone specimens and compared to results from a crystal plasticity finite element model. Results: Two sets of microstructure clones were tested in this study and displayed very consistent deformation responses within each clone set. Small observed differences in deformation invite investigation into microstructure stochasticity and the effect of small microstructural and loading differences. Conclusions: Microstructure clones represent a significant shift in understanding structure–property relationships. This work reshapes experimental crystal plasticity to allow for experiments that control for specific variables, quantification of microstructural stochasticity (and other sources of stochasticity), and opportunities for replicating experiments.

Crystal Plasticity↗
DOE OSTI2023

Extrusion parameter control optimization for DIW 3D printing using image analysis techniques

Material extrusion is a well-recognized facet of additive manufacturing that involves the fabrication of parts through the deposition of structural material from an extrusion head from a bulk supply. In the subdivision of Direct Ink Writing (DIW) additive manufacturing, challenges arise when the structural material is flowable, synchronous extrusion control and tool movement becomes critical for achieving high-quality parts with low defect populations. DIW techniques are most used in laboratory settings using expensive custom instruments and may require specialized 3D slicing software. Here, in this study, the fabrication of an inexpensive, consumer-friendly progressive cavity pump dispensing system is detailed, in which can create high-quality parts by executing G-code commands produced from a commercial slicing software. The precision and repeatability of the movement-synchronized material extrusion is demonstrated through a series of optimization schemes, entailing the alteration of various control parameters, which directly affect the extrusion properties demonstrated during a print. In situ diagnostics were implemented to evaluate the results of the established optimization experiment. Using a machine vision technique, images of the optimization prints are processed. Following this, a supervised machine learning model was trained to autonomously judge whether or not the extrusion parameters produced a passing or failing result. The machine learning scheme serves as a preliminary benchmark for future layer-by-layer evaluation of more complex DIW parts. The construction of the printer and development of in situ characterization capabilities demonstrates the ability for this printer to create high-fidelity DIW parts for a fraction of the price of other systems.

42 ENGINEERING↗
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2024

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystem function and soil quality. Salinity influences carbon cycling through direct effects on microbial activity and indirect alterations to soil physicochemical properties including cation exchange, pH, and soil organic carbon availability. Current models inadequately represent these complexities, relying on linear reduction functions that overlook specific physicochemical changes induced by salinity. To address this gap, we propose an integrated model framework, AquaMEND, that combines microbial-explicit carbon decomposition and geochemical models. This model allows cation exchange and surface complexation processes to capture solute chemistry and nutrient availability in soils upon saltwater intrusion. Using response functions that capture salinity impacts on both salt-sensitive and salt-resistant microbial processes, AquaMEND simulates how the abiotic and biotic mechanisms work individually and collectively to regulate organic and inorganic pools and fluxes. Here, the parallel structure of aqueous and non-aqueous phases, together with microbial functions, result in a versatile model for solving dynamic coupling of organics, minerals and microbes under various environmental settings.

54 ENVIRONMENTAL SCIENCES↗