A Newton type iterative inverse method for heat-conduction problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

双相介质波动方程孔隙率反演的同伦方法

从材料响应的理论合成应与实际测量数据相拟合这一出发点，将双相介质波劝方程参数的反演问题转化为非线性算子方程的零点求解问题，从而应用一种大范围收敛的同伦方尘土注来解非线性算子方程，并把这种方法用于Simon(1984)给出的具有解析的一维双相介质模型的数值模拟，最后的数值结果表明，给出的算法是十分有效的。     

A Newton's method scheme for solving free-surface flow problems

A Newton's method scheme is described for solving the system of non-linear algebraic equations arising when finite difference approximations are applied to the Navier–Stokes equations and their associated boundary conditions. The problem studied here is the steady, buoyancy-driven motion of a deformable bubble, assumed to consist of an inviscid, incompressible gas. The linear Newton system is solved using both direct and iterative equation solvers. The numerical results are in excellent agreement with previous work, and the method achieves quadratic convergence.     

二维边界层方程的迭代求解

针对二维边界层方程,提出了分析积分迭代法.首先将该方法应用于Blasius方程和Falkner-Skan方程的求解,数值计算结果稳定,计算精度高;然后对外部有势流不能达到自相似要求、必须二维求解的二维层流边界层问题,在分析积分迭代法中加上计算力学的松弛迭代法,形成了一套有效的算法.数值结果表明该方法用于层流边界层的计算是很有效的.     

索网天线找形分析的参变量变分有限元法

针对柔性空间索网天线的非线性特点,建立了基于参变量变分描述索网拉压非线性和共旋列式描述几何非线性的有限元控制方程,应用Lernke与改进牛顿法进行求解.通过对索网预张力平衡计算,证明改进牛顿法比Newton-Raphson法具有更强的收敛能力.进一步将力密度法迭代原理与有限元法结合应用于索网天线的非线性找形分析中,获得了理想的索网构型.本文的索网找形方法可广泛应用于空间索网天线结构的设计.     

Fully coupled finite volume solutions of the incompressible Navier–Stokes and energy equations using an inexact Newton method

An inexact Newton method is used to solve the steady, incompressible Navier–Stokes and energy equation. Finite volume differencing is employed on a staggered grid using the power law scheme of Patankar. Natural convection in an enclosed cavity is studied as the model problem. Two conjugate-gradient -like algorithms based upon the Lanczos biorthogonalization procedure are used to solve the linear systems arising on each Newton iteration. The first conjugate-gradient-like algorithm is the transpose-free quasi-minimal residual algorithm (TFQMR) and the second is the conjugate gradients squared algorithm (CGS). Incomplete lower-upper (ILU) factorization of the Jacobian matrix is used as a right preconditioner. The performance of the Newton- TFQMR algorithm is studied with regard to different choices for the TFQMR convergence criteria and the amount of fill-in allowed in the ILU factorization. Performance data are compared with results using the Newton-CGS algorithm and previous results using LINPACK banded Gaussian elimination (direct-Newton). The inexact Newton algorithms were found to be CPU competetive with the direct-Newton algorithm for the model problem considered. Among the inexact Newton algorithms, Newton-CGS outperformed Newton- TFQMR with regard to CPU time but was less robust because of the sometimes erratic CGS convergence behaviour.     

高阶谱元区域分解算法求解定常方腔驱动流

主要利用Jacobian-free的Newton-Krylov方法求解定常不可压缩Navier-Stokes方程,将基于高阶谱元法的区域分解Stokes算法的非定常时间推进步作为Newton迭代的预处理,回避了传统Newton方法Jacobian矩阵的显式装配,节省了程序内存,同时降低了Newton迭代线性系统的条件数,且没有非线性对流项的隐式求解,大大加快了收敛速度。对有分析解的Kovasznay流动的计算结果表明,本高阶谱元法在空间上有指数收敛的谱精度,且对定常解的Newton迭代是二次收敛的。本文模拟了二维方腔顶盖一致速度驱动流,同基准解符合得很好,表明本文方法是准确可靠的。本文还考虑了Re=800时方腔顶盖正弦速度驱动流,除得到已知的一个稳定对称解和一对稳定非对称解外,还获得了一对新的不稳定的非对称解。     

Several new types of finite-difference schemes for shallow-water equation with initial-boundary value and their numerical experiment

This paper proposed several new types of finite-difference methods for the shallow water equation in absolute coordinate system and put forward an effective two-step predictorcorrector method, a compact and iterative algorithm for five diagonal matrix. Then the iterative method was used for a multigrid procedure for shallow water equation. At last, an initial-boundary value problem was considered, and the numerical results show that the linear sinusoidal wave would successively evolve into conoidal wave.     

半空间SH波方程非线性井间双参数反演

以半空间的ＳＨ波方程出发,采用Ｂｏｒｎ迭代法求解半空间弹性介质中密度和剪切模量分布的非线性反演问题。首先,采用矩量法和正则化方法,给出井间反演积分方程的离散形式,然后应用Ｂｒｏｎ迭代法求解非线性反演问题。     

The auto-adjustable damping method for solving nonlinear equations

I.IntroductionTherearealotofnonlinearproblemsinengineering.Forexamples,theheattrallsfel'with'nonlinearradiation-convectionboundaryalldlimitedheatcapacity,thenonlinearhydrodynamicsandaerodynamics,thenonlillearstrtlcturaldynamicresponsesandtheIihe.Itisalmostimpossibletosolvethe11onlinearequationsbydirectmethodexceptfi)rsomespecialnonlinearequations.Thegeneralapproachislinearizingtheequationsandformingvariousiterationprocedures.Forthestronglynonlillearproblems,thesolutionobtainedintheitcrativep…     

Fast stochastic analysis for non-proportionally damped system

Fundamental principles from structural dynamics, pseudo excitation method and pseudo force method are used to develop a new fast stochastic method for seismic response analysis of non-proportionally damped system. In the method, the mathematical equation of non-proportionally damped system is expressed in the iterative form, based on which the inverse operation of the matrices is avoided. Moreover, the new method does not need the solution of any complex eigenvalue problem, in contrast to other methods found in the literature.     

Riemann solution for ideal isentropic magnetogasdynamics

In the present paper, we study the Riemann problem for quasilinear hyperbolic system of partial differential equations governing the one dimensional ideal isentropic magnetogasdynamics with transverse magnetic field. We discuss the properties of rarefaction waves, shocks and contact discontinuities. Differently from single equation methods rooted in the ideal gasdynamics, the new approach is based on the system of two nonlinear equations imposing the equality of total pressure and velocity, assuming as unknowns the two values of densities, on both sides of the contact discontinuity. Newton iterative method is used to obtain densities. The resulting exact solver is implemented with the examples of general applicability of the proposed approach. For comparisons with exact solution we also shown numerical results obtained by the total variation diminishing slope limiter centre scheme. It is shown that both analytical and numerical results demonstrate the broad applicability and robustness of the new Riemann solver.     

基于半光滑牛顿法的润滑液膜有限元空化算法

针对润滑液膜中空化问题,引入Fischer-Burmeister函数,提出一种求解满足质量守恒雷诺方程的半光滑牛顿迭代算法.该算法将空化问题的非线性互补关系转化为等式约束方程,避免了迭代计算中的不等式约束识别问题.算法可将空化约束方程与雷诺方程、力平衡方程、变形方程等同时纳入牛顿迭代方程组,有效解决了传统松弛迭代算法需要多重嵌套循环带来的效率低下问题及压力与膜厚的强耦合性带来的收敛困难问题.计算实例表明,该算法计算效率高、收敛性好,且易应用于弹流润滑分析中,在滑动轴承和机械端面密封等多种物理模型下均有良好的适用性.     

A Finite Element Method for 3-D Eddy Current Problems Using an Iterative Approach

An iterative method is proposed for a finite element approximation of three-dimensional eddy current problems. The method is based on an iterative method derived from a perturbation problem of magnetostatic problems. The TEAM model and a transformer model are considered as numerical examples. In both examples, BiConjugate Gradient (BiCG) method is applicable for the complex symmetric linear system arising in each step of the iterative procedure for a rather wide range of the perturbation parameter, and the present results seem to be suitable.     

An iterative \mathit{GL}(n,\mathbb{R}) method for solving non-linear inverse vibration problems

For the inverse vibration problem, a differential-algebraic equation (DAE) method is proposed to simultaneously estimate the time-dependent damping and stiffness coefficients by using two sets of displacement and velocity as input data. We combine the equations of motion and the supplemental data into a set of DAEs. We develop an implicit $\mathit{GL}(n,\mathbb{R})$ scheme and a Newton iterative algorithm to stably solve the DAEs to find the unknown structural coefficients. The unknown force is also recovered by the present method. A linear oscillator and a non-linear Duffing oscillator are used as testing examples. The estimated results are rather accurate and robust against random noise; hence, the new method can be used in the solutions of non-linear inverse vibration problems.     

Algorithm for phase contrast X-ray tomography based on nonlinear phase retrieval

A new algorithm for phase contrast X-ray tomography under holographic measurement was proposed in this paper. The main idea of the algorithm was to solve the nonlinear phase retrieval problem using the Newton iterative method. The linear equations for the Newton directions were proved to be ill-posed and the regularized solutions were obtained by the conjugate gradient method. Some numerical experiments with computer simulated data were presented. The efficiency, feasibility and the numerical stability of the algorithm were illustrated by the numerical experiments. Compared with the results produced by the linearized phase retrieval algorithm, we can see that the new algorithm is not limited to be only efficient for the data measured in the near-field of the Fresnel region and thus it has a broader validity range.     

关于Orr—Sommerfeld方程的Chebyshev谱方法的讨论

本文讨论了Ｏｒｒ－Ｓｏｍｍｅｒｅｌｄ方程的各种Ｃｈｅｂｙｓｈｅｖ谱离散方法，数值证明了Ｃｈｅｂｙｓｈｅｖ配置法离散Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程没有伪谱，并以此构造了适于任意平面平行速度剖面情形，对时间和时空稳定性模式一致有效的无伪谱的离离散方法，其中，对时空稳定性问题本文给出了一种新的迭代法可以快速有效地求出复频率的鞍点，对平面Ｐｏｉｓｅｕｉｌｌｅ流，Ｂｌａｓｉｕｓ边界层流和Ｇａｕｓｓ模型尾迹     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Shape reconstruction of an inverse boundary value problem of two‐dimensional Navier–Stokes equations

This paper is concerned with the problem of the shape reconstruction of two‐dimensional flows governed by the Navier–Stokes equations. Our objective is to derive a regularized Gauss–Newton method using the corresponding operator equation in which the unknown is the geometric domain. The theoretical foundation for the Gauss–Newton method is given by establishing the differentiability of the initial boundary value problem with respect to the boundary curve in the sense of a domain derivative. The numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible. Copyright © 2009 John Wiley & Sons, Ltd.