Multicomponent Fluid Research Articles

SOLEDGE3X integrated core-edge modelling of tungsten sources, migration, and radiation in WEST plasmas

The W impurity species was introduced in the multi-component edge fluid code SOLEDGE3X, in its 2D transport version. To do so, several physical ingredients have been implemented. Firstly, W gross erosion is calculated using the Garcia-Rosales formula, and W re-deposition is estimated with either a Groth–Tskhakaya formulation or a Neural Network trained on re-deposition data from a kinetic model. This fluid approach for W transport is compared with kinetic results obtained with ERO2.0 in edge and SOL domains. Moreover, as W radiates mainly in the core region, which was thus far not resolved in SOLEDGE3X, a 1D model has been implemented for the core plasma region assuming poloidal symmetry. The impact of Oxygen impurities on W content is also addressed (including in particular the impact of W oxide layers on the PFCs). A few applications of the model implementing a Deuterium plasma with Oxygen and Tungsten, taking into account the sputtering of Tungsten by Oxygen are shown.

Mass transport modelling of two partially miscible, multicomponent fluids in nanoporous media

High-pressure fluid transport in nanoporous media such as shale formations requires further understanding because conventional continuum approaches become inadequate due to their ultralow permeability and complexity of transport mechanisms. We propose a species-based approach for modelling two partially miscible, multicomponent fluids in nanoporous media – one that does not rely on conventional bulk fluid transport frameworks but on species movement. We develop a numerical model for species transport of partially miscible, non-ideal fluid mixtures using the chemical potential gradient as the driving force. The model considers the binary friction concept to include the friction between fluid molecules as well as between fluid molecules and pore walls, and incorporates the key multicomponent transport mechanisms – Knudsen, viscous and molecular diffusion. Under single-phase conditions, the system under consideration is quantified by introducing multicomponent Sherwood number (Sh), Péclet number (Pe) and fluid–solid friction modulus (φ). Despite the complexity of fluid transport in nanopores, the steady-state single-phase transport results reveal the contribution of diffusion in nanopores, where all parameters collapse on a set of master curves for the multicomponent Sh with a dependence on multicomponent Pe and φ. Unsteady state, two-phase transport modelling of the codiffusion process shows that light and intermediate alkanes are produced much higher than heavy alkanes when the vapour phase appears. We demonstrate that the pressure gradient is also crucial in promoting CO2 and alkane mixing during counterdiffusion processes. These results stress the need for a paradigm shift from classical bulk flow modelling to species-based transport modelling in nanoporous media.

Influence of Self-Heating on Landfill Leachate Migration

The hydrodynamic processes of landfill leachate migration in the base of a solid waste landfill can have a critical impact on the natural environment. In the case of improper operation of municipal solid waste placement facilities, highly contaminated leachate may penetrate into groundwater and subsequently into surface water. This work addresses fundamental issues of multicomponent fluid propagation in a multilayer porous medium, taking into account temperature inhomogeneities caused by waste decomposition with heat release. The regimes of diffusion and convection of leachate water penetrating into soil layers in the base of municipal solid waste facilities are numerically studied. Archival data from a set of field and laboratory measurements in the area of the operating landfill are used to model the features of pollutant propagation and determine migration parameters. The process of pollutant propagation and migration is described by quantitative values of dry residue content in leachate. Factors that have a significant impact on the migration of leachate are considered. The main ones are convective transfer, diffusion, and properties of the geological composition of the landfill base, which are taken into account in the mathematical formulation of the problem. The calculations show that leachate self-heating can substantially intensify leachate filtration and has to be taken into account in the assessment of leachate migration rate.

Computer modeling of employing binary/ternary organic blends in integrated HP-assisted HDH desalination systems

The global demand for potable water is rising, prompting the development of various energy systems for distilled water production. However, the significance of utilizing multicomponent working fluids in these systems has been largely overlooked. This study presents computer modeling of three HDH-based distillation units powered by two conventional heat pump cycles, namely, the simple and vapor injection heat pumps (first and second models) and an innovative heat pump cycle with an ejector expander (third model). The significance of the research lies in its pioneering investigation of utilizing two- and three-component mixtures in heat pump-based desalination units, which has not been previously explored. The study's primary aim is to determine whether using multicomponent working fluids instead of pure fluids or introducing structural modifications can more effectively enhance the performance of heat pump-based distillation units. The proposed models are simulated using EES and MATLAB software, with the study focusing on energetic, exergetic, exergoeconomic, and heat exchanger modeling to evaluate the feasibility of the configurations. The findings revealed that the structural modifications in the third scenario using R134a resulted in the highest GOR, with improvements of 44 % and 26.03 %, respectively, compared to the other scenarios. However, utilizing binary blend R22/R142b with different compositions improved the GOR of the first to three scenarios by 41.26 %, 29.06 %, and 11.87 %, respectively. Furthermore, in this case, the unit cost of distilled water for the first, second, and third scenarios increased by 12.87 %, 14.32 %, and 12.70 %, respectively. Finally, the first scenario has the highest NPV of 4.40 M$ and the shortest PP of 8.13 years. Therefore, utilizing the blend in a simple heat pump proves to be more efficient and cost-effective than implementing structural modifications.

Interactions of localized wave and dynamics analysis in the new generalized stochastic fractional potential-KdV equation.

In this paper, we investigate the new generalized stochastic fractional potential-Korteweg-de Vries equation, which describes nonlinear optical solitons and photon propagation in circuits and multicomponent plasmas. Inspired by Kolmogorov-Arnold network and our earlier work, we enhance the improved bilinear neural network method by using a large number of activation functions instead of neurons. This method incorporates the concept of simulating more complicated activation functions with fewer parameters, with more diverse activation functions to generate more complex and rare analytical solutions. On this basis, constraints are introduced into the method, reducing a significant amount of computational workload. We also construct neural network architectures, such as "2-3-1," "2-2-3-1," "2-3-3-1," and "2-3-2-1" using this method. Maple software is employed to obtain many exact analytical solutions by selecting appropriate parameters, such as the superposition of double-period lump solutions, lump-rogue wave solutions, and three interaction solutions. The results show that these solutions exhibit more complex waveforms than those obtained by conventional methods, which is of great significance for the electrical systems and multicomponent fluids to which the equation is applied. This novel method shows significant advantages when applied to fractional-order equations and is expected to be increasingly widely used in the study of nonlinear partial differential equations.

Comparison of the Reaction Characteristics of Different Fuels in the Supercritical Multicomponent Thermal Fluid Generation Process

Supercritical multicomponent thermal fluid technology is a new technology with obvious advantages in offshore heavy oil recovery. However, there is currently insufficient understanding of the generation characteristics of the supercritical multicomponent thermal fluid, which is not conducive to the promotion and application of this technology. In order to improve the economic benefits and applicability of the supercritical multicomponent thermal fluid thermal recovery technology, this article reports on indoor supercritical multicomponent thermal fluid generation experiments and compares the reaction characteristics of different fuels in the supercritical multicomponent thermal fluid generation process. The research results indicate that the main components of the products obtained from the supercritical water–crude oil/diesel reaction are similar. Compared to the supercritical water–crude oil reaction, the total enthalpy value of the supercritical multicomponent thermal fluid generated by the supercritical water–diesel reaction is higher, and the specific enthalpy is lower. When the thermal efficiency of the boiler is the same, the energy equilibrium concentration of crude oil is lower than that of diesel. The feasibility of using crude oil instead of diesel to prepare supercritical multicomponent thermal fluids is analyzed from three aspects: reaction mechanism, economic benefits, and technical conditions. It is believed that using crude oil instead of diesel to prepare supercritical multicomponent thermal fluids has good feasibility.

Composition and Injection Rate Co-Optimization Method of Supercritical Multicomponent Thermal Fluid Used for Offshore Heavy Oil Thermal Recovery

Supercritical multicomponent thermal fluid injection is a new technology with great potential for offshore heavy oil thermal recovery. In the process of thermal fluid generation, the reaction conditions including temperature, pressure, and the organic mass concentration in the reaction material will significantly affect its composition and injection rate and will further affect the thermal recovery and development quality of heavy oil. However, there is a lack of relevant research on the variation rules and control methods of the composition and injection rate of supercritical multicomponent thermal fluids, resulting in a lack of technical mechanisms for effective optimization. To fill this gap, a reaction molecular dynamics simulation method was used to simulate thermal fluid generation under different temperatures, pressures, and organic mass concentrations. The changes in thermal fluid composition and yield with reaction conditions were studied, and a control model of thermal fluid composition and yield was established. The proportional relationship between the thermal fluid generation scale of an offshore heavy oil platform and the simulated thermal fluid generation scale is analyzed, and a collaborative optimization method of thermal fluid composition and injection rate in field applications is proposed. The results show the following: (1) The higher the mass concentration of organic matter, the higher the content of supercritical carbon dioxide and supercritical nitrogen in thermal fluids, and the lower the content of supercritical water. (2) The higher the temperature and pressure, the higher the thermal fluid yield, and the higher the organic mass concentration, the lower the thermal fluid yield. (3) The component fitting model conforms to the power function relationship, and the coefficient of determination R2 is greater than 0.9; the yield fitting model conforms to the modified inverse linear logarithmic function relationship, the determination coefficient R2 is greater than 0.8, and the fitting degree is high. (4) The ratio between the actual injection rate of thermal fluids in the mine field and the molecular simulated thermal fluid yield is the ratio of organic matter mass in the platform thermal fluid generator and organic matter mass in the simulated box. (5) Based on the composition and yield control model, combined with the simulation of the ratio relationship between yield and injection rate in the field, a collaborative optimization method of thermal fluid composition and injection rate was established. The research results can provide an effective technical method for predicting, controlling, and optimizing the composition and injection rate of supercritical multicomponent thermal fluids.

Lattice Boltzmann modeling of forced imbibition dynamics in dual-wetted porous media

Forced imbibition dynamics are critical for enhancing recovery rates in reservoirs, as efficient fluid displacement directly impacts resource extraction. This study employs the Zou-He velocity boundary condition and a modified convective boundary condition (mCBC) within the multi-component Shan-Chen lattice Boltzmann method (LBM), thus facilitating unobstructed flow of multi-component fluids at the outlet while maintaining constant outlet pressure. The model's validity was established through outflow tests of immiscible droplets in channels. Its accuracy was further confirmed by comparing results from forced imbibition dynamics tests in dual-wetted pores with theoretical predictions. A symmetrical porous medium was constructed using the stacking method, achieving dual wettability by fixing the contact angle in the upper region and varying it in the lower region. We performed 20 simulation sets under unfavorable viscosity ratios and varying capillary numbers, focusing on overall displacement efficiency before and after breakthrough, and employing energy balance equations to evaluate the dominant forces. Results reveal that capillary forces predominantly dictate forced imbibition dynamics in low capillary number scenarios. In strongly wetted regions, the invading fluid fully occupies pore spaces; however, in weakly wetted regions, displacement efficiency significantly declines as wettability approaches neutrality, even nearing zero. As capillary numbers increase, viscous forces become more prominent, controlling dynamics and leading to fingering and trapping of defending fluids. In weakly wetted areas, the increasing influence of viscous forces enhances fluid displacement, resulting in significant improvements in overall displacement efficiency compared to conditions dominated by capillary forces.

Influence of Reaction Conditions on the Yield of Supercritical Multicomponent Thermal Fluids

Supercritical multicomponent thermal fluid (scMCTF) is a novel medium with great potential for heavy oil thermal recovery. The production rate of scMCTF will affect the injection efficiency of thermal fluid, and then affect the development effect of thermal recovery. However, at present, there are few reports on the production rate of each component of scMCTF, and their understanding is not clear. According to the existing production rate data of supercritical water (scH2O) gasification products, based on the generation mechanism of scMCTF, the production rate of thermal fluid generation products under different generation conditions was calculated, and its influencing factors were identified. The results show the following: (1) The factors affecting the production rate of scMCTF generation products can be divided into three categories: reaction raw material factors, reaction condition factors, and catalytic factors. (2) The hydrocarbon number of raw material, reaction temperature, reaction time, and catalyst concentration were positively correlated with the production rate of the product. (3) The concentration of the reaction raw material is negatively correlated with the production rate of the product. The higher the concentration of the raw material is, the lower the concentration of H2O is, and the steam reforming reaction is inhibited, which leads to the decrease in the production rate. (4) The effect of reaction pressure and catalyst load on the product is not significant. (5) The reaction product production rate increased first and then decreased with the ratio of H2O to oil in the raw material emulsion and the ratio of preheated H2O to raw material discharge. (6) The effect of metal salt catalysts is relatively stable, and the catalytic effect of simple metal catalysts is significantly different under the action of different types of accelerators, so it is necessary to study the degree of synergization of different accelerators on the catalytic effect. The results can lay a foundation for the subsequent experimental and theoretical research design.

Dissolution, phase behavior and mass transport of amorphous solid dispersions in aspirated human intestinal fluids

Amorphous solid dispersions (ASDs) typically show improved dissolution and generate supersaturated solutions, enhancing the oral bioavailability of poorly soluble drugs. To gain insights into intraluminal ASD behavior, we utilized two poorly soluble drugs with different crystallization tendencies, atazanavir and posaconazole, prepared as ASDs at a 10% drug loading with hydroxypropyl methylcellulose acetyl succinate (HPMCAS). We evaluated their release in aspirated fasted-state human intestinal fluid (FaHIF), and multi-component fasted-state simulated intestinal fluid (composite-FaSSIF), characterizing the supersaturation profiles and drug-rich nanodroplets that formed. Complete release was observed for atazanavir ASDs over a 90 min period. Flux for dissolved atazanavir ASDs remained high over the experimental time period of 3 h. In contrast, posaconazole solution concentrations were initially high and then decreased. Likewise, flux was initially high and then decreased where these changes are attributed to crystallization of the drug. Generation of spherical nano-sized amorphous droplets of ∼100-150 nm was found to occur in ex vivo FaHIF media for both ASDs, maximizing the diffusive flux during the supersaturation window. Moreover, buffer capacity differences were postulated to influence release rates of ASDs in simulated vs aspirated fluids. Importantly, the solution phase phenomena observed during ASD release in simulated fluids, namely amorphous nanodroplet formation and drug crystallization, were also found to occur in aspirated luminal fluids.

Research on the Phase Behavior of Multi-Component Thermal-Fluid-Heavy Oil Systems

Multi-component thermal luid technology optimizes development effects and has a strong adaptability, providing a new choice for the efficient development of heavy oil reservoirs. However, due to the significant differences between the phase behavior of multi-component thermal-fluid-heavy oil systems and conventional systems, and the lack of targeted and large-scale research, key issues such as the phase behavior of these systems are unclear. This research studies the phase behavior and influencing factors of emulsions and foamy oil in a multi-component thermal-fluid-heavy oil system through high-temperature and high-pressure PVT experiments, revealing the characteristics of the system’s special phase behavior. In the heavy oil emulsion system, the water content directly affects changes in the system’s phase state. The higher the temperature, the larger the phase transition point, and the two are positively correlated. As the stirring speed increases, the phase transition point first increases and then decreases. The amount of dissolved gas is negatively correlated with the size of the phase transition point, and dissolution can form foamy oil. In the heavy oil–foamy oil system, the dissolution capacity of CO2 is greater than that of multi-component gases, which is greater than that of N2. A high water content and high temperature are not conducive to the dissolution of multi-component gases. While an increase in stirring speed is beneficial for the dissolution of gases, there are limitations to its enhancement ability. Therefore, the development of multi-component thermal fluids should avoid the phase transition point of emulsions and promote the dissolution of multi-component gases.

Study on the adsorption and deformation laws of multi-components in shale oil with nanopores——insights from the molecular simulation

The development of shale oil reservoirs is significant for the petroleum industry. Fluids in shale oil play exist in nanoscale pores, and the interaction between fluids and shale rock surface walls can lead to the expansion and deformation of the entire pore medium thus affecting the development and production of oil reservoirs. We used an organic matter (graphene) pore model to represent the nanopore structure of shale reservoir. Methane, n-hexane and n-dodecane multi-component fluids were used to characterize shale oil. The adsorption and deformation law of multi-component shale oil in the nanopore was systematically studied by using the molecular simulation method. The results showed that the adsorption of multi-component shale oil in organic nanopores was positively correlated with pressure. The total adsorption amount increased with the increase of pressure. Firstly, when the ratio of CO2 was higher, it was more favorable for the development of crude oil and the large amount of CO2 sequestration. The adsorption amount was negatively correlated with temperature. Secondly, the adsorption amount of multicomponent shale oil is proportional to the pore size; The adsorption amount of large-size pore model is more sensitive to pressure changes than that of small-size pore model. Finally, the trend of adsorption deformation volume is similar to that of adsorption volume, which decreases with the increase of temperature. The shale deformation caused by CO2 injection at this time is much smaller than the other component ratios, indicating that the magnitude of deformation is positively related to the adsorption volume.

Comparison of Extraction Methods for the Detection of Avian Influenza Virus RNA in Cattle Milk.

Since early 2024, a multistate outbreak of highly pathogenic avian influenza H5N1 has been affecting dairy cattle in the USA. The influenza viral RNA concentrations in milk make it an ideal matrix for surveillance purposes. However, viral RNA detection in multi-component fluids such as milk can be complex, and optimization of influenza detection methods is thus required. Raw bulk tank milk and mastitis milk samples were artificially contaminated with an avian influenza strain and subjected to five extraction methods. HCoV-229E and synthetic RNA were included as exogenous internal process controls. Given the high viral load usually observed in individual raw milk samples, four out of five tested methods would enable influenza detection in milk with normal texture, over a time window of at least 2 weeks post-onset of clinical signs. Nevertheless, sample dilution 1:3 in molecular transport medium prior to RNA extraction provided the best results for dilution of inhibitory substances and a good recovery rate of influenza RNA, that reached 12.5 ± 1.2% and 10.4 ± 3.8% in two independent experiments in bulk milk and 11.2 ± 3.6% and 10.0 ± 2.9% on two cohorts of mastitis milk samples. We have also shown compatibility of an influenza RT-qPCR system with synthetic RNA detection for simultaneous validation of the RNA extraction and RT-qPCR processes.

High-order finite volume method for solving compressible multicomponent flows with Mie–Grüneisen equation of state

In this paper, we propose a new high-order finite volume method for solving the multicomponent fluids problem with Mie–Grüneisen EOS. Firstly, based on the cell averages of conservative variables, we develop a procedure to reconstruct the cell averages of the primitive variables in a high-order manner. Secondly, the high-order reconstructions employed in computing numerical fluxes are implemented in a characteristic-wise manner to reduce numerical oscillations as much as possible and obtain high-resolution results. Thirdly, advection equation within the governing system is rewritten in a conservative form with a source term to enhance the scheme’s performance. We utilize integration by parts and high-order numerical integration techniques to handle the source terms. Finally, all variables are evolved by using Runge–Kutta time discretization. All steps are carefully designed to maintain the equilibrium of pressure and velocity for the interface-only problem, which is crucial in designing a high-resolution scheme and adapting to more complex multicomponent problems. We have performed extensive numerical tests for both one- and two-dimensional problems to verify our scheme’s high resolution and accuracy.

Subducted carbon weakens the forearc mantle wedge in a warm subduction zone

Subducting oceanic plates carry large amounts of carbon into the Earth’s interior. The subducted carbon is mobilized by fluid and encounters ultramafic rocks in the mantle wedge, resulting in changes to the mineral assemblage and mechanical properties of the mantle. Here, we use thermodynamic modeling of interactions between carbon-bearing multi-component fluids and mantle rocks to investigate the down-dip variation in mineral assemblage in the forearc mantle along subduction megathrusts. We found that fluids rich in aqueous carbon are preferentially generated in a warm subduction zone (e.g., Nankai, SW Japan), causing a change in mineral assemblage from serpentine-rich at the mantle wedge corner to talc + carbonate-rich at greater depths. The transition caused by the infiltration of aqueous carbon may influence the depth of the boundary between the seismogenic and aseismic zones, and the down-dip limit of episodic tremor and slip.

Thermal magnetization transport analysis featuring darcian and multi-physical rheological characteristics in chemically reactive dual convective tangent-hyperbolic nanomaterial confined by vertical cone

The Soret-Dufour characteristics elaborate the phenomena of cross-diffusion where solutal gradients yield heat flux (Dufour aspect) and thermal gradients engender mass flux (Soret effect). Such consideration finds significance in multi-component fluid systems, influencing both heat-mass transportation processes. This study evaluates the porous medium characteristics in convectively heated tangent-hyperbolic nanomaterial driven by a magnetized convected cone. The analysis features Buongiorno nanomaterial model which accounts Brownian motion together with thermophoresis. Thermal transport expression which reports heat transfer characteristics is modeled by considering thermal radiation, convective thermal conditions and heat generation while mass transfer considers the chemical reaction effects. Dimensional expressions are converted into dimensionless forms deploying similarity transformations ensuing in a set of ODEs (ordinary differential expressions) which are computed by deploying bvp4c algorithm. Graphical and tabular behavior of velocity, nanoparticles concentration and temperature fields is inspected corresponding to related variables. This research highlights key findings, revealing a decrease in Nusselt number with increasing heat generation, Dufour parameter and thermophoresis parameter values while opposite trends are verified for radiation parameter values. The velocity function declines for escalating Weissenberg number, magnetic field and permeability parameters but it upsurges with material and buoyancy parameters.

The mineralogical evidence for the formation of gold-bearing pyrite through the metasomatism of siderite by the Au-As-S multicomponent fluids

The mineralogical evidence for the formation of gold-bearing pyrite through the metasomatism of siderite by the Au-As-S multicomponent fluids

Mobility-induced kinetic effects in multicomponent mixtures

We give an overview exploring the role of kinetics in multicomponent mixtures. Compared to the most commonly studied binary (single species plus solvent) case, multicomponent fluids show a rich interplay between kinetics and thermodynamics due to the possibility of fractionation, interdiffusion of mixture components and collective motion. This leads to a competition between multiple timescales that change depending on the underlying kinetics. At high densities, crowding effects are relevant and non-equilibrium structures can become long-lived. We present the main approaches for the study of kinetic effects in multicomponents mixtures, including the role of crowding, and explore their consequences for equilibrium and non-equilibrium scenarios. We conclude by identifying the main challenges in the field.

A compositional numerical study of vapor–liquid-adsorbed three-phase equilibrium calculation in a hydraulically fractured shale oil reservoir

The development of carbon capture, utilization, and storage technologies has notably advanced CO2-enhanced oil recovery (EOR) in shale oil reservoirs, which are characterized by abundant nanopores. These nanopores induce unique phase behaviors in hydrocarbons, challenging traditional phase equilibrium calculation methods. This paper presents a novel three-phase thermodynamic model (vapor–liquid-adsorbed three-phase equilibrium calculation) that addresses these challenges by considering the nanopore capillary pressure, critical parameter transitions, and material exchange between the adsorbed and bulk phases. Grounded in the multicomponent Langmuir–Freundlich adsorption equation and the Peng Robinson equation of state, this model is integrated into the MATLAB Reservoir Simulation Toolbox using an embedded discrete fracture model framework, enabling detailed study of CO2 and hydrocarbon phase behaviors within shale oil nanopores. The results reveal that there are significant nano-constrained effects on multicomponent fluid phase behavior, particularly in pores smaller than 20 nm, leading to notable changes in bubble and dew point pressures, as well as critical condensation pressures and temperatures. CO2 injection further complicates the system, enhancing interactions and expanding the coexistence region of the liquid and gas phases on the pressure–temperature diagram, especially across varying pore sizes. Optimization research on CO2 huff and puff technical parameters for shale oil reservoirs suggests the following optimal settings: a CO2 injection rate of 100 t/day, a shut-in time of 30 days, and six huff and puff cycles. The results of this study offer critical insights into CO2-EOR mechanisms in shale oil reservoirs and emphasize the importance of nanopore properties in EOR.

Giant Effect of CO2 Injection on Multiphase Fluid Adsorption and Shale Gas Production: Evidence from Molecular Dynamics.

Carbon dioxide (CO2) injection in unconventional gas-bearing shale reservoirs is a promising method for enhancing methane recovery efficiency and mitigating greenhouse gas emissions. The majority of methane is adsorbed within the micropores and nanopores (≤50 nm) of shale, which possess extensive surface areas and abundant adsorption sites for the sequestration system. To comprehensively discover the underlying mechanism of enhanced gas recovery (EGR) through CO2 injection, molecular dynamics (MD) provides a promising way for establishing the shale models to address the multiphase, multicomponent fluid flow behaviors in shale nanopores. This study proposes an innovative method for building a more practical shale matrix model that approaches natural underground environments. The grand canonical Monte Carlo (GCMC) method elucidates gas adsorption and sequestration processes in shale gas reservoirs under various subsurface conditions. The findings reveal that previously overlooked pore slits have a significant impact on both gas adsorption and recovery efficiency. Based on the simulation comparisons of absolute and excess uptakes inside the kerogen matrix and the shale slits, it demonstrates that nanopores within the kerogen matrix dominate the gas adsorption while slits dominate the gas storage. Regarding multiphase, multicomponent fluid flow in shale nanopores, moisture negatively influences gas adsorption and carbon storage while promoting methane recovery efficiency by CO2 injection. Additionally, saline solution and ethane further impede gas adsorption while facilitating displacement. Overall, this work elucidates the substantial effect of CO2 injection on fluid transport in shale formations and advances the comprehensive understanding of microscopic gas flow and recovery mechanisms with atomic precision for low-carbon energy development.