A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite...A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.展开更多
To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature...To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.展开更多
The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and nu...The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and numerical calculations.It is found that the longitudinal magnetic field can affect the heating by changing the level of the pressure heating along the longitudinal direction and that of the Ohmic heating along the direction which is perpendicular to both driving electric field and the applied transverse magnetic field,and a continuously increased longitudinal magnetic field can induce pressure heating to become dominant.Moreover,the electron temperature as well as proportion of some low energy electrons will increase if a small longitudinal magnetic field is introduced,which is attributed to the increased average electron energy.We believe that the research will provide guidance for optimizing the magnetic field configuration of some discharge systems having both transverse and longitudinal magnetic field.展开更多
The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separat...The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separate energy equations for the gas and solid phases coupled with an experimental kinetic model. A lumped kinetic model containing four main species CH4, O2, COx （CO2, CO）, and C2 （C2H4 and C2H6） was used with a plug flow reactor model as well. The results from the model agreed with the experimental data. The model was used to analyze the influence of temperature and feed gas composition on the conversion and selectivity of the reactor performance. The analytical results indicate that the conversion decreases, whereas, C2 selectivity increases by increasing gas hourly space velocity （GHSV） and the methane conversion also decreases by increasing the methane to oxygen ratio.展开更多
Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal...Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal around the roadway was divided into the seepage open zone, seepage orientation zone, seepage decay zone and original seepage zone along the radial direction of the roadway. The loaded gassy coal was treated as a viscoelastic and plastic softened medium, and the mechanical behaviors of the viscoelastic zone, plastic softened zone and broken zone around the roadway were analyzed with the consideration of the loading creep, softening and expansion effect of the gassy coal. According to the law of conservation of mass and the Darcy law, the flow-solid coupled model for the gas transportation of the coal around the roadway was established considering the dynamic evolution of the adsorption characteristics, porosity and permeability of the coal, and the simulation software COMSOL was utilized to numerically simulate the stress state and gas flow regularity around the coal, which provided meaningful reference for investigating the stability of the coal and rock around the roadway.展开更多
Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface...Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.展开更多
In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing appr...In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.展开更多
Based on the theory of continuum mechanics of multi-pbase media, a mathematical model and non-linear FEM equation of the coupling instability problem of solid-fluid biphase media for coal-methane outburst under finite...Based on the theory of continuum mechanics of multi-pbase media, a mathematical model and non-linear FEM equation of the coupling instability problem of solid-fluid biphase media for coal-methane outburst under finite deformation are established. The critical conditions of the surface instability are presented as the singularity of the total stiffness matrices of the coal body for coal-methaue outburst. That means the deformtion or the coal body emerges bifurcatiou phenomena. The numerical simulation of a typical outburst is made.展开更多
Enhancing plasma uniformity can be achieved by modifying coil and chamber structures in radio frequency inductively coupled plasma(ICP)to meet the demand for large-area and uniformly distributed plasma in industrial m...Enhancing plasma uniformity can be achieved by modifying coil and chamber structures in radio frequency inductively coupled plasma(ICP)to meet the demand for large-area and uniformly distributed plasma in industrial manufacturing.This study utilized a two-dimensional self-consistent fluid model to investigate how different coil configurations and chamber aspect ratios affect the radial uniformity of plasma in radio frequency ICP.The findings indicate that optimizing the radial spacing of the coil enhances plasma uniformity but with a reduction in electron density.Furthermore,optimizing the coil within the ICP reactor,using the interior point method in the Interior Point Optimizer significantly enhances plasma uniformity,elevating it from 56%to 96%within the range of the model sizes.Additionally,when the chamber aspect ratio k changes from 2.8 to 4.7,the plasma distribution changes from a center-high to a saddleshaped distribution.Moreover,the plasma uniformity becomes worse.Finally,adjusting process parameters,such as increasing source power and gas pressure,can enhance plasma uniformity.These findings contribute to optimizing the etching process by improving plasma radial uniformity.展开更多
The technology of pantograph sinking in the cavity is generally adopted in the new generation of high-speed trains in China for aerodynamic noise reduction in this region. This study takes a high-speed train with a 4-...The technology of pantograph sinking in the cavity is generally adopted in the new generation of high-speed trains in China for aerodynamic noise reduction in this region. This study takes a high-speed train with a 4-car formation and a pantograph as the research object and compares the aerodynamic acoustic performance of two scale models, 1/8 and 1/1, using large eddy simulation and Ffowcs Williams–Hawkings integral equation. It is found that there is no direct scale similarity between their aeroacoustic performance. The 1/1 model airflow is separated at the leading edge of the panhead and reattached to the panhead, and its vortex shedding Strouhal number(St) is 0.17. However, the 1/8 model airflow is separated directly at the leading edge of the panhead, and its St is 0.13. The cavity's vortex shedding frequency is in agreement with that calculated by the Rooster empirical formula. The two scale models exhibit some similar characteristics in distribution of sound source energy, but the energy distribution of the 1/8 model is more concentrated in the middle and lower regions. The contribution rates of their middle and lower regions to the radiated noise in the two models are 27.3% and 87.2%, respectively. The peak frequencies of the radiated noise from the 1/1 model are 307 and 571 Hz. The 307 Hz is consistent with the frequency of panhead vortex shedding, and the 571 Hz is more likely to be the result of the superposition of various components. In contrast, the peak frequencies of the radiated noise from the 1/8 scale model are 280 and 1970 Hz. The 280 Hz comes from the shear layer oscillation between the cavity and the bottom frame, and the 1970 Hz is close to the frequency at which the panhead vortex sheds. This shows that the scaled model results need to be corrected before applying to the full-scale model.展开更多
A three-dimensional geometric model of the oxidative coupling of methane （OCM） packed-bed reactor loaded with Na2WO4-Mn/SiO2 partic- ulate catalyst was set up, and an improved Stansch kinetic model was established t...A three-dimensional geometric model of the oxidative coupling of methane （OCM） packed-bed reactor loaded with Na2WO4-Mn/SiO2 partic- ulate catalyst was set up, and an improved Stansch kinetic model was established to calculate the OCM reactions using the computational fluid dynamics method and Fluent software. The simulation conditions were completely the same with the experimental conditions that the volume velocity of the reactant was 80 mL/min under standard state, the ratio of CH4/O2 was 3, the temperature and pressure were 800 ℃ and 1 atm, respectively. The contour of the characteristics parameters in the catalyst bed was analyzed, such as the species mass fractions, temperature, the heat flux on side wall surface, pressure, fluid density and velocity. The results showed that the calculated values matched well with the experimental values on the conversion of CH4 and the selectivity to products （C2H6, C2H4, CO2, CO） in the reactor outlet with an error range of 4-2%. The mass fractions of CH4 and O2 decreased from 0.6 and 0.4 in the catalyst bed inlet to 0.436 and 0.142 in the outlet, where the mass fractions of C2H6, C2H4, CO and CO2 were 0.035, 0.061, 0.032 and 0.106, respectively. Due to the existence of laminar boundary layer, the contours of each component bent upwards in the vicinity of the boundary layer. This OCM reaction was volume increase reaction and the total moles of products were greater than those of reactants. The flow field in the catalyst bed maintained constant temperature and pressure. The fluid density decreased gradually from 2.28 kg/m3 in the inlet of the catalyst bed to 2.22 kg/m3 in the outlet of the catalyst bed, while the velocity increased from 0.108 m/s to 0.115 m/s.展开更多
We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numeri...We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numerical simulation and experimental validation. The distribution characteristics and interaction mechanism of the flow field and electromagnetic field of the ICP wind tunnel are investigated using the multi-field coupling method of flow, electromagnetic, chemical, and thermodynamic field. The accuracy of the numerical simulation is validated by comparing the experimental results with the simulation results. Thereafter, the wind tunnel pressure, air velocity, electron density, Joule heating rate, Lorentz force, and electric field intensity obtained using the simulation are analyzed and discussed. The results indicate that for the 1.2-MW ICP wind tunnel, the maximum values of temperature, pressure, electron number density, and other parameters are observed during coil heating. The influence of the radial Lorentz force on the momentum transfer is stronger than that of the axial Lorentz force. The electron number density at the central axis and the amplitude and position of the Joule heating rate are affected by the radial Lorentz force. Moreover, the plasma in the wind tunnel is constantly in the subsonic flow state, and a strong eddy flow is easily generated at the inlet of the wind tunnel.展开更多
The virtual prototyping models of the mechanical, hydraulic and control system of the ITER tractor were built with CATIA, ADAMS and MATLAB/Simulink respectively according to its heavy load and high precision character...The virtual prototyping models of the mechanical, hydraulic and control system of the ITER tractor were built with CATIA, ADAMS and MATLAB/Simulink respectively according to its heavy load and high precision characteristics, and the data transfer between the different models was accomplished by the integration interface between different software. Consequently the virtual experimental platform for the multi-disciplinary co-simulation was established. A co-simulation study of the mechanical-hydraulic-control coupling system of the ITER tractor was carried out. The synchronization servo control of parallel hydraulic cylinders was implemented, and the tracking control of the preconcerted trajectory of the hydraulic cylinders was realized on the established experimental platform. This paper presents the optimization design and technology rebuilding for the complicated coupling system with its theoretic foundation and co-simulation virtual experimental platform.展开更多
The fluid-solid coupling theory, an interdisciplinary science between hydrodynamics and solid mechanics, is an important tool for response analysis and direct design of structures in naval architecture and ocean engin...The fluid-solid coupling theory, an interdisciplinary science between hydrodynamics and solid mechanics, is an important tool for response analysis and direct design of structures in naval architecture and ocean engineering. By applying the corresponding relations between generalized forces and generalized displacements, convolutions were performed between the basic equations of elasto-dynamics in the primary space and corresponding virtual quantities. The results were integrated and then added algebraically. In light of the fact that body forces and surface forces are both follower forces, the generalized quasi-complementary energy principle with two kinds of variables for an initial value problem is established in non-conservative systems. Using the generalized quasi-complementary energy principle to deal with the fluid-solid coupling problem and to analyze the dynamic response of structures, a method for using two kinds of variables simultaneously for calculation of force and displacement was derived.展开更多
In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The el...In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The electric current density, temperature and velocity distributions in the torch are obtained through the coupled iterative calculation about the electromagnetic equations described in a magnetic vector potential format and the modified fluid dynamics equations. The fluid-solid coupled calculation method is applied to guarantee the continuity of the electric current and heat transfer at the interface between the electrodes and fluid. The predicted location of the anodic arc root attachment and the arc voltage of the torch are consistent with corresponding experimental results. Through a specific analysis of the influence of mass flow rates and electric current on the torch outlet parameters, the total thermal efficiency, thermal loss of each part, and the laws of the variation of outlet parameters with the variation of mass flow rates and electric current was obtained. It is found that operation under a laminar condition with a limited area of the anode could increase the total thermal efficiency of the torch.展开更多
Numerical simulation studies on blasting in ore-pass with charge both 250 kg and 500 kg are carried out. Adopting fluid-solid coupling method, the whole process from detonating to shockwave transference in the ore-pas...Numerical simulation studies on blasting in ore-pass with charge both 250 kg and 500 kg are carried out. Adopting fluid-solid coupling method, the whole process from detonating to shockwave transference in the ore-pass reappears. The air shockwave＇s attenuation, chocking-up surface pressure and ore-wall press are analyzed. The results show that the biggest pressure of ore-wall happened at 0.43 m height, and the peak pressure are 167.0 MPa （250 kg） and 341.1 MPa （500 kg） individually. The choking-up zone＇s stress state from center to edge represents transition from press to strain, and the strain/press stress peak values are 9.2/ 4.8 MPa （250 kg） and 17.2/12.5 MPa （500 kg）. The shockwave＇s attenuation values got from both experiential formula and simulation are compared, the experiential formula values are 2.54 - 3.49 times of the simulation results.展开更多
A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fractu...A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.展开更多
The coupling of lower hybrid wave to the plasma is a crucial issue for efficient current drive in tokamaks. This paper establishes a new coupling model which assumes the antenna to be a curved face and the plasma to b...The coupling of lower hybrid wave to the plasma is a crucial issue for efficient current drive in tokamaks. This paper establishes a new coupling model which assumes the antenna to be a curved face and the plasma to be a cylinder. Power spectrum considering the coupling between wave-guides in both poloidal and toroidal direction is simply estimated and discussed. The effect of the poloidal wave vector on wave propagation, power deposition and driven current is also investigated with the help of lower hybrid current drive code. Results show that the poloidal wave vector affects the ray tracing, and also has effect on power deposition and driven current. The effect of the poloidal wave vector on power deposition and driven current profile depends on plasma parameters. Preliminary studies suggest that it seems possible to control the current profile by adjusting the poloidal phase difference between the waveguide in poloidal direction.展开更多
The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, ...The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, error function and normal distribution function, and a series of FEM equations of coupled fields of temperature and seepage were deduced and put forward. With the example of shaft ground freezing, the formation of freezing wall in seepage ground was simulated.展开更多
基金Supported by the National Science and Technology Major Project (2017ZX05013-005)。
文摘A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.
基金financially supported by National Natural Science Foundation of China(NSFC)(Nos.12105035 and 11935005)Guangdong Basic and Applied Basic Research Foundation(No.2021B1515120018)+1 种基金the Fundamental Research Funds for the Central Universities(No.DUT21TD104)the Advanced Space Propulsion Laboratory of BICE and Beijing Engineering Research Center of Efficient and Green Aerospace Propulsion Technology(No.Lab ASP-2020-01)。
文摘To maintain the high-density plasma source in inductively coupled plasma(ICP),very high radiofrequency power is often delivered to the antenna,which can heat the dielectric windows near the antenna to high temperature.This high temperature can modulate the plasma characteristics to a large degree.We thus study the effect of dielectric window temperature on plasma parameters in two different ICP structures based on COMSOL software.The distributions of various plasma species are examined at different dielectric window temperatures.The concentration of neutral gas is found to be largely modulated at high dielectric window temperature,which further affects the electron collision probability with neutrals and the electron temperature.However,the electron density profiles are barely affected by the dielectric window temperature,which is mainly concentrated at the center of the reactor due to the fixed power input and pressure.
基金supported by National Natural Science Foundation of China (Nos. 11975163 and 12175160)a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)
文摘The electron heating characteristics of magnetic enhancement capacitively coupled argon plasmas in presence of both longitudinal and transverse uniform magnetic field have been explored through both theoretical and numerical calculations.It is found that the longitudinal magnetic field can affect the heating by changing the level of the pressure heating along the longitudinal direction and that of the Ohmic heating along the direction which is perpendicular to both driving electric field and the applied transverse magnetic field,and a continuously increased longitudinal magnetic field can induce pressure heating to become dominant.Moreover,the electron temperature as well as proportion of some low energy electrons will increase if a small longitudinal magnetic field is introduced,which is attributed to the increased average electron energy.We believe that the research will provide guidance for optimizing the magnetic field configuration of some discharge systems having both transverse and longitudinal magnetic field.
文摘The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separate energy equations for the gas and solid phases coupled with an experimental kinetic model. A lumped kinetic model containing four main species CH4, O2, COx （CO2, CO）, and C2 （C2H4 and C2H6） was used with a plug flow reactor model as well. The results from the model agreed with the experimental data. The model was used to analyze the influence of temperature and feed gas composition on the conversion and selectivity of the reactor performance. The analytical results indicate that the conversion decreases, whereas, C2 selectivity increases by increasing gas hourly space velocity （GHSV） and the methane conversion also decreases by increasing the methane to oxygen ratio.
基金the financial support from the National Natural Science Foundation for Young Scientists of China (Nos.51604116 and 51604096)Natural Science Foundation ofHenbei Province (No.E2016508036)+1 种基金Hebei State Key Laboratory of Mine Disaster Prevention (No.KJZH2017K08)Basic and Frontier Technology Research Project of Henan Province in 2016 (No.162300410031)
文摘Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal around the roadway was divided into the seepage open zone, seepage orientation zone, seepage decay zone and original seepage zone along the radial direction of the roadway. The loaded gassy coal was treated as a viscoelastic and plastic softened medium, and the mechanical behaviors of the viscoelastic zone, plastic softened zone and broken zone around the roadway were analyzed with the consideration of the loading creep, softening and expansion effect of the gassy coal. According to the law of conservation of mass and the Darcy law, the flow-solid coupled model for the gas transportation of the coal around the roadway was established considering the dynamic evolution of the adsorption characteristics, porosity and permeability of the coal, and the simulation software COMSOL was utilized to numerically simulate the stress state and gas flow regularity around the coal, which provided meaningful reference for investigating the stability of the coal and rock around the roadway.
基金supported in part by R & D Department of China National Petroleum Corporation (2022DQ0604-01)National Natural Science Foundation of China (42204132)+3 种基金the China Postdoctoral Science Foundations (2020M680667, 2021T140661)Harvard-CUP Joint Laboratory on Petroleum Science“111” project (B13010)the financial support from the CAS Special Research Assistant Project。
文摘Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.
文摘In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.
文摘Based on the theory of continuum mechanics of multi-pbase media, a mathematical model and non-linear FEM equation of the coupling instability problem of solid-fluid biphase media for coal-methane outburst under finite deformation are established. The critical conditions of the surface instability are presented as the singularity of the total stiffness matrices of the coal body for coal-methaue outburst. That means the deformtion or the coal body emerges bifurcatiou phenomena. The numerical simulation of a typical outburst is made.
基金supported by the Scientific Research Foundation of Xijing University,China(No.XJ19T03)the Opening Project of Science and Technology on Reliability Physics and Application Technology of Electronic Component Laboratory(No.ZHD201701)the Natural Science Basic Research Plan in Shaanxi Province of China(No.2024JC-YBMS-342).
文摘Enhancing plasma uniformity can be achieved by modifying coil and chamber structures in radio frequency inductively coupled plasma(ICP)to meet the demand for large-area and uniformly distributed plasma in industrial manufacturing.This study utilized a two-dimensional self-consistent fluid model to investigate how different coil configurations and chamber aspect ratios affect the radial uniformity of plasma in radio frequency ICP.The findings indicate that optimizing the radial spacing of the coil enhances plasma uniformity but with a reduction in electron density.Furthermore,optimizing the coil within the ICP reactor,using the interior point method in the Interior Point Optimizer significantly enhances plasma uniformity,elevating it from 56%to 96%within the range of the model sizes.Additionally,when the chamber aspect ratio k changes from 2.8 to 4.7,the plasma distribution changes from a center-high to a saddleshaped distribution.Moreover,the plasma uniformity becomes worse.Finally,adjusting process parameters,such as increasing source power and gas pressure,can enhance plasma uniformity.These findings contribute to optimizing the etching process by improving plasma radial uniformity.
基金supported by the National Natural Science Foundation of China (No. 52272363)the Key Laboratory of Aerodynamic Noise Control (No. ANCL20200302),Aid Program for Science and Technology Innovative Research Team in Higher Educational Institutions of Hunan Province。
文摘The technology of pantograph sinking in the cavity is generally adopted in the new generation of high-speed trains in China for aerodynamic noise reduction in this region. This study takes a high-speed train with a 4-car formation and a pantograph as the research object and compares the aerodynamic acoustic performance of two scale models, 1/8 and 1/1, using large eddy simulation and Ffowcs Williams–Hawkings integral equation. It is found that there is no direct scale similarity between their aeroacoustic performance. The 1/1 model airflow is separated at the leading edge of the panhead and reattached to the panhead, and its vortex shedding Strouhal number(St) is 0.17. However, the 1/8 model airflow is separated directly at the leading edge of the panhead, and its St is 0.13. The cavity's vortex shedding frequency is in agreement with that calculated by the Rooster empirical formula. The two scale models exhibit some similar characteristics in distribution of sound source energy, but the energy distribution of the 1/8 model is more concentrated in the middle and lower regions. The contribution rates of their middle and lower regions to the radiated noise in the two models are 27.3% and 87.2%, respectively. The peak frequencies of the radiated noise from the 1/1 model are 307 and 571 Hz. The 307 Hz is consistent with the frequency of panhead vortex shedding, and the 571 Hz is more likely to be the result of the superposition of various components. In contrast, the peak frequencies of the radiated noise from the 1/8 scale model are 280 and 1970 Hz. The 280 Hz comes from the shear layer oscillation between the cavity and the bottom frame, and the 1970 Hz is close to the frequency at which the panhead vortex sheds. This shows that the scaled model results need to be corrected before applying to the full-scale model.
文摘A three-dimensional geometric model of the oxidative coupling of methane （OCM） packed-bed reactor loaded with Na2WO4-Mn/SiO2 partic- ulate catalyst was set up, and an improved Stansch kinetic model was established to calculate the OCM reactions using the computational fluid dynamics method and Fluent software. The simulation conditions were completely the same with the experimental conditions that the volume velocity of the reactant was 80 mL/min under standard state, the ratio of CH4/O2 was 3, the temperature and pressure were 800 ℃ and 1 atm, respectively. The contour of the characteristics parameters in the catalyst bed was analyzed, such as the species mass fractions, temperature, the heat flux on side wall surface, pressure, fluid density and velocity. The results showed that the calculated values matched well with the experimental values on the conversion of CH4 and the selectivity to products （C2H6, C2H4, CO2, CO） in the reactor outlet with an error range of 4-2%. The mass fractions of CH4 and O2 decreased from 0.6 and 0.4 in the catalyst bed inlet to 0.436 and 0.142 in the outlet, where the mass fractions of C2H6, C2H4, CO and CO2 were 0.035, 0.061, 0.032 and 0.106, respectively. Due to the existence of laminar boundary layer, the contours of each component bent upwards in the vicinity of the boundary layer. This OCM reaction was volume increase reaction and the total moles of products were greater than those of reactants. The flow field in the catalyst bed maintained constant temperature and pressure. The fluid density decreased gradually from 2.28 kg/m3 in the inlet of the catalyst bed to 2.22 kg/m3 in the outlet of the catalyst bed, while the velocity increased from 0.108 m/s to 0.115 m/s.
基金supported by the National Natural Science Foundation of China (Grant No. 11705143)the Open Foundation for Key Laboratories of National Defense Science and Technology of China (Grant No. 6142202031901)the Foundation for Research and Development of Applied Technology in Beilin District of Xi’an,China (Grant No. GX2047)。
文摘We take the established inductively coupled plasma(ICP) wind tunnel as a research object to investigate the thermal protection system of re-entry vehicles. A 1.2-MW high power ICP wind tunnel is studied through numerical simulation and experimental validation. The distribution characteristics and interaction mechanism of the flow field and electromagnetic field of the ICP wind tunnel are investigated using the multi-field coupling method of flow, electromagnetic, chemical, and thermodynamic field. The accuracy of the numerical simulation is validated by comparing the experimental results with the simulation results. Thereafter, the wind tunnel pressure, air velocity, electron density, Joule heating rate, Lorentz force, and electric field intensity obtained using the simulation are analyzed and discussed. The results indicate that for the 1.2-MW ICP wind tunnel, the maximum values of temperature, pressure, electron number density, and other parameters are observed during coil heating. The influence of the radial Lorentz force on the momentum transfer is stronger than that of the axial Lorentz force. The electron number density at the central axis and the amplitude and position of the Joule heating rate are affected by the radial Lorentz force. Moreover, the plasma in the wind tunnel is constantly in the subsonic flow state, and a strong eddy flow is easily generated at the inlet of the wind tunnel.
基金supported by design of the ITER transfer casks system (ITER International Team) ITA 23-01-CNthe Key Laboratory of Biomimetic Sensing and Advanced Robot Technology,Anhui Province,China
文摘The virtual prototyping models of the mechanical, hydraulic and control system of the ITER tractor were built with CATIA, ADAMS and MATLAB/Simulink respectively according to its heavy load and high precision characteristics, and the data transfer between the different models was accomplished by the integration interface between different software. Consequently the virtual experimental platform for the multi-disciplinary co-simulation was established. A co-simulation study of the mechanical-hydraulic-control coupling system of the ITER tractor was carried out. The synchronization servo control of parallel hydraulic cylinders was implemented, and the tracking control of the preconcerted trajectory of the hydraulic cylinders was realized on the established experimental platform. This paper presents the optimization design and technology rebuilding for the complicated coupling system with its theoretic foundation and co-simulation virtual experimental platform.
基金Supported by the National Natural Science Foundation under Grant No.10272034the Doctoral Education Foundation under Grant No.20060217020
文摘The fluid-solid coupling theory, an interdisciplinary science between hydrodynamics and solid mechanics, is an important tool for response analysis and direct design of structures in naval architecture and ocean engineering. By applying the corresponding relations between generalized forces and generalized displacements, convolutions were performed between the basic equations of elasto-dynamics in the primary space and corresponding virtual quantities. The results were integrated and then added algebraically. In light of the fact that body forces and surface forces are both follower forces, the generalized quasi-complementary energy principle with two kinds of variables for an initial value problem is established in non-conservative systems. Using the generalized quasi-complementary energy principle to deal with the fluid-solid coupling problem and to analyze the dynamic response of structures, a method for using two kinds of variables simultaneously for calculation of force and displacement was derived.
基金supported by National Natural Science Foundation of China (No. 50476081)
文摘In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The electric current density, temperature and velocity distributions in the torch are obtained through the coupled iterative calculation about the electromagnetic equations described in a magnetic vector potential format and the modified fluid dynamics equations. The fluid-solid coupled calculation method is applied to guarantee the continuity of the electric current and heat transfer at the interface between the electrodes and fluid. The predicted location of the anodic arc root attachment and the arc voltage of the torch are consistent with corresponding experimental results. Through a specific analysis of the influence of mass flow rates and electric current on the torch outlet parameters, the total thermal efficiency, thermal loss of each part, and the laws of the variation of outlet parameters with the variation of mass flow rates and electric current was obtained. It is found that operation under a laminar condition with a limited area of the anode could increase the total thermal efficiency of the torch.
基金Sponsored by the National Natural Science Foundation of China(50874109)
文摘Numerical simulation studies on blasting in ore-pass with charge both 250 kg and 500 kg are carried out. Adopting fluid-solid coupling method, the whole process from detonating to shockwave transference in the ore-pass reappears. The air shockwave＇s attenuation, chocking-up surface pressure and ore-wall press are analyzed. The results show that the biggest pressure of ore-wall happened at 0.43 m height, and the peak pressure are 167.0 MPa （250 kg） and 341.1 MPa （500 kg） individually. The choking-up zone＇s stress state from center to edge represents transition from press to strain, and the strain/press stress peak values are 9.2/ 4.8 MPa （250 kg） and 17.2/12.5 MPa （500 kg）. The shockwave＇s attenuation values got from both experiential formula and simulation are compared, the experiential formula values are 2.54 - 3.49 times of the simulation results.
基金Supported by the National Natural Science Foundation of China(51827804)CNPC Strategic Cooperation Science and Technology Major Project(ZLZX2020-01-05)Open Fund of State Key Laboratory of Rock Mechanics and Engineering(SKLGME021024).
文摘A fracture propagation model of radial well fracturing is established based on the finite element-meshless method.The model considers the coupling effect of fracturing fluid flow and rock matrix deformation.The fracture geometries of radial well fracturing are simulated,the induction effect of radial well on the fracture is quantitatively characterized,and the influences of azimuth,horizontal principle stress difference,and reservoir matrix permeability on the fracture geometries are revealed.The radial wells can induce the fractures to extend parallel to their axes when two radial wells in the same layer are fractured.When the radial wells are symmetrically distributed along the direction of the minimum horizontal principle stress with the azimuth greater than 15,the extrusion effect reduces the fracture length of radial wells.When the radial wells are symmetrically distributed along the direction of the maximum horizontal principal stress,the extrusion increases the fracture length of the radial wells.The fracture geometries are controlled by the rectification of radial borehole,the extrusion between radial wells in the same layer,and the deflection of the maximum horizontal principal stress.When the radial wells are distributed along the minimum horizontal principal stress symmetrically,the fracture length induced by the radial well decreases with the increase of azimuth;in contrast,when the radial wells are distributed along the maximum horizontal principal stress symmetrically,the fracture length induced by the radial well first decreases and then increases with the increase of azimuth.The fracture length induced by the radial well decreases with the increase of horizontal principal stress difference.The increase of rock matrix permeability and pore pressure of the matrix around radial wells makes the inducing effect of the radial well on fractures increase.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.10575104 and 10875149)Dean Foundation of Hefei Institute of Physical Science,Chinese Academy of Sciences
文摘The coupling of lower hybrid wave to the plasma is a crucial issue for efficient current drive in tokamaks. This paper establishes a new coupling model which assumes the antenna to be a curved face and the plasma to be a cylinder. Power spectrum considering the coupling between wave-guides in both poloidal and toroidal direction is simply estimated and discussed. The effect of the poloidal wave vector on wave propagation, power deposition and driven current is also investigated with the help of lower hybrid current drive code. Results show that the poloidal wave vector affects the ray tracing, and also has effect on power deposition and driven current. The effect of the poloidal wave vector on power deposition and driven current profile depends on plasma parameters. Preliminary studies suggest that it seems possible to control the current profile by adjusting the poloidal phase difference between the waveguide in poloidal direction.
文摘The coupling mechanism in freezing process of seepage ground was studied and a simplified coupling math model was proposed. The nonlinear and coupling problems of PDEs were well solved using the exponential function, error function and normal distribution function, and a series of FEM equations of coupled fields of temperature and seepage were deduced and put forward. With the example of shaft ground freezing, the formation of freezing wall in seepage ground was simulated.
