Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass diffusivity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
DOE OSTI2022

Microscale THMC Modeling of Pressure Solution in Salt Rock: Impacts of Geometry and Temperature

Abstract Pressure solution, a mechanism that involves tight coupling between the geometry and thermal-hydro-mechanical-chemical (THMC) processes, plays an important role in diagenesis. In this study, we make the first attempt to conduct microscale THMC modeling to understand and quantify the impacts of geometry and temperature on pressure solution, taking natural salt rock as an example. This modeling capability is achieved by expanding a novel MC code that we developed previously (Hu et al. J Geophys Res: Solid Earth 126:e2021JB023112, 2021) to include temperature effects. We first conduct a simulation of an example that involves a single brine inclusion within a single halite grain and find that the temperature impact is limited for that case. We then extract geometry from an image of a natural salt rock and conduct simulations with different cases: (A) only temperature and no stress, (B) only stress and no temperature, and (C) with both stress and temperature. These different cases result in quite different phenomena. In case A, dissolution and precipitation occur across the entire system due to isolated pore space reaching a localized mass balance between dissolution, precipitation, and diffusion. In case B, intense geometric features (e.g., major asperities, inclusions) in one area undergo stress concentration, thus dominating pressure solution in that area. In case C, pressure solution is spread out at contacting highly stressed geometric features close to the hotter side. We conclude that geometric features dominate stress distribution, thus dominating pressure solution in a natural salt rock that may be affected by the temperature if a sufficient temperature gradient is applied.

58 GEOSCIENCES↗
DOE OSTI2025

Enhancing water and oxygen transport through electrode engineering for AEM water electrolyzers

Anion-exchange membrane water electrolyzers (AEMWEs) can accelerate the deployment of more efficient and affordable hydrogen production solutions. Here, electrode structure is shown to affect water back-diffusion and oxygen transport, which, in return, governs overpotential behaviors in AEMWEs. Measurements indicate that electrode with copious catalytic sites produces water close to the AEM, creating a higher water gradient and driving water back-diffusion, which improves membrane hydration and mass transport. In situ measurement reveals a high pH gradient near the anode surface, which affects anode kinetics. Operando measurement shows reduced oxygen accumulation when decoupling oxygen production and transport on anode. Catalyst ink rheology and stability are tuned with additives to realize scalable fabrication of electrodes with enhanced transport features, allowing AEMWE to operate at 2 A cm−2 for over 1,000+ h at a 2.3 μV h−1 degradation rate. Analysis during and post-durability provides insights into degradation mechanisms. This work demonstrates an electrode design strategy for efficient and durable AEMWEs.

AEM↗
DOE OSTI2025

Regional-Scale Modeling Parameterizations for Secondary Organic Aerosol Formation from Isoprene Epoxydiols: Experimentally Based Evaluation and Optimization

Isoprene is an abundant volatile organic compound emitted from broadleaf forests. Under low nitric oxide concentrations, isoprene is photochemically oxidized to form gas-phase isoprene epoxydiols (IEPOX). In the presence of acidified sulfate aerosols, IEPOX enhances the secondary organic aerosol (SOA) formation. Predictions of IEPOX-SOA in regional-scale models, e.g., the Community Multiscale Air Quality Model (CMAQ), are uncertain due to homogeneous aerosol assumptions, underpredictions of water uptake (hygroscopicity), and aerosol surface area. Here, we used experimental measurements of IEPOX-SOA tracers, 2-methyltetrols (2-MT) and 2-methyltetrol sulfates (2-MTS), formed at initial IEPOX-to-inorganic sulfate ratios ranging from 1–10.5, at ∼50% relative humidity to constrain key IEPOX-SOA parameters: phase separation, organic shell diffusivity (D org ), acidity, hygroscopic growth, mass accommodation, and kinetics. The base CMAQ parametrization overpredicted experimental IEPOX-SOA with an average normalized mean bias (NMB average ) of 1.63. CMAQ with phase separation underpredicted IEPOX-SOA (NMB average = −0.71). Using the phase-separated model, CMAQ model performance was optimized (NMB average = 0.077) with an increased D org = 2 × 10 –16 m 2 s –1 and increased rate constants (k 2-MT = 1 × 10 –3 M 2 s –1 , k 2-MTS = 8.83 × 10 –3 M 2 s –1 ). The optimized model explicitly accounted for hygroscopic growth by utilizing experimentally derived growth rates, improving aerosol surface area predictions. Our model highlights the importance of the aerosol mixing state (homogeneous versus phase-separated), aerosol size dynamics, and hygroscopic growth in modeling heterogeneous reactive uptake of IEPOX.

aerosols↗
DOE OSTI2021

Enhancing carbon dioxide gas-diffusion electrolysis by creating a hydrophobic catalyst microenvironment

Electroreduction of carbon dioxide (CO 2 ) over copper-based catalysts provides an attractive approach for sustainable fuel production. While efforts are focused on developing catalytic materials, it is also critical to understand and control the microenvironment around catalytic sites, which can mediate the transport of reaction species and influence reaction pathways. Here, we show that a hydrophobic microenvironment can significantly enhance CO 2 gas-diffusion electrolysis. For proof-of-concept, we use commercial copper nanoparticles and disperse hydrophobic polytetrafluoroethylene (PTFE) nanoparticles inside the catalyst layer. Consequently, the PTFE-added electrode achieves a greatly improved activity and Faradaic efficiency for CO 2 reduction, with a partial current density >250 mA cm -2 and a single-pass conversion of 14% at moderate potentials, which are around twice that of a regular electrode without added PTFE. The improvement is attributed to a balanced gas/liquid microenvironment that reduces the diffusion layer thickness, accelerates CO 2 mass transport, and increases CO 2 local concentration for the electrolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Gas diffusion electrodes, reactor designs and key metrics of low-temperature CO 2 electrolysers

We report that CO 2 emissions can be recycled via low-temperature CO 2 electrolysis to generate products such as carbon monoxide, ethanol, ethylene, acetic acid, formic acid and propanol. In recent years, progress has been made towards an industrially relevant performance by leveraging the development of gas diffusion electrodes (GDEs), which enhance the mass transport of reactant gases (for example, CO 2 ) to the active electrocatalyst. Innovations in GDE design have thus set new benchmarks for CO 2 conversion activity. In this Review, we discuss GDE-based CO 2 electrolysers, in terms of reactor designs, GDE composition and failure modes, to identify the key advances and remaining shortfalls of the technology. This is combined with an overview of the partial current densities, efficiencies and stabilities currently achieved and an outlook on how phenomena such as carbonate formation could influence the future direction of the field. Our aim is to capture insights that can accelerate the development of industrially relevant CO 2 electrolysers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Effects of catalyst droplets on wire growth and the resulting branched structures during VLS growth

Vapor-liquid-solid (VLS) method is vastly employed to grow hierarchical structures with unique properties. However, key questions remain, such as what controls the branched structures and what the roles of catalyst droplet size are during the growth. Here, an in-depth understanding of the kinetics of the nucleation, growth, and subsequent coalescence processes of Bi liquid catalyst droplets is provided by direct observation of PbSe branched wire growth in an environmental transmission electron microscope. This brings a kinetic control of the branch density by varying the parameters, such as temperature. In addition, the dependence of wire growth rate on the catalyst droplet size is revealed, i.e., the smaller the catalyst size the larger the wire growth rate, unlike the wire growth controlled by Gibbs-Thomson effect, possibly due to different mass transport pathways and atomic surface diffusion. These results extend the fundamental understanding of the VLS growth mechanism of branched structures and benefit the structure design of hierarchical materials with tailored properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Defect-mediated diffusion of implanted Mg in GaN: Suppressing dopant redistribution by sequential thermal and microwave annealing

The diffusion behavior of Mg in Mg/N co-implanted GaN is investigated in response to a set of annealing conditions and methodologies, namely, 1000 °C/30 min thermal anneal, by high-temperature pulsed gyrotron microwave annealing at 1420 or 1500 °C, or by thermal and microwave annealing, sequentially. After 1000 °C annealing, the diffusion of Mg in GaN is found to be negligible, as measured by secondary ion mass spectrometry. Annealing by gyrotron microwave annealing alone induces the diffusion of Mg at a rate on the order of 10 −12 cm 2 /s. However, the use of a thermal anneal before microwave gyrotron annealing reduces this rate by an order of magnitude to 10 −13 cm 2 /s. We find that a model that considers Mg diffusion from an inhomogeneous medium that contains a defect-rich implanted region near-surface to a relatively pristine region below the implant range better explains the observed diffusion behavior than a conventional model that assumes a homogeneous medium. By analyzing the diffusion behavior using the Boltzmann–Matano method, we present a discussion of reduction in [V Ga ] by thermal annealing at 1000 °C, leading to a suppressed diffusion coefficient during subsequent high-temperature annealing relative to diffusion after 1420/1500 °C annealing alone. This effect holds potential for improvement in the precision of selectively doped regions for future applications based on the (Al)GaN material system. An improved field profile control in real devices can increase the breakdown and current-handling capabilities in power electronic applications.

Meyers, V.↗
DOE OSTI2022

Magnetic damping in ferromagnetic/heavy-metal systems: The role of interfaces and the relation to proximity-induced magnetism

Damping and spin transport in spintronic multilayered systems continues to be a topic of active research. The enhancement of damping in ferromagnet (FM)/spacer layer (SL)/heavy-metal (HM) thin-film systems was studied for Co 25 Fe 75 /SL/Pt with a nonmagnetic (NM) SL of either Au or Cu with variable thickness, in order to understand the correlation with proximity-induced magnetism (PIM) in the HM. Structural, PIM and magnetic damping measurements were undertaken on the same samples. Specifically, secondary ion mass spectroscopy, element specific x-ray magnetic reflectivity and x-ray magnetic circular dichroism at the Pt and Au L 3 edges, and ferromagnetic resonance methods were used. With increasing thickness of a Cu or Au SL directly between the FM and the Pt layer, the Pt PIM and the damping both fall rapidly, with a relationship between damping and PIM that depends on the SL material. The PIM observed in the Au layer showed a complex dependence on the layer thickness, suggesting some hybridization with the Pt. The role of the number and location of interfaces on the damping was demonstrated with the addition of a SL within the Pt layer, which showed that the specific details of the NM/HM interface also affects the damping. The insertion of a Cu SL within the Pt showed a measurable increase in the overall enhancement of the damping while the insertion of a Au SL into Pt had almost no effect on the damping. Finally, together these results demonstrate the role of both PIM and of additional interfaces in the enhancement of damping in FM/HM systems, which is not fully accounted for by existing theory

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

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↗
DOE OSTI2023

Aggregating in vitro-grown adipocytes to produce macroscale cell-cultured fat tissue with tunable lipid compositions for food applications

We present a method of producing bulk cell-cultured fat tissue for food applications. Mass transport limitations (nutrients, oxygen, waste diffusion) of macroscale 3D tissue culture are circumvented by initially culturing murine or porcine adipocytes in 2D, after which bulk fat tissue is produced by mechanically harvesting and aggregating the lipid-filled adipocytes into 3D constructs using alginate or transglutaminase binders. The 3D fat tissues were visually similar to fat tissue harvested from animals, with matching textures based on uniaxial compression tests. The mechanical properties of cultured fat tissues were based on binder choice and concentration, and changes in the fatty acid compositions of cellular triacylglyceride and phospholipids were observed after lipid supplementation (soybean oil) during in vitro culture. This approach of aggregating individual adipocytes into a bulk 3D tissue provides a scalable and versatile strategy to produce cultured fat tissue for food-related applications, thereby addressing a key obstacle in cultivated meat production.

Yuen Jr, John Se Kit (ORCID:0000000198541654)↗
DOE OSTI2022

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2024

Local Electric Field Effects on Water Dissociation in Bipolar Membranes Studied Using Core–Shell Catalysts

The local electric field strength is thought to affect the rate of water dissociation (WD) in bipolar membranes (BPMs) at the catalyst–nanoparticle surfaces. Here, we study core–shell nanoparticles, where the core is metallic, semiconducting, or insulating, to understand this effect. The nanoparticle cores were coated with a WD catalyst layer (TiO2 or HfO2) via atomic layer deposition (ALD), and the morphology was imaged with transmission electron microscopy. Irrespective of the core material, these core–shell catalysts displayed comparable WD overpotentials at optimal mass loading, despite the hypothesized differences in the electric field strength across the catalyst particle suggested by continuum electrostatic simulations. Substantial atomic interdiffusion between the core and shell was ruled out by X-ray absorption spectroscopy, X-ray photoelectron spectroscopy, and diffuse reflectance optical measurements. However, the optimal mass loading of catalyst was roughly 1 order of magnitude higher for the conductive and high dielectric core materials than for the low dielectric insulating cores. These findings are consistent with the hypothesis that electric field screening within the core material focuses the electric field drop between particles such that larger film thicknesses can be tolerated. Collectively, these data support the idea that it is the local electric field at the molecular level that controls proton-transfer rates and that the metal core/dielectric-shell constructs introduced here modulate that field. Further materials and synthetic design may enable optimization of the electric field strength across the proton-transfer trajectory at the material surface.

Sarma, Prasad V↗
DOE OSTI2023

Diffusion of Hydrophilic to Hydrophobic Forms of Nile Red in Aqueous C 12 EO 10 Gels by Variable Area Fluorescence Correlation Spectroscopy

In this work, solute diffusion within lyotropic liquid crystal gels prepared from a series of water and decaethylene glycol monododecyl ether (C 12 EO 10 ) mixtures was explored by variable area fluorescence correlation spectroscopy. Aqueous C 12 EO 10 gels were prepared in concentrations ranging from 55:45 to 70:30 wt% of surfactant and water. Small angle X-ray scattering revealed that these gels comprise hexagonal mesophases of cylindrical micelles. Micelle spacing was found to decrease with increasing C 12 EO 10 concentration. Three different Nile red (NR) dyes were employed as model solutes and were separately doped into the gels at nanomolar levels. These include a hydrophilic form of NR incorporating an anionic sulfonate group (NRSO 3 - ), a hydrophobic form incorporating a fourteen-carbon alkane tail (NRC 14 ), and commercial NR as an intermediate case. FCS data acquired from the gels revealed that NRSO 3 - diffused primarily in 3D. Its diffusion coefficient exhibited a monotonic decrease with increasing gel concentration and micelle packing density, consistent with confinement of its motions as a result of its exclusion from the micelle cores. NRC 14 exhibited the smallest diffusion coefficient, most likely due to its larger size and enhanced interactions with the micelle cores. NR yielded an intermediate diffusion coefficient and the most anomalous behavior of the three dyes, attributable to its facile partitioning between core and corona regions, and greater participation by 1D diffusion. The results of these studies afford an improved understanding of molecular mass transport through soft-matter nanomaterials like those being developed for use in drug delivery and membrane based chemical separations.

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

Product selectivity and mass transport in levulinic acid transfer hydrogenation by monolithic MIL-100, MIL-88B and ZIF-8@Pd MOFs

The diffusion processes between adsorbent and adsorbate naturally play a significant role in the efficiency and selectivity of the heterogenous catalytic process. This paper considers the importance of diffusion processes in the transfer hydrogenation reaction of levulinic acid to γ-valerolactone by MIL-88B, MIL-100, and ZIF8@Pd monolithic catalysts. Over a period of five catalytic cycles, it was shown that the Fe-based catalysts can achieve similar conversions to the ZIF-8 supported Pd, with the only current limitation being the lower aqueous stability of these MOFs. Diffusion constants were calculated using the ZLC method, with micropore diffusion limitation found for ZIF-8 and MIL-100 monolithic frameworks at 2.7 x 10 -8 and 4.6 x 10 -8 cm 2 s -1 respectively. This diffusion limitation was also confirmed by IR spectroscopy with an increasing concentration of C—H bands on the MOF substrate post-reaction. Mass transfer coefficients, also calculated by ZLC, revealed increased mass transport for the hydrophobic ZIF-8 framework, which perhaps aids in the γ-valerolactone selectivity over side products that are produced in the absence of catalytic material, as seen for MIL-88B and MIL-100 after multiple uses.

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

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗