Mixed Boundary Value Problem Research Articles

An efficient fourth-order convergent scheme based on half-step spline function for two-point mixed boundary value problems

The present work aims to introduce a novel numerical method for solving second-order two-point mixed boundary value problems. Mixed boundary value problems occur in various scientific fields, including quantum mechanics, fluid dynamics, and chemical reactor theory. The proposed method provides a highly accurate, fourth-order convergent numerical solution, achieved by implementing a compact finite difference method based on non-polynomial spline in tension approximations. Notably, the discretisation involves the use of half-step grid points, eliminating the need to modify the method at singularities while dealing with singular boundary value problems. The article also includes a comprehensive convergence analysis that demonstrates the theoretical order of convergence of the method. Additionally, the results obtained from solving six diverse mixed boundary value problems, including Burger's equation and Bratu-type equation, are compared to showcase the superiority and effectiveness of the proposed method.

Numerical solution of the mixed boundary value problem for the heat equation in two-dimensional domains of complex shape

A finite-difference computational algorithm is proposed for solving a mixed boundary value problem for heat equation given in a two-dimensional domains of complex shape. To solve the problem, generalized curvilinear coordinates are used. The physical domain is mapped to the computational domain (unit square) in the space of generalized coordinates. The original problem is written in curvilinear coordinates and approximated on a uniform grid in the computational domain. The obtained results are mapped on a non-uniform boundary-fitted difference grid in the physical domain. The second-order approximations of mixed Neumann – Dirichlet boundary conditions are constructed. The results of computational experiments are presented. The second order of accuracy of the presented computational algorithm is confirmed.

Investigation of an Exact Solution of a Mixed Boundary Value Problem Using the Residue Method

Investigation of an Exact Solution of a Mixed Boundary Value Problem Using the Residue Method

Thermoelastic fracture of two-dimensional hexagonal quasicrystal media weakened by a penny-shaped crack subjected to uniformly antisymmetric heat fluxes

This study focuses on the thermoelastic fracture behavior of quasicrystals, a class of solids that exhibit unique properties between traditional crystals and amorphous materials. The complex atomic structure of quasicrystals poses challenges in understanding their fracture mechanisms under thermoelastic loading conditions. Central to our investigation is an analytical examination of a penny-shaped crack in an infinite three-dimensional body composed of two-dimensional hexagonal quasicrystals, where the upper and lower crack surfaces are applied to a pair of uniformly antisymmetric heat fluxes. This crack problem requires simultaneous consideration of thermal-phason-phonon multiple field coupling. According to the symmetry of field variables with respect to the crack plane, the thermal-phason-phonon coupled crack problem is transformed into a mixed boundary value problem in the upper half space. The extended displacement discontinuities, encompassing both phonon and phason displacement discontinuities, as well as the temperature discontinuity, are chosen as the basic unknown variables to construct the boundary integral-differential equations governing the mixed boundary value problem. Based on these boundary governing equations and Fabrikant's potential theory method, the problem with the crack surface subjected to uniform antisymmetric heat fluxes is solved. The solutions of thermal-phason-phonon fields on the crack plane, and in the full space are given in closed-form. Numerical results are employed to validate the obtained analytical solutions and visually illustrate the spatial distribution of thermal-phonon-phason coupling fields in the vicinity of the crack. The study provides fundamental insights into the behavior of cracks in quasicrystals under thermal loading, with potential implications for the design of new materials and structures.

Existence result for a nonlinear mixed boundary value problem for the heat equation

In this paper we study the existence of solutions in parabolic Schauder space of a nonlinear mixed boundary value problem for the heat equation in a perforated domain. From a given regular open set Ω⊆Rn we remove a cavity ω⊆Ω. On the exterior boundary of Ω∖ω‾ we prescribe a Neumann boundary condition, while on the interior boundary we set a nonlinear Robin-type condition. Under suitable assumptions on the data and by means of Leray Schauder Fixed-Point Theorem, we prove the existence of (at least) one solution u∈C01+α2;1+α([0,T]×(Ω‾∖ω)).

Boundary Domain Integral Equations for Variable Coefficient Mixed BVP in 2D Unbounded Domain

In this paper, the direct segregated Boundary Domain Integral Equations (BDIEs) for the Mixed Boundary Value Problems (MBVPs) for a scalar second order elliptic Partial Differential Equation (PDE) with variable coefficient in unbounded (exterior) 2D domain is considered. Otar Chkadua, Sergey Mikhailov and David Natroshvili formulated both the interior and exterior 3D domain of the direct segregated systems of BDIEs for the MBVPs for a scalar second order elliptic PDE with a variable coefficients. On the other hand Sergey Mikhailov and Tamirat Temesgen formulated only the interior 2D domain of the direct segregated systems of BDIEs for the MBVPs for a scalar second order divergent elliptic PDE with a variable coefficients. However, in this paper we formulated the exterior 2D domain of the direct segregated systems of BDIEs for the MBVPs for a scalar second order divergent elliptic PDE with a variable coefficients. The aim of this work is to reduce the MBVPs to some direct segregated BDIEs with the use of an appropriate parametrix (Levi function). We examine the characteristics of corresponding parametrix-based integral volume and layer potentials in some weighted Sobolev spaces, as well as the unique solvability of BDIEs and their equivalence to the original MBVPs. This analysis is based on the corresponding properties of the MBVPs in weighted Sobolev spaces that are proved as well.

Analytical solutions to Mode I penny-shaped crack problems in two-dimensional hexagonal quasicrystals with piezoelectric effect

The paper studies a penny-shaped crack in an infinite three-dimensional body of two-dimensional hexagonal quasicrystal media with piezoelectric effect. The crack surfaces are applied combined electric and normal phonon loadings. Such a Model I crack problem is transformed into a mixed boundary value problem in the upper half-space, which is analytically solved using Fabrikant's potential theory method. The boundary integral-differential equations governing Model I crack problems are presented for two-dimensional hexagonal piezoelectric quasicrystals. The normal phonon displacement discontinuity and electric potential discontinuity across crack surfaces are taken as the unknown variables of boundary governing equations. Analytical solutions of all field variables are derived not only for the crack plane but also for the full space. Solutions in integral form are provided for the penny-shaped crack under arbitrarily distributed electric and normal phonon loadings. Closed-form solutions in terms of elementary functions are given for concentrated point loadings and uniformly distributed loadings, respectively. Key fracture mechanics parameters, such as crack surface extended displacements (i.e., normal phonon displacement, electric potential), crack tip extended stresses (i.e., normal phonon stress, electric displacement) distribution, and corresponding extended stress intensity factors, are clearly derived. Numerical results are utilized to verify the present analytical solutions and graphically illustrate the distribution of phonon-phason-electric coupling fields around the crack. The present solution can serve as a benchmark for both experimental and numerical investigations.

Some basic theorems on the modified Green–Lindsay model of linear generalized thermoelasticity

The present paper is concerned with a theory of thermoelasticity based upon a newly developed temperature-rate-dependent thermoelastic model that was proposed by Yu et al. [Meccanica. 2018;53:2543–2554], which has eliminated the discontinuity observed in the displacement field during transient motion. The basic governing equations for a homogeneous and isotropic medium are formulated to deduce the significant theorems. Firstly, we prove the uniqueness theorem concerning a mixed initial and boundary value problem within the framework of linear thermoelasticity. Furthermore, we establish a reciprocity theorem utilizing the Laplace transform, which provides insight into the mutual influence between different aspects of the system under consideration. Consequently, we proved the Somigliana and Green theorems. This generalized theorem extends the applicability and versatility of the classical theorems and also assists the issue of thermal and mechanical fields as well as is beneficial for geophysics and mining.

Numerical solution for the interior electric displacement of the moving DBY‐PS model for semi‐permeable cracked piezoelectric material

AbstractIn this study, we analysed a moving crack at the interface of an infinitely long piezoelectric bilayer using the Dugdale–Barenblatt yield (DBY) model and the polarisation saturation (PS) model. To model the moving crack problem, a Yoffe‐type crack moves at a constant subsonic speed on the interface of an infinitely long piezoelectric bilayer. The crack faces are assumed to be semi‐permeable, and at the boundary of the bilayer, in‐plane electrical and out‐of‐plane mechanical stresses are applied. Due to the application of electro‐mechanical loads, cracks propagate, mechanical yielding zones and electric saturation zones are developed. To arrest the crack from further propagation, mechanical yield stress and saturation electric displacement are applied at the developed zones. To address this problem analytically and numerically, the mixed boundary value problem is transformed into a set of coupled Fredholm integral equations (FIEs) of the second kind using the Fourier transform and the Copson method. The closed‐form analytical expressions for the length of the electrical saturation zone (ESZ), whether longer, shorter or equal to the mechanical yielding zone (MYZ), show dependence on external electro‐mechanical loads under semi‐permeable crack conditions. The algorithm to solve the electric crack condition parameter (ECCP) has been defined using numerical discretization and the bisection method. Illustrative examples demonstrate the proposed technique's effectiveness and suitability for Yoffe‐type moving cracks. The numerical results show the convergence of the ECCP. Furthermore, the numerical results show how mechanical and electrical zone lengths and energy release rate (ERR) are affected by electrical and mechanical loads, strip thickness and crack velocity. In addition, the size of the mechanical yielding zone is consistently promoted by electrical load, while the promotion or prevention of the electrical saturation zone by mechanical load depends on the relative sizes of the nonlinear zones.

Mixed Boundary Value Problem for a Monotone Equation with a Lower Order Term and Point Sources on the Right Side

Mixed Boundary Value Problem for a Monotone Equation with a Lower Order Term and Point Sources on the Right Side

Soap bubbles and convex cones: optimal quantitative rigidity

We consider a class of recent rigidity results in a convex cone Σ ⊆ R N \Sigma \subseteq \mathbb {R}^N . These include overdetermined Serrin-type problems for a mixed boundary value problem relative to Σ \Sigma , Alexandrov’s soap bubble-type results relative to Σ \Sigma , and Heintze-Karcher’s inequality relative to Σ \Sigma . Each rigidity result is obtained here by means of a single integral identity and holds true under weak integral overdeterminations in possibly non-smooth cones. Optimal quantitative stability estimates are obtained in terms of an L 2 L^2 -pseudodistance. In particular, the optimal stability estimate for Heintze-Karcher’s inequality is new even in the classical case Σ = R N \Sigma = \mathbb {R}^N . Stability bounds in terms of the Hausdorff distance are also provided. Several new results are established and exploited, including a new Poincaré-type inequality for vector fields whose normal components vanish on a portion of the boundary and an explicit (possibly weighted) trace theory – relative to the cone Σ \Sigma – for harmonic functions satisfying a homogeneous Neumann condition on the portion of the boundary contained in ∂ Σ \partial \Sigma . We also introduce new notions of uniform interior and exterior sphere conditions relative to the cone Σ ⊆ R N \Sigma \subseteq \mathbb {R}^N , which allow to obtain (via barrier arguments) uniform lower and upper bounds for the gradient in the mixed boundary value-setting. In the particular case Σ = R N \Sigma =\mathbb {R}^N , these conditions return the classical uniform interior and exterior sphere conditions (together with the associated classical gradient bounds of the Dirichlet setting).

Tilted spacetime positive mass theorem with arbitrary ends

In this paper, we prove the spacetime positive mass theorem for asymptotically flat spin initial data sets with arbitrary ends and a non-compact boundary. Moreover, we demonstrate a quantitative shielding theorem, subject to the tilted boundary dominant energy condition. Our results are established by solving a mixed boundary value problem for the Dirac-Witten operator with a Callias potential.

A generalized approach to solving the mixed boundary value problem governing self-diffusiophoresis

A generalized approach to solving the mixed boundary value problem governing self-diffusiophoresis

On the construction of the Stokes flow in a domain with cylindrical ends

Based on existence results for the Stokes operator and its solution properties in manifolds with cylindrical ends by Große et al. and Kohr et al., the Stokes flow in a three‐dimensional compact domain with circular openings through which the fluid enters and leaves through unbounded cylindrical pipes the Stokes flow is modeled as a mixed boundary value problem whereas in the cylindrical ends, the velocities and pressures are constant on every straight line in the cylindrical directions with initial values from the openings of . These values equal the velocities and pressures which are obtained from the mixed boundary values' solution in at the openings .

A High-Resolution Exponential Spline Method and Its Convergence Analysis for Two-Point Mixed Boundary Value Problems

In this paper, we present a new numerical method for effectively solving general second-order ordinary differential equations with mixed boundary conditions. Our approach utilizes a quasi-variable mesh enabling us the flexibility to adapt the mesh density according to different boundary layer problems. By employing a discretization technique that incorporates the construction of exponential spline, we achieve third-order accuracy at internal grid points, while boundary points exhibit fourth-order accuracy. Since the method is based on off-step grid points, it eliminates the need to modify the method while employing it to singular problems. We provide computational results to various numerical problems which arise in different physical phenomena like Burger’s equation and Sturm–Liouville equation. A comparison with recent findings underscores the superior performance of our method.

BOUNDARY-DOMAIN INTEGRAL EQUATIONS FOR VARIABLE-COEFFICIENT HELMHOLTZ BVPs IN 2D

In this paper, we construct boundary-domain integral equations (BDIEs) of the Dirichlet and mixed boundary value problems for a two-dimensional variable-coefficient Helmholtz equation. Using an appropriate parametrix, these problems are reduced to several BDIE systems. It is shown that the BVPs and the formulated BDIE systems are equivalent. Fredholm properties and unique solvability and invertibility of BDIE systems are investigated in appropriate Sobolev spaces.

Схлопывание присоединенной каверны при малых числах Фруда после отрывного удара плавающего эллиптического цилиндра

The plane problem of the collapse of an attached cavity formed as a result of a separation impact of an elliptical cylinder under the free surface of a heavy fluid is considered. It is assumed that after the impact the cylinder moves horizontally with constant velocity. At small Froude numbers, which correspond to small velocities of the cylinder, the free boundaries of the liquid are slightly perturbed and the process of cavity collapse is mainly reduced to the study of dynamics of breakaway points. The solution to this problem is constructed using asymptotic decompositions on a small parameter, which is the Froude number. In the main asymptotic approximation, the mixed boundary value problem of potential theory with one-sided restrictions on the body surface is formulated. On its basis, the dynamics of detachment points, the shape of a thin cavities, and the reaction of the environment to the body are determined. The results obtained can be used for solving practical problems of marine and ship hydrodynamics.

Mixed boundary value problems involving Sturm–Liouville differential equations with possibly negative coefficients

This paper is devoted to the study of a mixed boundary value problem for a complete Sturm–Liouville equation, where the coefficients can also be negative. In particular, the existence of infinitely many distinct positive solutions to the given problem is obtained by using critical point theory.

An alternative potential method for mixed steady‐state elastic oscillation problems

We consider an alternative approach to investigate three‐dimensional exterior mixed boundary value problems (BVPs) for the steady‐state oscillation equations of the elasticity theory for isotropic bodies. The unbounded domain occupied by an elastic body, , has a compact boundary surface , which is divided into two disjoint parts, the Dirichlet part and the Neumann part , where the displacement vector (the Dirichlet‐type condition) and the stress vector (the Neumann‐type condition) are prescribed, respectively. Our new approach is based on the classical potential method and has several essential advantages compared with the existing approaches. We look for a solution to the mixed BVP in the form of a linear combination of the single‐layer and double‐layer potentials with densities supported on the Dirichlet and Neumann parts of the boundary, respectively. This approach reduces the mixed BVP under consideration to a system of boundary integral equations, which contain neither extensions of the Dirichlet or Neumann data nor the Steklov–Poincaré‐type operator involving the inverse of a special boundary integral operator, which is not available explicitly for arbitrary boundary surface. Moreover, the right‐hand sides of the resulting boundary integral equations system are vector functions coinciding with the given Dirichlet and Neumann data of the problem in question. We show that the corresponding matrix integral operator is bounded and coercive in the appropriate ‐based Bessel potential spaces. Consequently, the operator is invertible, which implies unconditional unique solvability of the mixed BVP in the class of vector functions belonging to the Sobolev space and satisfying the Sommerfeld–Kupradze radiation conditions at infinity. We also show that the obtained matrix boundary integral operator is invertible in the ‐based Besov spaces and prove that under appropriate boundary data a solution to the mixed BVP possesses ‐Hölder continuity property in the closed domain with , where is an arbitrarily small number.

On one method for solving a mixed boundary value problem for a parabolic type equation using operators AT <sub>ƛ, J</sub>

A new method for obtaining a generalized solution of a mixed boundary value problem for a parabolic equation with boundary conditions of the third kind and a continuous initial condition is proposed. Generalized functions are understood in the sense of a sequential approach. The representative of the class of sequences, which is a generalized function, is obtained using the function interpolation operator, constructed using solutions of the Cauchy problem. The solution is obtained in the form of a series that converges uniformly inside the domain of the solution.