Integral Equation Formulation Research Articles

Semi-implicit Solver for the Heat Equation with Stefan–Boltzmann Law Boundary Condition

The surface energy balance on an atmosphereless body consists of solar irradiance, subsurface heat conduction, and thermal radiation to space by the Stefan–Boltzmann law. Here we extend the semi-implicit Crank–Nicolson method to this specific nonlinear boundary condition and validate its accuracy. A rapid change in incoming solar flux can cause a numerical instability, and several approaches to dampen this instability are analyzed. A predictor based on the Volterra integral equation formulation for the heat equation is also derived and can be used to improve accuracy and stability. The publicly available implementation provides a fast and robust thermophysical model that has been applied to lunar, Martian, and asteroidal surfaces, on occasion to millions of surface facets or parameter combinations.

The Numerical Solutions of weakly singular Fredholm integral equations of the Second kind Using Chebyshev Polynomials of the Second Kind

In this study, the second kind Chebyshev Polynomials were utilized to acquire interpolated solutions for the second kind Fredholm integral equations with weakly singular kernel. To accomplish this, the data, unknown, and kernel functions were converted into matrix form, and consequently we completely isolated the singularity of the kernel. The primary benefit of this method is the ability to change the form of integral equation to an equivalent algebraic system, which is easier to solve. The effectiveness of our technique was evaluated by applying it to three illustrated examples, and it was observed that the solutions obtained exhibit strong convergence towards the exact solutions.

Numerical simulation of an extensible capsule using regularized Stokes kernels and overset finite differences

In this paper, we present a novel numerical scheme for simulating deformable and extensible capsules suspended in a Stokesian fluid. The main feature of our scheme is a partition-of-unity (POU) based representation of the surface that enables asymptotically faster computations compared to spherical-harmonics based representations. We use a boundary integral equation formulation to represent and discretize hydrodynamic interactions. The boundary integrals are weakly singular. We use the quadrature scheme based on the regularized Stokes kernels by Tlupova and Beale 2019 (given in [34]). We also use partition-of unity based finite differences that are required for the computation of interfacial forces. Given an N-point surface discretization, our numerical scheme has fourth-order accuracy and O(N) asymptotic complexity, which is an improvement over the O(N2log⁡N) complexity of a spherical harmonics based spectral scheme that uses product-rule quadratures by Veerapaneni et al. 2011 [36]. We use GPU acceleration and demonstrate the ability of our code to simulate the complex shapes with high resolution. We study capsules that resist shear and tension and their dynamics in shear and Poiseuille flows. We demonstrate the convergence of the scheme and compare with the state of the art.

Tricomi problem and integral equations

Formulas for inverting integral equations that arise when studying the Tricomi problem for the Lavrentyev–Bitsadze equation were derived. Solvability conditions of an auxiliary overdetermined problem in the elliptic part of the mixed domain were found using the Green function method. A connection was established between the Green functions of the Dirichlet problem and problem N for the Laplace equation in the form of integral equations mutually inverting each other. Various integral equations were considered, including explicitly solvable ones, to which the Tricomi problem can be reduced. An explicit solution of the characteristic singular equation with a Cauchy kernel was obtained without involving the theory of boundary value problems for analytic functions.

Stress-driven nonlocal integral model with discontinuity for size-dependent buckling and bending of cracked nanobeams using Laplace transform

Due to the frequent occurrence of cracks as structures age, examining the size-dependent responses of cracked micro/nanobeams is critical for safely designing and optimizing high-tech devices such as those used in MEMS/NEMS. With this understanding, we present the size-dependent bending and buckling studies of cracked nanobeams with different boundary edges adopting the well-posed stress-driven nonlocal integral model with discontinuity. Since the presence of a crack, the beam is divided into two segments connected by a rotary spring and a linear spring, whose elastic coefficients can be determined through compliance calibration methods. We propose an efficient method to solve the original integral model directly by converting the integral form of the nonlocal constitutive equation into a general equation through the Laplace transform technique, equipped with three extra constraint conditions. Subsequently, the governing equations of two beam segments are solved by considering the standard boundary and compatibility conditions as well as the extra constraint conditions, while the bending deformation and buckling loads of various bounded nanobeams are presented analytically in the presence of a crack. After verifying the derived formulation, the influence of crack length and location as well as nonlocality on the beam deflections and buckling loads of the first several modes are investigated in detail. Communicated by Rossana Dimitri.

On the Domain Decomposition Method Preconditioning of Surface Integral Equation Formulations Solved by GMRES

On the Domain Decomposition Method Preconditioning of Surface Integral Equation Formulations Solved by GMRES

DETERMINATION OF THE DROP SIZE AND DISTRIBUTIONS IN SWIRL INJECTION IN CROSS FLOWS, IMPINGING, AND EFFERVESCENT INJECTORS

We have extended the primary atomization analysis to swirl injection in cross flows, impinging, and effervescent injectors. Using the integral form of the conservation equations, the drop size can be expressed in terms of injection and fluid parameters, the main variable being the liquid and gas velocities. Using the measured velocities as inputs to this <i>D</i><sub>32</sub>-equation, good agreements with experimental data are found for the drop size in the three spray geometries. Underlying physical mechanisms for the drop formation are also revealed from the analysis. The aerodynamic interaction between the swirl spray and cross flow results in reduction in momentum, with a corresponding decrease in kinetic energy that appears as surface tension of energy of many small droplets. Similarly, cancellation of the lateral momentum in impinging jets and internal deceleration in effervescent injectors are the key primary atomization routes. The use of the analytical drop size-velocity correlation has also been demonstrated for swirl sprays in cross flows. Therefore, this approach can be used to predict the drop size and distributions in different spray geometries, with appropriate changes in the velocity input terms and fluid properties.

Internal Resonant Characteristics Analysis of Arbitrarily Shaped Metallic Cavity with Holes Based on A Surface Integral Equation Formulation

Internal Resonant Characteristics Analysis of Arbitrarily Shaped Metallic Cavity with Holes Based on A Surface Integral Equation Formulation

Stochastic Modeling in Immunology Based On a Stage-Dependent Framework with Non-Markov Constraints for Individual Cell and Pathogen Dynamics

We present a systematic approach to modelling the responses of the immune system to virus infections. Two continuous-discrete stochastic models arising in mathematical immunology are developed and computationally implemented. The variables of the models are integer random variables that denote the quantity of individuals (cells and viral particles), and sets of unique types of individuals that take into account the current state and history of stay of individuals in some stages of their development. The distribution laws of the durations of the mentioned stages are different from exponential or geometric. A probabilistic description of a one-stage stochastic model of population dynamics is presented. A stochastic model of the development of HIV-1 infection in the lymph node in the initial period after infection of a healthy person is formulated. A computational algorithm based on the Monte Carlo method is given. Each of the stochastic models is complemented by a deterministic analogue in the form of integral and delay differential equations. The results of numerical simulation are presented.

Mathematical micro–macro modeling of fully coupled nonlinear magneto-elastic reinforced composites

In this paper, integral equation formulations and field-dependent micromechanical modeling of global nonlinear magneto-mechanical behavior of magneto-elastic composites under high strain and magnetic field are elaborated. The modeling is obtained based on an extension of micro–macro transition inclusion problem to nonlinear behavior using field-dependent and highly nonlinear localization tensors. A methodological procedure is elaborated based on newly introduced strain and magnetic field-dependent Green tensors, strain and magnetic field-dependent integral equations linked to tensors of Eshelby and micromechanical approaches. The field-dependent global moduli are predicted based on the Mori–Tanaka approach and self-consistent method. Iterative incremental schemes based on the Newton–Raphson and fixed-point algorithms are examined and established precise semi-analytic equations of the effective magneto-elastic properties of composites that are dependent on the magnetic-strain field for different types of inclusions. A numerical code is elaborated for numerical predictions and the obtained field-dependent effective magneto-elastic coefficients are obtained and presented for various volume fractions of inclusion, types and shapes of the reinforced nonlinear composites.

A Generalized Source Integral Equation for Enhanced Compression in Three Dimensions

An integral equation formulation, using generalized directional sources, for three-dimensional scattering by impenetrable and essentially convex bodies, is presented. This formulation increases the rank-deficiency of moment matrix blocks representing interactions between non-overlapping source and observer regions. Thereby, it enables enhanced low-rank approximation based matrix compression and the development of corresponding fast direct solvers. The directional sources are constructed by augmenting the conventional basis functions with spherical absorbing three-dimensional shields, on which auxiliary source distributions are defined. The bottleneck arising from the need to integrate over these distributions upon computing the modified Green’s function is removed by using efficient non-uniform sampling and tabulation of the modified Green’s function or its components, in a region-dependent manner. The formulation is studied and its favorable compressibility is demonstrated, for two fundamental types of compression strategies. The non-uniform sampling approach is employed also to facilitate the rank-revealing analysis of large off-diagonal blocks.

Efficient Calculation of the 3-D Rectangular Waveguide Green’s Functions Derivatives by the Ewald Method

In this contribution, the Ewald method has efficiently been applied to accelerate the computation of the rectangular waveguide Green’s functions (GFs) derivatives. Based on previous works, we have outlined new approximation formulae that avoid the evaluation of computationally expensive complementary error functions of complex argument, needed by the Ewald method. This is possible when the internal medium of the rectangular waveguide is homogeneous and lossless. On the other hand, different convergence numerical studies have been carried out, showing a similar convergence rate for computing the original GFs and their derivatives. Moreover, we have checked that the computational time is only slightly increased for obtaining the derivatives as compared to the original GFs, after the application of these new techniques. The newly derived expressions are useful for the evaluation of electromagnetic fields, the characterization of dielectric materials and step discontinuities between rectangular waveguides, and the analysis of rectangular cavities using integral equation formulations. For validation, the electric field produced by a surface electric current density with a rectangular pulse distribution has been evaluated, using the new proposed expressions. These results have been compared to simulations provided by full-wave finite elements commercial software to verify their correctness, exhibiting a good agreement.

Existence and stability of solution for a nonlinear Volterra integral equation with binary relation via fixed point results

The real-world problems that can be resolved by a dynamic mathematical model is an analogous form of a complex integral equation. In this article, we intend to exhibit the solution to a non-homogeneous, nonlinear Volterra integral equation endowed with an arbitrary binary relation utilizing the fixed point solution that is originally proved for a multivalued map and then reduced to a single valued self map. The results proved and the application demonstrated is an added asset to the present literature and can significantly contribute in determining the existence of a solution for a dynamic mathematical model in the form of a Volterra integral equation equipped with an arbitrary binary relation via fixed point findings. Furthermore, the stability of the solution of the Volterra integral equation is verified using the Hyers–Ulam stability concept.

Size-dependent nonlinear forced vibration of microbeam using two-phase local/nonlocal viscoelastic integral model with thermal effects

Combination strain-driven (εD) and stress-driven (σD) two-phase local/nonlocal integral models (TPNIM) with classic Kelvin-Voigt viscoelastic model, εD- and σD-two-phase local/nonlocal viscoelastic integral models (TPNVIM) are proposed to study the viscoelastic nonlinear forced vibration behaviors of microbeams in thermal environment. Based on von-Karman large deformation theory and the Hamilton’s principle, the differential governing equations of motion and boundary conditions are derived for nonlinear forced vibration of microbeam. Several non-dimensional variables are introduced to simplify the mathematical formulation. Nonlocal viscoelastic equations in integral form are transferred unitedly into equivalent differential form with constitutive constraints. The axial force of nonlocal microbeam with immovable ends is derived explicitly. Vibration frequency and corresponding linear vibration mode shape are derived through Laplace transformation for linear thermo-elastic vibration of microbeam using εD- and σD-TPNIMs. The Ritz-Galerkin approach is utilized to derive the reduced differential governing equation for nonlinear vibration of microbeam, which has exact same form as single-degree-of-freedom oscillator with cubic nonlinearity when material viscosity is neglected. The method of multiple scales (MMS) based on trigonometric functions is utilized to solve the reduced nonlinear differential equation approximately. The frequency-amplitude responses are addressed for primary and super-harmonic resonances of microbeam. The effects of viscous coefficient, amplitude of excitation force and environmental temperature variation as well as nonlocal length parameter on the nonlinear vibration behavior of microbeam are investigated numerically for different boundary conditions.

Nyström discretizations of boundary integral equations for the solution of 2D elastic scattering problems

We present three high-order Nyström discretization strategies of various boundary integral equation formulations of the impenetrable time-harmonic Navier equations in two dimensions. One class of such formulations is based on the four classical Boundary Integral Operators (BIOs) associated with the Green’s function of the Navier operator. We consider two types of Nyström discretizations of these operators, one that relies on Kussmaul–Martensen logarithmic splittings (Chapko et al., 2000; Domínguez and Turc, 2000), and the other on Alpert quadratures (Alpert, 1999). In addition, we consider an alternative formulation of Navier scattering problems based on Helmholtz decompositions of the elastic fields (Dong et al., 2021), which can be solved via a system of boundary integral equations that feature integral operators associated with the Helmholtz equation. Owing to the fact that some of the BIOs that are featured in those formulations are non-standard, we use Quadrature by Expansion (QBX) methods for their high order Nyström discretization. Alternatively, we use Maue integration by parts techniques to recast those non-standard operators in terms of single and double layer Helmholtz BIOs whose Nyström discretizations is amenable to the Kussmaul–Martensen methodology. We present a variety of numerical results concerning the high order accuracy that our Nyström discretization elastic scattering solvers achieve for both smooth and Lipschitz boundaries. We also present extensive comparisons regarding the iterative behavior of solvers based on different integral equations in the high frequency regime. Finally, we illustrate how some of the Nyström discretizations we considered can be incorporated seamlessly into the Convolution Quadrature (CQ) methodology to deliver high-order solutions of the time domain elastic scattering problems.

Spectral boundary integral equation method for simulation of 2D and 3D slip ruptures at bi‐material interfaces

AbstractThis paper presents a 3D spectral numerical scheme that can accurately model planarcracks of any shape residing on an interfacial plane between two elastic solids. The cracks propagate dynamically under the action of interfacial loads and interfacial constitutive laws. The scheme is a spectral form of the boundary integral equations (BIE) that relate the displacement discontinuity fields along the interfacial plane to the stress fields at the plane. In the BIE, the displacement discontinuities are expressed as a spatio‐temporal convolution of the traction components at the interface. This is in contrast with most previous studies which use a spatio‐temporal convolution of the displacement discontinuities or of the displacements at the interface. Due to the continuity of tractions across the interface, the present formulation is simpler, resulting in convolution kernels that can be written in closed‐form. Previous studies have relied on numerically obtained convolution kernels. The accuracy of the proposed scheme is validated with the known analytical solution of the 3D Lamb problem. Furthermore, new algorithms are developed for coupling the BIE with a slip‐weakening friction law, both for homogeneous materials as well as for bi‐material interfaces. This results in a widely applicable formulation for studying spontaneous rupture propagation. The algorithms are validated by comparing rupture simulation exercises to benchmark problems from Southern California Earthquake Center (SCEC). Finally, the methodology is used to simulate 3D frictional rupture propagation along a bi‐material interface.

Fractional approach for Dirac operator involving M-truncated derivative

In this study, we examine the basic spectral information for systems governed by the Dirac equation with distinct boundary conditions, utilizing a modified form of local derivatives known as M-truncated derivative (MTD). The spectral information discussed includes the representation of solutions in the form of integral equations, the asymptotics vector-valued eigenfunctions and eigenvalues, and their normalized forms, all within the context of the MTD method that incorporates truncated Mittag-Leffler functions. This type of MTD provides the features of integer-order operator theory. Also, by virtue of the parameters $\alpha $ and $\gamma$, we analyze and compare the solutions with graphs in terms of different potentials, different eigenvalues and different orders. Thus, the aim of this article is to consider spectral structure of Dirac system in frame of M-truncated derivative by proping with visual analysis.

An SIE Formulation With Triangular Discretization and Loop Analysis for Parameter Extraction of Arbitrarily Shaped Interconnects

A surface integral equation (SIE) formulation under the magneto-quasi-static assumption is proposed to efficiently and accurately model arbitrarily shaped interconnects in packages. Through decently transferring all electromagnetic quantities into circuit elements, the loop analysis is used to carefully construct matrix equations with an independent and complete set of unknowns based on graph theory. In addition, an efficient preconditioner is developed, and the proposed formulation is accelerated by the pre-corrected fast Fourier transform (pFFT). Four practical examples, including a rectangular metallic interconnect, bonding wire arrays, interconnects in a real-life circuit, and the power distribution network (PDN) used in packages, are carried out to validate its accuracy, efficiency, and scalability. Results show that the proposed formulation is accurate, efficient, and flexible to model complex interconnects in packages.

The measurement of true initial rates is not always absolutely necessary to estimate enzyme kinetic parameters

In the chapters dealing with enzyme reactions, the authors of all Biochemistry textbooks and of even more specialized texts consider that the characteristic parameters (kcat and Km) must be determined under initial or steady-state rate conditions. This implies the transformation of a very limited proportion of substrate (at most 10–20%) or a continuous recording of the product or substrate concentration vs. time. Both options can present practical difficulties. Is it possible to get around these very stringent conditions? Here we show that in the most favourable cases up to 70% of the substrate can be converted resulting in systematic errors on the parameters (that can easily be taken account of) if the simple Henri-Michaelis–Menten equation is utilised. Alternatively, the integrated form of the same equation directly yields excellent estimates of the same parameters. Our observations should greatly facilitate the task of researchers who study systems in which measurements of the reaction progress are painstaking or when substrate concentrations close to the detection limit must be used. The general conclusion is that it is not always absolutely necessary to determine initial or steady-state rates to obtain reliable estimations of the enzyme kinetic parameters..

Critical local well‐posedness for the fully nonlinear Peskin problem

AbstractWe study the problem where a one‐dimensional elastic string is immersed in a two‐dimensional steady Stokes fluid. This is known as the Stokes immersed boundary problem and also as the Peskin problem. We consider the case with equal viscosities and with a fully non‐linear tension law; this model has been called the fully nonlinear Peskin problem. In this case we prove local in time wellposedness for arbitrary initial data in the scaling critical Besov space . We additionally prove the optimal higher order smoothing effects for the solution. To prove this result we derive a new formulation of the boundary integral equation that describes the parametrization of the string, and we crucially utilize a new cancelation structure.