Approximation of the solution of a fourth order boundary value problem with nonsmooth coefficient

Summary A class of generalized finite element methods for the approximate solution of fourth order two point boundary value problem with nonsmooth coefficient is presented. The methods are based on the use of problem dependentL-splines incorporating the nonsmoothness of the coefficient. Stability is proved and optimal error estimates in theH ² norm are derived for the solution and postprocessed solution, under the assumption that the coefficient is of bounded variation. The relation of these methods to mixed methods is discussed.This research was sponsored by the Senate Research Committee of Syracuse University, Syracuse, NY 13210     

On the Stokes and Navier-Stokes System for Domains with Noncompact Boundary in Lq-spaces

We consider the L^q-theory of weak solutions of the Stokes and Navier-Stokes equations in two classes of unbounded domains with noncompact boundary, namely in perturbed half spaces which are obtained by a perturbation of the half space IRⁿ, and in aperture domains consisting of two disjoint half spaces separated by a wall but connected by a hole (aperture) through this wall. The proofs rest on the cut-off procedure and a new multiplier approach to the half space problem. In an aperture domain we additionally prescribe either the flux through the wall or the pressure drop at infinity to single out a unique solution. The nonlinear problem is solved for sufficiently small data and requires q =n/2, n ≥ 3, to estimate the nonlinearity.     

The Saffman-Taylor Fingers in the Limit of Very Large Surface Tension

In long rectangular Hele-Shaw cells the Saffman-Taylor instability generates steady advancing fingers of the less viscous fluid. The structure of these fingers is determined by the solution of a system of nonlinear integrodifferential equations for two unknown functions of the curvilinear coordinate defined along the boundary of the finger. The surface tension appears in these equations through a dimensionless coefficient called k. I analyze the solution when this coefficient tends to infinity. In this limit, the inviscid fluid tends to fill almost completely the cell, except for a thin layer of viscous fluid left on the sides of thickness of order k^?½.     

A Study of the Direct Spectral Transform for the Defocusing Davey-Stewartson II Equation the Semiclassical Limit

The defocusing Davey-Stewartson II equation has been shown in numerical experiments to exhibit behavior in the semiclassical limit that qualitatively resembles that of its one-dimensional reduction, the defocusing nonlinear Schrödinger equation, namely the generation from smooth initial data of regular rapid oscillations occupying domains of space-time that become well-defined in the limit. As a first step to studying this problem analytically using the inverse scattering transform, we consider the direct spectral transform for the defocusing Davey-Stewartson II equation for smooth initial data in the semiclassical limit. The direct spectral transform involves a singularly perturbed elliptic Dirac system in two dimensions. We introduce a WKB-type method for this problem, proving that it makes sense formally for sufficiently large values of the spectral parameter k by controlling the solution of an associated nonlinear eikonal problem, and we give numerical evidence that the method is accurate for such k in the semiclassical limit. Producing this evidence requires both the numerical solution of the singularly perturbed Dirac system and the numerical solution of the eikonal problem. The former is carried out using a method previously developed by two of the authors, and we give in this paper a new method for the numerical solution of the eikonal problem valid for sufficiently large k. For a particular potential we are able to solve the eikonal problem in closed form for all k, a calculation that yields some insight into the failure of the WKB method for smaller values of k. Informed by numerical calculations of the direct spectral transform, we then begin a study of the singularly perturbed Dirac system for values of k so small that there is no global solution of the eikonal problem. We provide a rigorous semiclassical analysis of the solution for real radial potentials at k=0, which yields an asymptotic formula for the reflection coefficient at k=0 and suggests an annular structure for the solution that may be exploited when k ≠ 0 is small. The numerics also suggest that for some potentials the reflection coefficient converges pointwise as ɛ↓ 0 to a limiting function that is supported in the domain of k-values on which the eikonal problem does not have a global solution. It is expected that singularities of the eikonal function play a role similar to that of turning points in the one-dimensional theory. © 2019 Wiley Periodicals, Inc.     

On the Cauchy Problem of the Camassa-Holm Equation

The purpose of this paper is to investigate the Cauchy problem of the Camassa-Holm equation. By using the abstract method proposed and studied by T. Kato and priori estimates, the existence and uniqueness of the global solution to the Cauchy problem of the Camassa-Holm equation in L ^p frame under certain conditions are obtained. In addition, the continuous dependence of the solution of this equation on the linear dispersive coefficient contained in the equation is obtained.     

Convergence of the interpolated coefficient finite element method for the two‐dimensional elliptic sine‐Gordon equations

An interpolated coefficient finite element method is presented and analyzed for the two‐dimensional elliptic sine‐Gordon equations with Dirichlet boundary conditions. It is proved that the discretization scheme admits at least one solution, and that a subsequence of the approximation solutions converges to an exact solution in L²‐norm as the mesh size tends to zero. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Note on the flux condition and pressure drop in the resolvent problem of the Stokes system

In [4] H. Sohr and the author considered theL ^q-theory of the resolvent problem of the generalized Stokes system in an aperture domain. This type of unbounded domain consists of two disjoint half spaces which are separated from each other by a wall but connected by a hole (aperture) in this wall. Due to this geometry the flux of the velocity field through the hole and the pressure drop at infinity are important physical and mathematical quantities. In this note we show that in order to single out a unique solution of the resolvent problem we must prescribe the flux for largeq, but that for smallq neither the flux nor the pressure drop can be prescribed. Only if the dimension is greater than two there is a certain range of values ofq where we must prescribe either the flux or the pressure drop. As a limit case we also investigate strong solutions of the Stokes system.     

Two‐Grid method for nonlinear parabolic equations by expanded mixed finite element methods

In this article, we develop a two‐grid algorithm for nonlinear reaction diffusion equation (with nonlinear compressibility coefficient) discretized by expanded mixed finite element method. The key point is to use two‐grid scheme to linearize the nonlinear term in the equations. The main procedure of the algorithm is solving a small‐scaled nonlinear equations on the coarse grid and dealing with a linearized system on the fine space using the Newton iteration with the coarse grid solution. Error estimation to the expanded mixed finite element solution is analyzed in detail. We also show that two‐grid solution achieves the same accuracy as long as the mesh sizes satisfy H = O(h^1/2). Two numerical experiments are given to verify the effectiveness of the algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Comparison of computational aeroacoustics with theory using the benchmark case of sound waves excited from a semi‐infinite duct

In this paper the benchmark problem of radiation from a semi infinite duct is discussed. The duct is modeled by an infinitesimal thin wall. Comparisons are made between the analytic solutions of R. M. Munt and the 2^nd CAAWorkshop with our CAA‐code. In all cases a 3D axisymmetric approach is chosen to describe the propagation of a single mode. The procedure consists of two steps. Step 1 concerns a numerical simulation of the near‐field‐region. Step 2 calculates the farfield‐solution by means of the acoustic‐analogy of FfowcsWilliams and Hawkings. Since the wall is of infinitesimal thickness, the edge needs to be treated in a special way. Two different approaches are persecuted and the results are compared. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Variational Principle for the Ackerberg-O'Malley Resonance Problem

A new variational principle is proposed for determining the asymptotic expansion of the solution of the Ackerberg-O'Malley resonance problem [Stud. Appl. Math. 49:277–295 (1970)] to any order in ε. The method used yields new higher-order results not permitted by the technique of Grasman and Matkowsky [SIAM J. Appl. Math. 32:588–597 (1977)]. Explicit results using the method are reported to O(ε) and confirmed with asymptotic expansions of the exact solution; the O(1) results agree with those reported in the literature. In the case where the coefficient functions are analytic, an exact solution is presented. It is not difficult to perform the higher-order calculations using the proposed variational approach, in contrast to the current methods in use.     

Problem of determining the nonstationary potential in a hyperbolic-type equation

We solve the problem of determining the hyperbolic equation coefficient depending on two variables. Some additional information is given by the trace of the direct problem solution on the hyperplane x = 0. We estimate the stability of the solution of the inverse problem under study and prove the uniqueness theorem. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 156, No. 2, pp. 220–225, August, 2008.     

The presentation of explicit analytical solutions of a class of nonlinear evolution equations

In this paper, we introduce a function set Ω_m. There is a conjecture that an arbitrary explicit travelling-wave analytical solution of a real constant coefficient nonlinear evolution equation is necessarily a linear (or nonlinear) combination of the product of some elements in Ω_m. A widespread applicable approach for solving a class of nonlinear evolution equations is established. The new analytical solutions to two kinds of nonlinear evolution equations are described with the aid of the guess.     

Unique solvability of the inverse problem of determination of the leading coefficient in a parabolic equation

We study existence and uniqueness of the solution for the inverse problem of determination of the unknown coefficient ϱ(t) multiplying u _t in a nondivergence parabolic equation. As additional information, the integral of the solution over the domain of space variables with some given weight function is specified. The coefficients of the equation depend both on time and on the space variables.     

Convergence analysis of an hp finite element method for singularly perturbed transmission problems in smooth domains

We consider a two‐dimensional singularly perturbed transmission problem with two different diffusion coefficients, in a domain with smooth (analytic) boundary. The solution will contain boundary layers only in the part of the domain where the diffusion coefficient is high and interface layers along the interface. Utilizing existing and newly derived regularity results for the exact solution, we prove the robustness of an hp finite element method for its approximation. Under the assumption of analytic input data, we show that the method converges at an “exponential” rate, provided the mesh and polynomial degree distribution are chosen appropriately. Numerical results illustrating our theoretical findings are also included. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

The Stokes equations with Fourier boundary conditions on a wall with asperities

We study the effect of the rugosity of a wall on the solution of the Stokes system complemented with Fourier boundary conditions. We consider the case of small periodic asperities of size ε. We prove that the velocity field, pressure and drag, respectively, converge to the velocity field, pressure and drag of a homogenized Stokes problem, where a different friction coefficient appears. This shows that, contrarily to the case of Dirichlet boundary conditions, rugosity is dominant here. Copyright © 2001 John Wiley & Sons, Ltd.     

Solving the Inverse Generalized Problem of Vertical Seismic Profiling

The dynamic inverse seismics problem is considered in a generalized setting. We investigate whether the wave propagation problem in a vertically nonhomogeneous medium is well-posed. We show that the regular part of the solution is an L ₂ function and the inverse problem, i.e., the determination of the reflection coefficient, is thus reducible to minimizing the error functional. The gradient of the functional is obtained in explicit form from the conjugate problem, and approximate formulas for its evaluation are derived. A regularization algorithm for the solution of the inverse problem is considered; simulation results using various excitation sources are reported.     

Beyond Zudilin's conjectured q-analogue of Schmidt's problem

Using the methodology of (rigorous) experimental mathematics, we give a simple and motivated solution to Zudilin's question concerning a q-analogue of a problem posed by Asmus Schmidt about a certain binomial coefficients sum. Our method is based on two simple identities that can be automatically proved using the Zeilberger and q-Zeilberger algorithms. We further illustrate our method by proving two further binomial coefficient sums.     

Stream-vorticity NS equation ADI scheme to investigate properties of impulsively started translatory flows around the rotational circular cylinder without wall vorticity conditions

The properties of flow around a circular cylinder impulsively started into translatory and, rotatory motion with rotational parameter a less than or equal to 8.0 and Reynolds number Re=100 and 200 are investigated in the present paper. The vorticity and stream function N-S equations are adopted here, with a 2nd-order spatial and temporal accuracy ADI (alternating direction implicit) scheme. Moreover the wall vorticity obtain through the principle of conservation of the total computational domain vorticity is determined by domain vorticity and stream function, therefore, through the wall vorticity iteration, the wall vorticity condition is not fixed during the time step. And the present model results indicate: (1) when α>4.0, vortex street suppression is obvious for the computational period (t<60) for all the Re numbers here studied; (2) the higher the αnumber for the same Reynolds number, the slower the upper main vortex proceeds; (3) the maximum instantaneous transverse coefficient exceeds the limitation 4π.     

Inverse problem of finding the coefficient of a lower derivative in a parabolic equation on the plane

We study the existence and uniqueness of the solution of the inverse problem of finding an unknown coefficient b(x) multiplying the lower derivative in the nondivergence parabolic equation on the plane. The integral of the solution with respect to time with some given weight function is given as additional information. The coefficients of the equation depend on the time variable as well as the space variable.     

Accelerate overrelaxation methods for rank deficient linear systems

In this paper, we apply accelerated overrelaxation (AOR) methods to find the least square solution of minimal norm to the linear system

Ay=b

where is a matrix of rank r and . We first augment the system to a block 4×4 consistent system, and then split the augmented coefficient matrix by AOR subproper splitting. Intervals for the two relaxation parameters where the AOR iteration matrix is semiconvergent are presented. Also, we provide a method to compute the least square solution of minimal norm to the system.