scholarly journals An integral equation–based numerical method for the forced heat equation on complex domains

2020 ◽  
Vol 46 (5) ◽  
Author(s):  
Fredrik Fryklund ◽  
Mary Catherine A. Kropinski ◽  
Anna-Karin Tornberg
Keyword(s):  
Integral Equation ◽  
Heat Equation ◽  
Helmholtz Equation ◽  
Singular Integrals ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Elliptic Pdes ◽  
Time Step ◽  
Modified Helmholtz Equation ◽  
High Regularity

Abstract Integral equation–based numerical methods are directly applicable to homogeneous elliptic PDEs and offer the ability to solve these with high accuracy and speed on complex domains. In this paper, such a method is extended to the heat equation with inhomogeneous source terms. First, the heat equation is discretised in time, then in each time step we solve a sequence of so-called modified Helmholtz equations with a parameter depending on the time step size. The modified Helmholtz equation is then split into two: a homogeneous part solved with a boundary integral method and a particular part, where the solution is obtained by evaluating a volume potential over the inhomogeneous source term over a simple domain. In this work, we introduce two components which are critical for the success of this approach: a method to efficiently compute a high-regularity extension of a function outside the domain where it is defined, and a special quadrature method to accurately evaluate singular and nearly singular integrals in the integral formulation of the modified Helmholtz equation for all time step sizes.

scholarly journals Convergence of the High-Accuracy Algorithm for Solving the Dirichlet Problem of the Modified Helmholtz Equation

Complexity ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-8
Author(s):  
Hu Li ◽  
Guang Zeng
Keyword(s):  
Dirichlet Problem ◽  
Helmholtz Equation ◽  
Singular Integrals ◽  
Quadrature Rule ◽  
High Accuracy ◽  
Boundary Integral ◽  
Modified Helmholtz Equation ◽  
Weakly Singular ◽  
Asymptotic Error ◽  
Error Expansion

In this paper, we derive the convergence for the high-accuracy algorithm in solving the Dirichlet problem of the modified Helmholtz equation. By the boundary element method, we transform the system to be a boundary integral equation. The high-accuracy algorithm using the specific quadrature rule is developed to deal with weakly singular integrals. The convergence of the algorithm is proved based on Anselone’s collective compact theory. Moreover, an asymptotic error expansion shows that the algorithm is of order Oh03. The numerical examples support the theoretical analysis.

scholarly journals Nonlinear two-dimensional free surface solutions of flow exiting a pipe and impacting a wedge

2021 ◽  
Vol 126 (1) ◽  
Author(s):  
Alex Doak ◽  
Jean-Marc Vanden-Broeck
Keyword(s):  
Surface Tension ◽  
Integral Equation ◽  
Free Surface ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Two Dimensional ◽  
Numerical Schemes ◽  
Additional Advantage ◽  
Infinite Wedge ◽  
Good Agreement

AbstractThis paper concerns the flow of fluid exiting a two-dimensional pipe and impacting an infinite wedge. Where the flow leaves the pipe there is a free surface between the fluid and a passive gas. The model is a generalisation of both plane bubbles and flow impacting a flat plate. In the absence of gravity and surface tension, an exact free streamline solution is derived. We also construct two numerical schemes to compute solutions with the inclusion of surface tension and gravity. The first method involves mapping the flow to the lower half-plane, where an integral equation concerning only boundary values is derived. This integral equation is solved numerically. The second method involves conformally mapping the flow domain onto a unit disc in the s-plane. The unknowns are then expressed as a power series in s. The series is truncated, and the coefficients are solved numerically. The boundary integral method has the additional advantage that it allows for solutions with waves in the far-field, as discussed later. Good agreement between the two numerical methods and the exact free streamline solution provides a check on the numerical schemes.

Analysis of Crack Propagation in Roller Bearings Using the Boundary Integral Equation Method—A Mixed-Mode Loading Problem

1988 ◽  
Vol 110 (3) ◽  
pp. 408-413 ◽  
Author(s):  
L. J. Ghosn
Keyword(s):  
Integral Equation ◽  
Stress Intensity ◽  
Crack Propagation ◽  
Boundary Integral Equation ◽  
Stress Intensity Factors ◽  
High Speed ◽  
Roller Bearing ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Intensity Factors

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

On modified Green functions in exterior problems for the Helmholtz equation

1982 ◽  
Vol 383 (1785) ◽  
pp. 313-332 ◽  
Keyword(s):  
Integral Equation ◽  
Green Function ◽  
Helmholtz Equation ◽  
Least Squares ◽  
Present Note ◽  
Boundary Integral ◽  
Green Functions ◽  
Exterior Problems ◽  
The Green Function ◽  
Definition Of

The question of non-uniqueness in boundary integral equation formu­lations of exterior problems for the Helmholtz equation has recently been resolved with the use of additional radiating multipoles in the definition of the Green function. The present note shows how this modification may be included in a rigorous formalism and presents an explicit choice of co­efficients of the added terms that is optimal in the sense of minimizing the least-squares difference between the modified and exact Green functions.

The boundary integral equation for curved solid/liquid interfaces propagating into a binary liquid with convection

2021 ◽  
Author(s):  
Ekaterina Titova ◽  
Dmitri Alexandrov
Keyword(s):  
Integral Equation ◽  
Boundary Integral Equation ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Curvilinear Coordinates ◽  
Forced Flow ◽  
Function Technique ◽  
Liquid Interfaces ◽  
Convective Boundary ◽  
Solid Liquid

Abstract The boundary integral method is developed for unsteady solid/liquid interfaces propagating into undercooled binary liquids with convection. A single integrodifferential equation for the interface function is derived using the Green function technique. In the limiting cases, the obtained unsteady convective boundary integral equation (CBIE) transforms into a previously developed theory. This integral is simplified for the steady-state growth in arbitrary curvilinear coordinates when the solid/liquid interface is isothermal (isoconcentration). Finally, we evaluate the boundary integral for a binary melt with a forced flow and analyze how the melt undercooling depends on P\'eclet and Reynolds numbers.

scholarly journals A boundary integral method with volume-changing objects for ultrasound-triggered margination of microbubbles

2017 ◽  
Vol 836 ◽  
pp. 952-997 ◽  
Author(s):  
Achim Guckenberger ◽  
Stephan Gekle
Keyword(s):  
Drug Delivery ◽  
Integral Equation ◽  
Boundary Integral Equation ◽  
Integral Method ◽  
Three Dimensional ◽  
Formal Proof ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Drug Uptake ◽  
Periodic Domains

A variety of numerical methods exist for the study of deformable particles in dense suspensions. None of the standard tools, however, currently include volume-changing objects such as oscillating microbubbles in three-dimensional periodic domains. In the first part of this work, we develop a novel method to include such entities based on the boundary integral method. We show that the well-known boundary integral equation must be amended with two additional terms containing the volume flux through the bubble surface. We rigorously prove the existence and uniqueness of the solution. Our proof contains as a subset the simpler boundary integral equation without volume-changing objects (such as red blood cell or capsule suspensions) which is widely used but for which a formal proof in periodic domains has not been published to date. In the second part, we apply our method to study microbubbles for targeted drug delivery. The ideal drug delivery agent should stay away from the biochemically active vessel walls during circulation. However, upon reaching its target it should attain a near-wall position for efficient drug uptake. Though seemingly contradictory, we show that lipid-coated microbubbles in conjunction with a localized ultrasound pulse possess precisely these two properties. This ultrasound-triggered margination is due to hydrodynamic interactions between the red blood cells and the oscillating lipid-coated microbubbles which alternate between a soft and a stiff state. We find that the effect is very robust, existing even if the duration in the stiff state is more than three times lower than the opposing time in the soft state.

Computations of Arbitrary Thin-Shell Vertical Cylinders in a Wave Field by a Hypersingular Integral-Equation Method

2014 ◽  
Author(s):  
Mohamed Hariri Nokob ◽  
Ronald W. Yeung
Keyword(s):  
Integral Equation ◽  
Thin Shell ◽  
Boundary Integral ◽  
The Body ◽  
Boundary Integral Method ◽  
Hypersingular Integral Equation ◽  
Wave Load ◽  
Hypersingular Integral ◽  
Method Of Solution ◽  
Vertical Cylinders

We present a numerical formulation and computational results for the hydrodynamic loads on bottom-mounted thin-shell vertical cylinders of arbitrary cross-sectional shapes, including open bodies. Such cylinders may undergo prescribed or free motion or may be subjected to a wave load. The formulation is based on linear theory and a hypersingular integral-equation results from sectional contours of zero thickness. The method reduces the fully three-dimensional problem into a number of two-dimensional ones in the horizontal plane and is therefore much faster than the usual boundary integral method used for water wave problems. This traditional method of solution is also known to become ill-conditioned as the body thickness decreases. As an example of the current method, radiation and diffraction loads are presented for the cases of a circular and square closed and opened cylinders with the effect of the increased opening size discussed at different frequencies. The free-surface elevation associated with this method of solution is presented for some cases as well.

scholarly journals On the Harmonic Problem with Nonlinear Boundary Integral Conditions

2014 ◽  
Vol 2014 ◽  
pp. 1-5
Author(s):  
Saker Hacene
Keyword(s):  
Integral Equation ◽  
Existence And Uniqueness ◽  
Integral Method ◽  
Nonlinear Boundary ◽  
Monotone Operators ◽  
Boundary Integral ◽  
Boundary Integral Method ◽  
Integral Conditions ◽  
Uniqueness Of The Solution ◽  
Harmonic Problem

In the present work, we deal with the harmonic problems in a bounded domain of ℝ2 with the nonlinear boundary integral conditions. After applying the Boundary integral method, a nonlinear boundary integral equation is obtained; the existence and uniqueness of the solution will be a consequence of applying theory of monotone operators.

Sign in / Sign up

Export Citation Format

Share Document