Non‐homogeneous Navier–Stokes systems with order‐parameter‐dependent stresses

We consider the Navier–Stokes system with variable density and variable viscosity coupled to a transport equation for an order‐parameter c. Moreover, an extra stress depending on c and ?c, which describes surface tension like effects, is included in the Navier–Stokes system. Such a system arises, e.g. for certain models of granular flows and as a diffuse interface model for a two‐phase flow of viscous incompressible fluids. The so‐called density‐dependent Navier–Stokes system is also a special case of our system. We prove short‐time existence of strong solution in L^q‐Sobolev spaces with q>d. We consider the case of a bounded domain and an asymptotically flat layer with a combination of a Dirichlet boundary condition and a free surface boundary condition. The result is based on a maximal regularity result for the linearized system. Copyright © 2010 John Wiley & Sons, Ltd.     

Strong solutions for the nonhomogeneous Navier–Stokes equations in unbounded domains

We show the existence of strong solutions for the nonhomogeneous Navier–Stokes equations in three‐dimensional domains with boundary uniformly of class C³. Under suitable assumptions, uniqueness is also proved. Copyright © 2009 John Wiley & Sons, Ltd.     

Generalized Stokes resolvent equations in an infinite layer with mixed boundary conditions

In this paper we prove unique solvability of the generalized Stokes resolvent equations in an infinite layer Ω₀ = ℝ^{n –1} × (–1, 1), n ≥ 2, in L^q ‐Sobolev spaces, 1 < q < ∞, with slip boundary condition of on the “upper boundary” ∂Ω⁺₀ = ℝ^{n –1} × {1} and non‐slip boundary condition on the “lower boundary” ∂Ω^–₀ = ℝ^{n –1} × {–1}. The solution operator to the Stokes system will be expressed with the aid of the solution operators of the Laplace resolvent equation and a Mikhlin multiplier operator acting on the boundary. The present result is the first step to establish an L^q ‐theory for the free boundary value problem studied by Beale [9] and Sylvester [22] in L ²‐spaces. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non‐standard Stokes and Navier–Stokes problems: existence and regularity in stationary case

This paper is devoted to Stokes and Navier–Stokes problems with non‐standard boundary conditions: we consider, in particular, the case where the pressure is given on a part of the boundary. These problems were studied by Bégue, Conca, Murat and Pironneau. They proved the existence of variational solutions, indicating that these were solutions of the initial non‐standard problems, if they are regular enough, but without specifying the conditions on the data which would imply this regularity. In this paper, first we show that the variational solutions, on supposing pressure on the boundary Γ₂ of regularity H^1/2 instead of H^?1/2, have their Laplacians in L² and, therefore, are solutions of non‐standard Stokes problem. Next, we give a result of regularity H², which we generalize, obtaining regularities W^{m, r}, m∈?, m?2, r?2. Finally, by a fixed‐point argument, we prove analogous results for the Navier–Stokes problem, in the case where the viscosity νis large compared to the data. Copyright © 2002 John Wiley & Sons, Ltd.     

On the asymptotic stability of steady solutions of the Navier–Stokes equations in unbounded domains

We consider the problem of the asymptotic behaviour in the L²‐norm of solutions of the Navier–Stokes equations. We consider perturbations to the rest state and to stationary motions. In both cases we study the initial‐boundary value problem in unbounded domains with non‐compact boundary. In particular, we deal with domains with varying and possibly divergent exits to infinity and aperture domains. Copyright © 2007 John Wiley & Sons, Ltd.     

Uncoupling evolutionary groundwater‐surface water flows using the Crank–Nicolson Leapfrog method

Consider an incompressible fluid in a region Ω_f flowing both ways across an interface into a porous media domain Ω_p saturated with the same fluid. The physical processes in each domain have been well studied and are described by the Stokes equations in the fluid region and the Darcy equations in the porous media region. Taking the interfacial conditions into account produces a system with an exactly skew symmetric coupling. Spatial discretization by finite element method and time discretization by Crank–Nicolson LeapFrog give a second‐order partitioned method requiring only one Stokes and one Darcy subphysics and subdomain solver per time step for the fully evolutionary Stokes‐Darcy problem. Analysis of this method leads to a time step condition sufficient for stability and convergence. Numerical tests verify predicted rates of convergence; however, stability tests reveal the problem of growth of numerical noise in unstable modes in some cases. In such instances, the addition of time filters adds stability. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

M ‐elliptic pseudo‐differential operators on L p (ℝn )

We construct the minimal and maximal extensions in L ^p (?ⁿ ), 1 < p < ∞, for M ‐elliptic pseudo‐differential operators initiated by Garello and Morando. We prove that they are equal and determine the domains of the minimal, and hence maximal, extensions of M ‐elliptic pseudo‐differential operators. For M ‐elliptic pseudodifferential operators with constant coefficients, the spectra and essential spectra are computed. An application to quantization is given. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Fast Computations with the Inverse to Harmonic Potential Operators via Domain Decomposition

In this paper a method for fast computations with the inverse to weakly singular, hypersingular and double layer potential boundary integral operators associated with the Laplacian on Lipschitz domains is proposed and analyzed. It is based on the representation formulae suggested for above-mentioned boundary operations in terms of the Poincare-Steklov interface mappings generated by the special decompositions of the interior and exterior domains. Computations with the discrete counterparts of these formulae can be efficiently performed by iterative substructuring algorithms provided some asymptotically optimal techniques for treatment of interface operators on subdomain boundaries. For both two- and three-dimensional cases the computation cost and memory needs are of the order O(N log^p N) and O(N log² N), respectively, with 1 ≤ p ≤ 3, where N is the number of degrees of freedom on the boundary under consideration (some kinds of polygons and polyhedra). The proposed algorithms are well suited for serial and parallel computations.     

A Crank‐Nicolson Petrov‐Galerkin method with quadrature for semi‐linear parabolic problems

We propose and analyze an application of a fully discrete C² spline quadrature Petrov‐Galerkin method for spatial discretization of semi‐linear parabolic initial‐boundary value problems on rectangular domains. We prove second order in time and optimal order H¹ norm convergence in space for the extrapolated Crank‐Nicolson quadrature Petrov‐Galerkin scheme. We demonstrate numerically both L² and H¹ norm optimal order convergence of the scheme even if the nonlinear source term is not smooth. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

A Factorization Property of R‐Bounded Sets of Operators on Lp‐Spaces

Let T be a bounded operator on L^p‐space, with 1 ≤ p < ∞. A theorem of W. B. Johnson and L. Jones asserts that after an appropriate change of density, T actually extends to a bounded operator on L². We show that if 𝒯 ⊂ B (L^p) is an R‐bounded set of operators, then the latter result holds for any T ∈ 𝒯 with a common change of density. Then we give applications including results on R‐sectorial operators.     

Existence of Weak Solutions Up to Collision for Viscous Fluid‐Solid Systems with Slip

We study in this paper the movement of a rigid solid inside an incompressible Navier‐Stokes flow within a bounded domain. We consider the case where slip is allowed at the fluid/solid interface through a Navier condition. Taking into account slip at the interface is very natural within this model, as classical no‐slip conditions lead to unrealistic collisional behavior between the solid and the domain boundary. We prove for this model existence of weak solutions of Leray type, up to collision, in three dimensions. The key point is that, due to the slip condition, the velocity field is discontinuous across the fluid/solid interface. This prevents obtaining global H¹ bounds on the velocity, which makes many aspects of the theory of weak solutions for Dirichlet conditions inappropriate. © 2014 Wiley Periodicals, Inc.     

On the well-posed coupling between free fluid and porous viscous flows

We present a well-posed model for the Stokes/Brinkman problem with a family of jump embedded boundary conditions (J.E.B.C.) on an immersed interface with weak regularity assumptions. It arises from a general framework recently proposed for fictitious domain problems. Our model is based on algebraic transmission conditions combining the stress and velocity jumps on the interface $\Sigma$ separating the fluid and porous domains. These conditions are well chosen to get the coercivity of the operator. Then, the general framework allows us to prove new results on the global solvability of some models with physically relevant stress or velocity jump boundary conditions for the momentum transport at a fluid–porous interface. The Stokes/Brinkman problem with Ochoa-Tapia and Whitaker (1995) [9], [10] interface conditions and the Stokes/Darcy problem with Beavers and Joseph (1967) [13] conditions are both proved to be well-posed, by an asymptotic analysis. Up to now, only the Stokes/Darcy problem with Saffman (1971) [15] approximate interface conditions with negligible tangential porous velocity was known to be well-posed.     

Initial‐to‐Interface Maps for the Heat Equation on Composite Domains

A map from the initial conditions to the function and its first spatial derivative evaluated at the interface is constructed for the heat equation on finite and infinite domains with n interfaces. The existence of this map allows changing the problem at hand from an interface problem to a boundary value problem which allows for an alternative to the approach of finding a closed‐form solution to the interface problem.     

First‐Order system least squares for generalized‐Newtonian coupled Stokes‐Darcy flow

The coupled problem for a generalized Newtonian Stokes flow in one domain and a generalized Newtonian Darcy flow in a porous medium is studied in this work. Both flows are treated as a first‐order system in a stress‐velocity formulation for the Stokes problem and a volumetric flux‐hydraulic potential formulation for the Darcy problem. The coupling along an interface is done using the well‐known Beavers–Joseph–Saffman interface condition. A least squares finite element method is used for the numerical approximation of the solution. It is shown that under some assumptions on the viscosity the error is bounded from above and below by the least squares functional. An adaptive refinement strategy is examined in several numerical examples where boundary singularities are present. Due to the nonlinearity of the problem a Gauss–Newton method is used to iteratively solve the problem. It is shown that the linear variational problems arising in the Gauss–Newton method are well posed. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1150–1173, 2015     

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.     

A remark on free boundary problem of 1‐D compressible Navier–Stokes equations with density‐dependent viscosity

In this paper, we prove the existence and uniqueness of the weak solution of the one‐dimensional compressible Navier–Stokes equations with density‐dependent viscosity µ(ρ)=ρ^θ with θ∈(0, γ?2], γ>1. The initial data are a perturbation of a corresponding steady solution and continuously contact with vacuum on the free boundary. The obtained results apply for the one‐dimensional Siant–Venant model of shallow water and generalize ones in (Arch. Rational Mech. Anal. 2006; 182: 223–253). Copyright © 2009 John Wiley & Sons, Ltd.     

A new frequency‐uniform coercive boundary integral equation for acoustic scattering

A new boundary integral operator is introduced for the solution of the soundsoft acoustic scattering problem, i.e., for the exterior problem for the Helmholtz equation with Dirichlet boundary conditions. We prove that this integral operator is coercive in L²(Γ) (where Γ is the surface of the scatterer) for all Lipschitz star‐shaped domains. Moreover, the coercivity is uniform in the wavenumber k = ω/c, where ω is the frequency and c is the speed of sound. The new boundary integral operator, which we call the “star‐combined” potential operator, is a slight modification of the standard combined potential operator, and is shown to be as easy to implement as the standard one. Additionally, to the authors' knowledge, it is the only second‐kind integral operator for which convergence of the Galerkin method in L²(Γ) is proved without smoothness assumptions on Γ except that it is Lipschitz. The coercivity of the star‐combined operator implies frequency‐explicit error bounds for the Galerkin method for any approximation space. In particular, these error estimates apply to several hybrid asymptoticnumerical methods developed recently that provide robust approximations in the high‐frequency case. The proof of coercivity of the star‐combined operator critically relies on an identity first introduced by Morawetz and Ludwig in 1968, supplemented further by more recent harmonic analysis techniques for Lipschitz domains. © 2011 Wiley Periodicals, Inc.     

Lipschitz Estimates in Almost‐Periodic Homogenization

We establish uniform Lipschitz estimates for second‐order elliptic systems in divergence form with rapidly oscillating, almost‐periodic coefficients. We give interior estimates as well as estimates up to the boundary in bounded C^1,α domains with either Dirichlet or Neumann data. The main results extend those in the periodic setting due to Avellaneda and Lin for interior and Dirichlet boundary estimates and later Kenig, Lin, and Shen for the Neumann boundary conditions. In contrast to these papers, our arguments are constructive (and thus the constants are in principle computable) and the results for the Neumann conditions are new even in the periodic setting, since we can treat nonsymmetric coefficients. We also obtain uniform W^1,p estimates.© 2016 Wiley Periodicals, Inc.     

Abstract wave equations with acoustic boundary conditions

We define an abstract setting to treat wave equations equipped with time‐dependent acoustic boundary conditions on bounded domains of R ⁿ . We prove a well‐posedness result and develop a spectral theory which also allows to prove a conjecture proposed in [13]. Concrete problems are also discussed. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions

In this paper, we consider low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions. The time discretization is based on the Euler implicit scheme, and the spatial discretization is based on the low‐order element (P₁−P₁ or P₁−P₀) for the approximation of the velocity and pressure. Moreover, some error estimates for the numerical solution of fully discrete stabilized finite element scheme are obtained. Finally, numerical experiments are performed to confirm our theoretical results.