Methane pyrolysis is a promising route for low-emission hydrogen (H₂) production, with solid carbon as a potentially valuable byproduct. Despite this potential, the economic feasibility of Catalytic Methane Pyrolysis (CMP) with fluidized bed reactors (FBR) has been insufficiently studied. This study develops a conceptual CMP plant using two novel isothermal reactor models—based on continuous stirred-tank reactor (CSTR) and plug-flow reactor (PFR) assumptions—to represent the operational extremes of FBRs. Our reactor models incorporate reaction and catalyst deactivation kinetics from experiments with nickel-supported catalysts, and the framework enables process simulations that account for the catalyst rate required to sustain reactor activity. These models address the lack of proper reduced-order FBR models and the reliance on oversimplified assumptions in the literature. In the baseline scenario, the conceptual plant yields an LCOH ranging from $3.89 to $4.79 per kilogram, defining the expected cost bounds for an FBR-based CMP plant. At a H₂ selling price of $5.00 per kilogram, the process achieves favorable payback time and net present value. Monte Carlo and sensitivity analyses indicate that CMP remains cost-competitive under economic uncertainties. Increased carbon sales could make CMP more economical than steam methane reforming, while unsold byproducts may incur costly sequestration. Reactor heating assessment shows methane combustion with carbon capture minimizes both cost and emissions. Overall, this work demonstrates the economic potential of CMP for H₂ production and provides a practical modeling framework for process evaluation.

Heat pumps, which recycle waste heat, are a promising technology for reducing CO₂ emissions. Efficiently using low-grade waste heat remains challenging due to the limitations of standard heat exchangers and the need for more effective working fluids. This work introduces a multi-scale methodology that combines force field-based Monte Carlo simulations, quantum mechanics, and equations of state to explore the potential of formic acid as a new reactive fluid in thermodynamic cycles. Formic acid exhibits dimerization behavior, forming cyclic dimers in the gas phase, which can enhance the thermodynamic efficiency of heat recovery systems. The dimerization reaction of formic acid is crucial because it integrates chemical energy into thermodynamic processes, potentially improving the performance of heat pumps and other energy systems. The study implements umbrella sampling in Monte Carlo simulations to compute the thermodynamic properties of HCOOH dimerization, including equilibrium constants, enthalpy, and entropy. Results from two different methods to study dimer formation, namely the dimer counter method and the potential of mean force method, show strong agreement with the enthalpy of dimerization of −60.46 kJ mol⁻¹ and −62.91 kJ mol⁻¹, and entropy of −137.36 J mol⁻¹K⁻¹ and −146.98 J mol⁻¹K⁻¹, respectively. A very good agreement of the Monte Carlo results with Quantum Mechanics and experimental data validates the accuracy of the simulations. For phase equilibrium properties, the Peng–Robinson equation of state, coupled with advanced mixing rules, was applied and compared to Monte Carlo simulations in the Gibbs ensemble. This approach enabled the determination of the Global Phase Equilibrium of the system, vaporization enthalpy, phase composition, vapor and liquid densities of the coexisting phases, and entropy as a function of temperature. The agreement between the thermodynamic model and Monte Carlo simulations confirms the reliability of the methodology in capturing the phase behavior of the system. The findings demonstrate a promising approach for discovering and characterizing new reactive fluids, contributing to more efficient and sustainable energy technologies.

Experimental screening of Metal Organic Frameworks (MOFs) for separation applications can be costly and time-consuming. Computational methods can provide many benefits in this process, as expensive compounds and a wide range of operating conditions can be tested while crucial mechanistic insights are gained. TAMOF-1, a recently developed MOF, stands out for its exceptional stability, robustness and cost-effective synthesis. Its good CO2 uptake capacity makes it a promising agent for flue gas separation applications. In this work, we combine experiments with simulations at the atomistic and numerical level to investigate the adsorption and separation of CO2 and N2. Using Monte Carlo simulations, we accurately reproduce experimental adsorption isotherms and elucidate the adsorption mechanisms. TAMOF-1 effectively separates CO2 from N2 because of preferential binding sites near Cu2+ atoms. To assess separation performance in equilibrium at different conditions along the entire isotherm pressure range, adsorbed mole fractions, selectivities, and the trade-off between selectivity and uptake (TSN) are calculated. The dynamic separation performance is assessed by breakthrough experiments and numerical simulations, demonstrating efficient dynamic separation of CO2 and N2, with CO2 being retained in the column.

Computation of the excess entropy (Formula presented.) from the second-order density expansion of the entropy holds strictly for infinite systems in the limit of small densities. For the reliable and efficient computation of (Formula presented.) it is important to understand finite-size effects. Here, expressions to compute (Formula presented.) and Kirkwood–Buff (KB) integrals by integrating the Radial Distribution Function (RDF) in a finite volume are derived, from which (Formula presented.) and KB integrals in the thermodynamic limit are obtained. The scaling of these integrals with system size is studied. We show that the integrals of (Formula presented.) converge faster than KB integrals. We compute (Formula presented.) from Monte Carlo simulations using the Wang–Ramírez–Dobnikar–Frenkel pair interaction potential by thermodynamic integration and by integration of the RDF. We show that (Formula presented.) computed by integrating the RDF is identical to that of (Formula presented.) computed from thermodynamic integration at low densities, provided the RDF is extrapolated to the thermodynamic limit. At higher densities, differences up to (Formula presented.) are observed.

One of the most promising energy carriers for transport applications are hydrogen-based energy carriers. NaBH₄ is a hydrogen energy carrier and produces hydrogen bubbles when it is dissolved in water. The formation of hydrogen bubbles hinders experimental measurements of the thermodynamic and transport properties of aqueous NaBH₄ solutions at elevated temperatures. Accurate knowledge of these properties is essential for the NaBH₄ hydrolysis reactor modeling and design. Molecular dynamics (MD) simulations provide the option to study the thermodynamic and transport properties of NaBH₄ aqueous solutions without hindering hydrogen bubble formation. In this work, a new force field is developed for BH₄^-, namely, the Delft force field of BH₄^- (DFF/BH₄^-), which, combined with additional force fields, can accurately describe experimental densities and viscosities of 0 to 5 m (mol salt/kg water) NaBH₄, 0 to 3 m NaB(OH)₄, and 1 m NaOH aqueous solutions at 295 K within 1.8% and 10.8% maximum deviation, respectively. Empirical fitting correlations are created for densities, viscosities, and self-diffusivities obtained from the MD simulations of 0 to 5 m NaBH₄, 0 to 5 m NaB(OH)₄, and 0 to 1 m NaOH aqueous solutions at 295-363 K for NaBH₄ hydrolysis reactor modeling and design purposes.

Nitric oxide, NO, is a free radical that forms dimers, (NO)2, at its vapor–liquid coexisting temperatures. In this work, we developed an all-atom force field for NO and (NO)2. To assess the performance of this force field, we computed the vapor–liquid equilibrium (VLE) properties of the reactive NO–(NO)2 system, as well as those of pure NO and pure (NO)2, using Continuous Fractional Component Monte Carlo (CFCMC) simulations. We then compared the results with the available experimental data and predictions from two previously developed force fields. For the reactive NO–(NO)2 system, we performed CFCMC simulations in the reactive Gibbs ensemble in which the formation of NO dimers, 2NO ⇌ (NO)2, is considered. The predicted coexistence vapor–liquid densities, dimer mole fractions in the liquid phase, saturated vapor pressures, and heats of vaporization using our force field in the temperature range 120 K to 170 K are in excellent agreement with experimental values. In addition, we conducted a systematic parameter study to analyze the sensitivity of the new force field parameters and the isolated molecule partition functions of (NO)2 on the VLE properties of the reactive NO–(NO)2 system. The results indicate that the VLE properties of the reactive NO–(NO)2 system are affected by both the force field parameters of the involved species as well as the isolated molecule partition functions of (NO)2.

Cyclopentyl methyl ether (CPME) is a promising green solvent due to its eco-friendly properties; it is produced by adding methanol (MeOH) to cyclopentene. Separation of the resulting product mixture containing CPME and MeOH is critical, and it requires vapor-liquid equilibrium (VLE) data. In this work, isobaric VLE data were measured experimentally using an ebulliometer in a 60.0–101.3 kPa pressure range for a binary system of CPME + MeOH. VLE data were modeled using excess Gibbs (G (Formula presented.)) energy-based models such as Wilson, NRTL, and UNIQUAC. The formation of an azeotrope was analyzed. Flash point, surface tension, Gibbs adsorption, and thickness of surface layer were estimated using the Wilson model, which can help in determining molecule interaction and overall behavior of the system. Atmospheric and high-pressure distillation columns were designed using Aspen Plus to study the separation of CPME + MeOH, and an economic evaluation of the same was carried out.

Accurate conductivity predictions of KOH(aq) are crucial for electrolysis applications. OH– is transferred in water by the Grotthuss transfer mechanism, thereby increasing its mobility compared to that of other ions. Classical and ab initio molecular dynamics struggle to capture this enhanced mobility due to limitations in computational costs or in capturing chemical reactions. Most studies to date have provided only qualitative descriptions of the structure during Grotthuss transfer, without quantitative results for the transfer rate and the resulting transport properties. Here, machine learning molecular dynamics is used to investigate 50,000 transfer events. Analysis confirmed earlier works that Grotthuss transfer requires a reduction in accepted and a slight increase in donated hydrogen bonds to the hydroxide, indicating that hydrogen-bond rearrangements are rate-limiting. The computed self-diffusion coefficients and electrical conductivities are consistent with experiments for a wide temperature range, outperforming classical interatomic force fields and earlier AIMD simulations.

Hydrate-based H2 storage is based on the mechanism of trapping H2 in water-based structures that are environmentally friendly and cost-efficient. Understanding the effects of common promoters on hydrate-based H2 storage at the molecular level is crucial for designing efficient storage systems, and for discovering novel promoters. Here, a series of μs-scale molecular dynamics simulations are performed to investigate the nucleation of binary H2 hydrates from gas-liquid two-phase solutions in the presence of various promoters, i.e., CH4, CO2, C2H6, C3H8, C5H10, and THF. The simulation results indicate that the H2 and promoter molecules first dissolve in water from the gas phase and then are absorbed on the cage-faces, promoting the nucleation and growth of binary H2 hydrates. THF is the most effective promoter for hydrate-based H2 storage, exhibiting high performance in converting H2 from the gas phase to hydrates. It is followed by CH4, C2H6, and CO2; C3H8 and C5H10 molecules are less effective H2 hydrate promoters. The presence of large promoter molecules enhances multi-occupied cage formation. The molecular insight into the nucleation of binary H2 hydrates with various promoters provided here not only contributes to a broader understanding of hydrate-based H2 storage but is expected to motivate further experimental and computational studies.

Adsorption simulations often assume a rigid framework, which can be exploited by replacing the expensive framework-adsorbate energy/force evaluation by interpolation of a precomputed energy grid. We present the implementation in RASPA3 of a triquintic interpolation algorithm by Boateng and Bradach and compare it to the tricubic algorithm of Lekien and Marsden. We extended the scheme to interpolation in fractional space to facilitate interpolation of non-rectangular frameworks and evaluated the accuracy. We find that the use of grids is advantageous for larger systems and/or large cutoffs, but generally the efficiency gains are modest (a factor of 2–5).

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.

We explore the impact of force field parameters and reaction equilibrium on the scaling behavior towards the critical point in reactive binary systems, focusing on NO ₂/N ₂O ₄. This system can be considered as a special single-component system since NO ₂ and N ₂O ₄ are in chemical equilibrium via the chemical reaction 2NO ₂⇌N ₂O ₄. We simplify the system by representing both components as single LJ particles, achieving excellent agreement with densities computed using molecular simulations in which all-atom force fields were used. We investigate the effect of force field parameters (ɛ and σ) on phase behavior and show that the critical exponent β remains constant, which means that intermolecular interactions do not affect the scaling to the critical point when the chemical reaction takes place. We also investigate the sensitivity of the reaction equilibrium constant and show that even small changes in isolated molecule partition functions lead to large differences in chemical equilibria. We show that the critical exponent β is different for systems with different reaction equilibrium constants, so a careful parameterization of β is needed for an accurate computation of critical temperatures of reactive mixtures. We perform a screening of reactive binary mixtures for a wide range of ideal gas reaction equilibrium constants, revealing key insights into the thermodynamic behavior and critical properties. Thereby we facilitate the efficient screening of reactive binary mixtures for various applications. Our results emphasize the importance of accurately parameterizing β and provide valuable insights into the critical scaling behavior of complex reactive systems.

Molecular-based equations of state for describing the thermodynamics of chain molecules are often based on mean-field like arguments that reduce the problem of describing the interactions between chains to a simpler one involving only nonbonded monomers. While for dense liquids such arguments are known to work well, at low density they are typically less appropriate due to an incomplete description of the effect of chain connectivity on the local environment of the chains’ monomer segments. To address this issue, we develop three semi-empirical approaches that significantly improve the thermodynamic description of chain molecules at low density. The approaches are developed for chain molecules with repulsive intermolecular forces; therefore, they could be used as reference models for developing equations of the state of real fluids based on perturbation theory. All three approaches are extensions of Wertheim’s first-order thermodynamic perturbation theory (TPT1) for polymerization. The first model, referred to as TPT1-v, incorporates a second-virial correction that is scaled to zero at liquid-like densities. The second model, referred to as TPT1-y, introduces a Helmholtz-energy contribution to account for correlations between next-nearest-neighbor segments within chain molecules. The third approach, called TPT-E, directly modifies TPT1 without utilizing an additional Helmholtz energy contribution. By employing TPT1 at the core of these approaches, we ensure an accurate description of mixtures and enable a seamless extension from chains of tangentially bonded hard-sphere segments of equal size to hetero-segmented chains, fused chains, and chains of soft repulsive segments (which are influenced by temperature). The low-density corrections implemented in TPT1 are designed to preserve these good characteristics, as confirmed through comparisons with novel molecular simulation results for the pressure of various chain fluids. TPT1-v exhibits excellent transferability across different chain types, but it relies on knowing the second virial coefficient of the chain molecules, which is non-trivial to obtain and determined here using Monte Carlo simulation. The TPT1-y model, on the other hand, achieves comparable accuracy to TPT1-v while being fully predictive, requiring no input besides the geometry of the chain molecules.

The confinement effect of porous materials on the thermodynamical equilibrium of the CO₂ hydrogenation reaction presents a cost-effective alternative to transition metal catalysts. In metal-organic frameworks, the type of metal center has a greater impact on the enhancement of formic acid production than the scale of confinement resulting from the pore size. The M-MOF-74 series enables a comprehensive study of how different metal centers affect HCOOH production, minimizing the effect of pore size. In this work, molecular simulations were used to analyze the adsorption of HCOOH and the CO₂ hydrogenation reaction in M-MOF-74, where M = Ni, Cu, Co, Fe, Mn, Zn. We combine classical simulations and density functional theory calculations to gain insights into the mechanisms that govern the low coverage adsorption of HCOOH in the surrounding of the metal centers of M-MOF-74. The impact of metal centers on the HCOOH yield was assessed by Monte Carlo simulations in the grand-canonical ensemble, using gas-phase compositions of CO₂, H₂, and HCOOH at chemical equilibrium at 298.15-800 K, 1-60 bar. The performance of M-MOF-74 in HCOOH production follows the same order as the uptake and the heat of HCOOH adsorption: Ni > Co > Fe > Mn > Zn > Cu. Ni-MOF-74 increases the mole fraction of HCOOH by ca. 10⁵ times compared to the gas phase at 298.15 K, 60 bar. Ni-MOF-74 has the potential to be more economically attractive for CO₂ conversion than transition metal catalysts, achieving HCOOH production at concentrations comparable to the highest formate levels reported for transition metal catalysts and offering a more valuable molecular form of the product.

Experimentally determining thermophysical properties for various compositions commonly found in CO₂ transportation systems is extremely challenging. To overcome this challenge, we performed Monte Carlo (MC) and Molecular Dynamics (MD) simulations of CO₂ rich mixtures to compute thermophysical properties such as densities, thermal expansion coefficients, isothermal compressibilities, heat capacities, Joule-Thomson coefficients, speed of sound, and viscosities at temperatures of (235-313) K and pressures of (20-200) bar. We computed thermophysical properties of pure CO₂ and CO₂ rich mixtures with N₂, Ar, H₂, and CH₄ as impurities of (1-10) mol % and showed good agreement with available Equations of State (EoS). We showed that impurities decrease the values of thermal expansion coefficients, isothermal compressibilities, heat capacities, and Joule-Thomson coefficients in the gas phase, while these values increase in the liquid and supercritical phases. In contrast, impurities increase the value of speed of sound in the gas phase and decrease it in the liquid and supercritical phases. We present an extensive data set of thermophysical properties for CO₂ rich mixtures with various impurities, which will help to design the safe and efficient operation of CO₂ transportation systems.

Linear regression (LR) is used to predict thermochemical properties of alkanes at temperatures (0–1000) K to study chemical reaction equilibria inside zeolites. The thermochemical properties of C1 until C10 isomers reported by Scott are used as training data sets in the LR model which is used to predict these properties for alkanes longer than C10 isomers. Second-order groups are used as independent variables which account for the interactions between the neighboring groups of atoms. This model accurately predicts Gibbs free energies, enthalpies, Gibbs free energies of formation, and enthalpies of formation for alkanes which exceeds the chemical accuracy of 1 kcal/mol and outperforms the group contribution methods developed by Benson et al., Joback and Reid, and Constantinou and Gani. Predictions from our model are used to compute the reaction equilibrium distribution of hydroisomerization of C10 and C14 isomers in MTW-type zeolite. Calculation of reaction equilibrium distribution inside zeolites also requires Henry coefficients of the isomers which can be computed using classical force field-based molecular simulations using the RASPA2 software for which we created an automated workflow. The reaction equilibrium distribution for C10 isomers obtained using the LR model and the training data set for this model are in very good agreement. The tools developed in this study will enable the computational study of hydroisomerization of long-chain alkanes (>C10).

H₂-CO₂ mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H₂-CO₂ mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H₂ due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO₂, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO₂ content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO₂-H₂-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO₂ and H₂ in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO₂ in NaCl brine decreased in the ternary system compared to the binary CO₂-NaCl brine system, the solubility of H₂ in NaCl brine increased less in the ternary system compared to the binary H₂-NaCl brine system. The cooperative effect of H₂-CO₂ enhances the H₂ solubility while suppressing the CO₂ solubility. The water content in the gas phase was found to be intermediate between H₂-NaCl brine and CO₂-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO₂-H₂ mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H₂-CO₂ mixtures.

Focused ultrasound has experimentally been found to enhance the diffusion of nanoparticles; our aim with this work is to study this effect closer using both experiments and non-equilibrium molecular dynamics. Measurements from single particle tracking of 40 nm polystyrene nanoparticles in an agarose hydrogel with and without focused ultrasound are presented and compared with a previous experimental study using 100 nm polystyrene nanoparticles. In both cases, we observed an increase in the mean square displacement during focused ultrasound treatment. We developed a coarse-grained non-equilibrium molecular dynamics model with an implicit solvent to investigate the increase in the mean square displacement and its frequency and amplitude dependencies. This model consists of polymer fibers and two sizes of nanoparticles, and the effect of the focused ultrasound was modeled as an external oscillating force field. A comparison between the simulation and experimental results shows similar mean square displacement trends, suggesting that the particle velocity is a significant contributor to the observed ultrasound-enhanced mean square displacement. The resulting diffusion coefficients from the model are compared to the diffusion equation for a two-time continuous time random walk. The model is found to have the same frequency dependency. At lower particle velocity amplitude values, the model has a quadratic relation with the particle velocity amplitude as described by the two-time continuous time random walk derived diffusion equation, but at higher amplitudes, the model deviates, and its diffusion coefficient reaches the non-hindered diffusion coefficient. This observation suggests that at higher ultrasound intensities in hydrogels, the non-hindered diffusion coefficient can be used.

The series of metal–organic frameworks M-MOF-74 gained popularity in the field of capture and separation of CO₂ due to the presence of numerous, highly reactive open-metal sites. The description of effective interactions between guest molecules and open-metal sites without accounting for polarization effects is challenging but it can significantly reduce the computational cost of simulations. In this study, we propose a non-polarizable force field for CO₂, and H₂ adsorption in M-MOF-74 (M = Ni, Cu, Co, Fe, Mn, Zn) by scaling the Coulombic interactions of M-MOF-74 atoms, and Lennard-Jones interaction potentials between the center of mass of H₂ and the open-metal centers. The presented force field is based on UFF and DREIDING parameters, characterized by high transferability and efficiency. The quantum behavior of H₂ at cryogenic temperatures is considered by incorporating Feynman–Hibbs quantum corrections. To validate the force field, the experimental isotherms of CO₂ at 298 K and 10⁻¹ – 10²kPa, the isotherms of H₂ at 77 K and 10⁻⁵ – 10²kPa, the corresponding enthalpy of adsorption, and the binding geometries in the M-MOF-74 series were reproduced using Monte Carlo simulations in the grand-canonical ensemble. The computed loadings, heats of CO₂ and H₂ adsorption, and binding geometries in M-MOF-74 are in very good agreement with the experimental values. The temperature transferability of the force field from 77 K to 87 K, and 298 K was shown for adsorption of H₂. The validated force field was used to study the adsorption and separation of CO₂/H₂ mixtures at 298 K. The adsorption of H₂ practically does not occur when CO₂ is present in the mixture. As indicated from simulated breakthrough curves, the breakthrough time of CO₂ in M-MOF-74 follows the same order as the uptake and the heat of CO₂ adsorption: Ni ¿ Co ¿ Fe ¿ Mn ¿ Zn ¿ Cu. Increasing the feed mole fraction of CO₂ in the breakthrough simulations from 0.1 to 0.9 speeds up the saturation of the adsorbent, leading to a faster exit of CO₂ with the column effluent. The application of the non-polarizable force field allows full investigation of the capture and separation of CO₂ in M-MOF-74, and can be expanded to study multi-component mixtures or industrial reactions in future research.

Knowledge on the kinetics of gas hydrate dissociation in clay pores at static and dynamic fluid conditions is a fundamental scientific issue for improving gas production efficiency from hydrate deposits using thermal stimulation and depressurization respectively. Here, molecular dynamics simulations were used to investigate poly- and mono-crystalline methane hydrates in Na-montmorillonite clay nanopores. Simulation results show that hydrate dissociation is highly sensitive to temperature and pressure gradients, but their effects differ. Temperature changes increase thermal instability of water and gas molecules, leading to layer-by-layer dissociation from the outer surface. Under flow conditions, laminar flow predominates in nano-pores, and non-Darcy flow occurs due to clay-fluid interactions. Viscous flow disrupts hydrogen bonding at the hydrate surface, enhancing kinetic instability of water. Grain boundaries of polycrystalline hydrates are less stable compared to bulk phases and preferentially decompose, forming new dissociation fronts. This accelerates dissociation compared to monocrystalline hydrates. Fracture occurs at the grain boundaries of polycrystalline hydrate in the fluid, resulting in separate hydrate crystal grains. This fracture process further accelerates hydrate dissociation. In flow systems, methane nanobubbles form in fluid and readily transport with fluid flow. Unlike surface nanobubbles at static conditions, these liquid nanobubbles exhibit mobility. The findings of this study can contribute to a better understanding of the complex phase transition behavior of hydrate in confined environment, and provide theoretical support for improving production control technology.