The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on mo...The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.展开更多
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.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution t...A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.展开更多
We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive ...We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.展开更多
Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulati...Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.展开更多
The Boussinesq equations,pivotal in the analysis of water wave dynamics,effectively model weakly nonlinear and long wave approximations.This study utilizes the complete discriminant system within a polynomial approach...The Boussinesq equations,pivotal in the analysis of water wave dynamics,effectively model weakly nonlinear and long wave approximations.This study utilizes the complete discriminant system within a polynomial approach to derive exact traveling wave solutions for the coupled Boussinesq equation.The solutions are articulated through soliton,trigonometric,rational,and Jacobi elliptic functions.Notably,the introduction of Jacobi elliptic function solutions for this model marks a pioneering advancement.Contour plots of the solutions obtained by assigning values to various parameters are generated and subsequently analyzed.The methodology proposed in this study offers a systematic means to tackle nonlinear partial differential equations in mathematical physics,thereby enhancing comprehension of the physical attributes and dynamics of water waves.展开更多
This paper introduces a bond-based peridynamics(BB-PD)algorithm for crack identification,integrating the Delaunay triangulation method to accurately identify the structural characteristics of threedimensional(3D)crack...This paper introduces a bond-based peridynamics(BB-PD)algorithm for crack identification,integrating the Delaunay triangulation method to accurately identify the structural characteristics of threedimensional(3D)cracks in rocks.A bond-based crack quantification standard is proposed to analyze the evolution of cracks of various sizes.A multi-attribute peridynamic model,developed using a multilayer algorithm,was employed to simulate the fracturing process of sandstone disks and semi-disks under varying temperatures,with the model calibrated and validated against experimental results.The simulation results show that temperature induces nonlinear degradation in the tensile strength and fracture toughness of sandstone,with 500℃ identified as the threshold temperature.Thermal cracks exhibit varying degrees of influence on Mode I cracks across different temperature ranges.Thermal damage significantly promotes the initiation and propagation of Mode I cracks in sandstone,thereby reducing its tensile strength and fracture toughness.Under applied loads,crack propagation in sandstone predominantly occurs during the failure stage,characterized by the rapid growth of longer cracks and a slow increase or reduction in shorter cracks.展开更多
Arsenic speciation in freshwater fish is crucial for providing meaningful consumption guidelines that allow the public to make informed decisions regarding its consumption.While marine fish have attractedmuch research...Arsenic speciation in freshwater fish is crucial for providing meaningful consumption guidelines that allow the public to make informed decisions regarding its consumption.While marine fish have attractedmuch research interest due to their higher arsenic content,research on freshwater fish is limited due to the challenges in quantifying and identifying arsenic species present at trace levels.We describe here a sensitivemethod and its application to the quantification of arsenic species in freshwater fish.Arsenic species from fish tissues were extracted using a methanol/water mixture(1:1 vol.ratio)and ultrasound sonication.Anion-exchange high-performance liquid chromatography(HPLC)enabled separation of arsenobetaine(AsB),inorganic arsenite(iAs^(Ⅲ)),dimethylarsinic acid(DMA),monomethylarsonic acid(MMA),inorganic arsenate(iAs^(Ⅴ)),and three new arsenic species.Inductively coupled plasma mass spectrometry(ICPMS)provided highly sensitive and specific detection of arsenic.A limit of detection of 0.25μg/kg(wet weight fish tissue)was achieved for the five target arsenic species:AsB,iAs^(Ⅲ),DMA,MMA,and iAs^(Ⅴ).A series of experimentswere conducted to ensure the accuracy and validity of the analytical method.The method was successfully applied to the determination of arsenic species in lakewhitefish,northern pike,and walleye,with AsB,DMA,and iAs^(Ⅴ) being frequently detected.Three new arsenic species were detected,but their chromatographic retention times did not match with those of any available arsenic standards.Future research is necessary to elucidate the identity of these new arsenic species detected in freshwater fish.展开更多
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.展开更多
In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SP...In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.展开更多
Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its co...Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.展开更多
Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading t...Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.展开更多
A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction pro...A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.展开更多
The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled ...The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.展开更多
42CrMo steel has the characteristics of high strength,high wear resistance,high impact resistance,and fatigue resistance.Therefore,drilling 42CrMo steel has always been a challenging task.The indexable drill bit has t...42CrMo steel has the characteristics of high strength,high wear resistance,high impact resistance,and fatigue resistance.Therefore,drilling 42CrMo steel has always been a challenging task.The indexable drill bit has the advantages of high processing efficiency and low processing cost and has been widely used in the field of aerospace hole processing.To better understand the machining mechanism of the indexable drill bit,this paper uses the Coupled EulerianLagrangian method(CEL)to simulate the three-dimensional drilling model for the first time.The simulation results of the drilling force obtained by the CEL method and Lagrangian method are compared with the experimental results.It is verified that the CEL method is easy to converge and can avoid the problem of program interruption caused by mesh distortion,and the CEL simulation value is more consistent with the actual value.Secondly,the simulation results of cutting force and blade cutting edge node temperature under different process parameters are extracted.The variation of time domain cutting force,frequency domain cutting force and tool temperature with process parameters are obtained.This study provides a new method for the prediction of cutting performance and the optimization of process parameters of indexable drills.展开更多
In this paper,two types of fractional nonlinear equations in Caputo sense,time-fractional Newell–Whitehead equation(FNWE)and time-fractional generalized Hirota–Satsuma coupled KdV system(HS-cKdVS),are investigated b...In this paper,two types of fractional nonlinear equations in Caputo sense,time-fractional Newell–Whitehead equation(FNWE)and time-fractional generalized Hirota–Satsuma coupled KdV system(HS-cKdVS),are investigated by means of the q-homotopy analysis method(q-HAM).The approximate solutions of the proposed equations are constructed in the form of a convergent series and are compared with the corresponding exact solutions.Due to the presence of the auxiliary parameter h in this method,just a few terms of the series solution are required in order to obtain better approximation.For the sake of visualization,the numerical results obtained in this paper are graphically displayed with the help of Maple.展开更多
The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when th...The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.展开更多
In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are c...In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are considered in terms of radial and non-radial basis functions.The time-dependent variable in the first radial scheme is directly considered as the normal space variables to formulate an"isotropic"space-time radial basis function.The second non-radial scheme considered relationship between time-dependent and spacedependent variables.Under such circumstance,we can get a one-step space-time meshless scheme.The numerical findings demonstrate that the proposed meshless schemes are precise,user-friendly,and effective in solving the coupled Burgers'equations.展开更多
This paper presents a new approach to synthesize admittance function polynomials and coupling matrices for coupled resonator filters. The N + 2 transversal network method is applied to study a coupled resonator f...This paper presents a new approach to synthesize admittance function polynomials and coupling matrices for coupled resonator filters. The N + 2 transversal network method is applied to study a coupled resonator filter. This method allowed us to determine the polynomials of the reflection and transmission coefficients. A study is made for a 4 poles filter with 2 transmission zeros between the N + 2 transversal network method and the one found in the literature. A MATLAB code was designed for the numerical simulation of these coefficients for the 6, 8, and 10 pole filter with 4 transmission zeros.展开更多
基金support of the National Natural Science Foundation of China(No.51675406)the Basic Research Project Group,China(No.514010106-205)。
文摘The implicit partition algorithm used to solve fluid–structure coupling problems has high accuracy,but it requires a long computation time.In this paper,a semi-implicit fluid–structure coupling algorithm based on modal force prediction-correction is proposed to improve the computational efficiency.In the pre-processing stage,the fluid domain is assumed to be a pseudo-elastic solid and merged with the solid domain to form a holistic system,and the normalized modal information of the holistic system is calculated and stored.During the sub-step cycle,the modal superposition method is used to obtain the response of the holistic system with the predicted modal force as the load,so that the deformation of the structure and the updating of the fluid mesh can be achieved simultaneously.After solving the Reynolds-averaged Navier-Stokes equations in the fluid domain,the predicted modal force is corrected and a new sub-step cycle is started until the converged result is obtained.In this method,the computation of the fluid equations and the updating of the dynamic mesh are done implicitly,while the deformation of the structure is done explicitly.Two numerical cases,vortex induced oscillation of an elastic beam and fluid–structure interaction of a final stage blade,are used to verify the efficiency and accuracy of the proposed algorithm.The results show that the proposed method achieves the same accuracy as the implicit method while the computational time is reduced.In the case of the vortex-induced oscillation problem,the computational time can be reduced to 18.6%.In the case of the final stage blade vibration,the computational time can be reduced to 53.8%.
基金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.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as ＂fictitious＂ fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure＇ s rigid body shape and behavior, a rigid body constraint is enforced on the ＂fictitious＂ fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain （DLM/ FD） method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A coupled thermal-hydro-mechanical cohesive phase-field model for hydraulic fracturing in deep coal seams is presented.Heat exchange between the cold fluid and the hot rock is considered,and the thermal contribution terms between the cold fluid and the hot rock are derived.Heat transfer obeys Fourier's law,and porosity is used to relate the thermodynamic parameters of the fracture and matrix domains.The net pressure difference between the fracture and the matrix is neglected,and thus the fluid flow is modeled by the unified fluid-governing equations.The evolution equations of porosity and Biot's coefficient during hydraulic fracturing are derived from their definitions.The effect of coal cleats is considered and modeled by Voronoi polygons,and this approach is shown to have high accuracy.The accuracy of the proposed model is verified by two sets of fracturing experiments in multilayer coal seams.Subsequently,the differences in fracture morphology,fluid pressure response,and fluid pressure distribution between direct fracturing of coal seams and indirect fracturing of shale interlayers are explored,and the effects of the cluster number and cluster spacing on fracture morphology for multi-cluster fracturing are also examined.The numerical results show that the proposed model is expected to be a powerful tool for the fracturing design and optimization of deep coalbed methane.
基金This work was supported by the National Natural Science Foundation of China(Grant No.12172160)Shenzhen Science and Technology Program(Grant No.JCYJ20220818100600002)+1 种基金South-ern University of Science and Technology(Grant No.Y01326127)the Department of Science and Technology of Guangdong Province(Grant Nos.2020B1212030001 and 2021QN020642).
文摘We propose a suite of strategies for the parallel solution of fully implicit monolithic fluid-structure interaction(FSI).The solver is based on a modeling approach that uses the velocity and pressure as the primitive variables,which offers a bridge between computational fluid dynamics(CFD)and computational structural dynamics.The spatiotemporal discretization leverages the variational multiscale formulation and the generalized-αmethod as a means of providing a robust discrete scheme.In particular,the time integration scheme does not suffer from the overshoot phenomenon and optimally dissipates high-frequency spurious modes in both subproblems of FSI.Based on the chosen fully implicit scheme,we systematically develop a combined suite of nonlinear and linear solver strategies.Invoking a block factorization of the Jacobian matrix,the Newton-Raphson procedure is reduced to solving two smaller linear systems in the multi-corrector stage.The first is of the elliptic type,indicating that the algebraic multigrid method serves as a well-suited option.The second exhibits a two-by-two block structure that is analogous to the system arising in CFD.Inspired by prior studies,the additive Schwarz domain decomposition method and the block-factorization-based preconditioners are invoked to address the linear problem.Since the number of unknowns matches in both subdomains,it is straightforward to balance loads when parallelizing the algorithm for distributed-memory architectures.We use two representative FSI benchmarks to demonstrate the robustness,efficiency,and scalability of the overall FSI solver framework.In particular,it is found that the developed FSI solver is comparable to the CFD solver in several aspects,including fixed-size and isogranular scalability as well as robustness.
基金financially supported by the Riset Pengembangan dan Penerapan(RPP),Diponegoro University 2023 research scheme with contract number 609-18/UN7.D2/PP/VIII/2023.
文摘Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.
基金supported by the National Natural Science Foundation of China(Grant No.11925204).
文摘The Boussinesq equations,pivotal in the analysis of water wave dynamics,effectively model weakly nonlinear and long wave approximations.This study utilizes the complete discriminant system within a polynomial approach to derive exact traveling wave solutions for the coupled Boussinesq equation.The solutions are articulated through soliton,trigonometric,rational,and Jacobi elliptic functions.Notably,the introduction of Jacobi elliptic function solutions for this model marks a pioneering advancement.Contour plots of the solutions obtained by assigning values to various parameters are generated and subsequently analyzed.The methodology proposed in this study offers a systematic means to tackle nonlinear partial differential equations in mathematical physics,thereby enhancing comprehension of the physical attributes and dynamics of water waves.
基金financially supported by the National Natural Science Foundation of China(Grant No.42077231).
文摘This paper introduces a bond-based peridynamics(BB-PD)algorithm for crack identification,integrating the Delaunay triangulation method to accurately identify the structural characteristics of threedimensional(3D)cracks in rocks.A bond-based crack quantification standard is proposed to analyze the evolution of cracks of various sizes.A multi-attribute peridynamic model,developed using a multilayer algorithm,was employed to simulate the fracturing process of sandstone disks and semi-disks under varying temperatures,with the model calibrated and validated against experimental results.The simulation results show that temperature induces nonlinear degradation in the tensile strength and fracture toughness of sandstone,with 500℃ identified as the threshold temperature.Thermal cracks exhibit varying degrees of influence on Mode I cracks across different temperature ranges.Thermal damage significantly promotes the initiation and propagation of Mode I cracks in sandstone,thereby reducing its tensile strength and fracture toughness.Under applied loads,crack propagation in sandstone predominantly occurs during the failure stage,characterized by the rapid growth of longer cracks and a slow increase or reduction in shorter cracks.
基金supported by Alberta Health,Alberta Innovates,the Canada Research Chairs Program,the Canadian Institutes of Health Research,and the Natural Sciences and Engineering Research Council of Canada。
文摘Arsenic speciation in freshwater fish is crucial for providing meaningful consumption guidelines that allow the public to make informed decisions regarding its consumption.While marine fish have attractedmuch research interest due to their higher arsenic content,research on freshwater fish is limited due to the challenges in quantifying and identifying arsenic species present at trace levels.We describe here a sensitivemethod and its application to the quantification of arsenic species in freshwater fish.Arsenic species from fish tissues were extracted using a methanol/water mixture(1:1 vol.ratio)and ultrasound sonication.Anion-exchange high-performance liquid chromatography(HPLC)enabled separation of arsenobetaine(AsB),inorganic arsenite(iAs^(Ⅲ)),dimethylarsinic acid(DMA),monomethylarsonic acid(MMA),inorganic arsenate(iAs^(Ⅴ)),and three new arsenic species.Inductively coupled plasma mass spectrometry(ICPMS)provided highly sensitive and specific detection of arsenic.A limit of detection of 0.25μg/kg(wet weight fish tissue)was achieved for the five target arsenic species:AsB,iAs^(Ⅲ),DMA,MMA,and iAs^(Ⅴ).A series of experimentswere conducted to ensure the accuracy and validity of the analytical method.The method was successfully applied to the determination of arsenic species in lakewhitefish,northern pike,and walleye,with AsB,DMA,and iAs^(Ⅴ) being frequently detected.Three new arsenic species were detected,but their chromatographic retention times did not match with those of any available arsenic standards.Future research is necessary to elucidate the identity of these new arsenic species detected in freshwater fish.
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金supported by the National Natural Science Foundation of China(Grant No.52201323).
文摘In this study,a common-node DEM-SPH coupling model based on the shared node method is proposed,and a fluid–structure coupling method using the common-node discrete element method-smoothed particle hydrodynamics(DS-SPH)method is developed using LS-DYNA software.The DEM and SPH are established on the same node to create common-node DEM-SPH particles,allowing for fluid–structure interactions.Numerical simulations of various scenarios,including water entry of a rigid sphere,dam-break propagation over wet beds,impact on an ice plate floating on water and ice accumulation on offshore structures,are conducted.The interaction between DS particles and SPH fluid and the crack generation mechanism and expansion characteristics of the ice plate under the interaction of structure and fluid are also studied.The results are compared with available data to verify the proposed coupling method.Notably,the simulation results demonstrated that controlling the cutoff pressure of internal SPH particles could effectively control particle splashing during ice crushing failure.
基金funded by the BeijingNatural Science Foundation of China(8222003)National Natural Science Foundation of China(41807180).
文摘Based on global initiatives such as the clean energy transition and the development of renewable energy,the pumped storage power station has become a new and significant way of energy storage and regulation,and its construction environment is more complex than that of a traditional reservoir.In particular,the stability of the rock strata in the underground reservoirs is affected by the seepage pressure and rock stress,which presents some challenges in achieving engineering safety and stability.Using the advantages of the numerical simulation method in dealing deal with nonlinear problems in engineering stability,in this study,the stability of the underground reservoir of the Shidangshan(SDS)pumped storage power station was numerically calculated and quantitatively analyzed based on fluid-structure coupling theory,providing an important reference for the safe operation and management of the underground reservoir.First,using the COMSOL software,a suitablemechanicalmodel was created in accordance with the geological structure and project characteristics of the underground reservoir.Next,the characteristics of the stress field,displacement field,and seepage field after excavation of the underground reservoir were simulated in light of the seepage effect of groundwater on the nearby rock of the underground reservoir.Finally,based on the construction specifications and Molar-Coulomb criterion,a thorough evaluation of the stability of the underground reservoir was performed through simulation of the filling and discharge conditions and anti-seepage strengthening measures.The findings demonstrate that the numerical simulation results have a certain level of reliability and are in accordance with the stress measured in the project area.The underground reservoir excavation resulted in a maximum displacement value of the rock mass around the caverns of 3.56 mm in a typical section,and the safety coefficient of the parts,as determined using the Molar-Coulomb criterion,was higher than 1,indicating that the project as a whole is in a stable state.
文摘Spurious forces are a significant challenge for multi-scale methods,e.g.,the coupled atomistic/discrete dislocation(CADD)method.The assumption of isotropic matter in the continuum domain is a critical factor leading to such forces.This study aims to minimize spurious forces,ensuring that atomic dislocations experience more precise forces from the continuum domain.The authors have already implemented this idea using a simplified and unrealistic slipping system.To create a comprehensive and realistic model,this paper considers all possible slip systems in the face center cubic(FCC)lattice structure,and derives the required relationships for the displacement fields.An anisotropic version of the three-dimensional CADD(CADD3D)method is presented,which generates the anisotropic displacement fields for the partial dislocations in all the twelve slip systems of the FCC lattice structure.These displacement fields are tested for the most probable slip systems of aluminum,nickel,and copper with different anisotropic levels.Implementing these anisotropic displacement fields significantly reduces the spurious forces on the slip systems of FCC materials.This improvement is particularly pronounced at greater distances from the interface and in more anisotropic materials.Furthermore,the anisotropic CADD3D method enhances the spurious stress difference between the slip systems,particularly for materials with higher anisotropy.
基金open foundation of the Hubei Key Laboratory of Theory and Application of Advanced Materials Mechanicsthe Open Foundation of Hubei Key Laboratory of Engineering Structural Analysis and Safety Assessment.
文摘A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.
基金the National Natural Science Foundation of China(No.51875062,No.52205336)the China Postdoctoral Science Foundation(No.2021M700567).
文摘The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.
基金supported by the Natural Science Foundation of Heilongjiang Province(No.TD2022E003)the National Natural Science Foundation of China(No.52175393)。
文摘42CrMo steel has the characteristics of high strength,high wear resistance,high impact resistance,and fatigue resistance.Therefore,drilling 42CrMo steel has always been a challenging task.The indexable drill bit has the advantages of high processing efficiency and low processing cost and has been widely used in the field of aerospace hole processing.To better understand the machining mechanism of the indexable drill bit,this paper uses the Coupled EulerianLagrangian method(CEL)to simulate the three-dimensional drilling model for the first time.The simulation results of the drilling force obtained by the CEL method and Lagrangian method are compared with the experimental results.It is verified that the CEL method is easy to converge and can avoid the problem of program interruption caused by mesh distortion,and the CEL simulation value is more consistent with the actual value.Secondly,the simulation results of cutting force and blade cutting edge node temperature under different process parameters are extracted.The variation of time domain cutting force,frequency domain cutting force and tool temperature with process parameters are obtained.This study provides a new method for the prediction of cutting performance and the optimization of process parameters of indexable drills.
基金supported by the National Natural Science Foundation of China(Grant No.12271433)。
文摘In this paper,two types of fractional nonlinear equations in Caputo sense,time-fractional Newell–Whitehead equation(FNWE)and time-fractional generalized Hirota–Satsuma coupled KdV system(HS-cKdVS),are investigated by means of the q-homotopy analysis method(q-HAM).The approximate solutions of the proposed equations are constructed in the form of a convergent series and are compared with the corresponding exact solutions.Due to the presence of the auxiliary parameter h in this method,just a few terms of the series solution are required in order to obtain better approximation.For the sake of visualization,the numerical results obtained in this paper are graphically displayed with the help of Maple.
基金Project supported by the National Natural Science Foundation of China(No.12102131)the Natural Science Foundation of Henan Province of China(No.242300420248)the International Science and Technology Cooperation Project of Henan Province of China(No.242102521010)。
文摘The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.
基金the Science and Technology Research Project of Henan Province (242102231052)the Key Scientific Research Plan of Colleges and Universities in Henan Province (23B140006)the Natural Science Foundation of Jiangxi Province (20224BAB201018)。
文摘In this paper,a simple direct space-time semi-analytical meshless scheme is proposed for the numerical approximation of the coupled Burgers'equations.During the whole solution procedure,two different schemes are considered in terms of radial and non-radial basis functions.The time-dependent variable in the first radial scheme is directly considered as the normal space variables to formulate an"isotropic"space-time radial basis function.The second non-radial scheme considered relationship between time-dependent and spacedependent variables.Under such circumstance,we can get a one-step space-time meshless scheme.The numerical findings demonstrate that the proposed meshless schemes are precise,user-friendly,and effective in solving the coupled Burgers'equations.
文摘This paper presents a new approach to synthesize admittance function polynomials and coupling matrices for coupled resonator filters. The N + 2 transversal network method is applied to study a coupled resonator filter. This method allowed us to determine the polynomials of the reflection and transmission coefficients. A study is made for a 4 poles filter with 2 transmission zeros between the N + 2 transversal network method and the one found in the literature. A MATLAB code was designed for the numerical simulation of these coefficients for the 6, 8, and 10 pole filter with 4 transmission zeros.
