The scattering of normally incident water waves by two surface-piercing inclined perforated barriers in water with a uniform finite depth is investigated within the framework of linear water wave theory.Considering th...The scattering of normally incident water waves by two surface-piercing inclined perforated barriers in water with a uniform finite depth is investigated within the framework of linear water wave theory.Considering that thin barriers are zero-thickness,a novel numerical method involving the the coupling of the dual boundary element method(DBEM)with damping layers is applied.In order to effectively damp out the reflected waves,two damping layers,instead of pseudoboundaries are implemented near the two side boundaries of the computational domain.Thus,the modified linearized free surface boundary conditions are formulated and used for solving both the ordinary boundary integral equation as well as the hypersingular boundary integral equation for degenerate boundaries.The newly developed numerical method is validated against analytical methods using the matched eigenfunction expansion method for the special case of two vertical barriers or the inclined angle to the vertical being zero.The influence of the length of the two damping layers has been discussed.Moreover,these findings are also validated against previous results for several cases.After validation,the numerical results for the reflection coefficient,transmission coefficient and dissipation coefficient are obtained by varying the inclination angle and porosity-effect parameter.The effects of both the inclination angle and the porosity on the amplitudes of wave forces acting on both the front and rear barriers are also investigated.It is found that the effect of the inclination angle mainly shifts the location of the extremal values of the reflection and the transmission coefficients.Additionally,a moderate value of the porosity-parameter is quite effective at dissipating wave energy and mitigating the wave loads on dual barriers.展开更多
With the aid of the properties of the hypersingular kernels, a geometric conversion approach was presented in this paper. The conversion leads to a general approach for the accurate and reliable numerical evaluation o...With the aid of the properties of the hypersingular kernels, a geometric conversion approach was presented in this paper. The conversion leads to a general approach for the accurate and reliable numerical evaluation of the hypersingular surface boundary integrals encountered in a variety of applications with boundary element method. Based on the conversion, the hypersingularity in the boundary integrals could be lowered by one order, resulting in the simplification of the computer code. Moreover, an integral transformation was introduced to damp out the nearly singular behavior of the kernels by the distance function defined in the local polar coordinate system for the nearly hypersingular case. The approach is simple to use, which can be inserted readily to computer code, thus getting rid of the dull routine deduction of formulae before the numerical implementations, as the expressions of these kernels are in general complicated. The numerical examples were given in three dimensional elasticity, verifying the effectiveness of the proposed approach, which makes it possible to observe numerically the behavior of the boundary integral values with hypersingular kernels across the boundary.展开更多
Based on the biological key-lock-principle common in various biological systems such as the human brain, this paper relates to a method and device for creating problem-solving complexes from individual elements that c...Based on the biological key-lock-principle common in various biological systems such as the human brain, this paper relates to a method and device for creating problem-solving complexes from individual elements that can be coupled with one another and that have different properties to solve problems. The problem solution can be carried out either serially with a large computer, or with several independent, hierarchically joined computers. In this system, an independent control unit that assumes a multitude of tasks and also acts as an interface with access to all participating computers, is assigned to each problem or object class according to the amount of potential problem-oriented solutions. Such a unit prepares the partial solutions found in its computer for the totality of the solutions computed in the associated computers, finally leading to a total problem solution.展开更多
应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝...应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。展开更多
In this paper, responses of a new dual-induction resistivity logging-while-drilling (LWD) tool in 3D inhomogeneous formation models are simulated by the vectorfinite element method (VFEM), the influences of the bo...In this paper, responses of a new dual-induction resistivity logging-while-drilling (LWD) tool in 3D inhomogeneous formation models are simulated by the vectorfinite element method (VFEM), the influences of the borehole, invaded zone, surroundingstrata, and tool eccentricity are analyzed, and calibration loop parameters and calibrationcoefficients of the LWD tool are discussed. The results show that the tool has a greater depthof investigation than that of the existing electromagnetic propagation LWD tools and is moresensitive to azimuthal conductivity. Both deep and medium induction responses have linearrelationships with the formation conductivity, considering optimal calibration loop parametersand calibration coefficients. Due to the different depths of investigation and resolution, deepinduction and medium induction are affected differently by the formation model parameters,thereby having different correction factors. The simulation results can provide theoreticalreferences for the research and interpretation of the dual-induction resistivity LWD tools.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51809209 and 11702244)the Open Fund of Key Laboratory of Water-Sediment Sciences and Water Disaster Prevention of Hunan Province(Grant No.2021SS04).
文摘The scattering of normally incident water waves by two surface-piercing inclined perforated barriers in water with a uniform finite depth is investigated within the framework of linear water wave theory.Considering that thin barriers are zero-thickness,a novel numerical method involving the the coupling of the dual boundary element method(DBEM)with damping layers is applied.In order to effectively damp out the reflected waves,two damping layers,instead of pseudoboundaries are implemented near the two side boundaries of the computational domain.Thus,the modified linearized free surface boundary conditions are formulated and used for solving both the ordinary boundary integral equation as well as the hypersingular boundary integral equation for degenerate boundaries.The newly developed numerical method is validated against analytical methods using the matched eigenfunction expansion method for the special case of two vertical barriers or the inclined angle to the vertical being zero.The influence of the length of the two damping layers has been discussed.Moreover,these findings are also validated against previous results for several cases.After validation,the numerical results for the reflection coefficient,transmission coefficient and dissipation coefficient are obtained by varying the inclination angle and porosity-effect parameter.The effects of both the inclination angle and the porosity on the amplitudes of wave forces acting on both the front and rear barriers are also investigated.It is found that the effect of the inclination angle mainly shifts the location of the extremal values of the reflection and the transmission coefficients.Additionally,a moderate value of the porosity-parameter is quite effective at dissipating wave energy and mitigating the wave loads on dual barriers.
文摘With the aid of the properties of the hypersingular kernels, a geometric conversion approach was presented in this paper. The conversion leads to a general approach for the accurate and reliable numerical evaluation of the hypersingular surface boundary integrals encountered in a variety of applications with boundary element method. Based on the conversion, the hypersingularity in the boundary integrals could be lowered by one order, resulting in the simplification of the computer code. Moreover, an integral transformation was introduced to damp out the nearly singular behavior of the kernels by the distance function defined in the local polar coordinate system for the nearly hypersingular case. The approach is simple to use, which can be inserted readily to computer code, thus getting rid of the dull routine deduction of formulae before the numerical implementations, as the expressions of these kernels are in general complicated. The numerical examples were given in three dimensional elasticity, verifying the effectiveness of the proposed approach, which makes it possible to observe numerically the behavior of the boundary integral values with hypersingular kernels across the boundary.
文摘Based on the biological key-lock-principle common in various biological systems such as the human brain, this paper relates to a method and device for creating problem-solving complexes from individual elements that can be coupled with one another and that have different properties to solve problems. The problem solution can be carried out either serially with a large computer, or with several independent, hierarchically joined computers. In this system, an independent control unit that assumes a multitude of tasks and also acts as an interface with access to all participating computers, is assigned to each problem or object class according to the amount of potential problem-oriented solutions. Such a unit prepares the partial solutions found in its computer for the totality of the solutions computed in the associated computers, finally leading to a total problem solution.
文摘应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。
基金supported by the National Oil and Gas Major Projects(No.2011ZX05020-002)
文摘In this paper, responses of a new dual-induction resistivity logging-while-drilling (LWD) tool in 3D inhomogeneous formation models are simulated by the vectorfinite element method (VFEM), the influences of the borehole, invaded zone, surroundingstrata, and tool eccentricity are analyzed, and calibration loop parameters and calibrationcoefficients of the LWD tool are discussed. The results show that the tool has a greater depthof investigation than that of the existing electromagnetic propagation LWD tools and is moresensitive to azimuthal conductivity. Both deep and medium induction responses have linearrelationships with the formation conductivity, considering optimal calibration loop parametersand calibration coefficients. Due to the different depths of investigation and resolution, deepinduction and medium induction are affected differently by the formation model parameters,thereby having different correction factors. The simulation results can provide theoreticalreferences for the research and interpretation of the dual-induction resistivity LWD tools.