An immersed interface method for anisotropic elliptic problems on irregular domains in 2D

This work is a generalization of the immersed interface method for discretization of a nondiagonal anisotropic Laplacian in 2D. This first‐order discretization scheme enforces weakly diagonal dominance of the numerical scheme whenever possible. A necessary and sufficient condition depending on the mesh size h for the existence of this scheme at an interior grid point is found in terms of the anisotropy matrix. A linear programming approach is introduced for finding the weights of the schemes. The method is tested with a parametrized family of anisotropic Poisson equations. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A one-sweep method for linear elliptic equations over irregular domains

The determination of the configuration of equilibrium in a number of problems in mechanics and structures such as torsion, deflection of elastic membranes,etc., involve the solution of variational problems defined over irregular regions. This problem, in turn, may be reduced to the solution of elliptic differential equations subject to boundary conditions. In this paper, we study a method for the solution of such a problem when the region is of irregular shape. The method consists in solving the problem over a larger, imbedding, rectangular domain subject to appropriate constraints such as to satisfy the conditions of the original problem at the boundary. In this paper, we introduce the constraints by considering appropriate factors on the Green's function of the auxiliary problem. A conveniently discretized version of the problem is then treated by invariant imbedding, yielding some earlier results plus some new ones, namely, a direct one-sweep procedure that minimizes storage requirements. In addition, the present solution appears to be very convenient when the solution is required at a limited number of points. The derivations are specialized to Laplace's equation, but the method can be applied readily to general systems of second-order elliptic equations with no essential modifications. Finally, the existence of the necessary matrices in the imbedding equations is established.     

Differential quadrature method for space-fractional diffusion equations on 2D irregular domains

In mathematical physics, the space-fractional diffusion equations are of particular interest in the studies of physical phenomena modelled by Lévy processes, which are sometimes called super-diffusion equations. In this article, we develop the differential quadrature (DQ) methods for solving the 2D space-fractional diffusion equations on irregular domains. The methods in presence reduce the original equation into a set of ordinary differential equations (ODEs) by introducing valid DQ formulations to fractional directional derivatives based on the functional values at scattered nodal points on problem domain. The required weighted coefficients are calculated by using radial basis functions (RBFs) as trial functions, and the resultant ODEs are discretized by the Crank-Nicolson scheme. The main advantages of our methods lie in their flexibility and applicability to arbitrary domains. A series of illustrated examples are finally provided to support these points.     

Integral equations in the wave scattering problem at an irregular interface of domains

Modeling the scattering of electromagnetic waves at an interface of media with different characteristics, one encounters the conjugation problem. Using the method of boundary integral equations and the theory of generalized potentials, we prove the classical resolvability of this problem. The boundary is assumed to be irregular. This means that the plane is divided into two domains by a curve which coincides with a straight line, except for a finite part, producing the irregularity. We propose algorithms for the approximate solution of the conjugation problem based on the spline methods for the solution of integral equations. We theoretically substantiate the computational scheme, namely, we prove the convergence and estimate the convergence rate.     

Global strong solutions of Navier-Stokes equations with interface boundary in three-dimensional thin domains

In this article, we study the spectrum of the Stokes operator in a 3D two layer domain with interface, obtain the asymptotic estimates on the spectrum of the Stokes operator as thickness ε goes to zero. Based on the spectral decomposition of the Stokes operator, a new average-like operator is introduced and applied to the study of Navier-Stokes equation in the two layer thin domains under interface boundary condition. We prove the global existence of strong solutions to the 3D Navier-Stokes equations when the initial data and external forces are in large sets as the thickness of the domain is small. This article is a continuation of our study on the Stokes operator under Navier friction boundary condition. Due to the viscosity distinction between the two layers, the Stokes operator displays radically different spectral structure from that under Navier friction boundary condition, then causes great difficulty to the analysis.     

An approximate layering method for the Navier-Stokes equations in bounded cylinders

Sunto Si sviluppa un metodo di costruzione di soluzioni approssimate per le equazioni di Navier-Stokes in cilindri limitati mediante un'adattamento del metodo detto di «layering». Si dimostrano alcuni teoremi di convergenza.     

An element-by-element BICGSTAB iterative method for three-dimensional steady Navier-Stokes equations

Construction of a stabilized Galerkin upwind finite element model for steady and incompressible Navier-Stokes equations in three dimensions is the main theme of this study. In the time-independent context, the weighted residuals statement is kept biased in favor of the upstream flow direction by adding an artificial damping term of physical plausibility to the Galerkin framework. This upwind approach has significant advantage of seeking solutions free from cross-stream diffusion error. Finite element solutions have been found by mixed formulation, implemented in quadratic cubic elements which are characterized as possessing the so-called LBB (Ladyzhenskaya-Babuška-Brezzi) condition. An element-by-element BICGSTAB solution solver is intended to alleviate difficulties regarding the asymmetry and indefiniteness arising from the use of a mixed formulation for incompressible fluid flows. The developed three-dimensional finite element code is first rectified by solving a problem amenable to analytic solution. A well-known lid-driven cavity flow problem in a cubical cavity is also studied.     

An analysis of a mixed finite element method for the Navier-Stokes equations

Summary A simple mixed finite element method is developed to solve the steady state, incompressible Navier-Stokes equations in a neighborhood of an isolated—but not necessarily unique—solution. Convergence is established under very mild restrictions on the triangulation, and, when the solution is sufficiently smooth, optimal error bounds are obtained.     

Clifford analysis and the Navier-Stokes equations over unbounded domains

We consider parabolic Dirac operators which do not involve fractional derivatives and use them to show the solvability of the in-stationary Navier-Stokes equations over time-varying domains. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Strong solutions of the Navier-Stokes equations in aperture domains

We consider the nonstationary Navier-Stokes equations in an aperture domain Ω⊂R³ consisting of two halfspaces separated by a wall, but connected by a hole in this wall. In this special domain one has to impose an auxiliary condition to single out a unique solution. This can be done by prescribing either the flux through the hole or the pressure drop between the two halfspaces. We construct suitable Stokes operators for both of the auxiliary conditions and show that they generate holomorphic semigroups. Then we prove the existence and uniqueness of solutions as well as a maximal regularity estimate for the Stokes equations subject to one of the auxiliary conditions. For the corresponding Navier-Stokes equations we prove existence and uniqueness of local in time solutions.

---

Sunto In questo lavoro consideriamo le equazioni di Navier-Stokes non stazionarie in un dominio con un’apertura, che consiste di due semispazi separati da una parete, ma collegati da un’apertura in quest’ultima. In questo dominio particolare è necessario imporre, per avere un’unica soluzione, una opportuna condizione ausiliaria. Questo può essere fatto sia assegnando il flusso attraverso l’apertura sia prescrivendo il salto di pressione tra i due semispazi. Qui costruiamo degli operatori di Stokes opportuni per ambedue i tipi di condizioni ausiliarie e mostriamo come essi generino semigruppi olomorfi. Dimostriamo, quindi, esistenza e unicità di soluzioni, assieme ad una stima di massima regolarità per le equazioni di Stokes soggette ad una delle condizioni ausiliarie. Per le corrispondenti equazioni di Navier-Stokes, dimostriamo esistenza e unicità di soluzioni locali nel tempo.

     

Nonlinear Galerkin method and two-step method for the Navier-Stokes equations

This article represents a new nonlinear Galerkin scheme for the Navier-Stokes equations. This scheme consists of a nonlinear Galerkin finite element method and a two-step difference method. Moreover, we also provide a Galerkin scheme. By convergence analysis, two numerical schemes have the same second-order convergence accuracy for the spatial discretization and time discretization if H is chosen such that H = O(h^2/3). However, the nonlinear Galerkin scheme is simpler than the Galerkin scheme, namely, this scheme can save a large amount of computational time. © 1996 John Wiley & Sons, Inc.     

Decay estimates for strong solutions of the Navier-Stokes equations in exterior domains

Based on theL _p−L_q-estimates for the Stokes semigroup on exterior domains, proven by Iwashita and refined by Maremonti-Solonnikov and Shibata, the precise time-asymptotic behaviour of solutions of the Navier-Stokes equations is derived. The estimates are improved, if some additional information on the initial value is available. The results apply also to weak solutions in three and four dimensions.

---

Sunto In questo lavoro forniamo l’andamento asintotico nel tempo preciso di soluzioni delle equazioni di Navier-Stokes in un dominio esterno, basandoci sulle stimeL _p−L _q per il semigruppo di Stokes fornite, per primo, da Iwashita e quindi rifinite da Maremonti-Solonnikov e Shibata. Tali stime, se sono note ulteriori informazioni sui dati iniziali, sono migliorate. Il risultato si applica anche a soluzioni deboli in dimensione tre e quattro.

     

Interface boundary value problem for the Navier-Stokes equations in thin two-layer domains

We study a system of 3D Navier-Stokes equations in a two-layer parallelepiped-like domain with an interface coupling of the velocities and mixed (free/periodic) boundary condition on the external boundary. The system under consideration can be viewed as a simplified model describing some features of the mesoscale interaction of the ocean and atmosphere. In case when our domain is thin (of order ε), we prove the global existence of the strong solutions corresponding to a large set of initial data and forcing terms (roughly, of order ε^−2/3). We also give some results concerning the large time dynamics of the solutions. In particular, we prove a spatial regularity of the global weak attractor.     

An approximate inertial manifolds approach to postprocessing the Galerkin method for the Navier-Stokes equations

In a recent paper we have introduced a postprocessing procedure for the Galerkin method for dissipative evolution partial differential equations with periodic boundary conditions. The postprocessing technique uses approximate inertial manifolds to approximate the high modes (the small scale components) in the exact solutions in terms of the Galerkin approximations, which in this case play the role of the lower modes (large scale components). This procedure can be seen as a defect-correction technique. But contrary to standard procedures, the correction is computed only when the time evolution is completed. Here we extend these results to more realistic boundary conditions. Specifically, we study in detail the two-dimensional Navier-Stokes equations subject to homogeneous (nonslip) Dirichlet boundary conditions. We also discuss other equations, such as reaction-diffusion systems and the Cahn-Hilliard equations.