scholarly journals The Numerical Analysis of the Flow on the Smooth and Nonsmooth Boundaries by IBEM/DBIEM

2019 ◽  
Vol 2019 ◽  
pp. 1-14
Author(s):  
Gang Xu ◽  
Guangwei Zhao ◽  
Jing Chen ◽  
Shuqi Wang ◽  
Weichao Shi
Keyword(s):  
Boundary Value ◽  
Three Dimensional ◽  
Boundary Surface ◽  
Tangential Velocity ◽  
Integral Equation Method ◽  
Boundary Integral ◽  
Surface Shape ◽  
Normal Vector ◽  
Computational Domain ◽  
Nonsmooth Boundary

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.

On mixed boundary value problems for the Helmholtz equation

1977 ◽  
Vol 77 (1-2) ◽  
pp. 65-77 ◽  
Author(s):  
R. Kress ◽  
G. F. Roach
Keyword(s):  
Helmholtz Equation ◽  
Boundary Value Problems ◽  
Boundary Integral Equations ◽  
Mixed Boundary ◽  
Boundary Value ◽  
Three Dimensional ◽  
Integral Equation Method ◽  
Boundary Integral ◽  
Theoretic Approach ◽  
Mixed Boundary Value Problems

SynopsisExistence and uniqueness theorems are obtained for a class of mixed boundary value problems associated with the three-dimensional Helmholtz equation. In this context the boundary of the region of interest is assumed to consist of the union of a finite number of disjoint, closed, bounded Lyapunov surfaces on some of which are imposed Dirichlet conditions whilst Neumann conditions are imposed on the remainder. An integral equation method is adopted throughout. The required boundary integral equations are generated by a modified layer theoretic approach which extends the work of Brakhage and Werner [1] and Leis [2, 3].

scholarly journals A unified boundary integral equation method for a class of second order elliptic boundary value problems

1984 ◽  
Vol 25 (4) ◽  
pp. 501-517 ◽  
Author(s):  
M. Rezayat ◽  
F. J. Rizzo ◽  
D. J. Shippy
Keyword(s):  
Integral Equation ◽  
Boundary Value Problems ◽  
Elliptic Boundary ◽  
Boundary Value ◽  
Three Dimensional ◽  
Integral Equation Method ◽  
Boundary Integral ◽  
Second Order ◽  
Elliptic Type ◽  
Integral Equation Formulation

AbstractA generalized integral equation formulation and a systematic numerical solution procedure are presented for a class of boundary value problems governed by a general second-order differential equation of elliptic type. Diverse numerical examples include problems of plane-wave scattering, three-dimensional fluid flow, and plane heat transfer for a body with a moving flame boundary. The last example employs certain representation functions useful to increase solution effectiveness in problems with an isolated integrable singularity.

Rigorous accounting diffraction on non-plane gratings irradiated by non-planar waves

2021 ◽  
Author(s):  
Leonid I. Goray
Keyword(s):  
Plane Wave ◽  
Wave Diffraction ◽  
Plane Waves ◽  
Three Dimensional ◽  
Integral Equation Method ◽  
Transverse Magnetic ◽  
Incident Field ◽  
Boundary Integral ◽  
Cylindrical Waves ◽  
Plane Wave Diffraction

Abstract The modified boundary integral equation method (MIM) is considered a rigorous theoretical application for the diffraction of cylindrical waves by arbitrary profiled plane gratings, as well as for the diffraction of plane/non-planar waves by concave/convex gratings. This study investigates two-dimensional (2D) diffraction problems of the filiform source electromagnetic field scattered by a plane lamellar grating and of plane waves scattered by a similar cylindrical-shaped grating. Unlike the problem of plane wave diffraction by a plane grating, the field of a localised source does not satisfy the quasi-periodicity requirement. Fourier transform is used to reduce the solution of the problem of localised source diffraction by the grating in the whole region to the solution of the problem of diffraction inside one Floquet channel. By considering the periodicity of the geometry structure, the problem of Floquet terms for the image can be formulated so that it enables the application of the MIM developed for plane wave diffraction problems. Accounting of the local structure of an incident field enables both the prediction of the corresponding efficiencies and the specification of the bounds within which the approximation of the incident field with plane waves is correct. For 2D diffraction problems of the high-conductive plane grating irradiated by cylindrical waves and the cylindrical high-conductive grating irradiated by plane waves, decompositions in sets of plane waves/sections are investigated. The application of such decomposition, including the dependence on the number of plane waves/sections and radii of the grating and wave front shape, was demonstrated for lamellar, sinusoidal and saw-tooth grating examples in the 0th & –1st orders as well as in the transverse electric and transverse magnetic polarisations. The primary effects of plane wave/section partitions of non-planar wave fronts and curved grating shapes on the exact solutions for 2D and three-dimensional (conical) diffraction problems are discussed.

Numerical simulations of nonlinear axisymmetric flows with a free surface

1987 ◽  
Vol 178 ◽  
pp. 195-219 ◽  
Author(s):  
Douglas G. Dommermuth ◽  
Dick K. P. Yue
Keyword(s):  
Free Surface ◽  
Numerical Simulations ◽  
Three Dimensional ◽  
Vertical Cylinder ◽  
Free Surface Flows ◽  
Boundary Integral ◽  
The Body ◽  
Special Treatment ◽  
Computational Domain ◽  
Vortex Sheets

A numerical method is developed for nonlinear three-dimensional but axisymmetric free-surface problems using a mixed Eulerian-Lagrangian scheme under the assumption of potential flow. Taking advantage of axisymmetry, Rankine ring sources are used in a Green's theorem boundary-integral formulation to solve the field equation; and the free surface is then updated in time following Lagrangian points. A special treatment of the free surface and body intersection points is generalized to this case which avoids the difficulties associated with the singularity there. To allow for long-time simulations, the nonlinear computational domain is matched to a transient linear wavefield outside. When the matching boundary is placed at a suitable distance (depending on wave amplitude), numerical simulations can, in principle, be continued indefinitely in time. Based on a simple stability argument, a regriding algorithm similar to that of Fink & Soh (1974) for vortex sheets is generalized to free-surface flows, which removes the instabilities experienced by earlier investigators and eliminates the need for artificial smoothing. The resulting scheme is very robust and stable.For illustration, three computational examples are presented: (i) the growth and collapse of a vapour cavity near the free surface; (ii) the heaving of a floating vertical cylinder starting from rest; and (iii) the heaving of an inverted vertical cone. For the cavity problem, there is excellent agreement with available experiments. For the wave-body interaction calculations, we are able to obtain and analyse steady-state (limit-cycle) results for the force and flow field in the vicinity of the body.

scholarly journals NUMERICAL ANALYSIS OF THE DYNAMICS OF THREE-DIMENSIONAL ANISOTROPIC BODIES BASED ON NON-CLASSICAL BOUNDARY INTEGRAL EQUATIONS

2021 ◽  
Vol 83 (1) ◽  
pp. 76-86
Author(s):  
A.A. Belov ◽  
A.N. Petrov
Keyword(s):  
Integral Equation ◽  
Boundary Element ◽  
Integral Equations ◽  
Boundary Integral Equations ◽  
Three Dimensional ◽  
Integral Equation Method ◽  
Boundary Integral ◽  
Classical Approach ◽  
Classical Boundary ◽  
Nonclassical Formulation

The application of non-classical approach of the boundary integral equation method in combination with the integral Laplace transform in time to anisotropic elastic wave modeling is considered. In contrast to the classical approach of the boundary integral equation method which is successfully implemented for solving three-dimensional isotropic problems of the dynamic theory of elasticity, viscoelasticity and poroelasticity, the alternative nonclassical formulation of the boundary integral equations method is presented that employs regular Fredholm integral equations of the first kind (integral equations on a plane wave). The construction of such boundary integral equations is based on the structure of the dynamic fundamental solution. The approach employs the explicit boundary integral equations. The inverse Laplace transform is constructed numerically by the Durbin method. A numerical solution of the dynamic problem of anisotropic elasticity theory based on the boundary integral equations method in a nonclassical formulation is presented. The boundary element scheme of the boundary integral equations method is built on the basis of a regular integral equation of the first kind. The problem is solved in anisotropic formulation for the load acting along the normal in the form of the Heaviside function on the cube face weakened by a cubic cavity. The obtained boundary element solutions are compared with finite element solutions. Numerical results prove the efficiency of using boundary integral equations on a single plane wave in solving three-dimensional anisotropic dynamic problems of elasticity theory. The convergence of boundary element solutions is studied on three schemes of surface discretization. The achieved calculation accuracy is not inferior to the accuracy of boundary element schemes for classical boundary integral equations. Boundary element analysis of solutions for a cube with and without a cavity is carried out.

Sign in / Sign up

Export Citation Format

Share Document