Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 469 records · Page 26

Model of Inclusion Evolution During Calcium Treatment in the Ladle Furnace

Calcium treatment of steel is typically employed to modify alumina inclusions to liquid calcium aluminates. However, injected calcium also reacts with the dissolved sulfur to form calcium sulfide. The current work aims to develop a kinetic model for the evolution of oxide and sulfide inclusions in Al-killed alloyed steel during Ca treatment in the ladle refining process. The model considers dissolution of the calcium from the calcium bubbles into the steel and reduction of calcium oxide in the slag to dissolved calcium. A steel–inclusion kinetic model is used for mass transfer to the inclusion interface and diffusion within the calcium aluminate phases formed on the inclusion. The inclusion–steel kinetic model is then coupled with a previously developed steel–slag kinetic model. The coupled inclusion–steel–slag kinetic model is applied to the chemical composition changes in molten steel, slag, and evolution of inclusions in the ladle. The result of calculations is found to agree well with an industrial heat for species in the steel as well as inclusions during Ca treatment.

36 MATERIALS SCIENCE↗

Modeling the Kinetics of Deactivation of Catalysts during the Upgrading of Bio-Oil

The fouling of catalysts for the upgrading of bio-oils appears to be very different from the fouling of catalysts for the hydroprocessing of petroleum-derived streams. There are two reasons for the differences: a) bio-oil contains polarizable components and phases that can stabilize reaction intermediates exhibiting charge separation and b) bio-oil components contain functional groups that contain O, notably carbonyls (>C=O). Aldol condensation of carbonyls affords very different pathways for the production of oligomeric, refractory deposits than does dehydrogenation/polymerization of petroleum-derived hydrocarbons. Colloquially, we refer to the bio-oil derived deposits as “gunk” to discriminate them from coke, the carbonaceous deposits encountered in petroleum refining. Classical gelation, appears to be a suitable model for the “gunking” reaction. Our work has helped explain the temperature range at which bio-oil should be pre-processed (“stabilized”) to confer longer lifetimes on the catalysts used for more severe processing. Stochastic modeling (kinetic Monte Carlo simulations) appears suitable to capture the rates of oligomerization of bio-oil. This work was supported by the US Department of Energy, Office of Energy Efficiency and Renewable Energy, Bioenergy Technologies Office. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle.

bio-oil upgrading↗

Kinetic and modeling studies of the mechanism of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2

Since its discovery over 15 years ago, the reversible dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 has remained one of the more intriguing hydrogen-cycling systems. While the mechanism of this reaction has been the subject of a good deal of speculation and computational studies, prior to this work it had not been probed through kinetic studies. Previous reports of the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 have not included kinetic studies. The present studies have shown that the dehydrogenation of Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 is suppressed by hydrogen pressure indicating that the rate-limiting step in this process involves hydrogen elimination. Computational modeling of kinetic data obtained from monitoring the hydrogen elimination from Mg(BH 4 ) 2 to Mg(B 3 H 8 ) 2 under static vacuum over a range of temperatures supports that the dehydrogenation occurs through a reversible three-step process in which the elimination of hydrogen from the [B 3 H 10 ] − intermediate is rate limiting. A mechanism involving the low energy transfer of neighboring BH3 groups is proposed to account for the formation of [B 3 H 8 ] − at relatively low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of Magnetite Oxidation Kinetics at the Particle Scale

The induration of magnetite pellets is a complex physico-chemical process that involves oxidation, sintering, and heat transfer. The thermal- and gas-composition profile that is experienced by the pellet in an induration reactor could result in the formation of a homogenous or heterogeneous pellet structure, which could affect the pellet quality. The oxidation kinetics of magnetite pellets from sintering studies have been studied at two levels, namely, the pellet scale and at the particle scale, and the findings of the latter are presented here. The rate of oxidation of the magnetite concentrate depends primarily on temperature, oxygen content in the oxidizing gas, and particle size. These factors are investigated in this study. It was found that the oxidation of the magnetite concentrate is comprised of two distinct stages, a primary stage with high rates followed by a secondary stage where rates decrease significantly. The isothermal oxidation behavior as analyzed by the Avrami kinetic model was found to fit better than the shrinking-core model. The partially oxidized particles were examined microstructurally to supplement the experimental and model results. The Avrami kinetic model for isothermal oxidation was extended to non-isothermal profiles using the superposition principle, and the model was validated experimentally.

36 MATERIALS SCIENCE↗

Homogenizing Advanced Alloys: Thermodynamic and Kinetic Simulations Followed by Experimental Results

Segregation of solute elements occurs in nearly all metal alloys during solidification. The resultant elemental partitioning can severely degrade as-cast material properties and lead to difficulties during post-processing (e.g., hot shorts and incipient melting). Many cast articles are subjected to a homogenization heat treatment in order to minimize segregation and improve their performance. Traditionally, homogenization heat treatments are based upon past practice or time-consuming trial and error experiments. Through the use of thermodynamic and kinetic modeling software, NETL has designed a systematic method to optimize homogenization heat treatments. Use of the method allows engineers and researchers to homogenize casting chemistries to levels appropriate for a given application. The method also allows for the adjustment of heat treatment schedules to fit limitations on in-house equipment (capability, reliability, etc.) while maintaining clear numeric targets for segregation reduction. In this approach, the Scheil module within Thermo-Calc is used to predict the as-cast segregation present within an alloy, and then diffusion controlled transformations is used to model homogenization kinetics as a function of time and temperature. Examples of computationally designed heat treatments and verification of their effects on segregation and properties of real castings are presented.

36 MATERIALS SCIENCE↗

The effect of model fidelity on prediction of char burnout for single-particle coal combustion

In this study, practical simulation of industrial-scale coal combustion relies on the ability to accurately capture the dynamics of coal subprocesses while also ensuring the computational cost remains reasonable. The majority of the residence time occurs post-devolatilization, so it is of great importance that a balance between the computational efficiency and accuracy of char combustion models is carefully considered. In this work, we consider the importance of model fidelity during char combustion by comparing combinations of simple and complex gas and particle-phase chemistry models. Detailed kinetics based on the GRI 3.0 mechanism and infinitely-fast chemistry are considered in the gas-phase. The Char Conversion Kinetics model and nth-Order Langmuir–Hinshelwood model are considered for char consumption. For devolatilization, the Chemical Percolation and Devolatilization and Kobayashi-Sarofim models are employed. The relative importance of gasification versus oxidation reactions in air and oxyfuel environments is also examined for various coal types. Results are compared to previously published experimental data collected under laminar, single-particle conditions. Calculated particle temperature histories are strongly dependent on the choice of gas phase and char chemistry models, but only weakly dependent on the chosen devolatilization model. Particle mass calculations were found to be very sensitive to the choice of devolatilization model, but only somewhat sensitive to the choice of gas chemistry and char chemistry models. High-fidelity models for devolatilization generally resulted in particle temperature and mass calculations that were closer to experimentally observed values.

01 COAL, LIGNITE, AND PEAT↗

Kinetic analyses for solid-state phase transition of metastable amorphous-AlO x (2.5 < x ≤ 3.0) nanostructures into crystalline alumina polymorphs

Solid-solid phase change materials (SS-PCMs) hold promise for energy storage/dissipation in batteries and energetic materials. Yet, phase change kinetics for SS-PCMs undergoing metastable to semi-stable/stable phase transformations remain relatively ill-studied because trapping metastable phases remain challenging. Recently, we demonstrated the kinetic entrapment and stabilization of a highly disordered and amorphous Al-oxide phase m-AlO x @C (x~2.5-3.0) via laser ablation synthesis in solution (LASiS). We report here, to our knowledge, the first chemical kinetics analysis for S-S phase transition of the m-AlO 3 @C nanocomposites (< 5–8 nm sizes) into semi-stable equilibrium alumina phases (θ/γ-Al 2 O 3 ) via disproportionation reaction, while releasing excess trapped gases. Our results indicate the atomic density of the AlO 3 structures to be ~5–10 times less than that of the final Al 2 O 3 phases, which led to the hypothesis of a volume shrinkage process during their phase transition. Temperature-dependent X-ray diffraction studies reveal the high-temperature phase transition for m-AlO 3 → θ/γ-Al 2 O 3 to follow contracting volume kinetics model, thereby validating our earlier hypothesis. Using the geometric volume contraction model, reaction kinetics analyses from Arrhenius plots reveal the activation energy barrier for the phase transition to be ~270±11 kJ/mol. This makes the activation energy barrier nearly identical to the oxidation of micron-sized Al particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurements of population densities of metastable and resonant levels of argon using laser induced fluorescence

We present a new approach to measure population densities of Ar I metastable and resonant excited states in low temperature Ar plasmas at pressures higher than 1 Torr. This approach combines the time resolved laser induced fluorescence technique with the kinetic model of Ar. The kinetic model of Ar is based on calculating the population rates of metastable and resonant levels by including contributions from the processes that affect population densities of Ar I excited states. In particular, we included collisional quenching processes between atoms in the ground state and excited states, since we are investigating plasma at higher pressures. We also determined time resolved population densities of Ar I 2 p excited states by employing optical emission spectroscopy technique. Time resolved Ar I excited state populations are presented for the case of the post-discharge of the supersonic flowing microwave discharge at pressures of 1.7 and 2.3 Torr. The experimental set-up consists of a pulsed tunable dye laser operating in the near infrared region and a cylindrical resonance cavity operating in TE{sub 111} mode at 2.45 GHz. Results show that time resolved population densities of Ar I metastable and resonant states oscillate with twice the frequency of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Simulation of Biogenic Fluid Catalytic Cracking (BFCC) Regenerators at Different Scales with MFIX-Exa

Catalytic Fast Pyrolysis (CFP) is a process that converts biomass into liquid intermediates suitable for transportation fuels by rapidly heating it in the presence of a catalyst, aiming to produce stable oils with reduced oxygen content. During CFP, the catalyst can become deactivated by the accumulation of coke, a carbon-rich deposit formed from the decomposition of biomass components. Unlike in petroleum refining, regenerating coked catalysts from biomass pyrolysis requires specific approaches due to the different chemical nature of the coke formed. An experimental technique, Temperature Programmed Oxidation (TPO), was used to study the de-coking process by gradually increasing temperature while monitoring the production of CO and CO2, which provides data for kinetic modeling. Utilizing data from TPO experiments, coke combustion kinetic model was developed to describe the rate of coke removal at different temperatures, allowing for simulation of regeneration processes. Then kinetic model is integrated into MFIX-Exa for the simulation of Biogenic Fluid Catalytic Cracker (BFCC) regenerator at different scales, enabling analysis of catalyst flow, temperature distribution, and regeneration efficiency under various operating conditions.

biogenic fluid catalytic cracking↗

Replicating HCCI-like autoignition behavior: What gasoline surrogate fidelity is needed?

This work seeks to characterize the fidelity needed in a gasoline surrogate with the intent to replicate the complex autoignition behavior exhibited within advanced combustion engines, and specifically Homogeneous Charge Compression Ignition (HCCI). A low-temperature gasoline combustion (LGTC) engine operating in HCCI mode and a rapid compression machine (RCM) are utilized to experimentally quantify fuel reactivity, through autoignition and preliminary heat release characteristics. Fuels considered include a research grade E10 U.S. gasoline (RD5-87), three multi-component surrogates (PACE-1, PACE-8, PACE-20), and a binary surrogate (PRF88.4). Each fuel was studied at lean/HCCI-like conditions covering a wide range of temperatures and pressures that are representative of naturally aspirated to high boost engine operation. Detailed chemical kinetic modeling is also undertaken using a recently updated gasoline surrogate kinetic model to simulate the RCM experiments and to provide chemical insight into surrogate-to-surrogate differences. The LGTC engine experiments demonstrate nearly identical reactivity between PACE-20 and RD5-87 across studied conditions, while faster phasing is seen for both PACE-1 and PACE-8 due to their stronger intermediate- and low-temperature heat release (ITHR/LTHR) at naturally aspirated and boosted conditions, respectively. The RCM experiments reveal typical low-temperature, negative temperature coefficient (NTC) and intermediate-temperature autoignition behaviors at all pressure conditions for RD5-87, which are qualitatively reproduced by all surrogates. Quantitative discrepancies in both autoignition and preliminary heat release are observed for all surrogates, while their ability to replicate RD5-87 autoignition behavior follows the order of PACE-20 > PACE-1 > PACE-8 > PRF88.4. Excellent mapping is obtained between the LGTC engine and the RCM, where the engine pressure-time trajectories can be characterized by the regimes represented by the RCM autoignition isopleths. The kinetic model performs commendably when simulating both autoignition and preliminary heat release of PACE-20, while typically overpredicting ignition delay times for PACE-1, PACE-8 and PRF88.4 at high-pressure and low-temperature/NTC conditions. Sensitivity and rate of production (ROP) analyses highlight surrogate-to-surrogate differences in the governing chemical kinetics where n-pentane initiates rapid OH branching at a faster rate and an earlier timing for PACE-20 than iso-pentane does for PACE-1 and PACE-8, making it computationally more reactive than the other surrogates. The current study highlights the need to include non-standardized properties, such as the lean/HCCI-like autoignition characteristics, in addition to ASTM properties (e.g., RON, MON) as metrics of fuel reactivity and targets to be matched when formulating high-fidelity surrogates that fully capture gasoline advanced combustion behavior such as HCCI-like autoignition.

42 ENGINEERING↗

Formic Acid Decomposition on Au catalysts: DFT, Microkinetic Modeling, and Reaction Kinetics Experiments

A combined theoretical and experimental approach is presented that uses a comprehensive mean-field microkinetic model, reaction kinetics experiments, and scanning transmission electron microscopy imaging to unravel the reaction mechanism and provide insights into the nature of active sites for formic acid (HCOOH) decomposition on Au/SiC catalysts. All input parameters for the microkinetic model are derived from periodic, self-consistent, generalized gradient approximation (GGA-PW91) density functional theory calculations on the Au(111), Au(100), and Au(211) surfaces and are subsequently adjusted to describe the experimental HCOOH decomposition rate and selectivity data. It is shown that the HCOOH decomposition follows the formate (HCOO) mediated path, with 100% selectivity toward the dehydrogenation products (CO21H2) under all reaction conditions. An analysis of the kinetic parameters suggests that an Au surface in which the coordination number of surface Au atoms is 4 may provide a better model for the active site of HCOOH decomposition on these specific supported Au catalysts.

active sites↗

A phase-field model coupled with lattice kinetics solver for modeling crystal growth in furnaces

In this study, we present a new numerical model for crystal growth in a vertical solidification system. This model takes into account the buoyancy induced convective flow and its effect on the crystal growth process. The evolution of the crystal growth interface is simulated using the phase-field method. Two novel phase-field models are developed to model the crystal growth interface in vertical gradient furnaces with two temperature profile setups: 1) fixed wall temperature profile setup and 2) time-dependent temperature profile setup. A semi-implicit lattice kinetics solver based on the Boltzmann equation is employed to model the unsteady incompressible flow. This model is used to investigate the effect of furnace operational conditions on crystal growth interface profiles and growth velocities. For a simple case of macroscopic radial growth, the phase-field model is validated against an analytical solution. Crystal growth in vertical gradient furnaces with two temperature profile setups have been also investigated using the developed model. The numerical simulations reveal that for a certain set of temperature boundary conditions, the heat transport in the melt near the phase interface is diffusion dominant and advection is suppressed.

Lin, Guang↗

Electro-Oxidation of Methane on Platinum under Ambient Conditions

Herein, we investigate the electrochemical conversion of methane to CO 2 on platinum electrodes under ambient conditions. Through a combination of experimentation, density functional theory (DFT) calculations, and ab initio kinetic modeling, we have developed an improved understanding of the reaction mechanism and the factors that determine catalyst activity. We hypothesized that the rate-determining methane activation step is thermochemical (i.e., CH 4 (g) → CH 3 * + H * ) as opposed to electrochemical based on a fitted barrier of approximately 0.96 eV that possesses minimal potential dependence. We developed a simple kinetic model based on the assumption of thermochemical methane activation as the rate-determining step, and the results match well with experimental data. Namely, the magnitude of the maximum current density and the electrode potential at which it is realized agree with our ab initio kinetic model. Finally, we expanded our kinetic model to include other transition metals via a descriptor-based analysis and found platinum to be the most active catalyst for the oxidation of methane, which is in line with previously published experimental observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A redox model for NO oxidation, NH 3 oxidation and high temperature standard SCR over Cu-SSZ-13

A kinetic model is developed to predict the influence of temperature and hydrothermal aging on the redox of active Cu sites under standard SCR, NO oxidation and NH 3 oxidation conditions over a practically relevant fully-formulated Cu-SSZ-13 catalyst. NO 2 /N 2 production during NO/NH 3 titration of Cu II sites is utilized to identify rate parameters associated with NO-only RHC (reduction half cycle) and NH 3 -only RHC respectively. Integral N 2 formation during subsequent NO + NH 3 titration is consistent with the production of one NO 2 per two CuII sites reduced during NO-only RHC and one N 2 per six Cu II sites reduced during NH 3 -only RHC. Decreased reduction of Cu II sites by NO/NH 3 upon hydrothermal aging, along with the production of one NO 2 per two CuII sites during NO-only RHC, is accordant with the involvement of proximal ZCuOH and oxygen-bridged dimeric Cu II sites. Oxidation of partially solvated and framework coordinated ZCu (Cu I ) sites occurs in presence of O 2 , does not produce N 2 and can lead to the consumption of Brønsted acid sites. A global OHC kinetic model is developed to predict transient and integral N 2 formation during exposure of CuI sites to a mixture of NO and O 2 . The resulting redox kinetic model quantitatively predicts NO and NH 3 consumption during isothermal transient response Cu redox (TRCR) protocols, along with temperature and age dependent steady-state standard SCR and oxidation conditions. The redox model presented in this work synthesizes recent kinetic, spectroscopic and computational findings to provide a foundational description of active site redox during standard SCR, NO oxidation and NH3 oxidation over Cu-SSZ-13.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A stochastic asymptotic-preserving scheme for a kinetic-fluid model for disperse two-phase flows with uncertainty

In this paper we consider a kinetic-fluid model for disperse two-phase flows with uncertainty. We propose a stochastic asymptotic-preserving (s-AP) scheme in the generalized polynomial chaos stochastic Galerkin (gPC-sG) framework, which allows the efficient computation of the problem in both kinetic and hydrodynamic regimes. The s-AP property is proved by deriving the equilibrium of the gPC version of the Fokker–Planck operator. The coefficient matrices that arise in a Helmholtz equation and a Poisson equation, essential ingredients of the algorithms, are proved to be positive definite under reasonable and mild assumptions. The computation of the gPC version of a translation operator that arises in the inversion of the Fokker–Planck operator is accelerated by a spectrally accurate splitting method. Numerical examples illustrate the s-AP property and the efficiency of the gPC-sG method in various asymptotic regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupled Kinetics of Ferrihydrite Transformation and As(V) Sequestration under the Effect of Humic Acids: A Mechanistic and Quantitative Study

In natural environments, kinetics of As(V) sequestration/release is usually coupled with the dynamic iron mineral transformation, which is further influenced by the presence of the natural organic matter (NOM). However, both quantitative and mechanistic understanding of this process is still lacking. In this study, we investigated the effect of humic acids (HA), a major component of NOM, on the coupled kinetics of ferrihydrite transformation into hematite/goethite and sequestration of As(V) in iron minerals. Time-resolved As(V) and HA interactions with iron minerals were studied in stirred-flow kinetic experiments, using aberration-corrected scanning transmission electron microscopy (STEM), and X-ray adsorption spectroscopy. Based on the experimental results, we developed a mechanistic kinetics model for As(V) reactions with iron minerals with varying iron mineral compositions and solution chemistry. Our results showed that HA slowed ferrihydrite transformation and the rates of As(V) speciation changes within iron minerals, while HA did not incorporate into iron minerals. Our kinetics model is able to account for the variations of iron mineral compositions, solution chemistry, and As(V) speciation, which has significant environmental implications on predicting the As(V) behavior in the environment.

ferrihydrite, humic acid↗