Boundary integral solutions to three-dimensional unconfined Darcy's Flow

1980 ◽  
Vol 16 (4) ◽  
pp. 651-658 ◽  
Author(s):  
Gerard P. Lennon ◽  
Philip L.-F. Liu ◽  
James A. Liggett
1978 ◽  
Vol 59 (2) ◽  
pp. 245-262 ◽  
Author(s):  
W.L. Meyer ◽  
W.A. Bell ◽  
B.T. Zinn ◽  
M.P. Stallybrass

Algorithms ◽  
2021 ◽  
Vol 14 (5) ◽  
pp. 129
Author(s):  
Yuan Li ◽  
Ni Zhang ◽  
Yuejiao Gong ◽  
Wentao Mao ◽  
Shiguang Zhang

Compared with continuous elements, discontinuous elements advance in processing the discontinuity of physical variables at corner points and discretized models with complex boundaries. However, the computational accuracy of discontinuous elements is sensitive to the positions of element nodes. To reduce the side effect of the node position on the results, this paper proposes employing partially discontinuous elements to compute the time-domain boundary integral equation of 3D elastodynamics. Using the partially discontinuous element, the nodes located at the corner points will be shrunk into the element, whereas the nodes at the non-corner points remain unchanged. As such, a discrete model that is continuous on surfaces and discontinuous between adjacent surfaces can be generated. First, we present a numerical integration scheme of the partially discontinuous element. For the singular integral, an improved element subdivision method is proposed to reduce the side effect of the time step on the integral accuracy. Then, the effectiveness of the proposed method is verified by two numerical examples. Meanwhile, we study the influence of the positions of the nodes on the stability and accuracy of the computation results by cases. Finally, the recommended value range of the inward shrink ratio of the element nodes is provided.


2019 ◽  
Vol 2019 ◽  
pp. 1-14
Author(s):  
Gang Xu ◽  
Guangwei Zhao ◽  
Jing Chen ◽  
Shuqi Wang ◽  
Weichao Shi

The value of the tangential velocity on the Boundary Value Problem (BVP) is inaccurate when comparing the results with analytical solutions by Indirect Boundary Element Method (IBEM), especially at the intersection region where the normal vector is changing rapidly (named nonsmooth boundary). In this study, the singularity of the BVP, which is directly arranged in the center of the surface of the fluid computing domain, is moved outside the computational domain by using the Desingularized Boundary Integral Equation Method (DBIEM). In order to analyze the accuracy of the IBEM/DBIEM and validate the above-mentioned problem, three-dimensional uniform flow over a sphere has been presented. The convergent study of the presented model has been investigated, including desingularized distance in the DBIEM. Then, the numerical results were compared with the analytical solution. It was found that the accuracy of velocity distribution in the flow field has been greatly improved at the intersection region, which has suddenly changed the boundary surface shape of the fluid domain. The conclusions can guide the study on the flow over nonsmooth boundaries by using boundary value method.


2003 ◽  
Vol 70 (2) ◽  
pp. 180-190 ◽  
Author(s):  
E. Pan

In this paper, three-dimensional Green’s functions in anisotropic elastic bimaterials with imperfect interface conditions are derived based on the extended Stroh formalism and the Mindlin’s superposition method. Four different interface models are considered: perfect-bond, smooth-bond, dislocation-like, and force-like. While the first one is for a perfect interface, other three models are for imperfect ones. By introducing certain modified eigenmatrices, it is shown that the bimaterial Green’s functions for the three imperfect interface conditions have mathematically similar concise expressions as those for the perfect-bond interface. That is, the physical-domain bimaterial Green’s functions can be obtained as a sum of a homogeneous full-space Green’s function in an explicit form and a complementary part in terms of simple line-integrals over [0,π] suitable for standard numerical integration. Furthermore, the corresponding two-dimensional bimaterial Green’s functions have been also derived analytically for the three imperfect interface conditions. Based on the bimaterial Green’s functions, the effects of different interface conditions on the displacement and stress fields are discussed. It is shown that only the complementary part of the solution contributes to the difference of the displacement and stress fields due to different interface conditions. Numerical examples are given for the Green’s functions in the bimaterials made of two anisotropic half-spaces. It is observed that different interface conditions can produce substantially different results for some Green’s stress components in the vicinity of the interface, which should be of great interest to the design of interface. Finally, we remark that these bimaterial Green’s functions can be implemented into the boundary integral formulation for the analysis of layered structures where imperfect bond may exist.


2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Sofia Sarraf ◽  
Ezequiel López ◽  
Laura Battaglia ◽  
Gustavo Ríos Rodríguez ◽  
Jorge D'Elía

In the boundary element method (BEM), the Galerkin weighting technique allows to obtain numerical solutions of a boundary integral equation (BIE), giving the Galerkin boundary element method (GBEM). In three-dimensional (3D) spatial domains, the nested double surface integration of GBEM leads to a significantly larger computational time for assembling the linear system than with the standard collocation method. In practice, the computational time is roughly an order of magnitude larger, thus limiting the use of GBEM in 3D engineering problems. The standard approach for reducing the computational time of the linear system assembling is to skip integrations whenever possible. In this work, a modified assembling algorithm for the element matrices in GBEM is proposed for solving integral kernels that depend on the exterior unit normal. This algorithm is based on kernels symmetries at the element level and not on the flow nor in the mesh. It is applied to a BIE that models external creeping flows around 3D closed bodies using second-order kernels, and it is implemented using OpenMP. For these BIEs, the modified algorithm is on average 32% faster than the original one.


Sign in / Sign up

Export Citation Format

Share Document