A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity poten...A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
Coastal wave energy resources have enormous exploitation potential due to shorter weather window,closer installation distance and lower maintenance cost.However,impact loads generated by depth variation from offshore ...Coastal wave energy resources have enormous exploitation potential due to shorter weather window,closer installation distance and lower maintenance cost.However,impact loads generated by depth variation from offshore to nearshore and wave-current interaction,may lead to a catastrophic damage or complete destruction to wave energy converters(WECs).This objective of this paper is to investigate slamming response of a coastal oscillating wave surge converter(OWSC)entering or leaving water freely.Based on fully nonlinear potential flow theory,a time-domain wave-current-structure interaction model combined with higher-order boundary element method(HOBEM),is developed to analyze the coupled hydrodynamic problem.The variable-depth seabed is considered in the model to illustrate the shallow water effect on impact loads and free surface profiles in coastal zone.A domain decomposition approach is utilized to simulate the overlapping phenomenon generated by a jet falling into water under gravity effect.Through a series of Lagrangian interpolation methods,the meshes on boundaries are rearranged to avoid the mismatch between element size on free surface and body surface.The present model is validated against the existing experimental and numerical results.Simulations are also provided for the effects of wave-current interaction and uneven local seabed on the slamming responses.It is found that the length of the splash jet increases for a following current and decreases for an opposing current,and that the slamming response of the OWSC device is sensitive to the geometric features of the uneven seabed.展开更多
In this paper,we discuss an algebraic multigrid(AMG)method for nearly incompressible elasticity problems in two-dimensions.First,a two-level method is proposed by analyzing the relationship between the linear finite e...In this paper,we discuss an algebraic multigrid(AMG)method for nearly incompressible elasticity problems in two-dimensions.First,a two-level method is proposed by analyzing the relationship between the linear finite element space and the quartic finite element space.By choosing different smoothers,we obtain two types of two-level methods,namely TL-GS and TL-BGS.The theoretical analysis and numerical results show that the convergence rates of TL-GS and TL-BGS are independent of the mesh size and the Young’s modulus,and the convergence of the latter is greatly improved on the order p.However the convergence of both methods still depends on the Poisson’s ratio.To fix this,we obtain a coarse level matrix with less rigidity based on selective reduced integration(SRI)method and get some types of two-level methods by combining different smoothers.With the existing AMG method used as a solver on the first coarse level,an AMG method can be finally obtained.Numerical results show that the resulting AMG method has better efficiency for nearly incompressible elasticity problems.展开更多
The present study concerns the modelization and numerical simulation for the heat and flow exchange characteristics in a novel configuration saturated with a nonNewtonian Ag-MgO hybrid nanofluid.The wavy shaped enclos...The present study concerns the modelization and numerical simulation for the heat and flow exchange characteristics in a novel configuration saturated with a nonNewtonian Ag-MgO hybrid nanofluid.The wavy shaped enclosure is equipped with onequarter of a conducting solid cylinder.The system of equations resulting from the mathematical modeling of the physical problem in its dimensionless form is discretized via the higher-order Galerkin-based finite element method(GFEM).The dependency of various factors and their interrelationships affecting the hydro-thermal behavior and heat exchange rate are delineated.The numerical experiments reveal that the best heat transfer rate is achieved for the pseudo-plastic hybrid nanoliquid with high Rayleigh number and thermal conductivity ratio and low Hartmann number.Besides,the power-law index has a major effect in deteriorating the heat convection at high Rayleigh number.展开更多
A novel,highly efficient and accurate adaptive higher-order finite element method(hp-FEM)is used to simulate a multi-frequency resistivity loggingwhile-drilling(LWD)tool response in a borehole environment.Presented in...A novel,highly efficient and accurate adaptive higher-order finite element method(hp-FEM)is used to simulate a multi-frequency resistivity loggingwhile-drilling(LWD)tool response in a borehole environment.Presented in this study are the vector expression of Maxwell’s equations,three kinds of boundary conditions,stability weak formulation of Maxwell’s equations,and automatic hpadaptivity strategy.The new hp-FEM can select optimal refinement and calculation strategies based on the practical formation model and error estimation.Numerical experiments show that the new hp-FEM has an exponential convergence rate in terms of relative error in a user-prescribed quantity of interest against the degrees of freedom,which provides more accurate results than those obtained using the adaptive h-FEM.The numerical results illustrate the high efficiency and accuracy of the method at a given LWD tool structure and parameters in different physical models,which further confirm the accuracy of the results using the Hermes library(http://hpfem.org/hermes)with a multi-frequency resistivity LWD tool response in a borehole environment.展开更多
Apertures generally exist in the sandwich structures attributing to mechanical connection and lightweight, which might induce failure of such structures. Thus, it is required to realize the impact of aperture on mecha...Apertures generally exist in the sandwich structures attributing to mechanical connection and lightweight, which might induce failure of such structures. Thus, it is required to realize the impact of aperture on mechanical behaviors of sandwich structures. If transverse shear deformations are unable to be described accurately, the reasonable prediction of dynamic behaviors of the form-core sandwich plates with apertures will meet severe challenges due to a large difference of transverse shear modulus at the adjacent layers. Thereby, such issue is less studied by using the efficient models and experimental testing, so an alternative sinusoidal-type finite element formulation is to be proposed to precisely predict dynamic response of the form-core sandwich structures with apertures. The proposed finite element formulation can meet beforehand compatible conditions of transverse shear stresses at the interfaces of adjacent laminates. In order to appraise strictly capability of the proposed model, experimental tests on natural frequencies of three groups of specimens with different apertures have been carried out. Moreover, four specimens in each group are tested to reduce the testing errors, which is less reported in the published literature. In addition,three-dimensional Finite Element Method(3-D FEM) is also selected to account for the good performance of the present model. Finally, the impact of aperture diameter on the natural frequencies of the sandwich structures is both experimentally and numerically investigated, which can serve as a reference for other researchers.展开更多
The body-fixed coordinate system is applied to the wave-body interaction problem of a small-depth elastic structure which has both rigid and elastic body motions in head waves.In the weakly non-linear assumption,the p...The body-fixed coordinate system is applied to the wave-body interaction problem of a small-depth elastic structure which has both rigid and elastic body motions in head waves.In the weakly non-linear assumption,the perturbation scheme is used and the expansion is conducted up to second-order to consider several non-linear quantities.To solve the boundary value problem,linearization is carried out based not on inertial coordinate but on body-fixed coordinate which could be accelerated by a motion of a body.At first,the main feature of the application of body-fixed coordinate system for a seakeeping problem is briefly described.After that the transformation of a coordinate system is extended to consider an elastic body motion and several physical variables are re-described in the generalized mode.It has been found that the deformation gradient could be used for the transformation of a coordinate system if several conditions are satisfied.Provided there are only vertical bending in elastic modes and the structure has relatively small depth,these conditions are generally satisfied.To calculate an elastic motion of a body,the generalized mode method is adopted and the mode shape is obtained by solving eigen-value problem of dynamic beam equation.In the boundary condition of the body-fixed coordinate system,the motion effect reflected to free-surface boundary is considered by extrapolating each mode shape to the horizontal direction from a body.At last,simple numerical tests are implemented as a validation process.The second-order hydrodynamic force of a freely floating hemisphere is first calculated in zero forward speed condition.Next,motion and added resistance of a ship with forward speed are considered at different flexibility to confirm the effect of an elastic body motion in body-fixed coordinate system.展开更多
基金financially supported by the National Natural Science Foundation of China (Grant Nos.52271276,52271319,and 52201364)the Natural Science Foundation of Jiangsu Province (Grant No.BK20201006)。
文摘A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52025112 and 51861130358)the State Key Laboratory of Ocean Engineering+1 种基金China(Shanghai Jiao Tong University)(Grant No.1905)the Newton Advanced Fellowships(Grant No.NAF\R1\180304)by the Royal Society。
文摘Coastal wave energy resources have enormous exploitation potential due to shorter weather window,closer installation distance and lower maintenance cost.However,impact loads generated by depth variation from offshore to nearshore and wave-current interaction,may lead to a catastrophic damage or complete destruction to wave energy converters(WECs).This objective of this paper is to investigate slamming response of a coastal oscillating wave surge converter(OWSC)entering or leaving water freely.Based on fully nonlinear potential flow theory,a time-domain wave-current-structure interaction model combined with higher-order boundary element method(HOBEM),is developed to analyze the coupled hydrodynamic problem.The variable-depth seabed is considered in the model to illustrate the shallow water effect on impact loads and free surface profiles in coastal zone.A domain decomposition approach is utilized to simulate the overlapping phenomenon generated by a jet falling into water under gravity effect.Through a series of Lagrangian interpolation methods,the meshes on boundaries are rearranged to avoid the mismatch between element size on free surface and body surface.The present model is validated against the existing experimental and numerical results.Simulations are also provided for the effects of wave-current interaction and uneven local seabed on the slamming responses.It is found that the length of the splash jet increases for a following current and decreases for an opposing current,and that the slamming response of the OWSC device is sensitive to the geometric features of the uneven seabed.
基金supported in part by NSF-10771178 and NSF-10672138 in Chinathe Basic Research Program of China under the grant 2005CB321702+1 种基金the Key Project of Chinese Ministry of Education and the Scientific Research Fund of Hunan Provincial Education Department(208093,07A068)the Provincial Natural Science Foundation of Hunan(07JJ6004)。
文摘In this paper,we discuss an algebraic multigrid(AMG)method for nearly incompressible elasticity problems in two-dimensions.First,a two-level method is proposed by analyzing the relationship between the linear finite element space and the quartic finite element space.By choosing different smoothers,we obtain two types of two-level methods,namely TL-GS and TL-BGS.The theoretical analysis and numerical results show that the convergence rates of TL-GS and TL-BGS are independent of the mesh size and the Young’s modulus,and the convergence of the latter is greatly improved on the order p.However the convergence of both methods still depends on the Poisson’s ratio.To fix this,we obtain a coarse level matrix with less rigidity based on selective reduced integration(SRI)method and get some types of two-level methods by combining different smoothers.With the existing AMG method used as a solver on the first coarse level,an AMG method can be finally obtained.Numerical results show that the resulting AMG method has better efficiency for nearly incompressible elasticity problems.
文摘The present study concerns the modelization and numerical simulation for the heat and flow exchange characteristics in a novel configuration saturated with a nonNewtonian Ag-MgO hybrid nanofluid.The wavy shaped enclosure is equipped with onequarter of a conducting solid cylinder.The system of equations resulting from the mathematical modeling of the physical problem in its dimensionless form is discretized via the higher-order Galerkin-based finite element method(GFEM).The dependency of various factors and their interrelationships affecting the hydro-thermal behavior and heat exchange rate are delineated.The numerical experiments reveal that the best heat transfer rate is achieved for the pseudo-plastic hybrid nanoliquid with high Rayleigh number and thermal conductivity ratio and low Hartmann number.Besides,the power-law index has a major effect in deteriorating the heat convection at high Rayleigh number.
基金The work for this paper was supported by the National Natural Science Foundation of China under Projects No.41074099。
文摘A novel,highly efficient and accurate adaptive higher-order finite element method(hp-FEM)is used to simulate a multi-frequency resistivity loggingwhile-drilling(LWD)tool response in a borehole environment.Presented in this study are the vector expression of Maxwell’s equations,three kinds of boundary conditions,stability weak formulation of Maxwell’s equations,and automatic hpadaptivity strategy.The new hp-FEM can select optimal refinement and calculation strategies based on the practical formation model and error estimation.Numerical experiments show that the new hp-FEM has an exponential convergence rate in terms of relative error in a user-prescribed quantity of interest against the degrees of freedom,which provides more accurate results than those obtained using the adaptive h-FEM.The numerical results illustrate the high efficiency and accuracy of the method at a given LWD tool structure and parameters in different physical models,which further confirm the accuracy of the results using the Hermes library(http://hpfem.org/hermes)with a multi-frequency resistivity LWD tool response in a borehole environment.
基金supported by SKLLIM1902the Natural Science Foundation in Shaanxi Province,China(No.2019JQ-909)。
文摘Apertures generally exist in the sandwich structures attributing to mechanical connection and lightweight, which might induce failure of such structures. Thus, it is required to realize the impact of aperture on mechanical behaviors of sandwich structures. If transverse shear deformations are unable to be described accurately, the reasonable prediction of dynamic behaviors of the form-core sandwich plates with apertures will meet severe challenges due to a large difference of transverse shear modulus at the adjacent layers. Thereby, such issue is less studied by using the efficient models and experimental testing, so an alternative sinusoidal-type finite element formulation is to be proposed to precisely predict dynamic response of the form-core sandwich structures with apertures. The proposed finite element formulation can meet beforehand compatible conditions of transverse shear stresses at the interfaces of adjacent laminates. In order to appraise strictly capability of the proposed model, experimental tests on natural frequencies of three groups of specimens with different apertures have been carried out. Moreover, four specimens in each group are tested to reduce the testing errors, which is less reported in the published literature. In addition,three-dimensional Finite Element Method(3-D FEM) is also selected to account for the good performance of the present model. Finally, the impact of aperture diameter on the natural frequencies of the sandwich structures is both experimentally and numerically investigated, which can serve as a reference for other researchers.
文摘The body-fixed coordinate system is applied to the wave-body interaction problem of a small-depth elastic structure which has both rigid and elastic body motions in head waves.In the weakly non-linear assumption,the perturbation scheme is used and the expansion is conducted up to second-order to consider several non-linear quantities.To solve the boundary value problem,linearization is carried out based not on inertial coordinate but on body-fixed coordinate which could be accelerated by a motion of a body.At first,the main feature of the application of body-fixed coordinate system for a seakeeping problem is briefly described.After that the transformation of a coordinate system is extended to consider an elastic body motion and several physical variables are re-described in the generalized mode.It has been found that the deformation gradient could be used for the transformation of a coordinate system if several conditions are satisfied.Provided there are only vertical bending in elastic modes and the structure has relatively small depth,these conditions are generally satisfied.To calculate an elastic motion of a body,the generalized mode method is adopted and the mode shape is obtained by solving eigen-value problem of dynamic beam equation.In the boundary condition of the body-fixed coordinate system,the motion effect reflected to free-surface boundary is considered by extrapolating each mode shape to the horizontal direction from a body.At last,simple numerical tests are implemented as a validation process.The second-order hydrodynamic force of a freely floating hemisphere is first calculated in zero forward speed condition.Next,motion and added resistance of a ship with forward speed are considered at different flexibility to confirm the effect of an elastic body motion in body-fixed coordinate system.
