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     

On non‐stationary viscous incompressible flow through a cascade of profiles

The paper deals with theoretical analysis of non‐stationary incompressible flow through a cascade of profiles. The initial‐boundary value problem for the Navier–Stokes system is formulated in a domain representing the exterior to an infinite row of profiles, periodically spaced in one direction. Then the problem is reformulated in a bounded domain of the form of one space period and completed by the Dirichlet boundary condition on the inlet and the profile, a suitable natural boundary condition on the outlet and periodic boundary conditions on artificial cuts. We present a weak formulation and prove the existence of a weak solution. Copyright © 2006 John Wiley & Sons, Ltd.     

Convergence of the solution for a domain decomposed,heterogeneously modeled flow past a concave profile problem

The free streamline problem investigated is that of fluid flow past a symmetric truncated concave‐shaped profile between walls. An open wake or cavity is formed behind the profile. Conformal mapping techniques are used to solve this problem. The problem formulated in the hodograph plane is decomposed into two nonoverlapping domains. Heterogeneous modeling is then used to describe the problems, i.e., a different governing differential equation in each domain. In one of these domains, a Baiocchi‐type transformation is used to obtain a fixed domain formulation for the part of the transformed problem containing an unknown boundary. In the other domain, the Baiocchi‐type transformation is extended across the boundary between the two domains, thus yielding a different problem formulation. This also assures that the dependent variables and their normal derivatives are continuous along this common boundary. The numerical solution scheme, a successive over‐relaxation approach, is applied over the whole problem domain with the use of a projection‐operation over only the fixed domain formulated part. Numerical results are obtained for the case of a truncated circular profile. These results are found to be in good agreement with another published result. The existence and uniqueness of the solution to the problem as a variational inequality is shown, and the convergence of the numerical solution using a domain decomposition method scheme is demonstrated by assuming some convergence property on the common interface of the two subdomains. © 2000 John Wiley & Sons, Inc. Numeer Methods Partial Differential Eq 16: 459–479, 2000     

Analysis of an augmented mixed‐primal formulation for the stationary Boussinesq problem

In this article, we propose and analyze a new mixed variational formulation for the stationary Boussinesq problem. Our method, which uses a technique previously applied to the Navier–Stokes equations, is based first on the introduction of a modified pseudostress tensor depending nonlinearly on the velocity through the respective convective term. Next, the pressure is eliminated, and an augmented approach for the fluid flow, which incorporates Galerkin‐type terms arising from the constitutive and equilibrium equations, and from the Dirichlet boundary condition, is coupled with a primal‐mixed scheme for the main equation modeling the temperature. In this way, the only unknowns of the resulting formulation are given by the aforementioned nonlinear pseudostress, the velocity, the temperature, and the normal derivative of the latter on the boundary. An equivalent fixed‐point setting is then introduced and the corresponding classical Banach Theorem, combined with the Lax–Milgram Theorem and the Babu?ka–Brezzi theory, are applied to prove the unique solvability of the continuous problem. In turn, the Brouwer and the Banach fixed‐point theorems are used to establish existence and uniqueness of solution, respectively, of the associated Galerkin scheme. In particular, Raviart–Thomas spaces of order k for the pseudostress, continuous piecewise polynomials of degree ≤ k+1 for the velocity and the temperature, and piecewise polynomials of degree ≤ k for the boundary unknown become feasible choices. Finally, we derive optimal a priori error estimates, and provide several numerical results illustrating the good performance of the augmented mixed‐primal finite element method and confirming the theoretical rates of convergence. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 445–478, 2016     

Stability analysis in the inverse Robin transmission problem

In this paper, we consider the conductivity problem with piecewise‐constant conductivity and Robin‐type boundary condition on the interface of discontinuity. When the quantity of interest is the jump of the conductivity, we perform a local stability estimate for a parameterized non‐monotone family of domains. We give also a quantitative stability result of local optimal solution with respect to a perturbation of the Robin parameter. In order to find an optimal solution, we propose a Kohn–Vogelius‐type cost functional over a class of admissible domains subject to two boundary values problems. The analysis of the stability involves the computation of first‐order and second‐order shape derivative of the proposed cost functional, which is performed rigorously by means of shape‐Lagrangian formulation without using the shape sensitivity of the states variables. © 2016 The Author. Mathematical Methods in the Applied Sciences Published by John Wiley & Sons Ltd.     

On an initial‐boundary value problem for a wide‐angle parabolic equation in a waveguide with a variable bottom

We consider the third‐order Claerbout‐type wide‐angle parabolic equation (PE) of underwater acoustics in a cylindrically symmetric medium consisting of water over a soft bottom B of range‐dependent topography. There is strong indication that the initial‐boundary value problem for this equation with just a homogeneous Dirichlet boundary condition posed on B may not be well‐posed, for example when B is downsloping. We impose, in addition to the above, another homogeneous, second‐order boundary condition, derived by assuming that the standard (narrow‐angle) PE holds on B, and establish a priori H² estimates for the solution of the resulting initial‐boundary value problem for any bottom topography. After a change of the depth variable that makes B horizontal, we discretize the transformed problem by a second‐order accurate finite difference scheme and show, in the case of upsloping and downsloping wedge‐type domains, that the new model gives stable and accurate results. We also present an alternative set of boundary conditions that make the problem exactly energy conserving; one of these conditions may be viewed as a generalization of the Abrahamsson–Kreiss boundary condition in the wide‐angle case. Copyright © 2008 John Wiley & Sons, Ltd.     

On a Non-symmetric Problem in Electrochemical Machining

We consider a free boundary value problem arising from a non-symmetric problem of electrochemical machining (ECM). After a conformal mapping of the unknown domain the problem is transformed to a non-smooth non-linear Riemann–Hilbert problem for holomorphic functions in the complex unit disk. A special technique allows to remove the singularities in the boundary condition. Utilizing existence results for smooth non-orientable Riemann–Hilbert problems, existence and uniqueness of a solution are shown. Finally, we propose an iterative method of Newton type for the effective numerical computation of the free boundary of the anode and present some test results. © 1997 B. G. Teubner Stuttgart–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     

Steady solutions of the Navier–Stokes equations with threshold slip boundary conditions

We establish the wellposedness of the time‐independent Navier–Stokes equations with threshold slip boundary conditions in bounded domains. The boundary condition is a generalization of Navier's slip condition and a restricted Coulomb‐type friction condition: for wall slip to occur the magnitude of the tangential traction must exceed a prescribed threshold, independent of the normal stress, and where slip occurs the tangential traction is equal to a prescribed, possibly nonlinear, function of the slip velocity. In addition, a Dirichlet condition is imposed on a component of the boundary if the domain is rotationally symmetric. We formulate the boundary‐value problem as a variational inequality and then use the Galerkin method and fixed point arguments to prove the existence of a weak solution under suitable regularity assumptions and restrictions on the size of the data. We also prove the uniqueness of the solution and its continuous dependence on the data. Copyright © 2006 John Wiley & Sons, Ltd.     

Integral equation methods for scattering by infinite rough surfaces

In this paper, we consider the Dirichlet and impedance boundary value problems for the Helmholtz equation in a non‐locally perturbed half‐plane. These boundary value problems arise in a study of time‐harmonic acoustic scattering of an incident field by a sound‐soft, infinite rough surface where the total field vanishes (the Dirichlet problem) or by an infinite, impedance rough surface where the total field satisfies a homogeneous impedance condition (the impedance problem). We propose a new boundary integral equation formulation for the Dirichlet problem, utilizing a combined double‐ and single‐layer potential and a Dirichlet half‐plane Green's function. For the impedance problem we propose two boundary integral equation formulations, both using a half‐plane impedance Green's function, the first derived from Green's representation theorem, and the second arising from seeking the solution as a single‐layer potential. We show that all the integral equations proposed are uniquely solvable in the space of bounded and continuous functions for all wavenumbers. As an important corollary we prove that, for a variety of incident fields including an incident plane wave, the impedance boundary value problem for the scattered field has a unique solution under certain constraints on the boundary impedance. Copyright © 2003 John Wiley & Sons, Ltd.     

Local and non‐local boundary quenching

The quenching problem is examined for a one‐dimensional heat equation with a non‐linear boundary condition that is of either local or non‐local type. Sufficient conditions are derived that establish both quenching and non‐quenching behaviour. The growth rate of the solution near quenching is also given for a power‐law non‐linearity. The analysis is conducted in the context of a nonlinear Volterra integral equation that is equivalent to the initial–boundary value problem. Copyright © 1999 John Wiley & Sons, Ltd.     

The solvability of a coupled thermoviscoelastic contact problem with small Coulomb friction and linearized growth of frictional heat

A coupled thermoviscoelastic frictional contact problem is investigated. The contact is modelled by the Signorini condition for the displacement velocities and the friction by the Coulomb law. The heat generated by friction is described by a non‐linear boundary condition with at most linear growth. The weak formulation of the problem consists of a variational inequality for the elasticity part and a variational equation for the heat conduction part. In order to prove the existence of a solution to this problem we first use an approximation of the Signorini condition by the penalty method. The existence of a solution for the approximate problem is shown using the fixed‐point theorem of Schauder. This theorem is applied to the composition of the solution operator for the contact problem with given temperature field and the solution operator for the heat equation problem with known displacement field. To obtain this proof, the unique solvability of both problems is necessary. Due to this reason it is necessary to introduce the penalty method. While the penalized contact problem has a unique solution, this is not clear for the original contact problem. The solvability of the original frictional contact problem is verified by an investigation of the limit for vanishing penalty parameter. Copyright © 1999 John Wiley & Sons, Ltd.     

On a time-splitting method for a scalar conservation law with a multiplicative stochastic perturbation and numerical experiments

In this paper, we present a numerical scheme for a first-order hyperbolic equation of nonlinear type perturbed by a multiplicative noise. The problem is set in a bounded domain D of ${\mathbb{R}^{d}}$ and with homogeneous Dirichlet boundary condition. Using a time-splitting method, we are able to show the existence of an approximate solution. The result of convergence of such a sequence is based on the work of Bauzet–Vallet–Wittbold (J Funct Anal, 2013), where the authors used the concept of measure-valued solution and Kruzhkov’s entropy formulation to show the existence and uniqueness of the stochastic weak entropy solution. Then, we propose numerical experiments by applying this scheme to the stochastic Burgers’ equation in the one-dimensional case.     

Indirect boundary integral formulation for the solution of shear‐thinning flow inside single rotor mixers

An efficient indirect boundary integral formulation for the evaluation of inelastic non‐Newtonian shear‐thinning flows at low Reynolds number is presented in this article. The formulation is based on the solution of a homogeneous Stokes flow field and the use of a particular solution for the nonlinear non‐Newtonian terms that yields the complete solution to the problem. Matrix multiplications are reduced in comparison to other means of handling nonlinear terms in boundary integral formulations such as the dual reciprocity method. The iterative solution of the nonlinear system of equations has been performed with a modified Newton‐Raphson method obtaining accurate results for values of the power law index as low as 0.4 without domain partitioning. Geometries such as Couette flow and a typical industrial polymer mixer have been analyzed with the proposed method obtaining good results with a reduction in computational cost compared with other equivalent formulations. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27:1610–1627, 2011     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

Mixed and hybrid formulations for the three‐dimensional magnetostatic problem

We propose mixed and hybrid formulations for the three‐dimensional magnetostatic problem. Such formulations are obtained by coupling finite element method inside the magnetic materials with a boundary element method. We present a formulation where the magnetic field is the state variable and the boundary approach uses a scalar Dirichlet‐Neumann map to describe the exterior domain. Also, we propose a second formulation where the magnetic induction is the state variable and a vectorial Dirichlet‐Neumann map is used to describe the outer field. Numerical discretizations with “edge” and “face” elements are proposed, and for each discrete problem we study an “inf‐sup” condition. © 2002 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 18: 85–104, 2002     

Multipole fast algorithm for the least‐squares approach of the method of fundamental solutions for three‐dimensional harmonic problems

In this article we describe an improvement in the speed of computation for the least‐squares method of fundamental solutions (MFS) by means of Greengard and Rokhlin's FMA. Iterative solution of the linear system of equations is performed for the equations given by the least‐squares formulation of the MFS. The results of applying the method to test problems from potential theory with a number of boundary points in the order of 80,000 show that the method can achieve fast solutions for the potential and its directional derivatives. The results show little loss of accuracy and a major reduction in the memory requirements compared to the direct solution method of the least squares problem with storage of the full MFS matrix. The method can be extended to the solution of overdetermined systems of equations arising from boundary integral methods with a large number of boundary integration points. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 828–845, 2003.     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

Scattering of sound from point sources by multiple circular cylinders using addition theorem and superposition technique

In this study, we use the addition theorem and superposition technique to solve the scattering problem with multiple circular cylinders arising from point sound sources. Using the superposition technique, the problem can be decomposed into two individual parts. One is the free‐space fundamental solution. The other is a typical boundary value problem (BVP) with specified boundary conditions derived from the addition theorem by translating the fundamental solution. Following the success of null‐field boundary integral formulation to solve the typical BVP of the Helmholtz equation with Fourier densities, the second‐part solution is easily obtained after collocating the observation point exactly on the real boundary and matching the boundary condition. The total solution is obtained by superimposing the two parts which are the fundamental solution and the semianalytical solution of the Helmholtz problem. An example was demonstrated to validate the present approach. The parameter study of size and spacing between cylinders are addressed. The results are well compared with the available theoretical solutions and experimental data. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2011