Casting microstructure evolution is difficult to describe quantitatively by only a separate simulation of dendrite scale or grain scale, and the numerical simulation of these two scales is difficult to render compatib...Casting microstructure evolution is difficult to describe quantitatively by only a separate simulation of dendrite scale or grain scale, and the numerical simulation of these two scales is difficult to render compatible. A three-dimensional cellular automaton model couplling both dendritic scale and grain scale is developed to simulate the microstructure evolution of the nickel-based single crystal superalloy DD406. Besides, a macro–mesoscopic/microscopic coupling solution algorithm is proposed to improve computational efficiency. The simulation results of dendrite growth and grain growth of the alloy are obtained and compared with the results given in previous reports. The results show that the primary dendritic arm spacing and secondary dendritic arm spacing of the dendritic growth are consistent with the theoretical and experimental results. The mesoscopic grain simulation can be used to obtain results similar to those of microscopic dendrites simulation. It is indicated that the developed model is feasible and effective.展开更多
In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and t...In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.展开更多
The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatu...The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatures,pressures and complex geological environments of deep strata frequently result in the coupling of multiple physical fields,including mechanical,thermal and hydraulic fields,during the fracturing of rocks.This review initially presents an overview of the coupling mechanisms of these physical fields,thereby elucidating the interaction processes ofmechanical,thermal,and hydraulic fields within rockmasses.Secondly,an in-depth analysis ofmulti-field coupling is conducted from both spatial and temporal perspectives,with the introduction of simulation methods for a range of scales.It emphasizes cross-scale coupling methodologies for the transfer of rock properties and physical field data,including homogenization techniques,nested coupling strategies and data-driven approaches.To address the discontinuous characteristics of the rock fracture process,the review provides a detailed explanation of continuousdiscontinuous couplingmethods,to elucidate the evolution of rock fracturing and deformationmore comprehensively.In conclusion,the review presents a summary of the principal points,challenges and future directions of multi-field coupling simulation research.It also puts forward the potential of integrating intelligent algorithms with multi-scale simulation techniques to enhance the accuracy and efficiency of multi-field coupling simulations.This offers novel insights into multi-field coupling simulation analysis in deep rock masses.展开更多
Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is sti...Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the key...A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.展开更多
The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. ...The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.展开更多
Ceramic Matrix Composite （CMC） turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress flu...Ceramic Matrix Composite （CMC） turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation, Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration （CVI） pro- cessed SiC/SiC-CMC are predicted by a Representative Volume Element （RVE） model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.展开更多
A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the clas...A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the classical multi-scale method(MSM)and the method of variation of parameters.Assuming that the forced excitation is a constant,one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM.Then,this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation.To certify the correctness and precision of the proposed analytical method,the van der Pol system with two scales subject to slowly periodic excitation is investigated;this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS),and quiescence(QS)responses.The approximate analytical expressions of the three types of responses are given by the MMSM,and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM).The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively,particularly when the excitation frequency is far less than the natural frequency of the system.展开更多
A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed cou...A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.展开更多
An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted F...An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.展开更多
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.展开更多
One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary condit...One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpola- tion can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin（EFG）-radial point interpola- tion method （RPIM） is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.展开更多
The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the c...The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.展开更多
A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(F...A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.展开更多
The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surfa...The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surface conditions to regional climate models( RCMs). There are two methods of downscaling: offline coupling and online coupling. The two kinds of coupling methods are described in detail by coupling the Weather Research and Forecasting model( WRF) with the Institute of Atmospheric Physics of Chinese Academy of Sciences Atmospheric General Circulation Model Version 4. 0( IAP AGCM4. 0) in the study. And the extreme precipitation event over Beijing on July 212012 is simulated by using the two coupling methods. Results show that online coupling method is of great value in improving the model simulation. Furthermore,the data exchange frequency of online coupling has some effect on simulation result.展开更多
A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radi...A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radial and axial magnetic pressure acting on the tube, the magnetic pressure is then used as boundary conditions to model the high velocity deformation of tube with DYNAFORM, The radial magnetic pressure on the tube decreases from the center to the tube end, axial magnetic pressure is greater near the location equal to the coil height and slight in the other region. The radial displacement of deformed workpicces is distributed uniformly near the tube center and decreases from the center to the end; Deformation from the location equal to coil height to the tube end is little. This distribution is consistent with the distribution of radial pressure; Effect of the axial magnetic pressure on deformation can be ignored, The calculated results show well agreements with the experimental results.展开更多
Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the be...Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the best catalytic performance among all samples. In addition, the effects of different addition sequences of Na, W, Mn, Zr, S and P on the catalytic performance were studied. The absence of Na before the addition of Mn and Zr in the catalysts preparation depressed the formation of the active phases of Mn2O3 and ZrO2 and decreased the activities of the catalysts significantly.展开更多
In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research gr...In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics.Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples.展开更多
Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applie...Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot＇s consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that （1） the effect of the stress path should be taken into account in the simulation of the soil consolidation question; （2） this decoupling method can predict the evolvement of pore water pressure; （3） the settlement using cam-clay model is less than that using numerical model because of dilatancy.展开更多
文摘Casting microstructure evolution is difficult to describe quantitatively by only a separate simulation of dendrite scale or grain scale, and the numerical simulation of these two scales is difficult to render compatible. A three-dimensional cellular automaton model couplling both dendritic scale and grain scale is developed to simulate the microstructure evolution of the nickel-based single crystal superalloy DD406. Besides, a macro–mesoscopic/microscopic coupling solution algorithm is proposed to improve computational efficiency. The simulation results of dendrite growth and grain growth of the alloy are obtained and compared with the results given in previous reports. The results show that the primary dendritic arm spacing and secondary dendritic arm spacing of the dendritic growth are consistent with the theoretical and experimental results. The mesoscopic grain simulation can be used to obtain results similar to those of microscopic dendrites simulation. It is indicated that the developed model is feasible and effective.
基金supported by Science and Technology on Reactor System Design Technology Laboratory,Chengdu,China(LRSDT2020106)
文摘In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.
基金supported by the National Natural Science Foundation of China(Grant Nos.42477185,41602308)the Zhejiang Provincial Natural Science Foundation of China(Grant No.LY20E080005)the Postgraduate Course Construction Project of Zhejiang University of Science and Technology(Grant No.2021yjskj05).
文摘The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatures,pressures and complex geological environments of deep strata frequently result in the coupling of multiple physical fields,including mechanical,thermal and hydraulic fields,during the fracturing of rocks.This review initially presents an overview of the coupling mechanisms of these physical fields,thereby elucidating the interaction processes ofmechanical,thermal,and hydraulic fields within rockmasses.Secondly,an in-depth analysis ofmulti-field coupling is conducted from both spatial and temporal perspectives,with the introduction of simulation methods for a range of scales.It emphasizes cross-scale coupling methodologies for the transfer of rock properties and physical field data,including homogenization techniques,nested coupling strategies and data-driven approaches.To address the discontinuous characteristics of the rock fracture process,the review provides a detailed explanation of continuousdiscontinuous couplingmethods,to elucidate the evolution of rock fracturing and deformationmore comprehensively.In conclusion,the review presents a summary of the principal points,challenges and future directions of multi-field coupling simulation research.It also puts forward the potential of integrating intelligent algorithms with multi-scale simulation techniques to enhance the accuracy and efficiency of multi-field coupling simulations.This offers novel insights into multi-field coupling simulation analysis in deep rock masses.
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
基金Projects (51175162, 50805045) supported by the National Natural Science Foundation of ChinaProject supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education,China
文摘A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.
基金Sponsored by National Science and Technology Support Plan of China (2009AA04Z143)Science and Technology Support Plan of Hebei Province of China (10212101D)Important Natural Science Foundation of Hebei Province of China (E2006001038)
文摘The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.
文摘Ceramic Matrix Composite （CMC） turbine guide vanes possess multi-scale stress and strain with inhomogeneity at the microscopic scale. Given that the macroscopic distribution cannot reflect the microscopic stress fluctuation, the macroscopic method fails to meet the requirements of stress and strain analysis of CMC turbine guide vanes. Furthermore, the complete thermodynamic properties of 2D woven SiC/SiC-CMC cannot be obtained through experimentation, Accordingly, a method to calculate the thermodynamic properties of CMC and analyze multi-scale stress and strain of the turbine guide vanes should be established. In this study, the multi-scale thermodynamic analysis is investigated. The thermodynamic properties of Chemical Vapor Infiltration （CVI） pro- cessed SiC/SiC-CMC are predicted by a Representative Volume Element （RVE） model with porosity, leading to the result that the relative error between the calculated in-plane tensile modulus and the experimental value is 4.2%. The macroscopic response of a guide vane under given conditions is predicted. The relative error between the predicted strain on the trailing edge and the experimental value is 9.7%. The calculation of the stress distribution of micro-scale RVE shows that the maximum value of microscopic stress, which is located in the interlayer matrix, is more than 1.5 times that of macroscopic stress in the same direction and the microscopic stress distribution of the interlayer matrix is related to the pore distribution of the composite.
基金the National Natural Science Foundation of China(Nos.11672191,11772206,and U1934201)the Hundred Excellent Innovative Talents Support Program in Hebei University(No.SLRC2017053)。
文摘A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the classical multi-scale method(MSM)and the method of variation of parameters.Assuming that the forced excitation is a constant,one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM.Then,this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation.To certify the correctness and precision of the proposed analytical method,the van der Pol system with two scales subject to slowly periodic excitation is investigated;this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS),and quiescence(QS)responses.The approximate analytical expressions of the three types of responses are given by the MMSM,and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM).The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively,particularly when the excitation frequency is far less than the natural frequency of the system.
基金Projects(51474251,51874351) supported by the National Natural Science Foundation of China
文摘A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.
基金Project (No. 50575143) supported by the National Natural Science Foundation of China
文摘An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.
基金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.
基金supported by the National Natural Science Foundation of China (No. 11172192)the College Postgraduate Research and Innovation Project of Jiangsu Province (No. CX10B 029Z)the Nominated Excellent Thesis for PHD Candidates Program of Soochow University (No. 23320957)
文摘One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpola- tion can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin（EFG）-radial point interpola- tion method （RPIM） is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.
基金supported by the National Science and Technology Major Project,China(No.2017-III-0010-0036).
文摘The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.
基金co-supported by the National Key R&D Program of China(Nos.2018YFF01012401-05)the National Natural Science Foundation of China(No.51975243)+2 种基金Jilin Provincial Department of Education(No.JJKH20180084KJ),Chinathe Fundamental Research Funds for the Central Universities and Jilin Provincial Department of Science&Technology Fund Project,China(Nos.20170101043JC and 20180520072JH)Graduate Innovation Fund of Jilin University,China(No.101832018C184).
文摘A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.
基金Supported by the National Natural Science Foundation of China(No.61602477)China Postdoctoral Science Foundation(No.2016M601158)National Key Research and Development Program of China(No.2016YFB0200804)
文摘The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surface conditions to regional climate models( RCMs). There are two methods of downscaling: offline coupling and online coupling. The two kinds of coupling methods are described in detail by coupling the Weather Research and Forecasting model( WRF) with the Institute of Atmospheric Physics of Chinese Academy of Sciences Atmospheric General Circulation Model Version 4. 0( IAP AGCM4. 0) in the study. And the extreme precipitation event over Beijing on July 212012 is simulated by using the two coupling methods. Results show that online coupling method is of great value in improving the model simulation. Furthermore,the data exchange frequency of online coupling has some effect on simulation result.
文摘A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radial and axial magnetic pressure acting on the tube, the magnetic pressure is then used as boundary conditions to model the high velocity deformation of tube with DYNAFORM, The radial magnetic pressure on the tube decreases from the center to the tube end, axial magnetic pressure is greater near the location equal to the coil height and slight in the other region. The radial displacement of deformed workpicces is distributed uniformly near the tube center and decreases from the center to the end; Deformation from the location equal to coil height to the tube end is little. This distribution is consistent with the distribution of radial pressure; Effect of the axial magnetic pressure on deformation can be ignored, The calculated results show well agreements with the experimental results.
基金supported by the financial support from National Natural Science Foundation of China (20676116)
文摘Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the best catalytic performance among all samples. In addition, the effects of different addition sequences of Na, W, Mn, Zr, S and P on the catalytic performance were studied. The absence of Na before the addition of Mn and Zr in the catalysts preparation depressed the formation of the active phases of Mn2O3 and ZrO2 and decreased the activities of the catalysts significantly.
基金The project supported by the 973 Program (G1999033108)the National Natural Science Foundation of China (10125211)
文摘In this paper,we present the applications of Boundary Element Method(BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems(MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics.Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples.
文摘Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot＇s consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that （1） the effect of the stress path should be taken into account in the simulation of the soil consolidation question; （2） this decoupling method can predict the evolvement of pore water pressure; （3） the settlement using cam-clay model is less than that using numerical model because of dilatancy.
