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.展开更多
The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges tha...The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges that affect its accuracy and efficiency.Firstly,the implicit handling of fluid-solid boundary conditions through parameter averaging strategy often results in low simulation accuracy.Secondly,surface topography can introduce staircase diffraction noise when grid spacing is large.To address these issues,this paper presents a novel approach.We derive an implicit expression for fluidsolid boundary conditions based on average medium theory,translating explicit boundary conditions into model parameter modification.This enables implicit handling of fluid-solid boundaries by modifying the parameters near the boundary.Furthermore,to mitigate staircase diffraction noise,we employ multiple interface discretization based on the superposition method.This effectively suppresses staircase diffraction noise without requiring grid refinement.The efficacy of our method in accurately modeling wave propagation phenomena in fluid-solid coupled media is demonstrated by numerical examples.Results align well with those obtained using the spectral element method(SEM),with significant reduction in staircase diffraction noise.展开更多
Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie conditio...Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.展开更多
In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin proj...In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin projectile.The platform's reliability is validated by reproducing the characteristic resonance instability of such projectiles.By coupling the solution of the Unsteady Reynolds-Averaged Navier-Stokes equations and the seven-degree-of-freedom RBD equations,the virtual flight simulations of fixed canard dual-spin projectiles at various curvature trajectories are achieved,and the dynamic mechanism of the trajectory following process is analyzed.The instability mechanism of the dynamic instability during trajectory following process of the fixed canard dual-spin projectile is elucidated by simulating the rolling/coning coupled forced motion,and subsequently validated through virtual flight simulations.The findings suggest that an appropriate yaw moment can drive the projectile axis to precession in the tangential direction of the trajectory,thereby enhancing the trajectory following stability.However,the damping of the projectile attains its minimum value when the forward body equilibrium rotational speed(-128 rad/s)is equal to the negative of the fast mode frequency of the projectile.Insufficient damping leads to the fixed canard dual-spin projectile exiting the dynamic stability domain during the trajectory following,resulting in weakly damped instability.Keeping the forward body not rotating or increasing the spin rates to-192 rad/s can enhance the projectile's damping,thereby improving its dynamic stability.展开更多
An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal str...An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.展开更多
Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression an...Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.展开更多
To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mat...To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.展开更多
This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating som...This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.展开更多
A coupled supersonic inlet-fan Navier–Stokes simulation method was developed by using COMSOL-CFD code. The flow turning, pressure rise and loss effects across blade rows of the fan and the inlet-fan interactions were...A coupled supersonic inlet-fan Navier–Stokes simulation method was developed by using COMSOL-CFD code. The flow turning, pressure rise and loss effects across blade rows of the fan and the inlet-fan interactions were taken into account as source terms of the governing equations without a blade geometry by a body force model. In this model, viscous effects in blade passages can also be calculated directly, which include the exchange of momentum between fluids and detailed viscous flow close to walls. NASA Rotor 37 compressor test rig was used to validate the ability of the body force model to estimate the real performance of blade rows. Calculated pressure ratio characteristics and the distribution of the total pressure, total temperature, and swirl angle in the span direction agreed well with experimental and numerical data. It is shown that the body force model is a promising approach for predicting the flow field of the turbomachinery. Then, coupled axisymmetric mixed compression supersonic inlet-fan simulations were conducted at Mach number 2.8 operating conditions. The analysis includes coupled steady-state performance, and effects of the fan on the inlet. The results indicate that the coupled simulation method is capable of simulating behavior of the supersonic inlet-fan system.展开更多
The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the bea...The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.展开更多
The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the con...The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.展开更多
The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction di...The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.展开更多
This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton）...This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.展开更多
Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynam...Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.展开更多
The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the ...The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.展开更多
A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite...A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.展开更多
In the present paper, we established a traveling wave solution by using modified Kudryashov method for the space-time fractional nonlinear partial differential equations. The method is used to obtain the exact solutio...In the present paper, we established a traveling wave solution by using modified Kudryashov method for the space-time fractional nonlinear partial differential equations. The method is used to obtain the exact solutions for different types of the space-time fractional nonlinear partial differential equations such as, the space-time fractional coupled equal width wave equation(CEWE) and the space-time fractional coupled modified equal width wave equation(CMEW), which are the important soliton equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform and properties of modified Riemann–Liouville derivative. We plot the exact solutions for these equations at different time levels.展开更多
The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a ...The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.展开更多
An improved analytical method to determine the content of 52 major, minor and trace elements in marine geological samples, using a HF-HCl-HNO_3 acid system with a high-pressure closed digestion method（HPCD）, is stud...An improved analytical method to determine the content of 52 major, minor and trace elements in marine geological samples, using a HF-HCl-HNO_3 acid system with a high-pressure closed digestion method（HPCD）, is studied by an inductively coupled plasma optical emission spectrometry（ICP-OES） and an inductively coupled plasma mass spectrometry（ICP-MS）. The operating parameters of the instruments are optimized, and the optimal analytical parameters are determined. The influences of optical spectrum and mass spectrum interferences, digestion methods and acid systems on the analytical results are investigated. The optimal spectral lines and isotopes are chosen, and internal standard element of rhodium is selected to compensate for matrix effects and analytical signals drifting. Compared with the methods of an electric heating plate digestion and a microwave digestion, a high-pressure closed digestion method is optimized with less acid, complete digestion,less damage for digestion process. The marine geological samples are dissolved completely by a HF-HCl-HNO_3 system, the relative error（RE） for the analytical results are all less than 6.0%. The method detection limits are 2–40μg/g by the ICP-OES, and 6–80 ng/g by ICP-MS. The methods are used to determine the marine sediment reference materials（GBW07309, GBW07311, GBW07313）, rock reference materials（GBW07103, GBW07104,GBW07105）, and cobalt-rich crust reference materials（GBW07337, GBW07338, GBW07339）, the obtained analytical results are in agreement with the certified values, and both of the relative standard deviation（RSD） and the relative error（RE） are less than 6.0%. The analytical method meets the requirements for determining 52 elements contents of bulk marine geological samples.展开更多
In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical e...In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.展开更多
基金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(Nos.42404134,U24B2031,42174160)the China Postdoctoral Science Foundation(No.2024M753204)the National Key R&D Program of China(Nos.2021YFA0716901,2022YFB3904601)。
文摘The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges that affect its accuracy and efficiency.Firstly,the implicit handling of fluid-solid boundary conditions through parameter averaging strategy often results in low simulation accuracy.Secondly,surface topography can introduce staircase diffraction noise when grid spacing is large.To address these issues,this paper presents a novel approach.We derive an implicit expression for fluidsolid boundary conditions based on average medium theory,translating explicit boundary conditions into model parameter modification.This enables implicit handling of fluid-solid boundaries by modifying the parameters near the boundary.Furthermore,to mitigate staircase diffraction noise,we employ multiple interface discretization based on the superposition method.This effectively suppresses staircase diffraction noise without requiring grid refinement.The efficacy of our method in accurately modeling wave propagation phenomena in fluid-solid coupled media is demonstrated by numerical examples.Results align well with those obtained using the spectral element method(SEM),with significant reduction in staircase diffraction noise.
基金a U.S. Federal Railroad Administration (FRA)BAA project,titled “Mitigation of Differential Movement at Railway Transitions for High-Speed Passenger Rail and Joint Passenger/Freight Corridors”the financial support provided by the China Scholarship Council (CSC),which funded Zhongyi Liu’s and Wenjing Li’s time and research efforts for this study
文摘Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.
基金supported by the National Natural Science Foundation of China(Grant Nos.U2141254 and U23B6009)。
文摘In this paper,a high-fidelity computational fluid dynamics(CFD)and rigid body dynamics(RBD)coupled platform for virtual flight simulation is developed to investigate the flight stability of fixed canard dual-spin projectile.The platform's reliability is validated by reproducing the characteristic resonance instability of such projectiles.By coupling the solution of the Unsteady Reynolds-Averaged Navier-Stokes equations and the seven-degree-of-freedom RBD equations,the virtual flight simulations of fixed canard dual-spin projectiles at various curvature trajectories are achieved,and the dynamic mechanism of the trajectory following process is analyzed.The instability mechanism of the dynamic instability during trajectory following process of the fixed canard dual-spin projectile is elucidated by simulating the rolling/coning coupled forced motion,and subsequently validated through virtual flight simulations.The findings suggest that an appropriate yaw moment can drive the projectile axis to precession in the tangential direction of the trajectory,thereby enhancing the trajectory following stability.However,the damping of the projectile attains its minimum value when the forward body equilibrium rotational speed(-128 rad/s)is equal to the negative of the fast mode frequency of the projectile.Insufficient damping leads to the fixed canard dual-spin projectile exiting the dynamic stability domain during the trajectory following,resulting in weakly damped instability.Keeping the forward body not rotating or increasing the spin rates to-192 rad/s can enhance the projectile's damping,thereby improving its dynamic stability.
基金the National Metal and Materials Technology Centerthe Thailand Research Fund+1 种基金the Office of Higher Education Commissionthe Chulalongkorn University for supporting the present research
文摘An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.
基金financially supported by the National Natural Science Foundation of China(Nos.42272153 and 42472195)the Research Fund of PetroChina Tarim Oilfield Company(No.671023060003)the Research Fund of China National Petroleum Corporation Limited(No.2023ZZ16YJ04).
文摘Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.
基金Project(2006AA03Z523) supported by the National High-Tech Research and Development Program of ChinaProject(08C26224302178) supported by the Innovation Foundation of Central South University,China
文摘To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.
基金supported by National Basic Research Program of China (Grant No 2006CB705500)Chang-Jiang Scholars and Innovative Research Team in University of China (Grant No IRT0605)the National Natural Science Foundation of China (Grant No70631001)
文摘This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.
基金support of National Natural Science Foundation of China (Nos. 51706008 and 51636001)China Postdoctoral Science Foundation (No. 2017M610742)Aeronautics Power Foundation of China (No. 6141B090315)
文摘A coupled supersonic inlet-fan Navier–Stokes simulation method was developed by using COMSOL-CFD code. The flow turning, pressure rise and loss effects across blade rows of the fan and the inlet-fan interactions were taken into account as source terms of the governing equations without a blade geometry by a body force model. In this model, viscous effects in blade passages can also be calculated directly, which include the exchange of momentum between fluids and detailed viscous flow close to walls. NASA Rotor 37 compressor test rig was used to validate the ability of the body force model to estimate the real performance of blade rows. Calculated pressure ratio characteristics and the distribution of the total pressure, total temperature, and swirl angle in the span direction agreed well with experimental and numerical data. It is shown that the body force model is a promising approach for predicting the flow field of the turbomachinery. Then, coupled axisymmetric mixed compression supersonic inlet-fan simulations were conducted at Mach number 2.8 operating conditions. The analysis includes coupled steady-state performance, and effects of the fan on the inlet. The results indicate that the coupled simulation method is capable of simulating behavior of the supersonic inlet-fan system.
基金supported by the National Natural Science Foundation of China(No.11605234)the Foundation of ASIPP(No.DSJJ-15-GC02)
文摘The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China(Grants 20113219110025,20133219110037)the National Natural Science Foundation of China(Grants 11102089,61304137)the Program for New Century Excellent Talents in University(NCET-10-0075)
文摘The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.
基金supported by the National Natural Science Foundation of China(Grant No.91130013)the Open Foundation of State Key Laboratory of HighPerformance Computing of China
文摘The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10871117 and 10571110)
文摘This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.
基金supported by the National Natural Science Foundation of China(51879159,51490675,11432009,and 51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for doctoral student
文摘Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.
基金supported by the National Natural Science Foundation of China(No.11472058)
文摘The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.
基金Supported by the National Science and Technology Major Project (2017ZX05013-005)。
文摘A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.
文摘In the present paper, we established a traveling wave solution by using modified Kudryashov method for the space-time fractional nonlinear partial differential equations. The method is used to obtain the exact solutions for different types of the space-time fractional nonlinear partial differential equations such as, the space-time fractional coupled equal width wave equation(CEWE) and the space-time fractional coupled modified equal width wave equation(CMEW), which are the important soliton equations. Both equations are reduced to ordinary differential equations by the use of fractional complex transform and properties of modified Riemann–Liouville derivative. We plot the exact solutions for these equations at different time levels.
基金supported by the National Basic Research Program of China (Grant 2013CB733004)
文摘The dynamic characteristics of a beam-cable coupled system are investigated using an improved Chebyshev spectral element method in order to observe the effects of adding cables on the beam. The system is modeled as a double Timoshenko beam system interconnected by discrete springs. Utilizing Chebyshev series expansion and meshing the system according to the locations of its connections, numerical results of the natural frequencies and mode shapes are obtained using only a few elements, and the results are validated by comparing them with the results of a finite-element method. Then the effects of the cable parameters and layout of connections on the natural frequencies and mode shapes of a fixed-pinned beam are studied. The results show that the modes of a beam-cable coupled system can be classified into two types, beam mode and cable mode, according to the dominant deformation. To avoid undesirable vibrations of the cable, its parameters should be controlled in a reasonable range, or the layout of the connections should be optimized.
基金The China Ocean Mineral Resources Research and Development Association Research Program of the State Oceanic Administration of China under contract No.DY125-13-R-07the National Natural Science Foundation of China under contract Nos 41322036 and 41230960+1 种基金the Shandong Provincial Natural Science Foundation of China under contract No.ZR2014DP009the Special Basic Research Funds for Central Public Research Institutes for The First Institute of Oceanography,State Oceanic Administration of China under contract Nos GY0213G06 and GY02-2012G35
文摘An improved analytical method to determine the content of 52 major, minor and trace elements in marine geological samples, using a HF-HCl-HNO_3 acid system with a high-pressure closed digestion method（HPCD）, is studied by an inductively coupled plasma optical emission spectrometry（ICP-OES） and an inductively coupled plasma mass spectrometry（ICP-MS）. The operating parameters of the instruments are optimized, and the optimal analytical parameters are determined. The influences of optical spectrum and mass spectrum interferences, digestion methods and acid systems on the analytical results are investigated. The optimal spectral lines and isotopes are chosen, and internal standard element of rhodium is selected to compensate for matrix effects and analytical signals drifting. Compared with the methods of an electric heating plate digestion and a microwave digestion, a high-pressure closed digestion method is optimized with less acid, complete digestion,less damage for digestion process. The marine geological samples are dissolved completely by a HF-HCl-HNO_3 system, the relative error（RE） for the analytical results are all less than 6.0%. The method detection limits are 2–40μg/g by the ICP-OES, and 6–80 ng/g by ICP-MS. The methods are used to determine the marine sediment reference materials（GBW07309, GBW07311, GBW07313）, rock reference materials（GBW07103, GBW07104,GBW07105）, and cobalt-rich crust reference materials（GBW07337, GBW07338, GBW07339）, the obtained analytical results are in agreement with the certified values, and both of the relative standard deviation（RSD） and the relative error（RE） are less than 6.0%. The analytical method meets the requirements for determining 52 elements contents of bulk marine geological samples.
基金supported by the National Natural Science Foundation of China(Grant No.11171038)
文摘In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.
