An Eulerian‐Lagrangian substructuring domain decomposition method for unsteady‐state advection‐diffusion equations

We develop an Eulerian‐Lagrangian substructuring domain decomposition method for the solution of unsteady‐state advection‐diffusion transport equations. This method reduces to an Eulerian‐Lagrangian scheme within each subdomain and to a type of Dirichlet‐Neumann algorithm at subdomain interfaces. The method generates accurate and stable solutions that are free of artifacts even if large time‐steps are used in the simulation. Numerical experiments are presented to show the strong potential of the method. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:565–583, 2001     

A Comparison of Preconditioners for the Steklov–Poincaré Formulation of the Fluid-Structure Coupling in Hemodynamics

A Fluid–Structure Interaction (FSI) problem can be reinterpreted as a heterogeneous problem with two subdomains. It is possible to describe the coupled problem at the interface between the fluid and the structure, yielding a nonlinear Steklov–Poincaré problem. The linear system can be linearized by Newton iterations on the interface and the resulting linear problem can be solved by the preconditioned GMRES method. In this work we investigate the behavior of preconditioners of Neumann–Neumann and Dirichlet–Neumann type. We find that, in the context of hemodynamics, the Dirichlet– Neumann, i.e., using Dirichlet boundary conditions on the fluid side and Neumann on the structure side, outperforms the Neumann–Neumann method, except when a weighting is used such that it basically reduces to the Dirichlet–Neumann method. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Well‐conditioned boundary integral formulations for high‐frequency elastic scattering problems in three dimensions

We construct and analyze a family of well‐conditioned boundary integral equations for the Krylov iterative solution of three‐dimensional elastic scattering problems by a bounded rigid obstacle. We develop a new potential theory using a rewriting of the Somigliana integral representation formula. From these results, we generalize to linear elasticity the well‐known Brakhage–Werner and combined field integral equation formulations. We use a suitable approximation of the Dirichlet‐to‐Neumann map as a regularizing operator in the proposed boundary integral equations. The construction of the approximate Dirichlet‐to‐Neumann map is inspired by the on‐surface radiation conditions method. We prove that the associated integral equations are uniquely solvable and possess very interesting spectral properties. Promising analytical and numerical investigations, in terms of spherical harmonics, with the elastic sphere are provided. Copyright © 2014 John Wiley & Sons, Ltd.     

Non‐local approach in mathematical problems of fluid–structure interaction

Three‐dimensional mathematical problems of interaction between elastic and scalar oscillation fields are investigated. An elastic field is to be defined in a bounded inhomogeneous anisotropic body occupying the domain Ω¯₁⊂ℝ³ while a physical (acoustic) scalar field is to be defined in the exterior domain Ω¯₂=ℝ³\Ω₁ which is filled up also by an anisotropic (fluid) medium. These two fields satisfy the governing equations of steady‐state oscillations in the corresponding domains together with special kinematic and dynamic transmission conditions on the interface ∂Ω₁. The problems are studied by the so‐called non‐local approach, which is the coupling of the boundary integral equation method (in the unbounded domain) and the functional‐variational method (in the bounded domain). The uniqueness and existence theorems are proved and the regularity of solutions are established with the help of the corresponding Steklov–Poincaré type operators and on the basis of the Gårding inequality and the Lax–Milgram theorem. In particular, it is shown that the physical fluid–solid acoustic interaction problem is solvable for arbitrary values of the frequency parameter. Copyright © 1999 John Wiley & Sons, Ltd.     

The null field equations for flexural oscillations of elastic plates

A null field method is constructed to solve the exterior Dirichlet, Neumann, and Robin boundary value problems associated with the high‐frequency harmonic oscillations of Mindlin‐type plates. The case of an infinite plate with a bounded elastic inclusion is also considered. Additionally, the completeness of certain sets of wavefunctions is investigated. Copyright © 2012 John Wiley & Sons, Ltd.     

Boundary layers for cellular flows at high Péclet numbers

We analyze the behavior of solutions of steady advection‐diffusion problems in bounded domains with prescribed Dirichlet data when the Péclet number Pe ? 1 is large. We show that the solution converges to a constant in each flow cell outside a boundary layer of width O(?^1/2), ? = Pe^?1, around the flow separatrices. We construct an ?‐dependent approximate “water pipe problem” purely inside the boundary layer that provides a good approximation of the solution of the full problem but has ?‐independent computational cost. We also define an asymptotic problem on the graph of streamline separatrices and show that solution of the water pipe problem itself may be approximated by an asymptotic, ?‐independent problem on this graph. Finally, we show that the Dirichlet‐to‐Neumann map of the water pipe problem approximates the Dirichlet‐to‐Neumann map of the separatrix problem with an error independent of the flow outside the boundary layers. © 2004 Wiley Periodicals, Inc.     

Semi‐analytical solution of MHD flow through boundary integrals on the pipe wall

A mathematical model is given for the magnetohydrodynamic (MHD) pipe flow as an inner Dirichlet problem in a 2D circular cross section of the pipe, coupled with an outer Dirichlet or Neumann magnetic problem. Inner Dirichlet problem is given as the coupled convection‐diffusion equations for the velocity and the induced current of the fluid coupling also to the outer problem, which is defined with the Laplace equation for the induced magnetic field of the exterior region with either Dirichlet or Neumann boundary condition. Unique solution of inner Dirichlet problem is obtained theoretically reducing it into two boundary integral equations defined on the boundary by using the corresponding fundamental solutions. Exterior solution is also given theoretically on the pipe wall with Poisson integral, and it is unique with Dirichlet boundary condition but exists with an additive constant obtained through coupled boundary and solvability conditions in Neumann wall condition. The collocation method is used to discretize these boundary integrals on the pipe wall. Thus, the proposed procedure is an improved theoretical analysis for combining the solution methods for the interior and exterior regions, which are consolidated numerically showing the flow behavior. The solution is simulated for several values of problem parameters, and the well‐known MHD characteristics are observed inside the pipe for increasing values of Hartmann number maintaining the continuity of induced currents on the pipe wall.     

Periodic Homogenization of Green and Neumann Functions

For a family of second‐order elliptic operators with rapidly oscillating periodic coefficients, we study the asymptotic behavior of the Green and Neumann functions, using Dirichlet and Neumann correctors. As a result we obtain asymptotic expansions of Poisson kernels and the Dirichlet‐to‐Neumann maps as well as optimal convergence rates in L^p and W^1,p for solutions with Dirichlet or Neumann boundary conditions. © 2014 Wiley Periodicals, Inc.     

Solvability,asymptotic analysis and regularity results for a mixed type interaction problem of acoustic waves and piezoelectric structures

In the paper, we investigate the mixed type transmission problem arising in the model of fluid–solid acoustic interaction when a piezoceramic elastic body (Ω⁺) is embedded in an unbounded fluid domain (Ω^?). The corresponding physical process is described by the boundary‐transmission problem for second‐order partial differential equations. In particular, in the bounded domain Ω⁺, we have a 4×4 dimensional matrix strongly elliptic second‐order partial differential equation, while in the unbounded complement domain Ω^?, we have a scalar Helmholtz equation describing acoustic wave propagation. The physical kinematic and dynamic relations mathematically are described by appropriate boundary and transmission conditions. With the help of the potential method and theory of pseudodifferential equations based on the Wiener–Hopf factorization method, the uniqueness and existence theorems are proved in Sobolev–Slobodetskii spaces. We derive asymptotic expansion of solutions, and on the basis of asymptotic analysis, we establish optimal Hölder smoothness results for solutions. Copyright © 2017 John Wiley & Sons, Ltd.     

Solutions to a model with Neumann boundary conditions for phase transitions driven by configurational forces

We study an initial boundary value problem of a model describing the evolution in time of diffusive phase interfaces in solid materials, in which martensitic phase transformations driven by configurational forces take place. The model was proposed earlier by the authors and consists of the partial differential equations of linear elasticity coupled to a nonlinear, degenerate parabolic equation of second order for an order parameter. In a previous paper global existence of weak solutions in one space dimension was proved under Dirichlet boundary conditions for the order parameter. Here we show that global solutions also exist for Neumann boundary conditions. Again, the method of proof is only valid in one space dimension.     

Lattice evolution solution for the nonlinear Poisson–Boltzmann equation in confined domains

The lattice evolution method for solving the nonlinear Poisson–Boltzmann equation in confined domain is developed by introducing the second-order accurate Dirichlet and Neumann boundary implements, which are consistent with the non-slip model in lattice Boltzmann method for fluid flows. The lattice evolution method is validated by comparing with various analytical solutions and shows superior to the classical numerical solvers of the nonlinear Poisson equations with Neumann boundary conditions. The accuracy and stability of the method are discussed. This lattice evolution nonlinear Poisson–Boltzmann solver is suitable for efficient parallel computing, complex geometry conditions, and easy extension to three-dimensional cases.     

Integral equations of the first kind in the theory of oscillating plates

A modified matrix of fundamental solutions is used to derive and solve first-kind integral equations for the problem of high-frequency harmonic oscillations of an infinite elastic plate with a hole when Dirichlet or Neumann conditions are prescribed on the boundary curve.     

Leray–Hopf solutions to a viscoelastoplastic fluid model with nonsmooth stress–strain relation

We consider a fluid model including viscoelastic and viscoplastic effects. The state is given by the fluid velocity and an internal stress tensor that is transported along the flow with the Zaremba–Jaumann derivative. Moreover, the stress tensor obeys a nonlinear and nonsmooth dissipation law as well as stress diffusion. We prove the existence of global-in-time weak solutions satisfying an energy inequality under general Dirichlet conditions for the velocity field and Neumann conditions for the stress tensor.     

On a resolvent estimate of the Stokes equation with Neumann–Dirichlet‐type boundary condition on an infinite layer

This paper is concerned with the standard Lp estimate of solutions to the resolvent problem for the Stokes operator on an infinite layer with ‘Neumann–Dirichlet‐type’ boundary condition. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

An inverse problem by eigenvalues of four spectra

Certain parts of the Dirichlet–Dirichlet, Neumann–Dirichlet, Dirichlet–Neumann and Neumann–Neumann spectra are used to find the potential of the Sturm–Liouville equation on a finite interval. This problem possesses a unique solution. Conditions are found necessary and sufficient for four sequences to be the corresponding parts of the four spectra.     

On the Existence of Trapped Modes in Channels of Arbitrary Cross-section

This article establishes the existence of trapped-mode solutions of a linearized water-wave problem. The fluid occupies a symmetric horizontal channel that is uniform everywhere apart from a confined region which either contains a thin vertical plate spanning the depth of the channel or has indentations in the channel walls. A trapped mode corresponds to an eigenvalue of a non-local Neumann–Dirichlet operator for an elliptic boundary-value problem in the fluid domain. The existence of such an eigenvalue is established by generalizing previous results concerning spectral theory for differential operators to this non-local operator. © 1997 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Analytic solutions for unsteady flows over a superhydrophobic surface

A class of unsteady analytic solutions for the evolution of viscous Newtonian fluids over a superhydrophobic surface is derived. The surface is represented by regular rectilinear riblets and the fluid is assumed to flow along the grooves’ direction. A mixed boundary condition is imposed on the surface, consisting in homogeneous Neumann conditions over the riblet voids and homogeneous Dirichlet conditions on the wall intervals. The transition between the above conditions is modelled through a Robin condition with non-constant smooth coefficients in a completely general manner. A global solution is derived and relevant examples, that can be fruitfully adopted as benchmark solutions for testing numerical solvers, are discussed.     

Boundary integral operators and boundary value problems for Laplace’s equation

In this paper, we define boundary single and double layer potentials for Laplace’s equation in certain bounded domains with d-Ahlfors regular boundary, considerably more general than Lipschitz domains. We show that these layer potentials are invertible as mappings between certain Besov spaces and thus obtain layer potential solutions to the regularity, Neumann, and Dirichlet problems with boundary data in these spaces.     

On the Diffraction by Biperiodic Anisotropic Structures

This article studies the scattering of electromagnetic waves by a nonmagnetic biperiodic structure. The structure consists of anisotropic optical materials and separates two regions with constant dielectric coefficients. The time harmonic Maxwell equations are transformed to an equivalent strongly elliptic variational problem for the magnetic field in a bounded biperiodic cell with nonlocal boundary conditions. This guarantees the existence of quasiperiodic magnetic fields in H 1 and electric fields in H (curl) solving Maxwell's equations. The uniqueness is proved for all frequencies excluding possibly a discrete set. The analytic dependence of these solutions on frequency and incident angles is studied.