Shape reconstruction for the time‐dependent Navier‐Stokes flow

This article deals with the shape reconstruction of a bounded domain with a viscous incompressible fluid driven by the time‐dependent Navier‐Stokes equations. For the approximate solution of the ill‐posed and nonlinear problem we propose a regularized Newton method. A theoretical foundation for the Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the interior boundary curve in the sense of the domain derivative. Numerical examples indicate the feasibility of our method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Newton方程周期解存在唯一性的新证明

本文首先将Newton方程周期边值问题转化为初值问题，然后在较弱的条件下利用微分连续法构造性地证明了该方程周期解的存在唯一性．证明方法同时也提供了一种计算该周期解的大范围收敛方法．     

一类奇摄动非线性边值问题激波解的Sinc—Galerkin方法

讨论了一类非线性奇摄动方程的激波问题．利用Sinc—Galerkin方法,构造出边值问题的激波解,并由Newton法得到其近似解．     

轴对称渗流中自由边界问题的近似解法

１．问题的提法 本文讨论各向异性非均匀介质的轴对称稳定渗流问题,我们延用［７］的记号。设有一水井（或油井）,其截面如图１所示,ｚ轴为对称轴,Ｄ为渗流区域,其边界为ＡＢＣＥＦＡ,Ｋ＝｛ｋｉｊ（ｒ,ｚ,ｈ,ｑ）｝为对称渗流张量,它依赖于柱坐标中的ｒ,ｚ,压头ｈ＝ｚ＋ｐ／ρｇ及渗流速率其中ｐ为点（ｒ,ｚ）处的压力,ρ为流体密度,ｇ为重力加速度．ｒ０为井的半径,Ｈ１为液面的高度,同时假设当ｒ≥Ｒ时,其渗流速度Ｖ＝０． 由渗流理论,有引入热函数,流函数 Ｄ中一点）,则满足下列一阶非线性方程组其中,若（ｉ＝ １…     

A globally convergent semi-smooth Newton method for control-state constrained DAE optimal control problems

We investigate a semi-smooth Newton method for the numerical solution of optimal control problems subject to differential-algebraic equations (DAEs) and mixed control-state constraints. The necessary conditions are stated in terms of a local minimum principle. By use of the Fischer-Burmeister function the local minimum principle is transformed into an equivalent nonlinear and semi-smooth equation in appropriate Banach spaces. This nonlinear and semi-smooth equation is solved by a semi-smooth Newton method. We extend known local and global convergence results for ODE optimal control problems to the DAE optimal control problems under consideration. Special emphasis is laid on the calculation of Newton steps which are given by a linear DAE boundary value problem. Regularity conditions which ensure the existence of solutions are provided. A regularization strategy for inconsistent boundary value problems is suggested. Numerical illustrations for the optimal control of a pendulum and for the optimal control of discretized Navier-Stokes equations conclude the article.     

On the Finite Element Analysis of Problems with Nonlinear Newton Boundary Conditions in Nonpolygonal Domains

The paper is concerned with the study of an elliptic boundary value problem with a nonlinear Newton boundary condition considered in a two-dimensional nonpolygonal domain with a curved boundary. The existence and uniqueness of the solution of the continuous problem is a consequence of the monotone operator theory. The main attention is paid to the effect of the basic finite element variational crimes: approximation of the curved boundary by a polygonal one and the evaluation of integrals by numerical quadratures. With the aid of some important properties of Zlamal's ideal triangulation and interpolation, the convergence of the method is analyzed.     

A boundary element method for the Dirichlet eigenvalue problem of the Laplace operator

The solution of eigenvalue problems for partial differential operators by using boundary integral equation methods usually involves some Newton potentials which may be resolved by using a multiple reciprocity approach. Here we propose an alternative approach which is in some sense equivalent to the above. Instead of a linear eigenvalue problem for the partial differential operator we consider a nonlinear eigenvalue problem for an associated boundary integral operator. This nonlinear eigenvalue problem can be solved by using some appropriate iterative scheme, here we will consider a Newton scheme. We will discuss the convergence and the boundary element discretization of this algorithm, and give some numerical results.     

Recovery of an interface from boundary measurement in an elliptic differential equation

We study the inverse problem of recovering an interior interface from a boundary measurement in an elliptic boundary value problem arising from a semiconductor transistor model. We set up a nonlinear least-squares formulation for solving the inverse problem, and establish the necessary derivatives with respect to the interface. We then propose both the Gauss–Newton iterative method and the conjugate gradient method for the least-squares problem, and present implementation of these methods using integral equations.     

A nonlinearity lagging for the solution of nonlinear steady state reaction diffusion problems

This paper concerns with the analysis of the iterative procedure for the solution of a nonlinear reaction diffusion equation at the steady state in a two dimensional bounded domain supplemented by suitable boundary conditions. This procedure, called Lagged Diffusivity Functional Iteration (LDFI)-procedure, computes the solution by “lagging” the diffusion term. A model problem is considered and a finite difference discretization for that model problem is described. Furthermore, properties of the finite difference operator are proved. Then, sufficient conditions for the convergence of the LDFI-procedure are given. At each stage of the LDFI-procedure a weakly nonlinear algebraic system has to be solved and the simplified Newton–Arithmetic Mean (Newton–AM) method is used. This method is particularly well suited for implementation on parallel computers. Numerical studies show the efficiency, for different test functions, of the LDFI-procedure combined with the simplified Newton–AM method. Better results are obtained when in the reaction diffusion equation also a convection term is present.     

基于遗传算法重建多个散射体的组合Newton法

本文研究由单个入射声波或电磁波及其远场数据反演多个柔性散射体边界的逆散射问题.通过建立边界到边界总场的非线性算子及其n6chet导数,本文首先给出了基于单层位势的组合Newton法.将组合Newton法转化为泛响优化问题,从而获得了该方法重建单个散射体的收敛性分析.然后,基于遗传算法和正则化参数选取的模型函数方法,给出...     

On an Augmented Lagrangian SQP Method for a Class of Optimal Control Problems in Banach Spaces

An augmented Lagrangian SQP method is discussed for a class of nonlinear optimal control problems in Banach spaces with constraints on the control. The convergence of the method is investigated by its equivalence with the generalized Newton method for the optimality system of the augmented optimal control problem. The method is shown to be quadratically convergent, if the optimality system of the standard non-augmented SQP method is strongly regular in the sense of Robinson. This result is applied to a test problem for the heat equation with Stefan-Boltzmann boundary condition. The numerical tests confirm the theoretical results.     

Numerical simulation for the shape reconstruction of a cavity

We consider the inverse problem of reconstructing the interior boundary curve of a cavity from the knowledge of the measurements on the exterior boundary. The domain derivative of the corresponding operator is presented, and this allows the investigation of the regularized Newton method for the solution of the ill‐posed and nonlinear problem. Numerical examples indicate the feasibility of our method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A fast and accurate algorithm for computing radial transonic flows

An efficient algorithm is described for calculating stationary one-dimensional transonic outflow solutions of the compressible Euler equations with gravity and heat source terms. The stationary equations are solved directly by exploiting their dynamical system form. Transonic expansions are the stable manifolds of saddle-point-type critical points, and can be obtained efficiently and accurately by adaptive integration outward from the critical points. The particular transonic solution and critical point that match the inflow boundary conditions are obtained by a two-by-two Newton iteration which allows the critical point to vary within the manifold of possible critical points. The proposed Newton Critical Point (NCP) method typically converges in a small number of Newton steps, and the adaptively calculated solution trajectories are highly accurate. A sample application area for this method is the calculation of transonic hydrodynamic escape flows from extrasolar planets and the early Earth. The method is also illustrated for an example flow problem that models accretion onto a black hole with a shock.     

混合边界条件下的三维定常 旋转Navier-Stokes方程的原始变量有限元计算

本文利用原始变量有限元法求解混合边界条件下的三维定常旋转Navier-Stokes方程,证明了离散问题解的存在唯一性,得到了有限元解的最优误差估计.给出了求解原始变量有限元逼近解的简单迭代算法,并证明了算法的收敛性.针对三维情况下计算资源的限制,采用压缩的行存储格式存储刚度矩阵的非零元素,并利用不完全的LU分解作预处理的GMRES方法求解线性方程组.最后分析了简单迭代和牛顿迭代的优劣对比,数值算例表明在同样精度下简单迭代更节约计算时间.     

Solving generalized equations via homotopies

Global Newton methods for computing solutions of nonlinear systems of equations have recently received a great deal of attention. By using the theory of generalized equations, a homotopy method is proposed to solve problems arising in complementarity and mathematical programming, as well as in variational inequalities. We introduce the concepts of generalized homotopies and regular values, characterize the solution sets of such generalized homotopies and prove, under boundary conditions similar to Smale’s [10], the existence of a homotopy path which contains an odd number of solutions to the problem. We related our homotopy path to the Newton method for generalized equations developed by Josephy [3]. An interpretation of our results for the nonlinear programming problem will be given.     

A Regularized Newton Method in Inverse Scattering Problem from Cavities

We consider the inverse problem to determine the shape of a open cavity embedded in the infinite ground plane from knowledge of the far-field pattern of the scattering of TM polarization.For its approximate solution we propose a regularized Newton iteration scheme.For a foundation of Newton type methods we establish the Fréchet differentiability of solution to the scattering problem with respect to the boundary of the cavity.Some numerical examples of the feasibility of the method are presented.     

Inverse scattering from an open arc

A Newton method is presented for the approximate solution of the inverse problem to determine the shape of a sound-soft or perfectly conducting arc from a knowledge of the far-field pattern for the scattering of time-harmonic plane waves. Fréchet differentiability with respect to the boundary is shown for the far-field operator, which for a fixed incident wave maps the boundary arc onto the far-field pattern of the scattered wave. For the sake of completeness, the first part of the paper gives a short outline on the corresponding direct problem via an integral equation method including the numerical solution.     

Regularized Newton iteration method for a penetrable cavity with internal measurements in inverse scattering problem

In this paper, we consider the inverse scattering problem of determining the shape of a cavity with a penetrable inhomogeneous medium of compact support from one source and a knowledge of measurements placed on a curve inside the cavity. First, the boundary value problem of the partial differential equations can be transformed into an equivalent system of nonlinear and ill-posed integral equations for the unknown boundary. Then, we apply the regularized Newton iterative method to reconstruct the boundary and prove the injectivity for the linearized system. Finally, we present some numerical examples to show the feasibility of our method.     

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