Viscous Flows in a Half Space Caused by Tangential Vibrations on Its Boundary

The paper is devoted to the studies of viscous flows caused by a vibrating boundary. The fluid domain is a half‐space, its boundary is a nondeformable plane that exhibits purely tangential vibrations. Such a simple geometrical setting allows us to study general boundary velocity fields and to obtain general results. From a practical viewpoint, such boundary conditions may be seen as the tangential vibrations of the material points of a stretchable plane membrane. In contrast to the classical boundary layer theory, we aim to build a global solution. To achieve this goal we employ the Vishik–Lyusternik approach, combined with two‐timing and averaging methods. Our main result is: we obtain a uniformly valid in the whole fluid domain approximation to the global solutions. This solution corresponds to general boundary conditions and to three different settings of the main small parameter. Our solution always include the inner part and outer part that both contain oscillating and non‐oscillating components. It is shown that the nonoscillating outer part of the solution is governed either by the full Navier–Stokes equations or the Stokes equations (both with the unit viscosity) and can be interpreted as a steady or unsteady streaming. In contrast to the existing theories of a steady streaming, our solutions do not contain any secular (infinitely growing with the inner normal coordinate) terms. The examples of the spatially periodic vibrations of the boundary and the angular torsional vibrations of an infinite rigid disc are considered. These examples are still brief and illustrative, while the core of the paper is devoted to the adaptation of the Vishik–Lyusternik method to the development of the general theory of vibrational boundary layers.     

A simple and robust boundary treatment for the forced Korteweg–de Vries equation

In this paper, we propose a simple and robust numerical method for the forced Korteweg–de Vries (fKdV) equation which models free surface waves of an incompressible and inviscid fluid flow over a bump. The fKdV equation is defined in an infinite domain. However, to solve the equation numerically we must truncate the infinite domain to a bounded domain by introducing an artificial boundary and imposing boundary conditions there. Due to unsuitable artificial boundary conditions, most wave propagation problems have numerical difficulties (e.g., the truncated computational domain must be large enough or the numerical simulation must be terminated before the wave approaches the artificial boundary for the quality of the numerical solution). To solve this boundary problem, we develop an absorbing non-reflecting boundary treatment which uses outward wave velocity. The basic idea of the proposing algorithm is that we first calculate an outward wave velocity from the solutions at the previous and present time steps and then we obtain a solution at the next time step on the artificial boundary by moving the solution at the present time step with the velocity. And then we update solutions at the next time step inside the domain using the calculated solution on the artificial boundary. Numerical experiments with various initial conditions for the KdV and fKdV equations are presented to illustrate the accuracy and efficiency of our method.     

Exponential spline solutions for a class of two point boundary value problems over a semi-infinite range

In this paper, we developed numerical methods of order O(h ²) and O(h ⁴) based on exponential spline function for the numerical solution of class of two point boundary value problems over a Semi-infinite range. The present approach gives better approximations over all the existing finite difference methods. Properties of the infinite linear system are established. Convergence analysis and a bound on the approximate solution are discussed. Test problem with various kinds of boundary conditions is included to illustrate the practical usefulness and superiority of our methods.     

Oseen and Stokes asymptotics for the problem on stationary motion of two immiscible liquids

The main result is an asymptotic formula for a solution to the conjugation problem for the Navier-Stokes equations describing the slow motion of two immiscible liquids such that one of them occupies a bounded domain Ω₁ ⊂ ℝ³, whereas the other occupies the exterior domain Ω₂=ℝ⁴∖Ω. Such a formula was obtained for a solution to the exterior problem with sticking conditions on the boundary in the works of Fischer, Hsiao, and Wendland. The result obtained is applied to the proof of the solvability of a free-boundary problem describing a uniform drop in an infinite liquid. Bibliography: 10 titles. Translated fromProblemy Matematicheskogo Analiza, No. 16. 1997, pp. 208–238.     

Fast solvers for 3D Poisson equations involving interfaces in a finite or the infinite domain

In this paper, numerical methods are proposed for Poisson equations defined in a finite or infinite domain in three dimensions. In the domain, there can exists an interface across which the source term, the flux, and therefore the solution may be discontinuous. The existence and uniqueness of the solution are also discussed. To deal with the discontinuity in the source term and in the flux, the original problem is transformed to a new one with a smooth solution. Such a transformation can be carried out easily through an extension of the jumps along the normal direction if the interface is expressed as the zero level set of a three-dimensional function. An auxiliary sphere is used to separate the infinite region into an interior and exterior domain. The Kelvin's inversion is used to map the exterior domain into an interior domain. The two Poisson equations defined in the interior and the exterior written in spherical coordinates are solved simultaneously. By choosing the mesh size carefully and exploiting the fast Fourier transform, the resulting finite difference equations can be solved efficiently. The approach in dealing with the interface has also been used with the artificial boundary condition technique which truncates the infinite domain. Numerical results demonstrate second order accuracy of our algorithms.     

Nonlocal Boundary Conditions for Higher–Order Parabolic Equations

This work deals with the efficient numerical solution of the two–dimensional one–way Helmholtz equation posed on an unbounded domain. In this case one has to introduce artificial boundary conditions to confine the computational domain. Here we construct with the Z –transformation so–called discrete transparent boundary conditions for higher–order parabolic equations schemes. These methods are Padé “Parabolic” approximations of the one–way Helmholtz equation and frequently used in integrated optics and (underwater) acoustics. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Steady transonic shocks and free boundary problems in infinite cylinders for the Euler equations

We establish the existence and stability of multidimensional transonic shocks (hyperbolic‐elliptic shocks) for the Euler equations for steady compressible potential fluids in infinite cylinders. The Euler equations, consisting of the conservation law of mass and the Bernoulli law for velocity, can be written as a second order nonlinear equation of mixed elliptic‐hyperbolic type for the velocity potential. The transonic shock problem in an infinite cylinder can be formulated into the following free boundary problem: The free boundary is the location of the multidimensional transonic shock which divides two regions of C^1,α flow in the infinite cylinder, and the equation is hyperbolic in the upstream region where the C^1,α perturbed flow is supersonic. We develop a nonlinear approach to deal with such a free boundary problem in order to solve the transonic shock problem in unbounded domains. Our results indicate that there exists a solution of the free boundary problem such that the equation is always elliptic in the unbounded downstream region, the uniform velocity state at infinity in the downstream direction is uniquely determined by the given hyperbolic phase, and the free boundary is C^1,α, provided that the hyperbolic phase is close in C^1,α to a uniform flow. We further prove that, if the steady perturbation of the hyperbolic phase is C^2,α, the free boundary is C^2,α and stable under the steady perturbation. © 2003 Wiley Periodicals Inc.     

The boundary integral method for the steady rotating Navier–Stokes equations in exterior domain (I): the existence of solution

In this paper, we apply the boundary integral method to the steady rotating Navier–Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and a infinite domain, we obtain a coupled problem by the steady rotating Navier–Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence of solution in a convex set.     

Maximum Principles for Boundary-Degenerate Second Order Linear Elliptic Differential Operators

We prove weak and strong maximum principles, including a Hopf lemma, for C ² subsolutions to equations defined by linear, second-order, linear, elliptic partial differential operators whose principal symbols vanish along a portion of the domain boundary. The boundary regularity property of the C ² subsolutions along this boundary vanishing locus ensures that these maximum principles hold irrespective of the sign of the Fichera function. Boundary conditions need only be prescribed on the complement in the domain boundary of the principal symbol's vanishing locus. We obtain uniqueness and a priori maximum principle estimates for C ² solutions to boundary value and obstacle problems defined by these boundary-degenerate elliptic operators with partial Dirichlet or Neumann boundary conditions. We also prove weak maximum principles and uniqueness for W ^{1, 2} solutions to the corresponding variational equations and inequalities defined with the aide of weighted Sobolev spaces. The domain is allowed to be unbounded when the operator coefficients and solutions obey certain growth conditions.     

Regularity of weak solutions of Maxwell's equations with mixed boundary‐conditions

In this paper global H^s‐ and L^p‐regularity results for the stationary and transient Maxwell equations with mixed boundary conditions in a bounded spatial domain are proved. First it is shown that certain elements belonging to the fractional‐order domain of the Maxwell operator belong to H^s(Ω) for sufficiently small s > 0. It follows from this regularity result that H^s(Ω) is an invariant subspace of the unitary group corresponding to the homogeneous Maxwell equations with mixed boundary conditions. In the case that a possibly non‐linear conductivity is present a L^p‐regularity theorem for the transient equations is proved. Copyright © 1999 John Wiley & Sons, Ltd.     

Superconvergence of Solution Derivatives of the Shortley–Weller Difference Approximation to Elliptic Equations with Singularities Involving the Mixed Type of Boundary Conditions

This paper presents a superconvergence analysis for the Shortley–Weller finite difference approximation of second-order self-adjoint elliptic equations with unbounded derivatives on a polygonal domain with the mixed type of boundary conditions. In this analysis, we first formulate the method as a special finite element/volume method. We then analyze the convergence of the method in a finite element framework. An O(h ^1.5)-order superconvergence of the solution derivatives in a discrete H ¹ norm is obtained. Finally, numerical experiments are provided to support the theoretical convergence rate obtained.     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

L2 ASYMPTOTIC FOR THE EQUILIBRIUM SOLUTION OF THE F-M EQUATIONS IN THE SPACE-PERIODIC CASE

The L² exponetial asymptotical stability for the equilibrium solution of the F-M equations in the space-periodic case (n = 2) is considered. Under some assumptions on the external force, it can be shown that the weak solution of F-M equations with initial and boundary conditions in space-periodic case approaches the stationary solution of the system exponetially when time t goes to infinite.     

Convergence of euler-stokes splitting of the navier-stokes equations

We consider approximation by partial time steps of a smooth solution of the Navier-Stokes equations in a smooth domain in two or three space dimensions with no-slip boundary condition. For small k > 0, we alternate the solution for time k of the inviscid Euler equations, with tangential boundary condition, and the solution of the linear Stokes equations for time k, with the no-slip condition imposed. We show that this approximation remains bounded in H^2,p and is accurate to order k in L^p for p > ∞. The principal difficulty is that the initial state for each Stokes step has tangential velocity at the boundary generated during the Euler step, and thus does not satisfy the boundary condition for the Stokes step. The validity of such a fractional step method or splitting is an underlying principle for some computational methods. © 1994 John Wiley & Sons, Inc.     

Global strong solution to the three dimensional nonhomogeneous incompressible magnetohydrodynamic equations with density‐dependent viscosity and resistivity

In this paper, we consider an initial boundary value problem for the 3‐dimensional nonhomogeneous incompressible magnetohydrodynamic equations with density‐dependent viscosity and resistivity coefficients over a bounded smooth domain. Global in time unique strong solution is proved to exist when the L² norms of initial vorticity and current density are both suitably small with arbitrary large initial density, and the vacuum of initial density is also allowed. Finally, we revisit the Navier‐Stokes model without electromagnetic effect. We find that this initial boundary problem also admits a unique global strong solution under other conditions. In particular, we prove small kinetic‐energy strong solution exists globally in time, which extends the recent result of Huang and Wang.     

On a grid-method solution of the Laplace equation in an infinite rectangular cylinder under periodic boundary conditions

We study the Dirichlet problem for the Laplace equation in an infinite rectangular cylinder. Under the assumption that the boundary values are continuous and bounded, we prove the existence and uniqueness of a solution to the Dirichlet problem in the class of bounded functions that are continuous on the closed infinite cylinder. Under an additional assumption that the boundary values are twice continuously differentiable on the faces of the infinite cylinder and are periodic in the direction of its edges, we establish that a periodic solution of the Dirichlet problem has continuous and bounded pure second-order derivatives on the closed infinite cylinder except its edges. We apply the grid method in order to find an approximate periodic solution of this Dirichlet problem. Under the same conditions providing a low smoothness of the exact solution, the convergence rate of the grid solution of the Dirichlet problem in the uniform metric is shown to be on the order of O(h ² ln h ⁻¹), where h is the step of a cubic grid.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite element solutions for three‐dimensional elliptic boundary value problems on unbounded domains

The finite element (FE) solutions of a general elliptic equation ?div([a_ij] ??u) + u = f in an exterior domain Ω, which is the complement of a bounded subset of R ³, is considered. The most common approach to deal with exterior domain problems is truncating an unbounded subdomain Ω_∞, so that the remaining part Ω_B = Ω\Ω _∞ is bounded, and imposing an artificial boundary condition on the resulted artificial boundary Γ_a = Ω _∞ ∩ Ω _B. In this article, instead of discarding an unbounded subdomain Ω_∞ and introducing an artificial boundary condition, the unbounded domain is mapped to a unit ball by an auxiliary mapping. Then, a similar technique to the method of auxiliary mapping, introduced by Babu?ka and Oh for handling the domain singularities, is applied to obtain an accurate FE solution of this problem at low cost. This method thus does have neither artificial boundary nor any restrictions on f. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Artificial boundary conditions for viscoelastic flows

The steady three-dimensional exterior flow of a viscoelastic non-Newtonian fluid is approximated by reducing the corresponding nonlinear elliptic–hyperbolic system to a bounded domain. On the truncation surface with a large radius R, nonlinear, local second-order artificial boundary conditions are constructed and a new concept of an artificial transport equation is introduced. Although the asymptotic structure of solutions at infinity is known, certain attributes cannot be found explicitly so that the artificial boundary conditions must be constructed with incomplete information on asymptotics. To show the existence of a solution to the approximation problem and to estimate the asymptotic precision, a general abstract scheme, adapted to the analysis of coupled systems of elliptic–hyperbolic type, is proposed. The error estimates, obtained in weighted Sobolev norms with arbitrarily large smoothness indices, prove an approximation of order O(R^−2+ε), with any ε>0. Our approach, in contrast to other papers on artificial boundary conditions, does not use the standard assumptions on compactly supported right-hand side f, leads, in particular, to pointwise estimates and provides error bounds with constants independent of both R and f. Copyright © 2007 John Wiley & Sons, Ltd.     

Free Boundary Value Problems for Abstract Elliptic Equations and Applications

The free boundary value problems for elliptic differential-operator equations are studied. Several conditions for the uniform maximal regularity with respect to boundary parameters and the Fredholmness in abstract L _p -spaces are given. In application, the nonlocal free boundary problems for finite or infinite systems of elliptic and anisotropic type equations are studied.