On the Solution of the Stokes Equation on Regions with Corners

The detailed behavior of solutions to Stokes equations on regions with corners has been historically difficult to characterize. The solutions to Stokes equations on regions with corners are known to develop singularities in the vicinity of corners; in particular, the solutions are known to have infinite oscillations along almost every ray that meet at the corner. While the nature of singularities for the differential equation have been analyzed in great detail, very little is known about the nature of singularities for the corresponding integral equations. In this paper, we observe that, when the Stokes equation is formulated as a boundary integral equation, the solutions are representable by rapidly convergent series of the form , where t is the distance from the corner and the parameters μ_j, β_j are real, and are determined via an explicit formula depending on the angle at the corner. In addition to being analytically perspicuous, these representations lend themselves to the construction of highly accurate and efficient numerical discretizations, significantly reducing the number of degrees of freedom required for the solution of the corresponding integral equations. The results are illustrated by several numerical examples. © 2020 Wiley Periodicals LLC     

Compatible wavelet-Galerkin method to solve stokes interior problem with circular boundary

This article addresses Neumann boundary value interior problem of Stokes equations with circular boundary. By using natural boundary element method, the Stokes interior problem is reduced into an equivalent natural integral equation with a hyper-singular kernel, which is viewed as Hadamard finite part. Based on trigonometric wavelet functions, the compatible wavelet space is constructed so that it can serve as Galerkin trial function space. In proposed compatible wavelet-Galerkin method, the simple and accurate computational formulae of the entries in stiffness matrix are obtained by singularity removal technique. It is also proved that the stiffness matrix is almost a block diagonal matrix, and its diagonal sub-blocks all are both symmetric and circulant submatrices. These good properties indicate that a 2 ^J+3 × 2 ^J+3 stiffness matrix can be determined only by its 2 ^J + 3J + 1 entries. It greatly decreases the computational complexity. Some error estimates for the compatible wavelet-Galerkin projection solutions are established. Finally, several numerical examples are given to demonstrate the validity of the proposed approach.     

The resolvent problem for the stokes system in exterior domains: An elementary approach

We consider the Stokes system with resolvent parameter in an exterior domain: under Dirichlet boundary conditions. Here Ω is a bounded domain with C² boundary, and [λ??\] ? [∞, 0], ν >0. Using the method of integral equations, we are able to construct solutions ( u , π) in L^p spaces. Our approach yields an integral representation of these solutions. By evaluating the corresponding integrals, we obtain L^p estimates that imply in particular that the Stokes operator on exterior domains generates an analytic semigroup in L^p.     

Exact solutions of the modified Korteweg-de Vries equation

We use the inverse scattering method to obtain a formula for certain exact solutions of the modified Korteweg-de Vries (mKdV) equation. Using matrix exponentials, we write the kernel of the relevant Marchenko integral equation as W( x + y;t ) = Ce ^{- ( x + y )A} e^{8A3 t} B\Omega \left( {x + y;t} \right) = Ce^{ - \left( {x + y} \right)A} e^{8A^3 t} BB, where the real matrix triplet (A,B,C) consists of a constant p×p matrix A with eigenvalues having positive real parts, a constant p×1 matrix B, and a constant 1× p matrix C for a positive integer p. Using separation of variables, we explicitly solve the Marchenko integral equation, yielding exact solutions of the mKdV equation. These solutions are constructed in terms of the unique solution P of the Sylvester equation AP + PA = BC or in terms of the unique solutions Q and N of the Lyapunov equations A^°Q + QA = C^°C and AN + NA^° = BB^°, where B^°denotes the conjugate transposed matrix. We consider two interesting examples.     

Singular perturbations for the exterior three-dimensional slow viscous flow problem

In this paper the rigorous justification of the formal asymptotic expansions constructed by the method of matched inner and outer expansions is established for the three-dimensional steady flow of a viscous, incompressible fluid past an arbitrary obstacle. The justification is based on the series representation of the solution to the Navier-Stokes equations due to Finn, and it involves the reductions of various exterior boundary value problems for the Stokes and Oseen equations to boundary integral equations of the first kind from which existence as well as asymptotic error estimates for the solutions are deduced. In particular, it is shown that the force exerted on the obstacle by the fluid admits the asymptotic representation F = A₀ + A₁Re + O(Re² ln Re⁻¹) as the Reynolds number Re → 0⁺, where the vectors A₀ and A₁ can be obtained from the method of matched inner and outer expansions.     

Periodic solutions of the Navier–Stokes equations in a perturbed half‐space and an aperture domain

We shall construct a periodic strong solution of the Navier–Stokes equations for some periodic external force in a perturbed half‐space and an aperture domain of the dimension n?3. Our proof is based on L^p–L^q estimates of the Stokes semigroup. We apply L^p–L^q estimates to the integral equation which is transformed from the original equation. As a result, we obtain the existence and uniqueness of periodic strong solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Boundary integral solution of the two-dimensional heat equation

Here we consider initial boundary value problems for the heat equation by using the heat potential representation for the solution. Depending on the choice of the representation we are led to a solution of the various boundary integral equations. We discuss the solvability of these equations in anisotropic Sobolev spaces. It turns out that the double-layer heat potential D and its spatial adjoint D′ have smoothing properties similar to the single-layer heat operator. This yields compactness of the operators D and D′. In addition, for any constant c ≠ 0, cI + D′ and cI + D′ are isomorphisms. Based on the coercivity of the single-layer heat operator and the above compactness we establish the coerciveness of the hypersingular heat operator. Moreover, we show an equivalence between the weak solution and the various boundary integral solutions. As a further application we describe a coupling procedure for an exterior initial boundary value problem for the non-homogeneous heat equation.     

Global existence for critical power ynag-mills-higgs equations in R3+1

We prove global existence for both smooth and finite energy solutions of the Yang-Mills-Higgs equations with a critical-power Higgs self-interaction in R ³⁺¹. The phrase “critical-power” denotes here a Higgs self-interaction behaving as φ ⁶ for large φ.     

Feedback boundary stabilization of the three-dimensional incompressible Navier–Stokes equations

We study the local stabilization of the three-dimensional Navier–Stokes equations around an unstable stationary solution w, by means of a feedback boundary control. We first determine a feedback law for the linearized system around w. Next, we show that this feedback provides a local stabilization of the Navier–Stokes equations. To deal with the nonlinear term, the solutions to the closed loop system must be in H^{3/2+ε,3/4+ε/2}(Q), with 0<ε. In [V. Barbu, I. Lasiecka, R. Triggiani, Boundary stabilization of Navier–Stokes equations, Mem. Amer. Math. Soc. 852 (2006); V. Barbu, I. Lasiecka, R. Triggiani, Abstract settings for tangential boundary stabilization of Navier–Stokes equations by high- and low-gain feedback controllers, Nonlinear Anal. 64 (2006) 2704–2746], such a regularity is achieved with a feedback obtained by minimizing a functional involving a norm of the state variable strong enough. In that case, the feedback controller cannot be determined by a well posed Riccati equation. Here, we choose a functional involving a very weak norm of the state variable. The compatibility condition between the initial state and the feedback controller at t=0, is achieved by choosing a time varying control operator in a neighbourhood of t=0.     

Onl p solutions of discrete Volterra equations

Summary In the paper a discrete analog to the Volterra nonlinear integral equation is discussed. Weighted norms are used to find sufficient conditions that all solutions of such equations are elements of anl ^p space.Some generalizations ofl ^p spaces are also considered and the corresponding sufficient conditions are established.     

On Some Second Order Differential Equations with Convergent Solutions

Via a special integral transformation, asymptotic integration results for ordinary differential equations are used to establish accurate asymptotic developments for radial solutions of the elliptic equation Δu + K(|x|)e ^u = 0, |x| > x ₀ > 0, in the bidimensional case.     

Solutions to the polynomial Dirac equations on unbounded domains in Clifford analysis

In this paper we study polynomial Dirac equation p(??)f = 0 including (?? ? λ)f = 0 with complex parameter λ and ??^kf = 0(k?1) as special cases over unbounded subdomains of ?^{n + 1}. Using the Clifford calculus, we obtain the integral representation theorems for solutions to the equations satisfying certain decay conditions at infinity over unbounded subdomains of ?^{n + 1}. Copyright © 2010 John Wiley & Sons, Ltd.     

The large time behavior of solutions to 3D Navier–Stokes equations with nonlinear damping

This paper studies the Cauchy problem of the 3D Navier–Stokes equations with nonlinear damping term | u | ^β?1u (β ≥ 1). For β ≥ 3, we derive a decay rate of the L²‐norm of the solutions. Then, the large time behavior is given by comparing the equation with the classic 3D Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.     

The Decoupled Potential Integral Equation for Time‐Harmonic Electromagnetic Scattering

We present a new formulation for the problem of electromagnetic scattering from perfect electric conductors. While our representation for the electric and magnetic fields is based on the standard vector and scalar potentials A ,φ in the Lorenz gauge, we establish boundary conditions on the potentials themselves rather than on the field quantities. This permits the development of a well‐conditioned second‐kind Fredholm integral equation that has no spurious resonances, avoids low‐frequency breakdown, and is insensitive to the genus of the scatterer. The equations for the vector and scalar potentials are decoupled. That is, the unknown scalar potential defining the scattered field, φ^scat, is determined entirely by the incident scalar potential φ^inc. Likewise, the unknown vector potential defining the scattered field, A ^scat is determined entirely by the incident vector potential A^inc. This decoupled formulation is valid not only in the static limit but for arbitrary ω ≥ 0$. © 2016 Wiley Periodicals, Inc.     

Stokes‐Fourier and acoustic limits for the Boltzmann equation: Convergence proofs

We establish a Stokes‐Fourier limit for the Boltzmann equation considered over any periodic spatial domain of dimension two or more. Appropriately scaled families of DiPerna‐Lions renormalized solutions are shown to have fluctuations that globally in time converge weakly to a unique limit governed by a solution of Stokes‐Fourier motion and heat equations provided that the fluid moments of their initial fluctuations converge to appropriate L² initial data of the Stokes‐Fourier equations. Both the motion and heat equations are both recovered in the limit by controlling the fluxes and the local conservation defects of the DiPerna‐Lions solutions with dissipation rate estimates. The scaling of the fluctuations with respect to Knudsen number is essentially optimal. The assumptions on the collision kernel are little more than those required for the DiPerna‐Lions theory and that the viscosity and heat conduction are finite. For the acoustic limit, these techniques also remove restrictions to bounded collision kernels and improve the scaling of the fluctuations. Both weak limits become strong when the initial fluctuations converge entropically to appropriate L² initial data. © 2001 John Wiley & Sons, Inc.     

Viscosity splitting method for three dimensional Navier-Stokes equations

Three dimensional initial boundary value problem of the Navier-Stokes equation is considered. The equation is split in an Euler equation and a non-stationary Stokes equation within each time step. Unlike the conventional approach, we apply a non-homogeneous Stokes equation instead of homogeneous one. Under the hypothesis that the original problem possesses a smooth solution, the estimate of theH ^s+1 norm, 0≦s<3/2, of the approximate solutions and the order of theL ² norm of the errors is obtained. This work was supported by the Science Foundation of Academia Sinica under grant (84)-103.     

The self-similar solutions of the Tricomi-type equations

This paper analyses the properties of the family of self-similar solutions of the generalized Tricomi equation u_tt - t^2k Du = 0 (2k ? \mathbbN)u_{{tt}} - t^{{2k}} \Delta u = 0\,(2k \in {{\mathbb{N}}}) in the domain \mathbbR ₊ ^{1 + n}{{\mathbb{R}}}_{ + }^{{1 + n}} by considering initial conditions on the functions and their derivatives, posed as the Cauchy problem with homogeneous initial data. For specific values of the power k ( = 1/2 or = 3/2) and n = 1 this problem has applications in the aerodynamics of airfoils operating in transonic flows of perfect or dense gases, respectively. An integral transformation is suggested and used to represent the solutions of the Cauchy problem with homogeneous initial functions in terms of fundamental solutions of the classical wave equation (the case k = 0). Then the Cauchy problem with homogeneous initial functions for the wave equation in \mathbbR^{1 + n}{{\mathbb{R}}}^{{1 + n}} is solved. These results are used to derive estimates of the upper bound for solutions’ size and to obtain the asymptotics for self-similar solutions of the wave equation and of the Tricomi-type equation in the neighbourhood of their light cones.     

On the relationship of some Fredholm integral equations of the first kind to a family of boundary value problems

It is shown that the solution of a class of forced convection heat transfer problems can be used as a vehicle for exhibiting a correspondence between certain boundary value problems and their associated integral equations. If the solution of the boundary value problem is known then the solution of the integral equation can be found by a simple calculation. This generalizes the relationship of certain solutions of Laplace's equation to integral equations having logarithmic or Cauchy kernels. When the boundary value problems have separable solutions, the solution of the integral equation can be reduced to the verification of a set of identities μ_n ∫_EP(t) φ_n(t) k(x ? t) dt = φ_n(x), x?E, where the {φ_n} form an orthonormal set on E.Once the method of solution has been derived from the physical approach it can be put on a firm mathematical basis.     

The exterior Dirichlet and the interior Neumann boundary value problems for the scalar Oseen equation

The scalar Oseen equation represents a linearized form of the Navier Stokes equations. We present an explicit potential theory for this equation and solve the exterior Dirichlet and interior Neumann boundary value problems via a boundary integral equations method in spaces of continuous functions on a C²-boundary, extending the classical approach for the isotropic selfadjoint Laplace operator to the anisotropic non-selfadjoint scalar Oseen operator. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Painlevé Classification of a Class of Differential Equations of the Second Order and Second Degree

In this paper we construct all Painlevé-type differential equations of the form (d²y/dx²)² = F(x,y,dy/dx), where F is rational in y and y′=dy/dx, locally analytic in x, and not a perfect square. No further simplifying assumptions are made, but it is found that the absence of a term linear in y″ in the class of equations under investigation forces F to be a polynomial in y and y′. We find exactly six distinct classes of second-degree Painlevé equations, denoted SD-I,??,SD-VI, some of which further subdivide into canonical subcases. Only the first three classes (or at least equations transformable to the first three classes) and part of the sixth have appeared previously in the literature, especially the work of Chazy and Bureau. The fourth and fifth classes are new. The unified treatment of SD-I, which we call the “master Painlevé equation,” is new. Complete solutions are given in terms of the classical Painlevé transcendents, elliptic functions, or solutions of linear equations. In an appendix, it is shown that a class of second-degree equations generalizing the Appell equation can always be reduced to a second-order linear equation.