Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering s...Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.展开更多
This work concerns numerical modeling of fluid-structure interaction（FSI） problems in a uniform smoothed particle hydrodynamics（SPH） framework. It combines a transport-velocity SPH scheme, advancing fluid motions,...This work concerns numerical modeling of fluid-structure interaction（FSI） problems in a uniform smoothed particle hydrodynamics（SPH） framework. It combines a transport-velocity SPH scheme, advancing fluid motions, with a total Lagrangian SPH formulation dealing with the structure deformations. Since both fluid and solid governing equations are solved in SPH framework, while coupling becomes straightforward, the momentum conservation of the FSI system is satisfied strictly. A well-known FSI benchmark test case has been performed to validate the modeling and to demonstrate its potential.展开更多
Pure commercial titanium was welded with two types of stainless steel,namely SUS 304 austenitic stainless steel and SUS 821L1 duplex stainless steel.The wavy interface of SUS 821L1 was smaller than that of SUS 304.The...Pure commercial titanium was welded with two types of stainless steel,namely SUS 304 austenitic stainless steel and SUS 821L1 duplex stainless steel.The wavy interface of SUS 821L1 was smaller than that of SUS 304.The vortex zone was observed from both longitudinal and transverse directions,and its composition was analyzed.The interface of Ti/SUS 821L11 was able to bear 401−431 MPa shear load while that of Ti/SUS 304 could withstand 352−387 MPa.The weldability window was used to analyze experimental phenomenon.Furthermore,the smoothed particle hydrodynamics(SPH)numerical simulation method was used to simulate the wavy interface.The trend of wavelength and amplitude change with strength and the stand-offs was consistent with the experimental results.展开更多
In this paper, the smoothed particle hydrodynamics（SPH） method is used to build a numerical wave-current tank（NWCT）. The wave is generated by using a piston-type wave generator and is absorbed by using a sponge la...In this paper, the smoothed particle hydrodynamics（SPH） method is used to build a numerical wave-current tank（NWCT）. The wave is generated by using a piston-type wave generator and is absorbed by using a sponge layer. The uniform current field is generated by simultaneously imposing the directional velocity and hydrostatic pressure in both inflow and outflow regions set below the NWCT. Particle cyclic boundaries are also implemented for recycling the Lagrangian fluid particles. Furthermore, to shorten the time to reach a steady state, a temporary rigid-lid treatment for the water surface is proposed. It turns out to be very effective for weakening the undesired oscillatory flow at the beginning stage of the current generation. The calculated water surface elevation and horizontal-velocity profile are validated against the available experimental data. Satisfactory agreements are obtained, demonstrating the good capability of the NWCT.展开更多
The violent water entry of flat plates is investigated using a Riemann-arbitrary Eulerian-Lagrangian（ALE） smoothed particle hydrodynamics（SPH） model. The test conditions are of interest for problems related to air...The violent water entry of flat plates is investigated using a Riemann-arbitrary Eulerian-Lagrangian（ALE） smoothed particle hydrodynamics（SPH） model. The test conditions are of interest for problems related to aircraft and helicopter emergency landing in water. Three main parameters are considered： the horizontal velocity, the approach angle（i.e., vertical to horizontal velocity ratio） and the pitch angle, a. Regarding the latter, small angles are considered in this study. As described in the theoretical work by Zhao and Faltinsen（1993）, for small a a very thin, high-speed jet of water is formed, and the time-spatial gradients of the pressure field are extremely high. These test conditions are very challenging for numerical solvers. In the present study an enhanced SPH model is firstly tested on a purely vertical impact with deadrise angle α=4°. An in-depth validation against analytical solutions and experimental results is carried out, highlighting the several critical aspects of the numerical modelling of this kind of flow, especially when pressure peaks are to be captured. A discussion on the main difficulties when comparing to model scale experiments is also provided. Then, the more realistic case of a plate with both horizontal and vertical velocity components is discussed and compared to ditching experiments recently carried out at CNR-INSEAN. In the latter case both 2-D and 3-D simulations are considered and the importance of 3-D effects on the pressure peak is discussed for α=4° and α=10°.展开更多
Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the ...Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.展开更多
This study aims to numerically and experimentally investigate the response of a medium strength rock material under unconfined compression loading up to failure. The unconfined compressive strength(UCS) is one of the ...This study aims to numerically and experimentally investigate the response of a medium strength rock material under unconfined compression loading up to failure. The unconfined compressive strength(UCS) is one of the most important parameters in characterising rock material behaviour. Hence the UCS is crucial in understanding the failure mechanism of fractured rocks. An effective approach to determine the UCS and to investigate the behaviours of rock materials under unconfined compression is essential in the majority of research fields of rock mechanics. The experimental configuration for the unconfined compression test, suggested by the protocols of the ASTM standard, has some limitations which affect the accuracy in determination of the real UCS. Among several alternative configurations proposed, the Mogi’s configuration seems to be the most appropriate one. Therefore, the ASTM and Mogi’s configurations were used to perform the tests on a medium strength rock material, i.e. Pietra Serena sandstone. The results using two configurations were discussed in terms of the differences. The tests were also replicated in LSDYNA using a finite element method(FEM) coupled smooth particle hydrodynamics(SPH) technique.This technique is employed in this study due to its capabilities to cope with large deformation issues related to the rocks. An advanced material model, called the Karagozian and Case Concrete(KCC) model,is implemented in the numerical simulations. The KCC model consists of three independent fixed failure surfaces and it can consider the damage accumulation based on the current state of stress among these failure surfaces. An equation-of-state(EOS) is used in conjunction with KCC material model for decoupling the volumetric and deviatoric responses. The numerical and experimental results were finally compared with the focus on the stress-strain diagram and the failure patterns. The comparison shows that the numerical results are in good agreement with the experimental results.展开更多
Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to...Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.展开更多
Over the past 13 years, ANDRITZ Hydro has developed an in-house tool based on the SPH-ALE method for applications in flow simulations in hydraulic turbines. The initial motivation is related to the challenging simulat...Over the past 13 years, ANDRITZ Hydro has developed an in-house tool based on the SPH-ALE method for applications in flow simulations in hydraulic turbines. The initial motivation is related to the challenging simulation of free surface flows in Pelton turbines, where highly dynamic water jets interact with rotating buckets, creating thin water jets traveling inside the housing and possibly causing disturbances on the runner. The present paper proposes an overview of industrial applications allowed by the developed tool, including design evaluation of Pelton runners and casings, transient operation of Pelton units and free surface flows in hydraulic structures.展开更多
With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during...With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during the transient flooding process. Therefore, an enhanced smoothed particle hydrodynamics（SPH） model is applied in this paper to investigate the response of a simplified cabin model under the condition of the transient water flooding. The enhanced SPH model is presented firstly including the governing equations, the diffusive terms, the boundary implementations and then an algorithm regarding the coupling motions of six degrees of freedom（6-DOF） between the structure and the fluid is described. In the numerical results, a non-damaged cabin floating under the rest condition is simulated. It is shown that a stable floating state can be reached and maintained by using the present SPH scheme. After that, three-dimensional（3-D） test cases of the damaged cabin with a hole at different locations are simulated. A series of model tests are also carried out for the validation. Fairly good agreements are achieved between the numerical results and the experimental data. Relevant conclusions are drawn with respect to the mechanism of the responses of the damaged cabin model under water flooding conditions.展开更多
Since its inception,the full Lagrangian meshless smoothed particle hydrodynamics(SPH)has experienced a tremendous enhancement in methodology and impacted a range of multi-physics applications in science and engineerin...Since its inception,the full Lagrangian meshless smoothed particle hydrodynamics(SPH)has experienced a tremendous enhancement in methodology and impacted a range of multi-physics applications in science and engineering.This review presents a concise survey on latest developments and achievements of the SPH method,including:(1)Brief review of theory and fundamental with kernel corrections,(2)The Riemann-based SPH method with dissipation limiting and high-order data reconstruction by using MUSCL,WENO and MOOD schemes,(3)Particle neighbor searching with particle sorting and efficient dual-criteria time stepping schemes,(4)Total Lagrangian formulation with stablized,dynamics relaxation and hourglass control schemes,(5)Fluid-structure interaction scheme with interface treatments and multi-resolution discretizations,(6)Novel applications of particle relaxation in SPH methodology for mesh and particle generations.Last but not least,benchmark tests for validating computational accuracy,convergence,robustness and efficiency are also supplied accordingly.展开更多
While Eulerian smoothed particle hydrodynamics(SPH)method has received increasing attention in scientific and industrial communities owing to its high spatial accuracy,it exhibits excessive numerical dissipation due t...While Eulerian smoothed particle hydrodynamics(SPH)method has received increasing attention in scientific and industrial communities owing to its high spatial accuracy,it exhibits excessive numerical dissipation due to the fact that the flux is derived in particle pair pattern.In this paper,we adopt a one-dimensional weighted essentially non-oscillatory(WENO)reconstruction to reduce the numerical dissipation and improve the overall accuracy particularly in capturing the contact discontinuity.The underlying principle is to construct a 4-point stencil along the interacting line of each particle pair and then the WENO scheme is applied to reconstruct the initial states of the Riemann problem which determines the flow flux.A set of benchmark tests for both compressible and incompressible flows are studied to investigate the accuracy,robustness and versatility of the proposed Eulerian SPH method with the WENO reconstruction(ESPH-WENO).展开更多
The mechanical responses and ultimate failure patterns of rocks are associated with the failure mechanism evolution.In this study,smoothed particle hydrodynamics(SPH)method with the mixed-mode failure model is propose...The mechanical responses and ultimate failure patterns of rocks are associated with the failure mechanism evolution.In this study,smoothed particle hydrodynamics(SPH)method with the mixed-mode failure model is proposed to probe into failure mechanism evolutions for disc specimens upon loading.The tensile damage model and the Drucker-Prager model are used to calculate the tensile failure and shear failure of the material,respectively.It is concluded that for flaw-unfilled disc specimens,the crack coalescence mechanism in the rock bridge area is affected by the flaw inclination angle and the material property.Considering disc specimens with filled flaws,the incremental rate of tensile damage grows more rapidly when the disc and filling material have a closer ratio of tensile strength to cohesion,which makes the entire specimen response greater brittleness.Furthermore,with the increasing non-uniformity of filling distribution,the incremental rate of tensile-activated damage decreases and the disc specimen performs more ductile.Besides,the influence of the fillings is greater when the flaw inclination angle is approaching 45°.It is proved that the proposed SPH method can be used to simulate the failure mechanism evolution of rocks,which lays a foundation for the study of more complex rock failure.展开更多
This paper presents a review on state-of-the-art of developments corresponding to fluid-structure interaction(FSI)solvers developed within the context of particle methods.The paper reviews and highlights the potential...This paper presents a review on state-of-the-art of developments corresponding to fluid-structure interaction(FSI)solvers developed within the context of particle methods.The paper reviews and highlights the potential robustness of entirely Lagrangian meshfree FSI solvers in reproducing FSI corresponding to extreme events and portrays the future perspectives for systematic developments towards reliable engineering applications with respect to rapid advances in technology and emergence of so-called advanced materials that can result in complex and highly non-linear structural responses.Accordingly,the paper highlights the necessity for reproduction of comprehensive structural responses,including viscoelastic,elastoplastic and progressive damages/failures,by the advanced FSI solvers developed within the context of particle methods.In this regard,extensions of the structure model are suggested to be conducted in a variationally consistent framework to ensure stability,accuracy and physical reliability including thermodynamic consistency.The paper reviews basics of mathematical and numerical modelling for entirely Lagrangian meshfree hydroelastic FSI solvers and presents a brief background on extensions of such solvers towards reproducing viscoelastic structural responses.Some preliminary numerical results on structural viscoelasticity achieved by an extended Hamiltonian SPH(HSPH)model are presented.This vision paper also concisely portrays the future perspectives for systematic development of particle-based FSI solvers.展开更多
This paper shows the results of the smooth particle hydrodynamics（SPH） modelling of the hydraulic jump at an abrupt drop,where the transition from supercritical to subcritical flow is characterised by several flow p...This paper shows the results of the smooth particle hydrodynamics（SPH） modelling of the hydraulic jump at an abrupt drop,where the transition from supercritical to subcritical flow is characterised by several flow patterns depending upon the inflow and tailwater conditions. SPH simulations are obtained by a pseudo-compressible XSPH scheme with pressure smoothing; turbulent stresses are represented either by an algebraic mixing-length model, or by a two-equation k-ε model. The numerical model is applied to analyse the occurrence of oscillatory flow conditions between two different jump types characterised by quasi-periodic oscillation,and the results are compared with experiments performed at the hydraulics laboratory of Bari Technical University. The purpose of this paper is to obtain a deeper understanding of the physical features of a flow which is in general difficult to be reproduced numerically,owing to its unstable character： in particular, vorticity and turbulent kinetic energy fields, velocity, water depth and pressure spectra downstream of the jump, and velocity and pressure cross-correlations can be computed and analysed.展开更多
Free surface flows are of significant interest in Computational Fluid Dynamics（CFD）. However, violent water wave impact simulation especially when free surface breaks or impacts on solid wall can be a big challenge ...Free surface flows are of significant interest in Computational Fluid Dynamics（CFD）. However, violent water wave impact simulation especially when free surface breaks or impacts on solid wall can be a big challenge for many CFD techniques. Smoothed Particle Hydrodynamics（SPH） has been reported as a robust and reliable method for simulating violent free surface flows. Weakly compressible SPH（WCSPH） uses an equation of state with a large sound speed, and the results of the WCSPH can induce a noisy pressure field and spurious oscillation of pressure in time history for wave impact problem simulation. As a remedy, the truly incompressible SPH（ISPH） technique was introduced, which uses a pressure Poisson equation to calculate the pressure. Although the pressure distribution in the whole field obtained by ISPH is smooth, the stability of the techniques is still an open discussion. In this paper, a new free surface identification scheme and solid boundary handling method are introduced to improve the accuracy of ISPH. This modified ISPH is used to study dam breaking flow and violent tank sloshing flows. On the comparative study of WCSPH and ISPH, the accuracy and efficiency are assessed and the results are compared with the experimental data.展开更多
The present work has been performed in the context of the European H2020 project increased SAfety and Robust certification for ditching of Aircrafts and Helicopters(SARAH)dedicated to improving the safety during aircr...The present work has been performed in the context of the European H2020 project increased SAfety and Robust certification for ditching of Aircrafts and Helicopters(SARAH)dedicated to improving the safety during aircraft ditching,together with a better understanding of the physics involved during those crucial events.Both numerical and experimental aspects are explored during this project.The present study focuses on the application of the smoothed particle hydrodynamics(SPH)method to the simulation of helicopter ditching,as this method has proved to be particularly adapted to free surface impact cases.Simulations are performed for three different impact configurations,for which the numerical solutions are compared with the experimental results(forces and kinematics)obtained at the wave basin of Ecole Centrale Nantes on a mock-up shape provided by Airbus Helicopters.Elements of sensitivity analysis are also provided when needed,to assess the role of some parameters involved in the helicopter behavior and the fluid pressure forces exerted during the impact.展开更多
Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element meth...Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.展开更多
As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the ...As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.展开更多
Microwave precondition has been highlighted as a promising technology for softening the rock mass prior to rock breakage by machine to reduce drill bit/cutter wear as well as inverse production rate.To numerically exp...Microwave precondition has been highlighted as a promising technology for softening the rock mass prior to rock breakage by machine to reduce drill bit/cutter wear as well as inverse production rate.To numerically explore the effect of numerical parameters on rock static strength simulation,and determine the numerical mechanical parameters of microwave-treated basalts for future drilling and cutting simulations,numerical models of uniaxial compression strength(UCS)and Brazilian tensile strength(BTS)were established with the coupling of smoothed particle hydrodynamics and finite element method(SPH-FEM).To eliminate the large rock strength errors caused by microwave-induced damage,the cohesion and internal friction angle of microwave-treated basalt specimens with the same microwave treatment parameters were calibrated based on a linear Mohr-Coulomb theory.Based on parametric sensitivity analysis of SPH simulation of UCS and BTS,experimental UCS and BTS values were simultaneously captured according to the same set of calibrated cohesion and internal friction angle data,and the UCS modeling results are in good agreement with experimental tests.Furthermore,the effect of microwave irradiation parameter on the basalt mechanical behaviors was evaluated.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.51779003)the National Key Research and Development Project of China(Grant No.2018YFB0704000)
文摘Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.
基金the financial support by Deutsche Forschungsgemeinschaft (Grant No. DFG HU1527/6-1)for the present work
文摘This work concerns numerical modeling of fluid-structure interaction（FSI） problems in a uniform smoothed particle hydrodynamics（SPH） framework. It combines a transport-velocity SPH scheme, advancing fluid motions, with a total Lagrangian SPH formulation dealing with the structure deformations. Since both fluid and solid governing equations are solved in SPH framework, while coupling becomes straightforward, the momentum conservation of the FSI system is satisfied strictly. A well-known FSI benchmark test case has been performed to validate the modeling and to demonstrate its potential.
文摘Pure commercial titanium was welded with two types of stainless steel,namely SUS 304 austenitic stainless steel and SUS 821L1 duplex stainless steel.The wavy interface of SUS 821L1 was smaller than that of SUS 304.The vortex zone was observed from both longitudinal and transverse directions,and its composition was analyzed.The interface of Ti/SUS 821L11 was able to bear 401−431 MPa shear load while that of Ti/SUS 304 could withstand 352−387 MPa.The weldability window was used to analyze experimental phenomenon.Furthermore,the smoothed particle hydrodynamics(SPH)numerical simulation method was used to simulate the wavy interface.The trend of wavelength and amplitude change with strength and the stand-offs was consistent with the experimental results.
基金Project supported by the National Natural Science Foun-dation of China(Grant Nos.51379144,51479135,51679167 and 51709201)the China Postdoctoral Science Foundation(Grant No.2017M621074)
文摘In this paper, the smoothed particle hydrodynamics（SPH） method is used to build a numerical wave-current tank（NWCT）. The wave is generated by using a piston-type wave generator and is absorbed by using a sponge layer. The uniform current field is generated by simultaneously imposing the directional velocity and hydrostatic pressure in both inflow and outflow regions set below the NWCT. Particle cyclic boundaries are also implemented for recycling the Lagrangian fluid particles. Furthermore, to shorten the time to reach a steady state, a temporary rigid-lid treatment for the water surface is proposed. It turns out to be very effective for weakening the undesired oscillatory flow at the beginning stage of the current generation. The calculated water surface elevation and horizontal-velocity profile are validated against the available experimental data. Satisfactory agreements are obtained, demonstrating the good capability of the NWCT.
基金partially received funding from the European Union's Horizon 2020 Research and Innovation Programme (Grant No. 724139)
文摘The violent water entry of flat plates is investigated using a Riemann-arbitrary Eulerian-Lagrangian（ALE） smoothed particle hydrodynamics（SPH） model. The test conditions are of interest for problems related to aircraft and helicopter emergency landing in water. Three main parameters are considered： the horizontal velocity, the approach angle（i.e., vertical to horizontal velocity ratio） and the pitch angle, a. Regarding the latter, small angles are considered in this study. As described in the theoretical work by Zhao and Faltinsen（1993）, for small a a very thin, high-speed jet of water is formed, and the time-spatial gradients of the pressure field are extremely high. These test conditions are very challenging for numerical solvers. In the present study an enhanced SPH model is firstly tested on a purely vertical impact with deadrise angle α=4°. An in-depth validation against analytical solutions and experimental results is carried out, highlighting the several critical aspects of the numerical modelling of this kind of flow, especially when pressure peaks are to be captured. A discussion on the main difficulties when comparing to model scale experiments is also provided. Then, the more realistic case of a plate with both horizontal and vertical velocity components is discussed and compared to ditching experiments recently carried out at CNR-INSEAN. In the latter case both 2-D and 3-D simulations are considered and the importance of 3-D effects on the pressure peak is discussed for α=4° and α=10°.
基金Project supported by the National Natural Science Foundation of China(No.52109068)the Water Conservancy Technology Project of Jiangsu Province of China(No.2022060)。
文摘Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.
文摘This study aims to numerically and experimentally investigate the response of a medium strength rock material under unconfined compression loading up to failure. The unconfined compressive strength(UCS) is one of the most important parameters in characterising rock material behaviour. Hence the UCS is crucial in understanding the failure mechanism of fractured rocks. An effective approach to determine the UCS and to investigate the behaviours of rock materials under unconfined compression is essential in the majority of research fields of rock mechanics. The experimental configuration for the unconfined compression test, suggested by the protocols of the ASTM standard, has some limitations which affect the accuracy in determination of the real UCS. Among several alternative configurations proposed, the Mogi’s configuration seems to be the most appropriate one. Therefore, the ASTM and Mogi’s configurations were used to perform the tests on a medium strength rock material, i.e. Pietra Serena sandstone. The results using two configurations were discussed in terms of the differences. The tests were also replicated in LSDYNA using a finite element method(FEM) coupled smooth particle hydrodynamics(SPH) technique.This technique is employed in this study due to its capabilities to cope with large deformation issues related to the rocks. An advanced material model, called the Karagozian and Case Concrete(KCC) model,is implemented in the numerical simulations. The KCC model consists of three independent fixed failure surfaces and it can consider the damage accumulation based on the current state of stress among these failure surfaces. An equation-of-state(EOS) is used in conjunction with KCC material model for decoupling the volumetric and deviatoric responses. The numerical and experimental results were finally compared with the focus on the stress-strain diagram and the failure patterns. The comparison shows that the numerical results are in good agreement with the experimental results.
基金the National Natural Science Foundation of China(No.52271316)the Guangdong Basic and Applied Basic Research Foundation(No.2023A1515030262).
文摘Predicting the response of liquefied natural gas(LNG)contained in vessels subjected to external waves is extremely important to ensure the safety of the transportation process.In this study,the coupled behavior due to ship motion and liquid tank sloshing has been simulated by the Smoothed-Particle Hydrodynamics(SPH)method.Firstly,the sloshing flow in a rectangular tank was simulated and the related loads were analyzed to verify and validate the accuracy of the present SPH solver.Then,a three-dimensional simplified LNG carrier model,including two prismatic liquid tanks and a wave tank,was introduced.Different conditions were examined corresponding to different wave lengths,wave heights,wave heading angles,and tank loading rates.Finally,the effects of liquid tank loading rate on LNG ship motions and sloshing loading were analyzed,thereby showing that the SPH method can effectively provide useful indications for the design of liquid cargo ships.
文摘Over the past 13 years, ANDRITZ Hydro has developed an in-house tool based on the SPH-ALE method for applications in flow simulations in hydraulic turbines. The initial motivation is related to the challenging simulation of free surface flows in Pelton turbines, where highly dynamic water jets interact with rotating buckets, creating thin water jets traveling inside the housing and possibly causing disturbances on the runner. The present paper proposes an overview of industrial applications allowed by the developed tool, including design evaluation of Pelton runners and casings, transient operation of Pelton units and free surface flows in hydraulic structures.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.U1430236,51609045)
文摘With the quasi-static analysis method, the terminal floating state of a damaged ship is usually evaluated for the risk assessment. But this is not enough since the ship has the possibility to lose its stability during the transient flooding process. Therefore, an enhanced smoothed particle hydrodynamics（SPH） model is applied in this paper to investigate the response of a simplified cabin model under the condition of the transient water flooding. The enhanced SPH model is presented firstly including the governing equations, the diffusive terms, the boundary implementations and then an algorithm regarding the coupling motions of six degrees of freedom（6-DOF） between the structure and the fluid is described. In the numerical results, a non-damaged cabin floating under the rest condition is simulated. It is shown that a stable floating state can be reached and maintained by using the present SPH scheme. After that, three-dimensional（3-D） test cases of the damaged cabin with a hole at different locations are simulated. A series of model tests are also carried out for the validation. Fairly good agreements are achieved between the numerical results and the experimental data. Relevant conclusions are drawn with respect to the mechanism of the responses of the damaged cabin model under water flooding conditions.
基金This work was supported by the Deutsche Forschungsgemeinschaft(DFG)(Grant No.HU1527/12-4)。
文摘Since its inception,the full Lagrangian meshless smoothed particle hydrodynamics(SPH)has experienced a tremendous enhancement in methodology and impacted a range of multi-physics applications in science and engineering.This review presents a concise survey on latest developments and achievements of the SPH method,including:(1)Brief review of theory and fundamental with kernel corrections,(2)The Riemann-based SPH method with dissipation limiting and high-order data reconstruction by using MUSCL,WENO and MOOD schemes,(3)Particle neighbor searching with particle sorting and efficient dual-criteria time stepping schemes,(4)Total Lagrangian formulation with stablized,dynamics relaxation and hourglass control schemes,(5)Fluid-structure interaction scheme with interface treatments and multi-resolution discretizations,(6)Novel applications of particle relaxation in SPH methodology for mesh and particle generations.Last but not least,benchmark tests for validating computational accuracy,convergence,robustness and efficiency are also supplied accordingly.
基金This work was supported by the China Scholarship Council(Grant No.No.201906120035)Chi Zhang and Xiangyu Hu would like to express their gratitude to Deutsche Forschungsgemeinschaft(DFG)for their sponsorship(Grant No.DFG HU1527/12-4).
文摘While Eulerian smoothed particle hydrodynamics(SPH)method has received increasing attention in scientific and industrial communities owing to its high spatial accuracy,it exhibits excessive numerical dissipation due to the fact that the flux is derived in particle pair pattern.In this paper,we adopt a one-dimensional weighted essentially non-oscillatory(WENO)reconstruction to reduce the numerical dissipation and improve the overall accuracy particularly in capturing the contact discontinuity.The underlying principle is to construct a 4-point stencil along the interacting line of each particle pair and then the WENO scheme is applied to reconstruct the initial states of the Riemann problem which determines the flow flux.A set of benchmark tests for both compressible and incompressible flows are studied to investigate the accuracy,robustness and versatility of the proposed Eulerian SPH method with the WENO reconstruction(ESPH-WENO).
基金This work is supported in part by the National Natural Science Foundation of China(No.51779084)。
文摘The mechanical responses and ultimate failure patterns of rocks are associated with the failure mechanism evolution.In this study,smoothed particle hydrodynamics(SPH)method with the mixed-mode failure model is proposed to probe into failure mechanism evolutions for disc specimens upon loading.The tensile damage model and the Drucker-Prager model are used to calculate the tensile failure and shear failure of the material,respectively.It is concluded that for flaw-unfilled disc specimens,the crack coalescence mechanism in the rock bridge area is affected by the flaw inclination angle and the material property.Considering disc specimens with filled flaws,the incremental rate of tensile damage grows more rapidly when the disc and filling material have a closer ratio of tensile strength to cohesion,which makes the entire specimen response greater brittleness.Furthermore,with the increasing non-uniformity of filling distribution,the incremental rate of tensile-activated damage decreases and the disc specimen performs more ductile.Besides,the influence of the fillings is greater when the flaw inclination angle is approaching 45°.It is proved that the proposed SPH method can be used to simulate the failure mechanism evolution of rocks,which lays a foundation for the study of more complex rock failure.
基金The authors would like to express their gratitude to Prof.De-cheng Wan at Shanghai Jiao Tong University for invitation for this vision paper.The first author,A.Khayyer,would like to express his sincere appreciation to Prof.Antonio J.Gil at Swansea University and Dr Chun Hean Lee at the University of Glasgow for discussions regarding viscoelastic and elastoplastic modelling.The authors appreciate the research grants by Japan Society for the Promotion of Science(JSPS)(Grant No.JP18K04368,JP21H01433 and JP21K14250).
文摘This paper presents a review on state-of-the-art of developments corresponding to fluid-structure interaction(FSI)solvers developed within the context of particle methods.The paper reviews and highlights the potential robustness of entirely Lagrangian meshfree FSI solvers in reproducing FSI corresponding to extreme events and portrays the future perspectives for systematic developments towards reliable engineering applications with respect to rapid advances in technology and emergence of so-called advanced materials that can result in complex and highly non-linear structural responses.Accordingly,the paper highlights the necessity for reproduction of comprehensive structural responses,including viscoelastic,elastoplastic and progressive damages/failures,by the advanced FSI solvers developed within the context of particle methods.In this regard,extensions of the structure model are suggested to be conducted in a variationally consistent framework to ensure stability,accuracy and physical reliability including thermodynamic consistency.The paper reviews basics of mathematical and numerical modelling for entirely Lagrangian meshfree hydroelastic FSI solvers and presents a brief background on extensions of such solvers towards reproducing viscoelastic structural responses.Some preliminary numerical results on structural viscoelasticity achieved by an extended Hamiltonian SPH(HSPH)model are presented.This vision paper also concisely portrays the future perspectives for systematic development of particle-based FSI solvers.
文摘This paper shows the results of the smooth particle hydrodynamics（SPH） modelling of the hydraulic jump at an abrupt drop,where the transition from supercritical to subcritical flow is characterised by several flow patterns depending upon the inflow and tailwater conditions. SPH simulations are obtained by a pseudo-compressible XSPH scheme with pressure smoothing; turbulent stresses are represented either by an algebraic mixing-length model, or by a two-equation k-ε model. The numerical model is applied to analyse the occurrence of oscillatory flow conditions between two different jump types characterised by quasi-periodic oscillation,and the results are compared with experiments performed at the hydraulics laboratory of Bari Technical University. The purpose of this paper is to obtain a deeper understanding of the physical features of a flow which is in general difficult to be reproduced numerically,owing to its unstable character： in particular, vorticity and turbulent kinetic energy fields, velocity, water depth and pressure spectra downstream of the jump, and velocity and pressure cross-correlations can be computed and analysed.
基金supported by the National Natural Science Foundations of China(Grant Nos.51009034 and 51279041)Fundamental Research Funds for the Central Universities(Grant Nos.HEUCDZ1202 and HEUCF120113)Pre-Research Foundation of General Armament Department of China(Grant No.9140A14020712CB01158)
文摘Free surface flows are of significant interest in Computational Fluid Dynamics（CFD）. However, violent water wave impact simulation especially when free surface breaks or impacts on solid wall can be a big challenge for many CFD techniques. Smoothed Particle Hydrodynamics（SPH） has been reported as a robust and reliable method for simulating violent free surface flows. Weakly compressible SPH（WCSPH） uses an equation of state with a large sound speed, and the results of the WCSPH can induce a noisy pressure field and spurious oscillation of pressure in time history for wave impact problem simulation. As a remedy, the truly incompressible SPH（ISPH） technique was introduced, which uses a pressure Poisson equation to calculate the pressure. Although the pressure distribution in the whole field obtained by ISPH is smooth, the stability of the techniques is still an open discussion. In this paper, a new free surface identification scheme and solid boundary handling method are introduced to improve the accuracy of ISPH. This modified ISPH is used to study dam breaking flow and violent tank sloshing flows. On the comparative study of WCSPH and ISPH, the accuracy and efficiency are assessed and the results are compared with the experimental data.
基金Supported by the European Union Horizon 2020 Research and Innovation Program(Grant No.724139)The authors thank the Pole de Calcul et de Donnees Marines(PCDM)for providing DATARMOR storage and computational resources(http://www.ifremer.fr/pcdm),and the Institut de Calcul Intensif(ICI)(Grant No.2018-P1804060).
文摘The present work has been performed in the context of the European H2020 project increased SAfety and Robust certification for ditching of Aircrafts and Helicopters(SARAH)dedicated to improving the safety during aircraft ditching,together with a better understanding of the physics involved during those crucial events.Both numerical and experimental aspects are explored during this project.The present study focuses on the application of the smoothed particle hydrodynamics(SPH)method to the simulation of helicopter ditching,as this method has proved to be particularly adapted to free surface impact cases.Simulations are performed for three different impact configurations,for which the numerical solutions are compared with the experimental results(forces and kinematics)obtained at the wave basin of Ecole Centrale Nantes on a mock-up shape provided by Airbus Helicopters.Elements of sensitivity analysis are also provided when needed,to assess the role of some parameters involved in the helicopter behavior and the fluid pressure forces exerted during the impact.
基金supported by the National Numerical Wind Tunnel Project (Grant No. NNW2019ZT2-B02)the National Natural Science Foundation of China (Grant Nos. 12032002,51779003,and 11902005)the SinoGerman Mobility Programme (Grant No. M-0210)。
文摘Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.
基金supported by the 111 China Project(Grant No.B16003)the National Natural Science Foundation of China(Grant Nos.11290151,11702022,and 11221202)
文摘As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics(SPH) method is used to model the compressible fluid, the natural coordinate formulation(NCF) and absolute nodal coordinate formulation(ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit(GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.
基金the National Natural Science Foundation of China (No. 51774323)the Natural Science Foundation of Hunan Province, China (No. 2020JJ4704)+1 种基金the Fundamental Research Funds for the Central Universities of Central South University, China (No. 2018zzts216) the financial support from the China Scholarship Councilthe support of the high-performance computer from Compute Canada
文摘Microwave precondition has been highlighted as a promising technology for softening the rock mass prior to rock breakage by machine to reduce drill bit/cutter wear as well as inverse production rate.To numerically explore the effect of numerical parameters on rock static strength simulation,and determine the numerical mechanical parameters of microwave-treated basalts for future drilling and cutting simulations,numerical models of uniaxial compression strength(UCS)and Brazilian tensile strength(BTS)were established with the coupling of smoothed particle hydrodynamics and finite element method(SPH-FEM).To eliminate the large rock strength errors caused by microwave-induced damage,the cohesion and internal friction angle of microwave-treated basalt specimens with the same microwave treatment parameters were calibrated based on a linear Mohr-Coulomb theory.Based on parametric sensitivity analysis of SPH simulation of UCS and BTS,experimental UCS and BTS values were simultaneously captured according to the same set of calibrated cohesion and internal friction angle data,and the UCS modeling results are in good agreement with experimental tests.Furthermore,the effect of microwave irradiation parameter on the basalt mechanical behaviors was evaluated.
