Search

A generalized boundary integral equation method for boundary value problems of two-dimensional isotropic lattice Laplacian is proposed in this paper. The proposed method is an extension of the classical boundary integral equation method with notable advantage. By utilizing the asymptotic expression of the fundamental solution at infinity, this method effectively addresses the challenge of numerical integration involving singular integral kernels. The introduction of Green’s formulas, Dirichlet and Neumann traces, and other tools which are parallel to the traditional integral equation method, form a solid foundation for the development of the generalized boundary integral equation method. The solvability of boundary integral equations and the solvability of lattice interface problem are important guarantees for the feasibility of this method, and these are emphasized in this paper. Subsequently, the generalized boundary integral equation method is applied to boundary value problems equipped with either Dirichlet or Neumann boundary conditions. Simple numerical examples demonstrate the accuracy and effectiveness of the generalized boundary integral equation method. [ABSTRACT FROM AUTHOR]

The paper is concerned with using boundary element methods (BEM) for the accurate evaluation of fields in structures with clusters. For large-scale problems, the BEM system is solved iteratively by speeding up matrix-to-vector multiplications by applying a kernel-independent fast multipole method. Multiplication starts with source-to-multipole (S2M) translations, whose accuracy predefines the overall accuracy. We aim to increase the accuracy of these translations. The intensities of sources are assembled into clusters by an algorithm suggested. Each of them is characterized by its representative source, whose intensity equals the sum of the intensities of cluster sources. Thus, with growing distance, its field tends toward the field of the cluster. The accuracy of S2M translations is increased by subtracting from and adding to the far field of the cluster the far field of its representative source, and by using the proposed modified kernel to evaluate the difference of the fields, which decreases faster than the field of the cluster itself. Numerical results for typical kernels show a notable increase in the accuracy provided by the modified S2M translations. Keeping in them merely the added field is acceptable for many practical applications. This simplifies the modified S2M translations by avoiding calculation and storing matrices specific to each of the clusters. The improved translations may be also used for multipole-to-multipole translations, performed on next, after leaves, levels in upward running a hierarchical tree. [ABSTRACT FROM AUTHOR]

The Navier equation is the governing equation of elastic waves, and computing its solution accurately and rapidly has a wide range of applications in geophysical exploration, materials science, etc. In this paper, we focus on the efficient and high-precision numerical algorithm for the time harmonic elastic wave scattering problems from cornered domains via the boundary integral equations in two dimensions. The approach is based on the combination of Nyström discretization, analytical singular integrals and kernel-splitting method, which results in a high-order solver for smooth boundaries. It is then combined with the recursively compressed inverse preconditioning (RCIP) method to solve elastic scattering problems from cornered domains. Numerical experiments demonstrate that the proposed approach achieves high accuracy, with stabilized errors close to machine precision in various geometric configurations. The algorithm is further applied to investigate the asymptotic behavior of density functions associated with boundary integral operators near corners, and the numerical results are highly consistent with the theoretical formulas. [ABSTRACT FROM AUTHOR]

Accurate evaluation of nearly singular integrals plays an important role in many boundary integral equation based numerical methods. In this paper, we propose a variant of singularity swapping method to accurately evaluate the layer potentials for arbitrarily close targets. Our method is based on the global trapezoidal rule and trigonometric interpolation, resulting in an explicit quadrature formula. The method achieves spectral accuracy for nearly singular integrals on closed analytic curves. In order to extract the singularity from the complexified distance function, an efficient root finding method is proposed based on contour integration. Through the change of variables, we also extend the quadrature method to integrals on the piecewise analytic curves. Numerical examples for Laplace and Helmholtz equations show that high-order accuracy can be achieved for arbitrarily close field evaluation. [ABSTRACT FROM AUTHOR]

This paper is concerned with the cavity scattering problem in an infinite thin plate, where the out-of-plane displacement is governed by the two-dimensional biharmonic wave equation. Based on an operator splitting, the scattering problem is recast into a coupled boundary value problem for the Helmholtz and modified Helmholtz equations. A novel boundary integral formulation is proposed for the coupled problem. By introducing an appropriate regularizer, the well-posedness is established for the system of boundary integral equations. Moreover, the convergence analysis is carried out for the semi- and full-discrete schemes of the boundary integral system by using the collocation method. Numerical results show that the proposed method is highly accurate for both smooth and nonsmooth examples. [ABSTRACT FROM AUTHOR]

In this paper, we find an equivalent boundary integral equation to the classical ∂ ̄ -Neumann problem. The new equation contains an equivalent regularity to the global regularity of the ∂ ̄ -Neumann problem. We also use the integral equations to observe a new-found rigidity property of the ∂ ̄ -Neumann problem on the unit ball in ℂ 2 . [ABSTRACT FROM AUTHOR]

This paper describes a locally corrected trapezoidal quadrature method for the discretization of singular and hypersingular boundary integral operators (BIOs) that arise in solving boundary value problems for elliptic partial differential equations. The quadrature is based on a uniform grid in parameter space coupled with the standard punctured trapezoidal rule. A key observation is that the error incurred by the singularity in the kernel can be expressed exactly using generalized Euler--Maclaurin formulas that involve the Riemann zeta function in 2 dimensions (2D) and the Epstein zeta functions in 3 dimensions (3D). These expansions are exploited to correct the errors via local stencils at the singular point using a novel systematic moment-fitting approach. This new method provides a unified treatment of all common BIOs (Laplace, Helmholtz, Stokes, etc.). We present numerical examples that show convergence of up to 32nd-order in 2D and 9th-order in 3D with respect to the mesh size. [ABSTRACT FROM AUTHOR]

This paper presents a novel boundary integral equation (BIE) formulation for the two-dimensional time-harmonic water-waves problem. It utilizes a complex-scaled Laplace free-space Green's function, resulting in a BIE posed on the infinite boundaries of the domain. The perfectly matched layer (PML) coordinate stretching that is used to render propagating waves exponentially decaying allows for the effective truncation and discretization of the BIE unbounded domain. We show through a variety of numerical examples that, despite the logarithmic growth of the complex-scaled Laplace free-space Green's function, the truncation errors are exponentially small with respect to the truncation length. Our formulation uses only simple function evaluations (e.g., complex logarithms and square roots), hence avoiding the need to compute the involved water-wave Green's function. Finally, we show that the proposed approach can also be used to find complex resonances through a linear eigenvalue problem since the Green's function is frequency-independent. [ABSTRACT FROM AUTHOR]

Let Ω be the multiply connected domain in the extended complex plane ℂ ¯ obtained by removing m non-overlapping rectilinear segments from the infinite strip S = { z : Im z < π / 2 } . In this paper, we present an iterative method for numerical computation of a conformally equivalent bounded multiply connected domain G in the interior of the unit disk D and the exterior of m non-overlapping smooth Jordan curves. We demonstrate the utility of the proposed method through two applications. First, we estimate the capacity of condensers of the form (S,E) where E ⊂ S is a union of disjoint segments. Second, we determine the streamlines associated with uniform incompressible, inviscid and irrotational flow past disjoint segments in the strip S. [ABSTRACT FROM AUTHOR]

In this paper, a parallel probabilistic method using the Feynman--Kac formula of killed Brownian motions is proposed to solve the mixed boundary value problems (BVPs) of 3D Laplace equations. To avoid using reflecting Brownian motions and the calculation of their local time L(t) in the Feynman--Kac representation of solutions for Neumann and Robin BVPs, the proposed method uses an iterative approach to approximate the solutions where each iteration will use the Feynman--Kac formula to solve a pure Dirichlet problem, thus only involving killed Brownian motions. First, the boundary of the domain is decomposed with overlapping local patches formed by the intersection of hemispheres superimposed on the domain boundary. The iteration starts with an arbitrary initial guess for the Dirichlet data on the Neumann and Robin boundaries; then, using the Feynman--Kac formula for a pure Dirichlet problem with the current available Dirichlet data on the whole boundary, the solution over the hemispheres can be obtained by the Feynman--Kac formula for the killed Brownian motions, sampled by a walk-on-spheres (WOS) algorithm. Second, by solving a local boundary integral equation (BIE) over each hemisphere and a local patch on the domain boundary, the Dirichlet data on the Neumann and Robin boundaries can be updated. By continuing this process, the proposed hybrid probabilistic and deterministic BIE-WOS method gives a highly parallel algorithm for the global solution of any mixed-type BVPs of the Laplace equations. Numerical results of various mixed interior and exterior BVPs demonstrate the parallel efficiency and accuracy of the proposed method. [ABSTRACT FROM AUTHOR]

Fast and high-order accurate algorithms for three-dimensional elastic scattering are of great importance when modeling physical phenomena in mechanics, seismic imaging, and many other fields of applied science. In this paper, we develop a novel boundary integral formulation for the three-dimensional elastic scattering based on the Helmholtz decomposition of elastic fields, which converts the Navier equation to a coupled system consisted of Helmholtz and Maxwell equations. An FFT-accelerated separation of variables solver is proposed to efficiently invert boundary integral formulations of the coupled system for elastic scattering from axisymmetric rigid bodies. In particular, by combining the regularization properties of the singular boundary integral operators and the FFT-based fast evaluation of modal Green’s functions, our numerical solver can rapidly solve the resulting integral equations with a high-order accuracy. Several numerical examples are provided to demonstrate the efficiency and accuracy of the proposed algorithm, including geometries with corners at different wavenumbers. [ABSTRACT FROM AUTHOR]

This paper concerns an inverse elastic scattering problem which is to determine a rigid obstacle from time domain scattered field data for a single incident plane wave. By using the Helmholtz decomposition, we reduce the initial-boundary value problem for the time domain Navier equation to a coupled initial-boundary value problem for wave equations, and prove the uniqueness of the solution for the coupled problem by employing the energy method. The retarded single layer potential is introduced to establish a set of coupled boundary integral equations, and uniqueness is discussed for the solution of these boundary integral equations. Based on the convolution quadrature method for time discretization, the coupled boundary integral equations are reformulated into a system of boundary integral equations in the s -domain, and then a convolution quadrature based nonlinear integral equation method is proposed for the inverse problem. Numerical experiments are presented to show the feasibility and effectiveness of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, we consider the reconstruction of the shape and the impedance function of an obstacle from measurements of the scattered field at a collection of receivers outside the object. The data is assumed to be generated by plane waves impinging on the unknown obstacle from multiple directions and at multiple frequencies. This inverse problem can be reformulated as an optimization problem: that of finding band-limited shape and impedance functions which minimize the L2 distance between the computed value of the scattered field at the receivers and the given measurement data. The optimization problem is highly non-linear, non-convex, and ill-posed. Moreover, the objective function is computationally expensive to evaluate (since a large number of Helmholtz boundary value problems need to be solved at every iteration in the optimization loop). The recursive linearization approach (RLA) proposed by Chen has been successful in addressing these issues in the context of recovering the sound speed of an inhomogeneous object or the shape of a sound-soft obstacle. We present an extension of the RLA for the recovery of both the shape and impedance functions of the object. The RLA is, in essence, a continuation method in frequency where a sequence of single frequency inverse problems is solved. At each higher frequency, one attempts to recover incrementally higher resolution features using a step assumed to be small enough to ensure that the initial guess obtained at the preceding frequency lies in the basin of attraction for Newton’s method at the new frequency. We demonstrate the effectiveness of this approach with several numerical examples. Surprisingly, we find that one can recover the shape with high accuracy even when the measurements are generated by sound-hard or sound-soft objects, eliminating the need to know the precise boundary conditions appropriate for modeling the object under consideration. While the method is effective in obtaining high-quality reconstructions for many complicated geometries and impedance functions, a number of interesting open questions remain regarding the convergence behavior of the approach. We present numerical experiments that suggest underlying mechanisms of success and failure, pointing out areas where improvements could help lead to robust and automatic tools for the solution of inverse obstacle scattering problems. [ABSTRACT FROM AUTHOR]

The hot spots conjecture is only known to be true for special geometries. This paper shows numerically that the hot spots conjecture can fail to be true for easy to construct bounded domains with one hole. The underlying eigenvalue problem for the Laplace equation with Neumann boundary condition is solved with boundary integral equations yielding a non-linear eigenvalue problem. Its discretization via the boundary element collocation method in combination with the algorithm by Beyn yields highly accurate results both for the first non-zero eigenvalue and its corresponding eigenfunction which is due to superconvergence. Additionally, it can be shown numerically that the ratio between the maximal/minimal value inside the domain and its maximal/minimal value on the boundary can be larger than 1 + 10− 3. Finally, numerical examples for easy to construct domains with up to five holes are provided which fail the hot spots conjecture as well. [ABSTRACT FROM AUTHOR]

In this paper, we investigate numerically the long-time dynamics of a two-dimensional dendritic precipitate. We focus our study on the self-similar scaling behavior of the primary dendritic arm with profile x ∼ t α¹ and y ∼ t α2, and explore the dependence of parameters α1 and α2 on applied driving forces of the system (e.g. applied far-field flux or strain). We consider two dendrite forming mechanisms: the dendritic growth driven by (i) an anisotropic surface tension and (ii) an applied strain at the far-field of the elastic matrix. We perform simulations using a spectrally accurate boundary integral method, together with a rescaling scheme to speed up the intrinsically slow evolution of the precipitate. The method enables us to accurately compute the dynamics far longer times than could previously be accomplished. Comparing with the original work on the scaling behavior α1 = 0.6 and α2 = 0.4 [Phys. Rev. Lett. 71(21) (1993) 3461–3464], where a constant flux was used in a diffusion only problem, we found at long times this scaling still serves a good estimation of the dynamics though it deviates from the asymptotic predictions due to slow retreats of the dendrite tip at later times. In particular, we find numerically that the tip grows self-similarly with α1 = 1/3 and α2 = 1/3 if the driving flux J ∼ 1/R(t) where R(t) is the equivalent size of the evolving precipitate. In the diffusive growth of precipitates in an elastic media, we examine the tip of the precipitate under elastic stress, under both isotropic and anisotropic surface tension, and find that the tip also follows a scaling law. [ABSTRACT FROM AUTHOR]

• Novel boundary-based ML method for 2D solid mechanics. • Requires only boundary discretization, enabling faster learning. • Replaces unstable differential operators by stable integral operators. • Higher precision and efficiency than traditional ML methods. In this paper, we make the first attempt to adopt the boundary integrated neural networks (BINNs) for the numerical solution of two-dimensional (2D) elastostatic and piezoelectric problems. The proposed BINNs combine the artificial neural networks with the exact boundary integral equations (BIEs) to effectively solve the boundary value problems based on the corresponding partial differential equations (PDEs). The BIEs are utilized to localize all the unknown physical quantities on the boundary, which are approximated by using artificial neural networks and resolved via a training process. In contrast to many traditional neural network methods based on a domain discretization, the present BINNs offer several distinct advantages. Firstly, by embedding the analytical BIEs into the learning procedure, the present BINNs only need to discretize the boundary of the problem domain, which reduces the number of the unknowns and can lead to a faster and more stable learning process. Secondly, the differential operators in the original PDEs are substituted by integral operators, which can effectively eliminate the need for additional differentiations of the neural networks (high-order derivatives of the neural networks may lead to instabilities in the learning process). Thirdly, the loss function of the present BINNs only contains the residuals of the BIEs, as all the boundary conditions have been inherently incorporated within the formulation. Therefore, there is no necessity for employing any weighting functions, which are commonly used in most traditional methods to balance the gradients among the different objective functions. Extensive numerical experiments show that the present BINNs are much easier to train and can usually provide more accurate solutions as compared to many traditional neural network methods. [Display omitted] [ABSTRACT FROM AUTHOR]

Problem of electromagnetic wave scattering by ideal conductors with thin dielectric coatings is considered. The problem is reduced to a system of volume and surface hypersingular integral equations. We solve the system by Galerkin's method, using piecewise linear RWG functions for the surface equation, whilst for the volume equation we propose a combined system of piecewise linear and piecewise constant functions with prismatic supports. The method is tested on model problems, which showed the applicability of the proposed algorithm to the solution of scattering problems. [ABSTRACT FROM AUTHOR]

In this paper, we present a novel technique to accurately model scattering phenomena at two-dimensional circular and rectangular structures consisting of arbitrary homogeneous materials, including magnetic media in particular. The proposed formalism utilizes a differential surface admittance operator, which invokes a single source equivalence theorem to replace the inside material by its surrounding medium, while introducing an equivalent surface current density. The arbitrary magnetic contrast can be combined with an arbitrary electrical conductivity. As such, the skin effect is rigorously taken into account, making our method ideally suited for broadband modeling of good conductors as well. It is demonstrated that an appropriate choice of the basis functions for the discretized problem is critical to obtain a convergent result when magnetic contrast is introduced. The method is analytically validated for the case of a circular cylinder and additional numerical results illustrate the correctness of the technique for (combinations of) rectangular cylinders, through comparison with a Poggio-Miller-Chan-Harrington-Wu-Tsai approach, a volume integral equation method and a commercial solver. [ABSTRACT FROM AUTHOR]

The paper is concerned with using boundary element methods (BEM) for the accurate evaluation of fields in structures with clusters. For large-scale problems, the BEM system is solved iteratively by speeding up matrix-to-vector multiplications by applying a kernel-independent fast multipole method. Multiplication starts with source-to-multipole (S2M) translations, whose accuracy predefines the overall accuracy. We aim to increase the accuracy of these translations. The intensities of sources are assembled into clusters by an algorithm suggested. Each of them is characterized by its representative source, whose intensity equals the sum of the intensities of cluster sources. Thus, with growing distance, its field tends toward the field of the cluster. The accuracy of S2M translations is increased by subtracting from and adding to the far field of the cluster the far field of its representative source, and by using the proposed modified kernel to evaluate the difference of the fields, which decreases faster than the field of the cluster itself. Numerical results for typical kernels show a notable increase in the accuracy provided by the modified S2M translations. Keeping in them merely the added field is acceptable for many practical applications. This simplifies the modified S2M translations by avoiding calculation and storing matrices specific to each of the clusters. The improved translations may be also used for multipole-to-multipole translations, performed on next, after leaves, levels in upward running a hierarchical tree.

A fast and spectrally accurate Ewald summation method for the evaluation of stokeslet, stresslet and rotlet potentials of three-dimensional Stokes flow is presented. This work extends the previously developed Spectral Ewald method for Stokes flow to periodic boundary conditions in any number (three, two, one, or none) of the spatial directions, in a unified framework. The periodic potential is split into a short-range and a long-range part, where the latter is treated in Fourier space using the fast Fourier transform. A crucial component of the method is the modified kernels used to treat singular integration. We derive new modified kernels, and new improved truncation error estimates for the stokeslet and stresslet. An automated procedure for selecting parameters based on a given error tolerance is designed and tested. Analytical formulas for validation in the doubly and singly periodic cases are presented. We show that the computational time of the method scales like O (N log ⁡ N) for N sources and targets, and investigate how the time depends on the error tolerance and window function, i.e. the function used to smoothly spread irregular point data to a uniform grid. The method is fastest in the fully periodic case, while the run time in the free-space case is around three times as large. Furthermore, the highest efficiency is reached when applying the method to a uniform source distribution in a primary cell with low aspect ratio. The work presented in this paper enables efficient and accurate simulations of three-dimensional Stokes flow with arbitrary periodicity using e.g. boundary integral and potential methods. • A unified framework for fast summation of Stokes potentials in arbitrary periodicity. • Improved error estimates permit accurate and automatic selection of parameters. • New modified kernels ensure the truncation error to be independent of periodicity. • New analytical formulas permit validation in doubly and singly periodic cases. • Our method scales like O (N log ⁡ N) for N sources and target points. [ABSTRACT FROM AUTHOR]

We propose and study a method for finding quasi-resonances for a linear acoustic transmission problem in frequency domain. Starting from an equivalent boundary-integral equation we perform Galerkin boundary element discretization and look for the minima of the smallest singular value of the resulting matrix as a function of the wave number k. We develop error estimates for the impact of Galerkin discretization on singular values and devise a heuristic adaptive algorithm for finding the minima in prescribed k-intervals. Our method exclusively relies on the solution of eigenvalue problems for real k, in contrast to alternative approaches that rely on extension to the complex plane. [ABSTRACT FROM AUTHOR]

The scattering of a time-harmonic electromagnetic wave by a penetrable chiral obstacle in an achiral environment is considered. The method of fundamental solutions is employed in order to obtain numerically the solution of the problem using fundamental solutions in dyadic form. Surface vector potentials in terms of dyadic fundamental solutions together with the associated boundary integral operators are defined and their regularity properties are presented. Based on the dependence of the solution to the boundary data, appropriate systems of functions containing elements of dyadic fundamental solutions on the surface of the scatterer are constructed. Completeness and linear independence for these systems are proved with the usage of surface vector potentials. Using the transmission conditions, the scattering problem is transformed into a linear algebraic system with a coefficient matrix which consists of chiral and achiral blocks. [ABSTRACT FROM AUTHOR]

We consider the computation of the transmission eigenvalue problem based on a boundary integral formulation. The problem is formulated as the eigenvalue problem of a holomorphic Fredholm operator function. A Fourier–Galerkin method is proposed for the integral equations. The approximation properties of the associated discrete operators are analyzed and some convergence results of the eigenvalues are obtained. We present the details of the implementation and employ the spectral projection method to compute the eigenvalues. Numerical examples validate the effectiveness and accuracy of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, the Cauchy problem for the Helmholtz equation is investigated. The objective is to recover the missing data on some part of the boundary of a bounded domain from overspecified data on the remaining part of the boundary. We propose a preconditioned Krylov algorithm to solve this ill-posed problem, based on the representation of the solution with surface integral equations and the Steklov–Poincaré operator. We give a theoretical and numerical validation of the proposed method conducted in the 3D setting. We show the fast convergence of the proposed algorithm tested on various synthetic examples. The numerical results show a high precision of the reconstruction obtained for different levels of noisy data. [ABSTRACT FROM AUTHOR]

We consider eigenvalue problems for dielectric cylindrical scatterers of arbitrary cross with generalized conditions at infinity that enable one to take into account complex eigenvalues. The existence of resonance (scattering) frequencies associated with these eigenvalues is proved. The technique involves the determination of characteristic numbers (CNs) of the Fredholm operator-valued functions of the problems constructed using Green's potentials. Separating principal parts in the form of meromorphic operator pencils, we apply the operator generalization of Rouché's theorem to verify the occurrence of CNs in close proximities of the pencil poles. The results are illustrated in detail using the case of a dielectric cylinder of circular cross section. [ABSTRACT FROM AUTHOR]

The first and second initial-boundary value problems for second-order parabolic systems with coefficients satisfying the Dini condition in a semibounded plane domain with a nonsmooth lateral boundary admitting cusps are considered. Theorems on the uniqueness of classical solutions of these problems in the class of functions that are continuous and bounded together with their first spatial derivatives in the closure of this domain are proved. [ABSTRACT FROM AUTHOR]

We are concerned with the numerical computation of electrostatic forces/torques in only piece-wise homogeneous materials using the boundary element method (BEM). Conventional force formulas based on the Maxwell stress tensor yield functionals that fail to be continuous on natural trace spaces. Thus their use in conjunction with BEM incurs slow convergence and low accuracy. We employ the remedy discovered in [P. Panchal and R. Hiptmair, Electrostatic force computation with boundary element methods, SMAI J. Comput. Math.8 (2022), 49–74]. Motivated by the virtual work principle which is interpreted using techniques of shape calculus, and using the adjoint method from shape optimization, we derive stable interface-based force functionals suitable for use with BEM. This is done in the framework of single-trace direct boundary integral equations for second-order transmission problems. Numerical tests confirm the fast asymptotic convergence and superior accuracy of the new formulas for the computation of total forces and torques. [ABSTRACT FROM AUTHOR]

We analyze the mathematical model that describes the heat generated by electromagnetic nanoparticles. We use the known optical properties of the nanoparticles to control the support and amount of the heat needed around a nanoparticle. Precisely, we show that the dominant part of the heat around the nanoparticle is the electric field multiplied by a constant which is explicitly and solely dependent on the permittivity and quantities related to the eigenvalues and eigenfunctions of the magnetization (or the Newtonian) operator, defined on the nanoparticle, and inversely proportional to the distance to the nanoparticle. The nanoparticles are described via the Lorentz model. If the used incident frequency is chosen related to the plasmonic frequency ω p (via the magnetization operator), then the nanoparticle behaves as a plasmonic one; if it is chosen related to the undamped resonance frequency ω 0 (via the Newtonian operator), then it behaves as a dielectric one. The two regimes exhibit different optical behaviors. In both cases, we estimate the generated heat and discuss advantages of each incident frequency regime. The analysis is based on time-domain integral equation techniques, avoiding the use of (formal) Fourier type transformations. [ABSTRACT FROM AUTHOR]

The determination of rigid-body motion caused by incident waves is a familiar problem in mechanics. Good examples are the motion of a ship in the presence of water waves and the motion of a rigid structure in the presence of seismic waves. The basic goals are to determine the motion of the rigid body and the effects of the motion on the wave field, assuming linear theory is adequate. Although the underlying mathematical problems are similar, several solution methods have evolved, depending on the physical problems of interest. For ship motions, the standard approach is to decompose the problem into seven subproblems, one for each of the six rigid-body modes and one to take account of the incident wave. This approach is reviewed and then adapted to problems in acoustics and to problems in elastodynamics, such as those that arise in simple examples of soil–structure interaction. It is argued that the resulting approach for elastodynamic problems has clear advantages over those currently in use. [ABSTRACT FROM AUTHOR]

We develop a numerical method for solving three-dimensional problems of fluid filtration and absorption in a piecewise homogeneous medium by means of boundary integral equations. This method is applied to a simulation of the lymph flow in a lymph node. The lymph node is considered as a piecewise homogeneous domain containing porous media. The lymph flow is described by Darcy's law. Taking into account the lymph absorption, we propose an integral representation for the velocity and pressure fields, where the lymph absorption imitates the lymph outflow from a lymph node through a system of capillaries. The original problem is reduced to a system of boundary integral equations, and a numerical algorithm for solving this system is provided. We simulate the lymph velocity and pressure as well as the total lymph flux. The method is verified by comparison with experimental data. [ABSTRACT FROM AUTHOR]

We consider the Helmholtz transmission problem with piecewise-constant material coefficients and the standard associated direct boundary integral equations. For certain coefficients and geometries, the norms of the inverses of the boundary integral operators grow rapidly through an increasing sequence of frequencies, even though this is not the case for the solution operator of the transmission problem; we call this phenomenon that of spurious quasi-resonances. We give a rigorous explanation of why and when spurious quasi-resonances occur and propose modified boundary integral equations that are not affected by them. [ABSTRACT FROM AUTHOR]

The Robin problem for the diffusion equation in non-homogeneous media partial differential equation is reduced to a system of direct segregated parametrix-based Boundary-Domain Integral Equations (BDIEs). We use a parametrix different from the one applied by Chkadua, Mikhailov, Natroshvili. in [4]. We prove the equivalence between the original BVP and the corresponding BDIE system. The invertibility and Fredholm properties of the boundary-domain integral operators are also analysed. [ABSTRACT FROM AUTHOR]

We provide a new analytical and computational study of the transmission eigenvalues with a conductive boundary condition. These eigenvalues are derived from the scalar inverse scattering problem for an inhomogeneous material with a conductive boundary condition. The goal is to study how these eigenvalues depend on the material parameters in order to estimate the refractive index. The analytical questions we study are: deriving Faber–Krahn type lower bounds, the discreteness and limiting behavior of the transmission eigenvalues as the conductivity tends to infinity for a sign changing contrast. We also provide a numerical study of a new boundary integral equation for computing the eigenvalues. Lastly, using the limiting behavior we will numerically estimate the refractive index from the eigenvalues provided the conductivity is sufficiently large but unknown. [ABSTRACT FROM AUTHOR]

The first and second initial-boundary value problems for Petrovskii parabolic systems of the second order with coefficients satisfying the Dini condition in plane domains with nonsmooth lateral boundaries admitting, in particular, cusps are considered. Theorems on the uniqueness of classical solutions of these problems in the class of functions that are continuous and bounded together with their first spatial derivatives in the closure of these domains are proved. [ABSTRACT FROM AUTHOR]

The second initial-boundary value problem for a second-order Petrovskii parabolic system with constant coefficients in a semibounded plane domain with a nonsmooth lateral boundary is considered. The uniqueness of a solution to this problem in the class is proved. The minimum condition on the boundary function under which the solution of the problem belongs to is investigated. A constructive solution is obtained by applying the boundary integral equation method. [ABSTRACT FROM AUTHOR]

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 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 and âˆ'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 3D (conical) diffraction problems are discussed. [ABSTRACT FROM AUTHOR]

This work is devoted to studying a direct and inverse scattering problem for a magnetoelastic layer having a defect, in the frame of the electromagnetic theory. In terms of the displacement field over the defect's contour, a coupled system of boundary integral equations is formulated, for magnetically permeable and impermeable defects. To identify the position and size of the defect, an efficient numerical algorithm is developed by using the quasi-Newton iterative method. In order to check the influence of the magnetic field upon the scattering waves from the layer, a series of numerical examples is presented with different noise levels. The results showed that the magnetic field has a sensitive effect on the identification process when the external magnetic field increases, especially for the materials having a high magnetic permeability factor μ r . Also, a special inverse problem for predicting the external applied magnetic field, upon a copper layer having a defect with various sizes, has been performed. [ABSTRACT FROM AUTHOR]

A two-dimensional pulsating flow of a viscous fluid in a plane channel whose wall has rectangular microcavities partially or completely filled with a compressible gas is investigated. This problem formulation can clarify the friction reduction mechanism in a laminar sublayer of a turbulent viscous boundary layer flow over a textured stripped superhydrophobic surface containing periodically arranged rectangular micro-cavities filled with gas. It is assumed that the dimensions of the cavities are much smaller than the channel thickness. On the macroscale, the problem of one-dimensional unsteady viscous flow in a plane channel with no-slip conditions on the walls and a harmonic variation of the pressure difference is solved. The solution obtained in this way is used for formulating non-stationary in time and periodic in space boundary conditions for the flow on the scale of a chosen cavity (microscale), with the instantaneous volume of the gas bubble in the cavity depending on the instantaneous pressure over the cavity. The flow on the microscale near a cavity with a gas bubble occurs in the Stokes regime. The numerical solution is obtained using an original version of the boundary element method. A parametric numerical study of the flow field in a pulsating shear flow over a cavity with a compressible gas bubble is performed. The averaged parameters characterizing the effective 'velocity slip' of viscous fluid and the friction reduction in a pulsating flow over a stripped superhydrophobic surface are calculated. [ABSTRACT FROM AUTHOR]

The lymph node (LN) represents a key structural component of the lymphatic system network responsible for the fluid balance in tissues and the immune system functioning. Playing an important role in providing the immune defense of the host organism, LNs can also contribute to the progression of pathological processes, e.g., the spreading of cancer cells. To gain a deeper understanding of the transport function of LNs, experimental approaches are used. Mathematical modeling of the fluid transport through the LN represents a complementary tool for studying the LN functioning under broadly varying physiological conditions. We developed an artificial neural network (NN) model to describe the lymph node drainage function. The NN model predicts the flow characteristics through the LN, including the exchange with the blood vascular systems in relation to the boundary and lymphodynamic conditions, such as the afferent lymph flow, Darcy's law constants and Starling's equation parameters. The model is formulated as a feedforward NN with one hidden layer. The NN complements the computational physics-based model of a stationary fluid flow through the LN and the fluid transport across the blood vessel system of the LN. The physical model is specified as a system of boundary integral equations (IEs) equivalent to the original partial differential equations (PDEs; Darcy's Law and Starling's equation) formulations. The IE model has been used to generate the training dataset for identifying the NN model architecture and parameters. The computation of the output LN drainage function characteristics (the fluid flow parameters and the exchange with blood) with the trained NN model required about 1000-fold less central processing unit (CPU) time than computationally tracing the flow characteristics of interest with the physics-based IE model. The use of the presented computational models will allow for a more realistic description and prediction of the immune cell circulation, cytokine distribution and drug pharmacokinetics in humans under various health and disease states as well as assisting in the development of artificial LN-on-a-chip technologies. [ABSTRACT FROM AUTHOR]

When using boundary integral equation methods, we represent solutions of a linear partial differential equation as layer potentials. It is well-known that the approximation of layer potentials using quadrature rules suffer from poor resolution when evaluated closed to (but not on) the boundary. To address this challenge, we provide modified representations of the problem's solution. Similar to Gauss's law used to modify Laplace's double-layer potential, we use modified representations of Laplace's single-layer potential and Helmholtz layer potentials that avoid the close evaluation problem. Some techniques have been developed in the context of the representation formula or using interpolation techniques. We provide alternative modified representations of the layer potentials directly (or when only one density is at stake). Several numerical examples illustrate the efficiency of the technique in two and three dimensions. [ABSTRACT FROM AUTHOR]