Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Chemical solutions”

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 487 records · Page 27

Support for ISU's 2003 Summer Session Program and the Theme Day on "Upcoming Technologies for Space"

Electric propulsion, whether powered by solar generators or by nuclear reactors, offers a valuable solution to the problems posed by chemical in-space propulsion. It maybe argued that solar electric propulsion is more applicable to near-Earth missions or to missions to the nearer planets, whereas nuclear electric propulsion remains operational even in the vicinity of the far outer planets. Taking into account considerations such as cost and safety it can be concluded that both types of electric propulsion are complementary having their specific niches of application. A third technology based on non-chemical sources for high efficiency propulsion is solar thermal propulsion, which may find its own specific niche in near-Earth missions or missions to the inner planets.

Finarelli, Margaret G.↗

Electrochemical Hydrogen Peroxide Generator

Two-electron reduction of oxygen to produce hydrogen peroxide is a much researched topic. Most of the work has been done in the production of hydrogen peroxide in basic media, in order to address the needs of the pulp and paper industry. However, peroxides under alkaline conditions show poor stabilities and are not useful in disinfection applications. There is a need to design electrocatalysts that are stable and provide good current and energy efficiencies to produce hydrogen peroxide under acidic conditions. The innovation focuses on the in situ generation of hydrogen peroxide using an electrochemical cell having a gas diffusion electrode as the cathode (electrode connected to the negative pole of the power supply) and a platinized titanium anode. The cathode and anode compartments are separated by a readily available cation-exchange membrane (Nafion 117). The anode compartment is fed with deionized water. Generation of oxygen is the anode reaction. Protons from the anode compartment are transferred across the cation-exchange membrane to the cathode compartment by electrostatic attraction towards the negatively charged electrode. The cathode compartment is fed with oxygen. Here, hydrogen peroxide is generated by the reduction of oxygen. Water may also be generated in the cathode. A small amount of water is also transported across the membrane along with hydrated protons transported across the membrane. Generally, each proton is hydrated with 3-5 molecules. The process is unique because hydrogen peroxide is formed as a high-purity aqueous solution. Since there are no hazardous chemicals or liquids used in the process, the disinfection product can be applied directly to water, before entering a water filtration unit to disinfect the incoming water and to prevent the build up of heterotrophic bacteria, for example, in carbon based filters. The competitive advantages of this process are: 1. No consumable chemicals are needed in the process. The only raw materials needed are water and oxygen or air. 2. The product is pure and can therefore be used in disinfection applications directly or after proper dilution with water. 3. Oxygen generated in the anode compartment is used in the electrochemical reduction process; in addition, external oxygen is used to establish a high flow rate in the cathode compartment to remove the desired product efficiently. Exiting oxygen can be recycled after separation of liquid hydrogen peroxide product, if so desired. 4. The process can be designed for peroxide generation under microgravity conditions. 5. High concentrations of the order of 6-7 wt% can be generated by this method. This method at the time of this reporting is superior to what other researchers have reported. 6. The cell design allows for stacking of cells to increase the hydrogen peroxide production. 7. The catalyst mix containing a diquaternary ammonium compound enabled not only higher concentration of hydrogen peroxide but also higher current efficiency, improved energy efficiency, and catalyst stability. 8. The activity of the catalyst is maintained even after repeated periods of system shutdown. 9. The catalyst system can be extended for fuel-cell cathodes with suitable modifications.

Tennakoon, Charles L. K.↗

Process for Encapsulating Protein Crystals

A process for growing protein crystals encapsulated within membranes has been invented. This process begins with the encapsulation of a nearly saturated aqueous protein solution inside semipermeable membranes to form microcapsules. The encapsulation is effected by use of special formulations of a dissolved protein and a surfactant in an aqueous first liquid phase, which is placed into contact with a second, immiscible liquid phase that contains one or more polymers that are insoluble in the first phase. The second phase becomes formed into the semipermeable membranes that surround microglobules of the first phase, thereby forming the microcapsules. Once formed, the microcapsules are then dehydrated osmotically by exposure to a concentrated salt or polymer solution. The dehydration forms supersaturated solutions inside the microcapsules, thereby enabling nucleation and growth of protein crystals inside the microcapsules. By suitable formulation of the polymer or salt solution and of other physical and chemical parameters, one can control the rate of transport of water out of the microcapsules through the membranes and thereby create physicochemical conditions that favor the growth, within each microcapsule, of one or a few large crystals suitable for analysis by x-ray diffraction. The membrane polymer can be formulated to consist of low-molecular-weight molecules that do not interfere with the x-ray diffraction analysis of the encapsulated crystals. During dehydration, an electrostatic field can be applied to exert additional control over the rate of dehydration. This protein-crystal-encapsulation process is expected to constitute the basis of protein-growth experiments to be performed on the space shuttle and the International Space Station. As envisioned, the experiments would involve the exposure of immiscible liquids to each other in sequences of steps under microgravitational conditions. The experiments are expected to contribute to knowledge of the precise conditions under which protein crystals form. By enhancing the ability to grow crystals suitable for x-ray diffraction analysis, this knowledge can be expected to benefit not only the space program but also medicine and the pharmaceutical industry.

Morrison, Dennis R.↗

Environmental Barrier Coatings for Ceramic Matrix Composites - An Overview

SiC/SiC Ceramic Matrix Composites (CMCs) are increasingly being considered as structural materials for advanced power generation equipment because of their light weight, higher temperature capability, and oxidation resistance. Limitations of SiC/SiC CMCs include surface recession and component cracking and associated chemical changes in the CMC. The solutions pursued to improve the life of SiC/SiC CMCs include the incorporation of coating systems that provide surface protection, which has become known as an Environmental Barrier Coating (EBC). The development of EBCs for the protection of gas turbine hot section CMC components was a continuation of coating development work for corrosion protection of silicon-based monolithics. Work on EBC development for SiC/SiC CMCs has been ongoing at several national laboratories and the original gas turbine equipment manufacturers. The work includes extensive laboratory, rig and engine testing, including testing of EBC coated SiC/SiC CMCs in actual field applications. Another EBC degradation issue which is especially critical for CMC components used in aircraft engines is the degradation from glassy deposits of calcium-magnesium-aluminosilicate (CMAS) with other minor oxides. This paper addresses the need for and properties of external coatings on SiC/SiC CMCs to extend their useful life in service and the retention of their properties.

EBC↗

Enabling Mission Flexibility to Battery Driven Deep Space Endeavors With Generalized Battery-Health-Monitoring Using Physics-Based and Data-Driven Reduced-Order Models

The needs and requirements for an electrochemical energy storage for deep space exploration is well explored. It is often understood that different mission sites and environmental conditions require different battery chemistries or technologies. Additionally, various engineering solutions are deployed to overcome specific chemical challenges. One often overlooked need is the “health” monitoring of an electrochemical storage system. The term generalized health monitoring, as envisioned in this work, refers to the monitoring of various aspects such as electrode health, electrolyte health, reaction pathway health, cooling system health, sensor health, and BMS health [1]. Generalized health monitoring allows mission leads, engineers, and scientists to incorporate flexibility in mission designs, make on-the-fly mission changes, and extend the duration of science missions. Moreover, it enables automation and data-driven decision-making without compromising safety and performance. Recently, our group developed a hierarchy of thermal reduced-order models (TROM) by combining a physics-based modeling approach and data-driven model reduction techniques applied to flight data [2]. The resulting TROMs were found to be not only accurate but also identifiable from the flight data. Consequently, the coefficient of variance of the model parameters is small over the course of hundreds of flights, allowing for monitoring the parameter evolution trajectories as the battery ages and degrades. These parameters constitute the metrics of the generalized health of a battery. Monitoring their evolution allows such models to be used for anomaly detection and prognostics, improving early detection of abnormal behavior and thus enabling timely maintenance, longer battery life, and enhanced battery safety. For this presentation, the practicality of the thermal model will be validated on a pack of 14cells under various topology configurations such as 1S14P, 2P7S, 7S2P, and 1P14S. It is well known that manufacturing and non-uniform aging lead to variability in the performance of a cell, which is exacerbated by cell balancing during active load. Additionally, in extreme scenarios, the paramount objective is to complete the mission, regardless of the stresses on the battery. Topology-induced balancing issues further stress the battery. The goal of this study is to determine if the noise (identifiability) in the reduced-order thermal model parameters is sensitive to topology, cell spacing, cooling strategy, and manufacturing or age variability. The variability in cells is considered by assuming a multimodal distribution for microscopic parameters of a cell (such as porosity, tortuosity, reaction kinetics, volumetric thermal conductivity, and volumetric heat capacity). The compounded effect of manufacturing variability, topological selection, cooling strategies, and cell balancing ages each cell in a battery differently. The study aims to clarify whether the challenge in extracting maximum information depends on the minimum number of sensors or models used for data extraction.

Automation↗

Compact Superconducting RF Accelerators for Water Treatment

Jefferson Lab is actively exploring the environmental applications of compact irradiation facilities utilizing Superconducting Radio Frequency (SRF) accelerator technology ? a well-established technology primarily used in large research machines. Recent advancements, such as the application of thin films of superconducting alloys in radio frequency (RF) resonators and the development of high-capacity cryocoolers, have made it feasible to employ this key technology in compact, standalone irradiation facilities. The main advantage over conventional systems based on normal conducting technology is the significantly enhanced efficiency in beam power generation, enabling exposure to substantially higher radiation doses. In the realm of addressing contaminants in wastewater treatment, particularly concerning persistent substances known as "forever chemicals," this technology presents an attractive solution, given the required considerable effort i.e. dose to break them down. Experiment Numbers:

Vennekate, John↗

Enhancement of chemical oscillations by self-generated convective flows

Chemical feedback loops in fluids can produce not only chemical oscillations, but also density variations that generate solutal buoyancy forces, which in turn initiate fluid flow. Using analytical and computational models, we herein examine how the reaction-induced flows alter the chemical oscillations in a fluid-filled chamber whose top and bottom walls are coated with different enzymes. Due to this chemo-fluidic coupling, the systems form oscillating flow patterns, which combine the characteristic size of the buoyancy-driven convection rolls with the frequency of the chemical oscillations. With changes in the distance between the enzyme-coated walls, the convective flows not only enhance or suppress the chemical oscillations, but also substantially increase the amplitude and frequency of the oscillations and extend the regime of the oscillatory behavior. These design principles can facilitate the development of artificial biochemical networks that act as chemical clocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Challenges and Opportunities in Decarbonizing the U.S. Energy System

The United States has pledged to develop a 100% carbon-free electric power system by 2035 and a net-zero-emissions economy by 2050. While important advancements have been made in the scale, performance, and economics of clean energy technologies, meeting the nation's ambitious goals will not only require their deployment at scale, but also additional innovation and effective integration of different solutions. Technological developments across the broad suite of low-carbon energy solutions are advancing rapidly, with ongoing innovations in renewable electricity generation, industrial processes, and energy-saving technologies and services, including LED lighting, induction heating, electric vehicles, energy storage solutions, and mobility as a service, plus smart devices, controls, and more efficient and smart buildings. Combining renewable electricity with biotic and abiotic pathways to produce chemicals, fuels, and materials promises to deliver new solutions. Grid-interactive buildings and communities, integrating transportation infrastructure and vehicles, are likely to be significant components of any zero-carbon energy strategy. Low-carbon industrial manufacturing will also make strong contributions to a net-zero economy. While the technical prospects appear promising, variations in the state of infrastructure, jurisdictional and social equity, pollution, economic and socio-cultural constraints, energy resource availability, and supply chain dynamics found in different locations present a range of challenges and demand customized solutions. This paper provides a critical review and offers new insights into the technical, infrastructure, analytic, political, and economic challenges faced in translating the nation's ambitious net-zero-emissions goals into feasible and reliable implementation action plans.

circular economy↗

Colloidal Behavior of Plutonium Oxide in Concentrated Electrolyte Solutions

The Hanford Site in Washington State manages legacy high-level radioactive waste streams that display major chemistry and engineering challenges, including the high salt levels and pH values that correspond to conditions under which many classical concepts describing chemical reactivity cannot be applied. One particular challenge that needs to be tackled at Hanford is that Pu concentrations, [Pu], in the soluble phases of the tank wastes are higher than expected based on the solubility of crystalline PuO2, which is widely accepted to be caused by the formation of PuO2 colloid, consisting of nano- to submicron-sized particles (PuO2 NPs). Fundamental research underpinning the behavior of PuO2 NPs under conditions not only relevant to the Hanford tank waste but at high ionic strength in general is needed to reliably predict the chemical reactivity of PuO2 NPs and develop engineering solutions to safely and efficiently process high-level radioactive waste into forms suitable for long-term storage. In this work, we study the behavior of PuO2 NPs (particle size ~100 nm) under high ionic strength conditions by reacting it with highly concentrated (up to 5 M) salt solutions. We explore different electrolyte compositions to elucidate the impact of different anions (NO3-, Cl-, ClO4-, SO42-, C2O42-, CO32-) on the stability of PuO2 NP in the acidic and alkaline pH regime. PuO2 NP aggregation and precipitation as function electrolyte concentration is tracked by a combination of liquid scintillation counting, dynamic light scattering for determination of particle size distributions, and zeta potentials as a proxy for particle charge. At acidic pH, electrolytes containing non-coordinating anions, such as NaNO3, NaCl, and NaClO4 mostly stabilize PuO2 NPs over a large electrolyte concentration range, showing only subtle differences in their reactivity. Other electrolyte anions show a more pronounced effect on the PuO2 NP stability: SO42-, binds directly to the particles’ surface, reverses the particle charge, and precipitates the PuO2 NPs efficiently even at intermediate sulfate concentrations (>0.1 M). In contrast, C2O42- is found to lead to high [Pu] in solution, in the milli-molar range, even at mildly acidic pH (~4). Thermodynamic modeling of the dissolved Pu concentrations using PHREEQC is unable to predict the observed [Pu] in the acidic pH regime, supporting the influence of colloids in maintaining elevated [Pu]. It is noteworthy that the current thermodynamic databases do not include constants for colloidal Pu phases and cannot accurately predict many of the high ionic strength solutions relevant to this work. The mechanisms and models responsible for these observations will need further investigation in the future. At high pH values (~12), PuO2 NPs exhibits classical sol-gel chemistry, meaning that upon destabilization of the colloidal sol, for example by addition of concentrated NaOH, highly porous and viscid PuO2 coagulates are formed that consist of a three-dimensional network likely held together by physical interactions. The PuO2 NP coagulate shows no significant reversibility of the aggregation when contacted with concentrated brines; however, PuO2 NPs can be efficiently resuspended in solution by addition of diluted electrolytes, alkaline solutions containing high amounts of carbonate, or simple addition of water. Especially carbonate is shown to stabilize PuO2 NPs in solution at high pH, characterized by stable colloidal suspensions that are resistant against sedimentation during centrifugation. Thermodynamic modeling of the carbonate system was able to predict an increasing dissolved Pu concentration with increasing carbonate concentration. However, the model was profoundly sensitive to the fixed redox potential and does not include any thermodynamic constants for colloidal Pu species.

Neumann, Julia↗

From suppressed void growth to significant void swelling in NiCoFeCr complex concentrated solid-solution alloy

Void swelling can result in dimensional instability and undermine the safe operation of nuclear reactors. Current strategies to inhibit void swelling mainly focus on enhancing defect absorption and recombination by introducing high-density defect sinks. Complex concentrated solid-solution alloys (CSAs), including high-entropy alloys, can withstand severe radiation damage due to their inherent chemical complexity without interfaces. However, the underlying mechanisms for void suppression in CSAs are far from clear. In this research, we studied the void evolution with respect to irradiation depths, doses, and temperatures in equiatomic NiCoFeCr under 3 MeV Ni ion irradiations. At relatively low doses (16 and 54 displacements per atom, dpa), voids form mainly outside of the ion-damaged region, and void formation in the peak damage region is suppressed, leading to negligible swelling. However, with further increase of dose (86 up to 250 dpa), significant void growth occurs in the peak damage region and extended dislocation lines dominate instead of short dislocation lines and loops formed at lower doses. From 500 to 700 °C, the dislocation density decreases while dislocations grow. Although the overall void swelling increases dramatically from 500 to 580 °C at 54 dpa, void growth in the peak damage region is still suppressed. The transition from suppressed void growth to significant void swelling is attributed to dislocation evolution and local chemical inhomogeneity (enrichment of Fe/Cr in the matrix) at higher doses. Our study shows that controlling element diffusion and defect evolution through tuning chemical complexity can further enhance the swelling resistance of CSAs.

36 MATERIALS SCIENCE↗

Astrochemical Bistability: Autocatalysis in Oxygen Chemistry

The origin of bistable solutions in the kinetic equations describing the chemistry of dense interstellar clouds is explained as being due to the autocatalysis and feedback of oxygen nuclei from the oxygen dimer (O2).We identify four autocatalytic processes that can operate in dense molecular clouds, driven respectively by reactions of H+, He+, C+, and S+ with O2. We show that these processes can produce the bistable solutions found in previous studies, as well as the dependence on various model parameters such as the helium ionization rate, the sulfur depletion and the H+3 electron recombination rate. We also show that iongrain neutralizations are unlikely to affect the occurrence of bistability in dense clouds. It is pointed out that many chemical models of astronomical sources should have the potential to show bistable solutions. It is pointed out that many chemical models of astronomical sources should have the potential to show bistable solutions.

Gwenaelle Dufour↗

Direct Optical Lithography of CsPbX 3 Nanocrystals via Photoinduced Ligand Cleavage with Postpatterning Chemical Modification and Electronic Coupling

Microscale patterning of solution-processed nanomaterials is important for integration in functional devices. Colloidal lead halide perovskite (LHP) nanocrystals (NCs) can be particularly challenging to pattern due to their incompatibility with polar solvents and lability of surface ligands. Here, we introduce a direct photopatterning approach for LHP NCs through the binding and subsequent cleavage of a photosensitive oxime sulfonate ester (-C=N-OSOO-). The photosensitizer binds to the NCs through its sulfonate group and is cleaved at the N-O bond during photoirradiation with 405 nm light. This bond cleavage decreases the solubility of the NCs, which allows patterns to emerge upon development with toluene. Postpatterning ligand exchange results in photoluminescence quantum yields of up to 79%, while anion exchange provides tunability in the emission wavelength. Finally, the patterned NC films show photoconductive behavior, demonstrating that good electrical contact between the NCs can be established.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

CAG12 - A CSCM based procedure for flow of an equilibrium chemically reacting gas

The Conservative Supra Characteristic Method (CSCM), an implicit upwind Navier-Stokes algorithm, is extended to the numerical simulation of flows in chemical equilibrium. The resulting computer code known as Chemistry and Gasdynamics Implicit - Version 2 (CAG12) is described. First-order accurate results are presented for inviscid and viscous Mach 20 flows of air past a hemisphere-cylinder. The solution procedure captures the bow shock in a chemically reacting gas, a technique that is needed for simulating high altitude, rarefied flows. In an initial effort to validate the code, the inviscid results are compared with published gasdynamic and chemistry solutions and satisfactorily agreement is obtained.

Green, M. J.↗

Direct Observations of Solute Dispersion in Rocks With Distinct Degree of Sub‐Micron Porosity

Abstract The transport of chemical species in rocks is affected by their structural heterogeneity to yield a wide spectrum of local solute concentrations. To quantify such imperfect mixing, advanced methodologies are needed that augment the traditional breakthrough curve analysis by probing solute concentration within the fluids locally. Here, we demonstrate the application of asynchronous, multimodality imaging by X‐ray computed tomography (XCT) and positron emission tomography (PET) to the study of passive tracer experiments in laboratory rock cores. The four‐dimensional concentration maps measured by PET reveal specific signatures of the transport process, which we have quantified using fundamental measures of mixing and spreading. We observe that the extent of solute spreading correlate strongly with the strength of subcore‐scale porosity heterogeneity measured by XCT, while dilution is enhanced in rocks containing substantial sub‐micron porosity. We observe that the analysis of different metrics is necessary, as they can differ in their sensitivity to the strength and forms of heterogeneity. The multimodality imaging approach is uniquely suited to probe the fundamental difference between spreading and mixing in heterogeneous media. We propose that when multi‐dimensional data is available, mixing and spreading can be independently quantified using the same metric. We also demonstrate that one‐dimensional transport models have limited predictive ability toward the internal evolution of the solute concentration, when the model is solely calibrated against the effluent breakthrough curves. The data set generated in this study can be used to build realistic digital rock models and to benchmark transport simulations that account deterministically for rock property heterogeneity.

Kurotori, Takeshi [Department of Chemical Engineer↗

Investigation of the Solution Chemistry of Hybrid Organic–Inorganic Indium Halides for New Material Discovery

Recently, metal halide perovskites (MHPs) have emerged as a new class of materials for optical and electronic applications such as solar cells and ionizing radiation detectors. Although the solution-processability of MHPs is among their greatest advantages, the solution chemistries of most metal halide systems and their relationship with the observed structural and chemical diversity are poorly understood. In this work, we study the solution chemistry of a model indium halide system, methylammonium (MA)–In–Br, using a combination of the UV–vis spectroscopy, electrospray ionization mass spectrometry (ESI-MS) measurements, small-angle X-ray scattering (SAXS), and density functional theory (DFT) calculations. Our results show that indium could form either octahedral [InBr 6 3– ] or tetrahedral [InBr 4 – ] anions in solution or a combination of both, depending on the loading ratios of MABr and InBr 3 reactants. Understanding the solution chemistry of this system and recognizing the optical fingerprints of these polyanions allow for targeted crystallization of two novel compounds: MAInBr 4 featuring tetrahedral [InBr 4 – ] anions and MA 2 InBr 5 containing both octahedral [InBr 6 3– ] and tetrahedral [InBr 4 – ] anions. Further increase of the MABr content leads to the formation of previously reported MA 4 InBr 7 , containing only octahedral [InBr 6 3– ] anions separated by Br – anions. Furthermore, our results suggest that understanding the solution chemistry of multinary metal halide systems could be a valuable tool for discovering functional materials for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical simulation of a low-emission gas turbine combustor using KIVA-II

A modified version of the KIVA-II code was used to obtain a multidimensional numerical solution for the turbulent two-phase chemically reacting flows inside a staged turbine combustor (STC). The STC under consideration is equipped with an advanced airblast fuel nozzle and encompasses a fuel nozzle (FN), a rich-burn (RB) zone, a converging connecting pipe, a quick-quench (QQ) zone, a diverging connecting pipe, and a lean-combustion (LC) zone. The STC was divided into two subsystems, namely, FN/RB zone and QQ/LC zones, and the numerical solutions were obtained separately for each subsystem. Preliminary data characterize the major features of the flow and temperature fields inside the STC. Information on velocity, temperature, and some critical species in the FN/RB zone is presented. In the QQ/LC zones, formation of the co- and counter-rotating bulk flow and the sandwiched-ring-shaped temperature field can be clearly seen. The calculations of the mass-weighted standard deviation and the pattern factor of temperature indicated that the mixing performance of the STC is very promising.

Yang, S. L.↗

A chemical kinetics theory for interpreting the non-monotonic temperature dependence of enzymatic reactions

Abstract. One notable observation of enzymatic chemical reactions is that, for a given abundance of enzymes and substrates, temperature increases cause reaction rates to first increase consistent with the Arrhenius relationship, then plateau, and finally fall off quickly to zero at high temperatures. While many mathematical functions have been used to describe this pattern, we here propose a chemical kinetics theory which successfully replicates this observation and provides insights into the processes responsible for these dynamics. The chemical kinetics theory combines the law of mass action, von Smoluchowski's diffusion-limited chemical reaction theory, and Eyring's transition state theory. This new theory reveals that the thermally reversible enzyme denaturation ensured by the ceaseless thermal motion of molecules and ions in an enzyme solution explains the plateau and subsequent decrease in chemical reaction rates with increasing temperature. The temperature-dependent affinity parameter (K) that relates enzymes and substrates through their binding also affects the shape of the emergent temperature response. We demonstrate that with an increase in substrate availability, K shifts the optimal temperature, where reaction rates plateau, towards higher values. Further, we show that the chemical kinetics theory accurately represents 12 sets of published enzyme assay data and includes the popular mechanistic model by Ratkowsky et al. (2005) as a special case. Given its good performance and solid theoretical underpinning, we believe this new theory will facilitate the construction of more mechanistic-based environmental biogeochemical models.

Tang, Jinyun↗