This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the ...This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.展开更多
Hydroformed parts are widely used in industrial automotive parts because of their higher stiffness and fatigue strength and reduced weight relative to their corresponding cast and welded parts.This paper reports a hyd...Hydroformed parts are widely used in industrial automotive parts because of their higher stiffness and fatigue strength and reduced weight relative to their corresponding cast and welded parts.This paper reports a hydraulicforming experimental platform for rectangular tube fittings that was constructed to conduct an experiment on the hydraulic forming of rectangular tube fittings.A finite element model was established on the basis of the fluid–solid coupling method and simulation analysis.The correctness of the simulation analysis and the feasibility of the fluid–solid coupling method for hydraulic forming simulation analysis were verified by comparing the experimental results with the simulation results.On the basis of the simulation analysis of the hydraulic process of the torsion beam using the fluid–solid coupling method,a sliding mold suitable for the hydroforming of torsion beams was designed for its structural characteristics.The effects of fluid characteristics,shaping pressure,axial feed rate,and friction coefficient on the wall thicknesses of torsions beams during formation were investigated.Fluid movement speed was related to tube deformation.Shaping pressure had a significant effect on rounded corners and straight edges.The axial feed speed was increased,and the uneven distribution of wall thicknesses was effectively improved.Although the friction coefficient had a nonsignificant effect on the wall thickness of the ladder-shaped region,it had a significant influence on a large deformation of wall thickness in the V-shaped area.In this paper,a method of fluid-solid coupling simulation analysis and sliding die is proposed to study the high pressure forming law in torsion beam.展开更多
Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved ...Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.展开更多
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.展开更多
Abstract: A joint solution model of variabk:-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated...Abstract: A joint solution model of variabk:-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method (FEM). The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/(m2. K). When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the: decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.展开更多
The seismic response characteristics of underground structures in saturated soils are investigated.A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to sim...The seismic response characteristics of underground structures in saturated soils are investigated.A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to simulate the dynamic behavior of saturated soils.The accuracy of the model is validated using a classic example in literature.The performance of the model is verified by its application on simulating the seismic response characteristics of a subway station built in saturated soils.The merits of the model are demonstrated by comparing the difference of the seismic response of an underground structure in saturated soils between using the fully coupling model and a single-phase medium model.The study finds that the fully coupling model developed herein can simulate the dynamic response characteristics of the underground structures in saturated soils with high accuracy.The seismic response of the underground structure tends to be underestimated by using the single-phase medium model compared with using the fully coupling model,which provides a weaker confining action to the underground structure.展开更多
The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contami...The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.展开更多
As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments (HBS) is almost inevitable. The variation in temperature and p...As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments (HBS) is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.展开更多
About 75% water-inrush accidents in China are caused by geological structure such as faults, therefore, it is necessary to investigate the water-inrush mechanism of faults to provide references for the mining activity...About 75% water-inrush accidents in China are caused by geological structure such as faults, therefore, it is necessary to investigate the water-inrush mechanism of faults to provide references for the mining activity above confined water. In this paper, based on the fluid-solid coupling theory, we built the stress-seepage coupling model for rock, then we combined with an example of water-inrush caused by fault, studied the water-inrush mechanism by using the numerical software COMSOL Mutiphysics, analyzed the change rule of shear stress, vertical stress, plastic area and water pressure for stope with a fault, and estimated the water-inrush risk at the different distances between working faces and the fault. The numerical simula- tion results indicate that: (1) the water-inrush risk will grow as the decrease of the distance between working face and the fault; (2) the failure mode of the rock in floor with fault is shear failure; (3) the rock between water-containing fault and working face failure is the reason for water-inrush.展开更多
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.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Latti...Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.展开更多
In the paper, 3-D analysis method with unitive schemes is set up, which is used to resolve the uplift with multiple moving boundaries and multiple nonlinear coupling for anchored liquid storage tanks. hi it, an algori...In the paper, 3-D analysis method with unitive schemes is set up, which is used to resolve the uplift with multiple moving boundaries and multiple nonlinear coupling for anchored liquid storage tanks. hi it, an algorithm of quasi-harmonious finite elements for arbitrary quadrilateral of thin plates and shells is built up to analyze the multiple coupling problems of general thin plates and shells structures with three dimensions, the complementary equations for analyzing uplifting moving boundary problems are deduced. The axial symmetry and presumption of beam type mode are not used. In it, an algorithm is put forward for analyzing the Navier-Stokes problems of unsteady, three-dimensional, and viscous liquid with sloshing of moving boundary surfaces in large amplitude under ALE frame by scheme of time-split-steps to which linear potential theory is not applied. The algorithms can be used to analyze the solid-liquid multiple nonlinear coupling problems with 3-D moving boundary with friction in multiple places.展开更多
The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift...The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute twophase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersedphase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.展开更多
A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two metho...A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.展开更多
The Helmholtz equation is sometimes supplemented by conditions that include the specification of the boundary value of the divergence of the unknown.In this paper, we study the vector Helmholtz problem in domains of b...The Helmholtz equation is sometimes supplemented by conditions that include the specification of the boundary value of the divergence of the unknown.In this paper, we study the vector Helmholtz problem in domains of both C^(1,1)and Lipschitz.We es- tablish a rigorous variational analysis such as equivalence,existence and uniqueness. And we propose finite element approximations based on the uncoupled solutions.Fi- nally we present a convergence analysis and error estimates.展开更多
This work deals with the analysis of soil-structure interaction modeling of pipeline problems in static behavior using the coupling between FEM (finite element method) and BEM (boundary element method). The repres...This work deals with the analysis of soil-structure interaction modeling of pipeline problems in static behavior using the coupling between FEM (finite element method) and BEM (boundary element method). The representation of the pipe is made by MEF using one fmite element in the cylindrical panel formulated from the theory of equivalent discrete layers (Layerwise theory), proposed by J. N. Reddy. The soil is represented by elastic continum infimite or semi-infinite and modeled using boundary elements with special curved surface, associated with cylindrical panel used to represent the soil-structure interaction within the soil, especially at the contact surface with the pipe.展开更多
A new calculation formula of THM coupling stress intensity factor was derived by the boundary collocation method, in which an additional constant stress function was successfully introduced for the cracked specimen wi...A new calculation formula of THM coupling stress intensity factor was derived by the boundary collocation method, in which an additional constant stress function was successfully introduced for the cracked specimen with hydraulic pressure applied on its crack surface. Based on the newly derived formula, THM coupling fracture modes (including tensile, shear and mixed fracture mode) can be predicted by a new fracture criterion of stress intensity factor ratio, where the maximum axial load was measured by self-designed THM coupling fracture test. SEM analyses of THM coupling fractured surface indicate that the higher the temperature and hydraulic pressure are and the lower the confining pressure is, the more easily the intergranular (tension) fracture occurs. The transgranular (shear) fracture occurs in the opposite case while the mixed-mode fracture occurs in the middle case. The tested THM coupling fracture mechanisms are in good agreement with the predicted THM coupling fracture modes, which can verify correction of the newly-derived THM coupling stress intensity factor formula.展开更多
In the context of global warming and intensified human activities,glacier instability in plateau regions has increased,and glacier debris flows have become active,which poses a significant threat to the lives and prop...In the context of global warming and intensified human activities,glacier instability in plateau regions has increased,and glacier debris flows have become active,which poses a significant threat to the lives and property of people and socioeconomic development.The mass movement process of glacier debris flows is extremely complex,so this paper uses the 2018 Sedongpu glacier debris flow event on the Qinghai-Tibet Plateau as an example and applies a numerical simulation method to invert the whole process of mass movement.In view of the interaction between phases in the process of motion,we use the fluid-solid coupling method to describe the mass movement.The granular-flow model and drift-flux model are employed in FLOW3D software to study the mass movement process of glacier debris flows and explore their dynamic characteristics.The results indicate that the glacier debris flow lasted for 700 s,and the movement process was roughly divided into four stages,including initiation,scraping,surging and deposition;the depositional characteristics calculated by the fluid-solid coupling model are consistent with the actual survey results and have good reliability;strong erosion occurs during the mass movement,the clear volume amplification effect,and the first wave climbs 17.8 m across the slope.The fluid-solid coupling method can better simulate glacier debris flows in plateau regions,which is helpful for the study of the mechanism and dynamic characteristics of such disasters.展开更多
Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two...Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.展开更多
基金supported by the Natural Science Foundation of Hebei Province,China (Grant No.A2021502004)the Fundamental Research Funds for the Central Universities (Grant No.2024MS126).
文摘This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.
基金Natural Science Foundation of China(Grant No.81960332)Guangxi Provincial Innovation driven Development Project(Grant No.GKAA17204062)+1 种基金Guangxi Provincial Natural Science Foundation(Grant No.2016GXNSFAA380211)Liuzhou Municipal Scientific Research and Technology Development Plan(Grant No.2016C050203)。
文摘Hydroformed parts are widely used in industrial automotive parts because of their higher stiffness and fatigue strength and reduced weight relative to their corresponding cast and welded parts.This paper reports a hydraulicforming experimental platform for rectangular tube fittings that was constructed to conduct an experiment on the hydraulic forming of rectangular tube fittings.A finite element model was established on the basis of the fluid–solid coupling method and simulation analysis.The correctness of the simulation analysis and the feasibility of the fluid–solid coupling method for hydraulic forming simulation analysis were verified by comparing the experimental results with the simulation results.On the basis of the simulation analysis of the hydraulic process of the torsion beam using the fluid–solid coupling method,a sliding mold suitable for the hydroforming of torsion beams was designed for its structural characteristics.The effects of fluid characteristics,shaping pressure,axial feed rate,and friction coefficient on the wall thicknesses of torsions beams during formation were investigated.Fluid movement speed was related to tube deformation.Shaping pressure had a significant effect on rounded corners and straight edges.The axial feed speed was increased,and the uneven distribution of wall thicknesses was effectively improved.Although the friction coefficient had a nonsignificant effect on the wall thickness of the ladder-shaped region,it had a significant influence on a large deformation of wall thickness in the V-shaped area.In this paper,a method of fluid-solid coupling simulation analysis and sliding die is proposed to study the high pressure forming law in torsion beam.
基金support of this work by the National Natural Science Foundation of China(No.51405096)the Fundamental Research Funds for the Central Universities(HEUCF210710).
文摘Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.
基金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.
基金Project(20080431380) supported by China Postdoctoral Science Foundation
文摘Abstract: A joint solution model of variabk:-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method (FEM). The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/(m2. K). When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the: decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.
基金National Natural Science Foundation of People’s Republic of China under Grant Nos.51178011 and 51778386the Key Fundamental Study Development Project of People’s Republic of China under Grant No.2011CB013602。
文摘The seismic response characteristics of underground structures in saturated soils are investigated.A fully fluid-solid coupling dynamic model is developed and implemented into ABAQUS with a user-defined element to simulate the dynamic behavior of saturated soils.The accuracy of the model is validated using a classic example in literature.The performance of the model is verified by its application on simulating the seismic response characteristics of a subway station built in saturated soils.The merits of the model are demonstrated by comparing the difference of the seismic response of an underground structure in saturated soils between using the fully coupling model and a single-phase medium model.The study finds that the fully coupling model developed herein can simulate the dynamic response characteristics of the underground structures in saturated soils with high accuracy.The seismic response of the underground structure tends to be underestimated by using the single-phase medium model compared with using the fully coupling model,which provides a weaker confining action to the underground structure.
文摘The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.
基金supported by the Major National Science and Technology Program(Nos.2008ZX05026-00411 and 2011ZX05026-004-08)the Program for Changjiang Scholars and Innovative Research Team in University(No.RT1086)
文摘As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments (HBS) is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.
基金Supported by the National Basic Research Program of China (2010CB226800) the National Natural Science Foundation of China (50904065) the Program for New Century Excellent Talents in University (NCET-09-0728)
文摘About 75% water-inrush accidents in China are caused by geological structure such as faults, therefore, it is necessary to investigate the water-inrush mechanism of faults to provide references for the mining activity above confined water. In this paper, based on the fluid-solid coupling theory, we built the stress-seepage coupling model for rock, then we combined with an example of water-inrush caused by fault, studied the water-inrush mechanism by using the numerical software COMSOL Mutiphysics, analyzed the change rule of shear stress, vertical stress, plastic area and water pressure for stope with a fault, and estimated the water-inrush risk at the different distances between working faces and the fault. The numerical simula- tion results indicate that: (1) the water-inrush risk will grow as the decrease of the distance between working face and the fault; (2) the failure mode of the rock in floor with fault is shear failure; (3) the rock between water-containing fault and working face failure is the reason for water-inrush.
基金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.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
基金WJD,JYZ,CLC,ZX,and ZGY were supported by the National Natural Science Foundation of China(Grant Number 51705143)the Education Department of Hunan Province(Grant Number 22B0464)the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant Number QL20230249).
文摘Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.
文摘In the paper, 3-D analysis method with unitive schemes is set up, which is used to resolve the uplift with multiple moving boundaries and multiple nonlinear coupling for anchored liquid storage tanks. hi it, an algorithm of quasi-harmonious finite elements for arbitrary quadrilateral of thin plates and shells is built up to analyze the multiple coupling problems of general thin plates and shells structures with three dimensions, the complementary equations for analyzing uplifting moving boundary problems are deduced. The axial symmetry and presumption of beam type mode are not used. In it, an algorithm is put forward for analyzing the Navier-Stokes problems of unsteady, three-dimensional, and viscous liquid with sloshing of moving boundary surfaces in large amplitude under ALE frame by scheme of time-split-steps to which linear potential theory is not applied. The algorithms can be used to analyze the solid-liquid multiple nonlinear coupling problems with 3-D moving boundary with friction in multiple places.
基金The project supported by the National Natural Science Foundation of China(90205024)Russian Foundation for Basic Research(RFBR and(RFBR-NSFC-39004)
文摘The present paper describes a numerical two-way coupling model for shock-induced laminar boundary-layer flows of a dust-laden gas and studies the transverse migration of fine particles under the action of Saffman lift force. The governing equations are formulated in the dilute twophase continuum framework with consideration of the finiteness of the particle Reynolds and Knudsen numbers. The full Lagrangian method is explored for calculating the dispersedphase flow fields (including the number density of particles) in the regions of intersecting particle trajectories. The computation results show a significant reaction of the particles on the two-phase boundary-layer structure when the mass loading ratio of particles takes finite values.
基金Project supported by China Postdoctoral Science Foundation (No.2004036145)
文摘A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.
文摘The Helmholtz equation is sometimes supplemented by conditions that include the specification of the boundary value of the divergence of the unknown.In this paper, we study the vector Helmholtz problem in domains of both C^(1,1)and Lipschitz.We es- tablish a rigorous variational analysis such as equivalence,existence and uniqueness. And we propose finite element approximations based on the uncoupled solutions.Fi- nally we present a convergence analysis and error estimates.
文摘This work deals with the analysis of soil-structure interaction modeling of pipeline problems in static behavior using the coupling between FEM (finite element method) and BEM (boundary element method). The representation of the pipe is made by MEF using one fmite element in the cylindrical panel formulated from the theory of equivalent discrete layers (Layerwise theory), proposed by J. N. Reddy. The soil is represented by elastic continum infimite or semi-infinite and modeled using boundary elements with special curved surface, associated with cylindrical panel used to represent the soil-structure interaction within the soil, especially at the contact surface with the pipe.
基金Project(11072269)supported by the National Natural Science Foundation of ChinaProject(20090162110066)supported by the Research Fund for the Doctoral Program of Higher Education of China
文摘A new calculation formula of THM coupling stress intensity factor was derived by the boundary collocation method, in which an additional constant stress function was successfully introduced for the cracked specimen with hydraulic pressure applied on its crack surface. Based on the newly derived formula, THM coupling fracture modes (including tensile, shear and mixed fracture mode) can be predicted by a new fracture criterion of stress intensity factor ratio, where the maximum axial load was measured by self-designed THM coupling fracture test. SEM analyses of THM coupling fractured surface indicate that the higher the temperature and hydraulic pressure are and the lower the confining pressure is, the more easily the intergranular (tension) fracture occurs. The transgranular (shear) fracture occurs in the opposite case while the mixed-mode fracture occurs in the middle case. The tested THM coupling fracture mechanisms are in good agreement with the predicted THM coupling fracture modes, which can verify correction of the newly-derived THM coupling stress intensity factor formula.
基金supported by the National Natural Science Foundation of China(Nos.U20A20111,41977229)the Sichuan Youth Science and Technology Innovation Research Team Project(No.2020JDTD0006).
文摘In the context of global warming and intensified human activities,glacier instability in plateau regions has increased,and glacier debris flows have become active,which poses a significant threat to the lives and property of people and socioeconomic development.The mass movement process of glacier debris flows is extremely complex,so this paper uses the 2018 Sedongpu glacier debris flow event on the Qinghai-Tibet Plateau as an example and applies a numerical simulation method to invert the whole process of mass movement.In view of the interaction between phases in the process of motion,we use the fluid-solid coupling method to describe the mass movement.The granular-flow model and drift-flux model are employed in FLOW3D software to study the mass movement process of glacier debris flows and explore their dynamic characteristics.The results indicate that the glacier debris flow lasted for 700 s,and the movement process was roughly divided into four stages,including initiation,scraping,surging and deposition;the depositional characteristics calculated by the fluid-solid coupling model are consistent with the actual survey results and have good reliability;strong erosion occurs during the mass movement,the clear volume amplification effect,and the first wave climbs 17.8 m across the slope.The fluid-solid coupling method can better simulate glacier debris flows in plateau regions,which is helpful for the study of the mechanism and dynamic characteristics of such disasters.
基金the Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period(Grant No.2015BAK17B01)Science Foundation of Institute of Engineering Mechanics,CEA under Grant No.2014B10+1 种基金Natural Science Foundation of Heilongjiang Province of China under Grant No.LC201403National Natural Science Foundation under Grant No.51378479 and No.51108431
文摘Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.