A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes...A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.展开更多
Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow...Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.展开更多
Based on turbulent theory, a 3D coupled model of fluid flow and solidification was built using finite difference method and used to study the influence of superheating degree and casting speed on fluid flow and solidi...Based on turbulent theory, a 3D coupled model of fluid flow and solidification was built using finite difference method and used to study the influence of superheating degree and casting speed on fluid flow and solidification, analyze the interaction between shell and molten steel, and compare the temperature distribution under different technological conditions. The results indicate that high superheating degree can lengthen the liquid-core depth and make the crack and breakout possible, so suitable superheating should be controlled within 35℃ according to the simulation results. Casting speed which is one of the most important technological parameters of improving production rate, should be controlled between 0. 85 m/min and 1.05 m/min and the caster has great potential in the improvement of blank quality.展开更多
Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the...Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.展开更多
A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The res...A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.展开更多
In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma propertie...In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.展开更多
Study on gas–liquid flow in stirred tank with two combinations of dual-impeller(six-bent-bladed turbine(6BT)+six-inclined-blade down-pumping turbine(6 ITD),the six-bent-bladed turbine(6BT)+six-inclinedblade up-pumpin...Study on gas–liquid flow in stirred tank with two combinations of dual-impeller(six-bent-bladed turbine(6BT)+six-inclined-blade down-pumping turbine(6 ITD),the six-bent-bladed turbine(6BT)+six-inclinedblade up-pumping turbine(6ITU))was conducted using computational fluid dynamics(CFD)and population balance model(PBM)(CFD-PBM)coupled model.The local bubble size was captured by particle image velocimetry(PIV)measurement.The gas holdup,bubble size distribution and gas–liquid interfacial area were explored at different conditions through numerical simulation.The results showed that the 4 mm bubbles accounted for the largest proportion of 33%at the gas flow rates Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1) for combined impeller of 6BT+6ITU,while the bubbles of 4.7 mm and 5.5 mm were the largest proportion for 6BT+6ITD combination,i.e.25%at Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1),respectively,which indicated that 6BT+6ITU could reduce bubble size effectively and promote gas dispersion.In addition,the gas holdup around impellers was increased obviously with the speed compared with gas flow rate.So it was concluded that 6ITU impeller could be more conductive to the bubble dispersion with more uniform bubble size,which embodied the advantages of 6BT+6ITU combination in gas–liquid mixing.展开更多
From the viewpoint of interaction mechanics for solid and gas, a coupled mathematical model was presented for solid coal/rock deformation and gas leak flow in parallel deformable coal seams. Numerical solutions using ...From the viewpoint of interaction mechanics for solid and gas, a coupled mathematical model was presented for solid coal/rock deformation and gas leak flow in parallel deformable coal seams. Numerical solutions using the SIP (Strong Implicit Proce- dure) method to the coupled mathematical model for double parallel coal seams were also developed in detail. Numerical simulations for the prediction of the safety range using protection layer mining were performed with experimental data from a mine with potential danger of coal/gas outbursts. Analyses show that the numerical simulation results are consistent with the measured data in situ.展开更多
This paper gives a numerical study on the flow and temperature fields in an induced plasma reactor, which worked in 0.5 ATM with air as a working gas. We employed a two-dimensional mode of an inductively coupled plas...This paper gives a numerical study on the flow and temperature fields in an induced plasma reactor, which worked in 0.5 ATM with air as a working gas. We employed a two-dimensional mode of an inductively coupled plasma to calculate the temperature and flow field of the reactor as well as the generator. The algorithm is based on the solutions of the two-dimensional continuity, momentum, and energy equations in term of vorticity, stream function and enthalpy. An upwind finite-difference scheme was adopted to solve those equations with appropriate boundary conditions. The computed results show that there is a flat region with little parameter change in the reactor, that the diameter of the region is not much larger than that of the generator and that a deep change of parameter exists in the outer side of the region.展开更多
In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)...In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)model that can simultaneously consider the fracture evolution and non-Darcy flow dynamics in studying the thermo-hydro-mechanical(THM)coupling processes for heat extraction in geothermal reservoir.We further employed the model on the Habanero enhanced geothermal systems(EGS)project located in Australia.First,our findings illustrate a clear spatial-temporal variation in the thermal stress and pressure perturbations,as well as uneven spatial distribution of shear failure in 3D fracture networks.Activated shear failure is mainly concentrated in the first fracture cluster.Secondly,channeling flow have also been observed in DFNs during heat extraction and are further intensified by the expansion of fractures driven by thermal stresses.Moreover,the combined effect of non-Darcy flow and fracture evolution triggers a rapid decline in the resulting heat rate and temperature.The NR-DFN model framework and the Habanero EGS's results illustrate the importance of both fracture evolution and non-Darcy flow on the efficiency of EGS production and have the potential to promote the development of more sustainable and efficient EGS operations for stakeholders.展开更多
Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recogni...Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.展开更多
A comprehensive mathematical model has been developed to describe the interaction of the multiple physics fields during the conventional DC casting and LFEC (low frequency electromagnetic casting) process. The model i...A comprehensive mathematical model has been developed to describe the interaction of the multiple physics fields during the conventional DC casting and LFEC (low frequency electromagnetic casting) process. The model is based on a combination of the commercial finite element package ANSYS and the commercial finite volume package FLUENT, with the former for the calculation of the electromagnetic field and the latter for the calculation of the magnetic driven fluid flow, heat transfer and solidification. Moreover, the model has been verified against the temperature measurements obtained from two 7XXX aluminum alloy billets of 200mm diameter, cast during the conventional DC casting and the LFEC casting processes. In addition, a measurement of the sump shape of the billets were carried out by using addition melting metal of Al-30%Cu alloy into the billets during casting process. There was a good agreement between the calculated results and the measured results. Further, comparison of the calculated results during the LFEC process with that during the conventional DC casting process indicated that velocity patterns, temperature profiles and the sump depth are strongly modified by the application of a low frequency electromagnetic field during the DC casting.展开更多
We present a workflow linking coupled fluid-flow and geomechanical simulation with seismic modelling to predict seismic anisotropy induced by non-hydrostatic stress changes. We generate seismic models from coupled sim...We present a workflow linking coupled fluid-flow and geomechanical simulation with seismic modelling to predict seismic anisotropy induced by non-hydrostatic stress changes. We generate seismic models from coupled simulations to examine the relationship between reservoir geometry, stress path and seismic anisotropy. The results indicate that geometry influences the evolution of stress,which leads to stress-induced seismic anisotropy. Although stress anisotropy is high for the small reservoir, the effect of stress arching and the ability of the side-burden to support the excess load limit the overall change in effective stress and hence seismic anisotropy. For the extensive reservoir, stress anisotropy and induced seismic anisotropy are high. The extensive and elongate reservoirs experience significant compaction, where the inefficiency of the developed stress arching in the side-burden cannot support the excess load.The elongate reservoir displays significant stress asymmetry,with seismic anisotropy developing predominantly along the long-edge of the reservoir. We show that the link betweenstress path parameters and seismic anisotropy is complex,where the anisotropic symmetry is controlled not only by model geometry but also the nonlinear rock physics model used. Nevertheless, a workflow has been developed to model seismic anisotropy induced by non-hydrostatic stress changes, allowing field observations of anisotropy to be linked with geomechanical models.展开更多
In light of previous work [Phys. Rev. E 60 4000 （1999）], a modified coupled-map car-following model is proposed by considering the headways of two successive vehicles in front of a considered vehicle described by th...In light of previous work [Phys. Rev. E 60 4000 （1999）], a modified coupled-map car-following model is proposed by considering the headways of two successive vehicles in front of a considered vehicle described by the optimal velocity function. The non-jam conditions are given on the basis of control theory. Through simulation, we find that our model can exhibit a better effect as p = 0.65, which is a parameter in the optimal velocity function. The control scheme, which was proposed by Zhao and Gao, is introduced into the modified model and the feedback gain range is determined. In addition, a modified control method is applied to a mixed traffic system that consists of two types of vehicle. The range of gains is also obtained by theoretical analysis. Comparisons between our method and that of Zhao and Gao are carried out, and the corresponding numerical simulation results demonstrate that the temporal behavior of traffic flow obtained using our method is better than that proposed by Zhao and Gao in mixed traffic systems.展开更多
Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformab...Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformable coal seam can be essentially considered to be compressible with time dependent and mixed permeation and diffusion through a pore cleat deformable heterogeneous and anisotropy medium. Based on this new viewpoint, a coupled mathematical model for coal seam deformation and gas leakage in a double coal seam was formulated and numerical simulations for gas emission from the coal seam are presented. It is found that coupled models might be closer to reality.展开更多
As the demand for more efficient and adaptable power distribution systems intensifies, especially in rural areas, innovative solutions like the Capacitor-Coupled Substation with a Controllable Network Transformer (CCS...As the demand for more efficient and adaptable power distribution systems intensifies, especially in rural areas, innovative solutions like the Capacitor-Coupled Substation with a Controllable Network Transformer (CCS-CNT) are becoming increasingly critical. Traditional power distribution networks, often limited by unidirectional flow capabilities and inflexibility, struggle to meet the complex demands of modern energy systems. The CCS-CNT system offers a transformative approach by enabling bidirectional power flow between high-voltage transmission lines and local distribution networks, a feature that is essential for integrating renewable energy sources and ensuring reliable electrification in underserved regions. This paper presents a detailed mathematical representation of power flow within the CCS-CNT system, emphasizing the control of both active and reactive power through the adjustment of voltage levels and phase angles. A control algorithm is developed to dynamically manage power flow, ensuring optimal performance by minimizing losses and maintaining voltage stability across the network. The proposed CCS-CNT system demonstrates significant potential in enhancing the efficiency and reliability of power distribution, making it particularly suited for rural electrification and other applications where traditional methods fall short. The findings underscore the system's capability to adapt to varying operational conditions, offering a robust solution for modern power distribution challenges.展开更多
A mathematical model for heap bioleaching is developed to analyze heat transfer, oxygen flow, target ion distribution and oxidation leaching rate in the heap. The model equations are solved with Comsol Multiphysics so...A mathematical model for heap bioleaching is developed to analyze heat transfer, oxygen flow, target ion distribution and oxidation leaching rate in the heap. The model equations are solved with Comsol Multiphysics software. Numerical simulation results show the following facts： Concentration of oxygen is relatively high along the boundary of the slope, and low in the center part where leaching rate is slow. Temper- ature is relatively low along the slope and reaches the highest along the bottom region near the slope, with difference being more than 6℃. Concentration of target mental ions is the highest in the bottom region near the slope. Oxidation leaching rate is relatively large in the bottom and slope part with a fast reaction rate, and small in the other part with low oxygen concentration.展开更多
In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a h...In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.展开更多
This paper presents the realization of two-way coupling of the unsaturated-saturated flow interactions of the SWAT2000 and MODFLOW96 models on the basis of the integrated surface/groundwater model SWATMOD99, and its a...This paper presents the realization of two-way coupling of the unsaturated-saturated flow interactions of the SWAT2000 and MODFLOW96 models on the basis of the integrated surface/groundwater model SWATMOD99, and its application in Hetao Irrigation District （HID）, Inner Mongolia, China. Major revisions and enhancements were made to the SWAT2000 and MODFLOW models for simulating the detailed hydrologic budget and coupled unsaturated and saturated interactions, and irrigation canal hydrology for the HID. The simulation results of seasonal groundwater recharge to and evaporate from the shallow groundwater, and the annual water budget over the district are presented and discussed. The results implied the necessity of two-way coupling of the unsaturated-saturated interactions when groundwater is shallow, and the feasibility of making comprehensive use of the information coming from both the surface water and groundwater models to make a more physically-based assessment of the coupled interactions.展开更多
Cemented tailings backfill(CTB) is made by mixing cement, tailings and water together, thus cement hydration and water seepage flow are the two crucial factors affecting the quality of CTB. Cement hydration process ca...Cemented tailings backfill(CTB) is made by mixing cement, tailings and water together, thus cement hydration and water seepage flow are the two crucial factors affecting the quality of CTB. Cement hydration process can release significant amount of heat to raise the temperature of CTB and in turn increase the rate of cement hydration. Meanwhile, the progress of cement hydration consumes water and produces hydration products to change the pore structures within CTB, which further influences the hydraulic behavior of CTB. In order to understand the hydraulic behavior of CTB, a numerical model was developed by coupling the hydraulic,thermal and hydration equations. This model was then implemented into COMSOL Multiphysics to simulate the evolutions of temperature and water seepage flow within CTB versus curing time. The predicted outcomes were compared with correspondent experimental results, proving the validity and availability of this model. By taking advantage of the validated model, effects of various initial CTB and curing temperatures, cement content, and CTB's geometric shapes on the hydraulic behavior of CTB were demonstrated numerically. The presented conclusions can contribute to preparing more environmentally friendly CTB structures.展开更多
基金This study was supported by the National Natural Science Foundation of China(U22B2075,52274056,51974356).
文摘A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.
基金This research was conducted within the DECOVALEX-2023 projectDECOVALEX is an international research project comprising participants from industry,government,and academia,focusing on development of understanding,models and codes in complex coupled problems in sub-surface geological and engineering applications.DECOVALEX-2023 is the current phase of the project.The authors appreciate the DECOVALEX-2023 Funding Organisations Andra,BASE,BGE,BGR,CAS,CNSC,COVRA,US DOE,ENRESA,ENSI,JAEA,KAERI,NWMO,RWM,SÚRAO,SSM and Taipower for their financial and technical support of the work described in this paper.The statements made in the paper are,however,solely those of the authors and do not necessarily reflect those of the Funding Organisations.This work was further supported by the German Federal Ministry for Economic Affairs and Climate Action(BMWK).
文摘Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.
文摘Based on turbulent theory, a 3D coupled model of fluid flow and solidification was built using finite difference method and used to study the influence of superheating degree and casting speed on fluid flow and solidification, analyze the interaction between shell and molten steel, and compare the temperature distribution under different technological conditions. The results indicate that high superheating degree can lengthen the liquid-core depth and make the crack and breakout possible, so suitable superheating should be controlled within 35℃ according to the simulation results. Casting speed which is one of the most important technological parameters of improving production rate, should be controlled between 0. 85 m/min and 1.05 m/min and the caster has great potential in the improvement of blank quality.
基金the National High Technical Reasearch and Development Programme of China （No. 2003AA327140） the National Natural Science Foundation of China （No. 50374081）.
文摘Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.
文摘A mathematical model for coupled multiphase fluid flow and sedi- mentation deformation is developed based on fluid-solid interaction mechanism.A finite difference-finite element numerical approach is presented.The results of an example show that the fluid-solid coupled effect has great influence on multiphase fluid flow and reservoir recovery performances,and the coupled model has practical significance for oilfield development.
基金Project supported by the National Natural Science Foundations of China(Grant No.11305023)
文摘In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.
基金supported by the National Natural Science Foundation of China(52176040)Shandong Provincial Natural Science Foundation of China(ZR2018LE015)。
文摘Study on gas–liquid flow in stirred tank with two combinations of dual-impeller(six-bent-bladed turbine(6BT)+six-inclined-blade down-pumping turbine(6 ITD),the six-bent-bladed turbine(6BT)+six-inclinedblade up-pumping turbine(6ITU))was conducted using computational fluid dynamics(CFD)and population balance model(PBM)(CFD-PBM)coupled model.The local bubble size was captured by particle image velocimetry(PIV)measurement.The gas holdup,bubble size distribution and gas–liquid interfacial area were explored at different conditions through numerical simulation.The results showed that the 4 mm bubbles accounted for the largest proportion of 33%at the gas flow rates Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1) for combined impeller of 6BT+6ITU,while the bubbles of 4.7 mm and 5.5 mm were the largest proportion for 6BT+6ITD combination,i.e.25%at Q=0.76 m^(3)·h^(-1) and 22%at Q=1.52 m^(3)·h^(-1),respectively,which indicated that 6BT+6ITU could reduce bubble size effectively and promote gas dispersion.In addition,the gas holdup around impellers was increased obviously with the speed compared with gas flow rate.So it was concluded that 6ITU impeller could be more conductive to the bubble dispersion with more uniform bubble size,which embodied the advantages of 6BT+6ITU combination in gas–liquid mixing.
文摘From the viewpoint of interaction mechanics for solid and gas, a coupled mathematical model was presented for solid coal/rock deformation and gas leak flow in parallel deformable coal seams. Numerical solutions using the SIP (Strong Implicit Proce- dure) method to the coupled mathematical model for double parallel coal seams were also developed in detail. Numerical simulations for the prediction of the safety range using protection layer mining were performed with experimental data from a mine with potential danger of coal/gas outbursts. Analyses show that the numerical simulation results are consistent with the measured data in situ.
文摘This paper gives a numerical study on the flow and temperature fields in an induced plasma reactor, which worked in 0.5 ATM with air as a working gas. We employed a two-dimensional mode of an inductively coupled plasma to calculate the temperature and flow field of the reactor as well as the generator. The algorithm is based on the solutions of the two-dimensional continuity, momentum, and energy equations in term of vorticity, stream function and enthalpy. An upwind finite-difference scheme was adopted to solve those equations with appropriate boundary conditions. The computed results show that there is a flat region with little parameter change in the reactor, that the diameter of the region is not much larger than that of the generator and that a deep change of parameter exists in the outer side of the region.
基金funded by the National Natural Science Foundation of China (No.U22A20166)Science and Technology Foundation of Guizhou Province (No.QKHJC-ZK[2023]YB074)+2 种基金Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical EngineeringInstitute of Rock and Soil MechanicsChinese Academy of Sciences (No.SKLGME022009)。
文摘In fractured geothermal reservoirs,the fracture networks and internal fluid flow behaviors can significantly impact the thermal performance.In this study,we proposed a non-Darcy rough discrete fracture network(NR-DFN)model that can simultaneously consider the fracture evolution and non-Darcy flow dynamics in studying the thermo-hydro-mechanical(THM)coupling processes for heat extraction in geothermal reservoir.We further employed the model on the Habanero enhanced geothermal systems(EGS)project located in Australia.First,our findings illustrate a clear spatial-temporal variation in the thermal stress and pressure perturbations,as well as uneven spatial distribution of shear failure in 3D fracture networks.Activated shear failure is mainly concentrated in the first fracture cluster.Secondly,channeling flow have also been observed in DFNs during heat extraction and are further intensified by the expansion of fractures driven by thermal stresses.Moreover,the combined effect of non-Darcy flow and fracture evolution triggers a rapid decline in the resulting heat rate and temperature.The NR-DFN model framework and the Habanero EGS's results illustrate the importance of both fracture evolution and non-Darcy flow on the efficiency of EGS production and have the potential to promote the development of more sustainable and efficient EGS operations for stakeholders.
基金supported by the National Natural Science Foundation of China(No.52174045)。
文摘Natural fractures(NFs)are common in shale and tight reservoirs,where staged multi-cluster fracturing of horizontal wells is a prevalent technique for reservoir stimulation.While NFs and stress interference are recognized as significant factors affecting hydraulic fracture(HF)propagation,the combined influence of these factors remains poorly understood.To address this knowledge gap,a novel coupled hydromechanical-damage(HMD)model based on the phase field method is developed to investigate the propagation of multi-cluster HFs in fractured reservoirs.The comprehensive energy functional and control functions are established,while incorporating dynamic fluid distribution between multiple perforation clusters and refined changes in rock mechanical parameters during hydraulic fracturing.The HMD coupled multi-cluster HF propagation model investigates various scenarios,including single HF and single NF,reservoir heterogeneity,single HF and NF clusters,and multi-cluster HFs with NF clusters.The results show that the HMD coupling model can accurately capture the impact of approach angle(θ),stress difference and cementation strength on the interaction of HF and NF.The criterion of the open and cross zones is not fixed.The NF angle(a)is not a decisive parameter to discriminate the interaction.According to the relationship between approach angle(θ)and NF angle(a),the contact relationship of HF can be divided into three categories(θ=a,θ<a,andθ>a).The connected NF can increase the complexity of HF by inducing it to form branch fracture,resulting in a fractal dimension of HF as high as2.1280 at angles of±45°.Inter-fracture interference from the heel to the toe of HF shows the phenomenon of no,strong and weak interference.Interestingly,under the influence of NFs,distant HFs from the injection can become dominant fractures.However,as a gradually increases,inter-fracture stress interference becomes the primary factor influencing HF propagation,gradually superseding the dominance of NF induced fractures.
文摘A comprehensive mathematical model has been developed to describe the interaction of the multiple physics fields during the conventional DC casting and LFEC (low frequency electromagnetic casting) process. The model is based on a combination of the commercial finite element package ANSYS and the commercial finite volume package FLUENT, with the former for the calculation of the electromagnetic field and the latter for the calculation of the magnetic driven fluid flow, heat transfer and solidification. Moreover, the model has been verified against the temperature measurements obtained from two 7XXX aluminum alloy billets of 200mm diameter, cast during the conventional DC casting and the LFEC casting processes. In addition, a measurement of the sump shape of the billets were carried out by using addition melting metal of Al-30%Cu alloy into the billets during casting process. There was a good agreement between the calculated results and the measured results. Further, comparison of the calculated results during the LFEC process with that during the conventional DC casting process indicated that velocity patterns, temperature profiles and the sump depth are strongly modified by the application of a low frequency electromagnetic field during the DC casting.
基金the sponsors of the IPEGG project, BG, BP, Statoilthe Research Council UK (EP/K035878/1+1 种基金 EP/K021869/1 NE/L000423/1) for financial support
文摘We present a workflow linking coupled fluid-flow and geomechanical simulation with seismic modelling to predict seismic anisotropy induced by non-hydrostatic stress changes. We generate seismic models from coupled simulations to examine the relationship between reservoir geometry, stress path and seismic anisotropy. The results indicate that geometry influences the evolution of stress,which leads to stress-induced seismic anisotropy. Although stress anisotropy is high for the small reservoir, the effect of stress arching and the ability of the side-burden to support the excess load limit the overall change in effective stress and hence seismic anisotropy. For the extensive reservoir, stress anisotropy and induced seismic anisotropy are high. The extensive and elongate reservoirs experience significant compaction, where the inefficiency of the developed stress arching in the side-burden cannot support the excess load.The elongate reservoir displays significant stress asymmetry,with seismic anisotropy developing predominantly along the long-edge of the reservoir. We show that the link betweenstress path parameters and seismic anisotropy is complex,where the anisotropic symmetry is controlled not only by model geometry but also the nonlinear rock physics model used. Nevertheless, a workflow has been developed to model seismic anisotropy induced by non-hydrostatic stress changes, allowing field observations of anisotropy to be linked with geomechanical models.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11372166,11372147,61074142,and 11072117)the Scientific Research Fund of Zhejiang Province,China(Grant No.LY13A010005)+1 种基金the Disciplinary Project of Ningbo City,China(Grant No.SZXL1067)the K.C.Wong Magna Fund in Ningbo University,China,and the Government of the Hong Kong Administrative Region,China(Grant No.119011)
文摘In light of previous work [Phys. Rev. E 60 4000 （1999）], a modified coupled-map car-following model is proposed by considering the headways of two successive vehicles in front of a considered vehicle described by the optimal velocity function. The non-jam conditions are given on the basis of control theory. Through simulation, we find that our model can exhibit a better effect as p = 0.65, which is a parameter in the optimal velocity function. The control scheme, which was proposed by Zhao and Gao, is introduced into the modified model and the feedback gain range is determined. In addition, a modified control method is applied to a mixed traffic system that consists of two types of vehicle. The range of gains is also obtained by theoretical analysis. Comparisons between our method and that of Zhao and Gao are carried out, and the corresponding numerical simulation results demonstrate that the temporal behavior of traffic flow obtained using our method is better than that proposed by Zhao and Gao in mixed traffic systems.
文摘Based on the new viewpoint of solid and gas interaction mechanics, gas leakage in a double deformable coal seam can be understood. That is, under the action of geophysical fields, the methane flow in a double deformable coal seam can be essentially considered to be compressible with time dependent and mixed permeation and diffusion through a pore cleat deformable heterogeneous and anisotropy medium. Based on this new viewpoint, a coupled mathematical model for coal seam deformation and gas leakage in a double coal seam was formulated and numerical simulations for gas emission from the coal seam are presented. It is found that coupled models might be closer to reality.
文摘As the demand for more efficient and adaptable power distribution systems intensifies, especially in rural areas, innovative solutions like the Capacitor-Coupled Substation with a Controllable Network Transformer (CCS-CNT) are becoming increasingly critical. Traditional power distribution networks, often limited by unidirectional flow capabilities and inflexibility, struggle to meet the complex demands of modern energy systems. The CCS-CNT system offers a transformative approach by enabling bidirectional power flow between high-voltage transmission lines and local distribution networks, a feature that is essential for integrating renewable energy sources and ensuring reliable electrification in underserved regions. This paper presents a detailed mathematical representation of power flow within the CCS-CNT system, emphasizing the control of both active and reactive power through the adjustment of voltage levels and phase angles. A control algorithm is developed to dynamically manage power flow, ensuring optimal performance by minimizing losses and maintaining voltage stability across the network. The proposed CCS-CNT system demonstrates significant potential in enhancing the efficiency and reliability of power distribution, making it particularly suited for rural electrification and other applications where traditional methods fall short. The findings underscore the system's capability to adapt to varying operational conditions, offering a robust solution for modern power distribution challenges.
基金Project supported by the National Natural Science Foundation of China (Nos. 50934002 and 50774011)the Postdoctoral Science Foundation of China (No. 20090450014)the DoctoralNatural Science Foundation of China (No. 20070008038)
文摘A mathematical model for heap bioleaching is developed to analyze heat transfer, oxygen flow, target ion distribution and oxidation leaching rate in the heap. The model equations are solved with Comsol Multiphysics software. Numerical simulation results show the following facts： Concentration of oxygen is relatively high along the boundary of the slope, and low in the center part where leaching rate is slow. Temper- ature is relatively low along the slope and reaches the highest along the bottom region near the slope, with difference being more than 6℃. Concentration of target mental ions is the highest in the bottom region near the slope. Oxidation leaching rate is relatively large in the bottom and slope part with a fast reaction rate, and small in the other part with low oxygen concentration.
基金supported by the Project of the National Key R&D Program(Grant No.2021YFA1000202)National Natural Science Foundation of China(Grant Nos.12120101001,12001326 and 12171283)+2 种基金Natural Science Foundation of Shandong Province(Grant Nos.ZR2021ZD03,ZR2020QA032 and ZR2019ZD42)China Postdoctoral Science Foundation(Grant Nos.BX20190191 and 2020M672038)the Startup Fund from Shandong University(Grant No.11140082063130)。
文摘In this paper,we first establish a new fractional magnetohydrodynamic(MHD)coupled flow and heat transfer model for a generalized second-grade fluid.This coupled model consists of a fractional momentum equation and a heat conduction equation with a generalized form of Fourier law.The second-order fractional backward difference formula is applied to the temporal discretization and the Legendre spectral method is used for the spatial discretization.The fully discrete scheme is proved to be stable and convergent with an accuracy of O(τ^(2)+N-r),whereτis the time step-size and N is the polynomial degree.To reduce the memory requirements and computational cost,a fast method is developed,which is based on a globally uniform approximation of the trapezoidal rule for integrals on the real line.The strict convergence of the numerical scheme with this fast method is proved.We present the results of several numerical experiments to verify the effectiveness of the proposed method.Finally,we simulate the unsteady fractional MHD flow and heat transfer of the generalized second-grade fluid through a porous medium.The effects of the relevant parameters on the velocity and temperature are presented and analyzed in detail.
基金supported by the National Basic Research Program of China (2010CB951002)the Knowledge Innovation Project of Chinese Academy of Sciences (KZXC2-YW-BR-12)
文摘This paper presents the realization of two-way coupling of the unsaturated-saturated flow interactions of the SWAT2000 and MODFLOW96 models on the basis of the integrated surface/groundwater model SWATMOD99, and its application in Hetao Irrigation District （HID）, Inner Mongolia, China. Major revisions and enhancements were made to the SWAT2000 and MODFLOW models for simulating the detailed hydrologic budget and coupled unsaturated and saturated interactions, and irrigation canal hydrology for the HID. The simulation results of seasonal groundwater recharge to and evaporate from the shallow groundwater, and the annual water budget over the district are presented and discussed. The results implied the necessity of two-way coupling of the unsaturated-saturated interactions when groundwater is shallow, and the feasibility of making comprehensive use of the information coming from both the surface water and groundwater models to make a more physically-based assessment of the coupled interactions.
基金Project(SKLCRSM13KFB05)supported by State Key Laboratory for Coal Resources and Safe Mining(China University of Mining&Technology)
文摘Cemented tailings backfill(CTB) is made by mixing cement, tailings and water together, thus cement hydration and water seepage flow are the two crucial factors affecting the quality of CTB. Cement hydration process can release significant amount of heat to raise the temperature of CTB and in turn increase the rate of cement hydration. Meanwhile, the progress of cement hydration consumes water and produces hydration products to change the pore structures within CTB, which further influences the hydraulic behavior of CTB. In order to understand the hydraulic behavior of CTB, a numerical model was developed by coupling the hydraulic,thermal and hydration equations. This model was then implemented into COMSOL Multiphysics to simulate the evolutions of temperature and water seepage flow within CTB versus curing time. The predicted outcomes were compared with correspondent experimental results, proving the validity and availability of this model. By taking advantage of the validated model, effects of various initial CTB and curing temperatures, cement content, and CTB's geometric shapes on the hydraulic behavior of CTB were demonstrated numerically. The presented conclusions can contribute to preparing more environmentally friendly CTB structures.
