Domain sensitivity analysis of the elastic far-field patterns in scattering from nonsmooth obstacles

We consider elastic scattering problems described by the Dirichlet or the Neumann boundary value problem for the elastodynamic equation in the exterior of a 2D bounded domain or in the exterior of a crack. The boundary of the domain is assumed to have a finite set of corner points where the scattered wave may have singular behaviour. The paper is concerned with the sensitivity of the far scattered field with respect to small perturbations of the shape of the scatterer. Using a modification of the method of adjoint problems (K. Dems, Z. Mróz, Internat. J. Solids Structures 20 (1984) 527-552) we obtain a representation for the shape derivative which is well suited for a numerical realization with boundary element methods and which shows in some cases directly the influence of the singularities of the solution on the sensitivity of the far-field patterns.     

Asymptotic analysis of Neumann periodic optimal boundary control problem

An optimal boundary control problem in a domain with oscillating boundary has been investigated in this paper. The controls are acting periodically on the oscillating boundary. The controls are applied with suitable scaling parameters. One of the major contribution is the representation of the optimal control using the unfolding operator. We then study the limiting analysis (homogenization) and obtain two limit problems according to the scaling parameters. Another notable observation is that the limit optimal control problem has three controls, namely, a distributed control, a boundary control, and an interface control. Copyright © 2016 John Wiley & Sons, Ltd.     

A fictitious domain approach to the numerical solution of PDEs in stochastic domains

We present an efficient method for the numerical realization of elliptic PDEs in domains depending on random variables. Domains are bounded, and have finite fluctuations. The key feature is the combination of a fictitious domain approach and a polynomial chaos expansion. The PDE is solved in a larger, fixed domain (the fictitious domain), with the original boundary condition enforced via a Lagrange multiplier acting on a random manifold inside the new domain. A (generalized) Wiener expansion is invoked to convert such a stochastic problem into a deterministic one, depending on an extra set of real variables (the stochastic variables). Discretization is accomplished by standard mixed finite elements in the physical variables and a Galerkin projection method with numerical integration (which coincides with a collocation scheme) in the stochastic variables. A stability and convergence analysis of the method, as well as numerical results, are provided. The convergence is “spectral” in the polynomial chaos order, in any subdomain which does not contain the random boundaries.     

A mixed formulation for the direct approximation of <Emphasis Type="Italic">L</Emphasis><Superscript>2</Superscript>-weighted controls for the linear heat equation

This paper deals with the numerical computation of null controls for the linear heat equation. The goal is to compute approximations of controls that drive the solution from a prescribed initial state to zero at a given positive time. In [Fernandez-Cara & Münch, Strong convergence approximations of null controls for the 1D heat equation, 2013], a so-called primal method is described leading to a strongly convergent approximation of distributed control: the controls minimize quadratic weighted functionals involving both the control and the state and are obtained by solving the corresponding optimality conditions. In this work, we adapt the method to approximate the control of minimal square integrable-weighted norm. The optimality conditions of the problem are reformulated as a mixed formulation involving both the state and its adjoint. We prove the well-posedeness of the mixed formulation (in particular the inf-sup condition) then discuss several numerical experiments. The approach covers both the boundary and the inner situation and is valid in any dimension.     

Numerical method for coupling the macro and meso scales in stochastic chemical kinetics

A numerical method is developed for simulation of stochastic chemical reactions. The system is modeled by the Fokker–Planck equation for the probability density of the molecular state. The dimension of the domain of the equation is reduced by assuming that most of the molecular species have a normal distribution with a small variance. The numerical approximation preserves properties of the analytical solution such as non-negativity and constant total probability. The method is applied to a nine dimensional problem modelling an oscillating molecular clock. The oscillations stop at a fixed point with a macroscopic model but they continue with our two dimensional, mixed macroscopic and mesoscopic model. Dedicated to the memory of Germund Dahlquist (1925–2005). AMS subject classification (2000)  65M20, 65M60     

无界区域Stokes 问题非重叠型区域分解算法及其收敛性

本文研究无界区域Stokes方程外问题的利用有限元法和自然边界归化的非蕈叠型区域分解算法,此方法对无界区域Stokes问题非常有效.给出连续和离散情形的D-N算法及其收敛性分析,得到算法收敛的充要条件及充分条件,并得到最优的松弛因子和压缩因子,最后给出数值算例予以验证.     

Highly accurate numerical solutions with repeated Richardson extrapolation for 2D laplace equation

A theoretical basis is presented for the repeated Richardson extrapolation (RRE) to reduce and estimate the discretization error of numerical solutions for heat conduction. An example application is described for the 2D Laplace equation using the finite difference method, a domain discretized with uniform grids, second-order accurate approximations, several variables of interest, Dirichlet boundary conditions, grids with up to 8,193 × 8,193 nodes, a multigrid method, single, double and quadruple precisions and up to twelve Richardson extrapolations. It was found that: (1) RRE significantly reduces the discretization error (for example, from 2.25E-07 to 3.19E-32 with nine extrapolations and a 1,025 × 1,025 grid, yielding an order of accuracy of 19.1); (2) the Richardson error estimator works for numerical results obtained with RRE; (3) a higher reduction of the discretization error with RRE is achieved by using higher calculation precision, a larger number of extrapolations, a larger number of grids and correct error orders; and (4) to obtain a given value error, much less CPU time and RAM memory are required for the solution with RRE than without it.     

Simulation of wake effects of windfarms using an Actuator Disk implementation

Wind turbines extract energy from the approaching flow field resulting in reduced wind speeds, increased turbulence and a wake downstream of the wind turbine. Wakes extending up to 20 km downstream of a windfarm are reported [1]. In [2] a multi-scale numerical method is proposed that addresses length scales of the boundary layer, windfarm and meso-scale. As a step towards full realization of the multi scale approach the implementation of an Actuator Disk (AD) model into OpenFOAM including a proper atmospheric boundary layer (ABL) is presented. For a meso-scale area of open sea, the standard ABL conditions of OpenFOAM can not be used. The pre implemented ABL boundary condition for OpenFOAM is restricted to roughness lengths greater than 0.001 m. The boundary condition are modified in order to achieve adequate homogeneous flow properties throughout the considered domain enabling almost constant velocity profiles. A new AD class is implemented in OpenFOAM to establish an arbitrary number and location of windturbines. The class uses the standard k-ϵ turbulence model. The wake effect of windfarms composed of windturbines with different heights and different lateral spacing is investigated. Specifically the windfarm Trianel Borkum in the North Sea of Germany is simulated. Within this work wakes extending up to 20 km in length could be observed, suggesting to study the ABL-windfarm interaction in more detail. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of a difference scheme for the heat equation in a long strip by artificial boundary conditions

The numerical solution of the heat equation on a strip in two dimensions is considered. An artificial boundary is introduced to make the computational domain finite. On the artificial boundary, an exact boundary condition is proposed to reduce the original problem to an initial‐boundary value problem in a finite computational domain. A difference scheme is constructed by the method of reduction of order to solve the problem in the finite computational domain. It is proved that the difference scheme is uniquely solvable, unconditionally stable and convergent with the convergence order 2 in space and order 3/2 in time in an energy norm. A numerical example demonstrates the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

On the isoperimetric problem for the higher eigenvalues of the Robin and Wentzell Laplacians

We consider the problem of minimising the kth eigenvalue, k ≥ 2, of the (p-)Laplacian with Robin boundary conditions with respect to all domains in \mathbbR^N{\mathbb{R}^N} of given volume. When k = 2, we prove that the second eigenvalue of the p-Laplacian is minimised by the domain consisting of the disjoint union of two balls of equal volume, and that this is the unique domain with this property. For p = 2 and k ≥ 3, we prove that in many cases a minimiser cannot be independent of the value of the constant in the boundary condition, or equivalently of the domain’s volume. We obtain similar results for the Laplacian with generalised Wentzell boundary conditions.     

On approximation in the fictitious domain method for a bidimensional transport equation

In this article a numerical method for solving a two‐dimensional transport equation in the stationary case is presented. Using the techniques of the variational calculus, we find the approximate solution for a homogeneous boundary‐value problem that corresponds to a square domain D₂. Then, using the method of the fictitious domain, we extend our algorithm to a boundary value problem for a set D that has an arbitrary shape. In this approach, the initial computation domain D (called physical domain) is immersed in a square domain D₂. We prove that the solution obtained by this method is a good approximation of the exact solution. The theoretical results are verified with the help of a numerical example. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A consistency analysis for the numerical solution of boundary integral equations

A method is presented for assessing the nature of the error incurred in the boundary integral equation (BIE) solution of both harmonic and biharmonic boundary value problems (BVPs). It is shown to what order of accuracy the governing partial differential equation is actually represented by the approximating numerical scheme, and how raising the order of the boundary ‘shape functions’ affects this representation. The effect of varying both the magnitude and the aspect ratio of the solution domain is investigated; it is found that the present technique may suggest an optimum nondimensional scaling for the BIE solution of a particular harmonic or biharmonic BVP.     

Coupled Galerkin and parameterization methods for optimal control of discretely connected parallel beams

Galerkin and wavelet methods for optimal boundary control of a couple of discretely connected parallel beams are proposed. First, the problem with boundary controls is converted into a problem with distributed controls. The problem is, then, reduced by a Galerkin-based approach into determining the optimal control of a linear time-invariant lumped parameter system, which will be solved by a wavelet-based method using Legendre wavelets. The integration-operational matrix and Kronecker product are utilized to significantly simplify the optimization problem into a system of linear equations. A numerical example is presented to demonstrate the applicability and the efficiency of the proposed method.     

Reproducing kernel method for singularly perturbed turning point problems having twin boundary layers

A numerical method is proposed for solving singularly perturbed turning point problems exhibiting twin boundary layers based on the reproducing kernel method (RKM). The original problem is reduced to two boundary layers problems and a regular domain problem. The regular domain problem is solved by using the RKM. Two boundary layers problems are treated by combining the method of stretching variable and the RKM. The boundary conditions at transition points are obtained by using the continuity of the approximate solution and its first derivatives at these points. Two numerical examples are provided to illustrate the effectiveness of the present method. The results compared with other methods show that the present method can provide very accurate approximate solutions.     

Global and domain-decomposed mixed methods for the solution of Maxwell's equations with application to magnetotellurics

A collection of global and domain decomposition mixed finite element schemes for the approximate solution of the harmonic Maxwell's equations on a bounded domain with absorbing boundary conditions at the artificial boundaries are presented. The numerical procedures allow us to solve efficiently the direct problem in magnetotellurics consisting of determining the electromagnetic scattered field in a two–dimensional earth model of arbitrary conductivity properties. Convergence results for the numerical procedures are derived. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 407–437, 1998     

Interface problem in holonomic elastoplasticity

The three-dimensional interface problem with the homogeneous Lamé system in an unbounded exterior domain and holonomic material behaviour in a bounded interior Lipschitz domain is considered. Existence and uniqueness of solutions of the interface problem are obtained rewriting the exterior problem in terms of boundary integral operators following the symmetric coupling procedure. The numerical approximation of the solutions consists in coupling of the boundary element method (BEM) and the finite element method (FEM). A Céa-like error estimate is presented for the discrete solutions of the numerical procedure proving its convergence.     

Stability of a finite element method for 3D exterior stationary Navier-Stokes flows

We consider numerical approximations of stationary incompressible Navier-Stokes flows in 3D exterior domains, with nonzero velocity at infinity. It is shown that a P1-P1 stabilized finite element method proposed by C. Rebollo: A term by term stabilization algorithm for finite element solution of incompressible flow problems, Numer. Math. 79 (1998), 283–319, is stable when applied to a Navier-Stokes flow in a truncated exterior domain with a pointwise boundary condition on the artificial boundary.     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

An efficient method for solving non-linear singularly perturbed two points boundary-value problems of fractional order

In this paper, we discuss a numerical solution of a class of non-linear fractional singularly perturbed two points boundary-value problem. The method of solution consists of solving reduced problem and boundary layer correction problem. A series method is used to solve the boundary layer correction problem, and then the series solutions is approximated by the Pade’ approximant of order [m, m]. Some theoretical results are established and proved. Two numerical examples are discussed to illustrate the efficiency of the present scheme.