Investigation and application of point implicit Runge–Kutta methods to inviscid flow problems

We derive and investigate point implicit Runge–Kutta methods to significantly improve the convergence rate to approximate steady‐state solutions of inviscid flows. It turns out that the point implicit Runge–Kutta can be interpreted as a preconditioned explicit Runge–Kutta method, where the preconditioner arises naturally as local derivative of the residual function. Moreover, many preconditioners suggested in the literature so far are identified as special case of our general ansatz. Conditions will be formulated such that explicit Runge–Kutta methods with local time stepping are equivalent to point implicit methods. In numerical examples, we will demonstrate the improved convergence rates. Copyright © 2011 John Wiley & Sons, Ltd.     

A dual‐weighted trust‐region adaptive POD 4‐D Var applied to a finite‐volume shallow water equations model on the sphere

In this paper we study solutions of an inverse problem for a global shallow water model controlling its initial conditions specified from the 40‐yr ECMWF Re‐analysis (ERA‐40) data sets, in the presence of full or incomplete observations being assimilated in a time interval (window of assimilation) with or without background error covariance terms. As an extension of the work by Chen et al. (Int. J. Numer. Meth. Fluids 2009), we attempt to obtain a reduced order model of the above inverse problem, based on proper orthogonal decomposition (POD), referred to as POD 4D‐Var for a finite volume global shallow water equation model based on the Lin–Rood flux‐form semi‐Lagrangian semi‐implicit time integration scheme. Different approaches of POD implementation for the reduced inverse problem are compared, including a dual‐weighted method for snapshot selection coupled with a trust‐region POD adaptivity approach. Numerical results with various observational densities and background error covariance operator are also presented. The POD 4‐D Var model results combined with the trust‐region adaptivity exhibit similarity in terms of various error metrics to the full 4D Var results, but are obtained using a significantly lesser number of minimization iterations and require lesser CPU time. Based on our previous and current work, we conclude that POD 4‐D Var certainly warrants further studies, with promising potential of its extension to operational 3‐D numerical weather prediction models. Copyright © 2011 John Wiley & Sons, Ltd.     

Reconstruction of Hp -Functions: Best Approximation,Regularization and Optimal Error Estimates

Let U be the open unit disc of the complex plane. We explicitly construct the best pointwise approximation for determining a function in the Hardy space H ^p (U) from measured data at a countable set of points in U whenever the data is exact. A regularization scheme is given to deal with the case of nonexact data. Moreover, optimal error estimates are also studied.     

Accuracies and conservation errors of various ghost fluid methods for multi-medium Riemann problem

Since the (original) ghost fluid method (OGFM) was proposed by Fedkiw et al. in 1999 [5], a series of other GFM-based methods such as the gas–water version GFM (GWGFM), the modified GFM (MGFM) and the real GFM (RGFM) have been developed subsequently. Systematic analysis, however, has yet to be carried out for the various GFMs on their accuracies and conservation errors. In this paper, we develop a technique to rigorously analyze the accuracies and conservation errors of these different GFMs when applied to the multi-medium Riemann problem with a general equation of state (EOS). By analyzing and comparing the interfacial state provided by each GFM to the exact one of the original multi-medium Riemann problem, we show that the accuracy of interfacial treatment can achieve “third-order accuracy” in the sense of comparing to the exact solution of the original mutli-medium Riemann problem for the MGFM and the RGFM, while it is of at most “first-order accuracy” for the OGFM and the GWGFM when the interface approach is actually near in balance. Similar conclusions are also obtained in association with the local conservation errors. A special test method is exploited to validate these theoretical conclusions from the numerical viewpoint.     

A penalized Fischer-Burmeister NCP-function

We introduce a new NCP-function in order to reformulate the nonlinear complementarity problem as a nonsmooth system of equations. This new NCP-function turns out to have stronger theoretical properties than the widely used Fischer-Burmeister function and other NCP-functions suggested previously. Moreover, numerical experience indicates that a semismooth Newton method based on this new NCP-function performs considerably better than the corresponding method based on the Fischer-Burmeister function. Received: March 10, 1997 / Accepted: February 15, 2000?Published online May 12, 2000     

A strong uniqueness theorem for planar vector fields

This work establishes a strong uniqueness property for a class of planar locally integrable vector fields. A result on pointwise convergence to the boundary value is also proved for bounded solutions.Dedicated to Constantine Dafermos on his 60^th birthdayWork supported in part by CNPq, FINEP, FAPESP and a Research Incentive Fund grant of Temple University.     

Peano type theorem for abstract parabolic equations

We consider parabolic problems with non-Lipschitz nonlinearity in the different scales of Banach spaces and prove local-in-time existence theorem.     

用次微分表达的对偶问题

本利用次微分建立了多目标规划的一个新的对偶问题，并给出其弱、强和逆对偶性，得到了一个新的次梯度的定义，并用其建立了一个新的对偶问题。     

Brachistochrone with Coulomb friction as the solution of an isoperimetrical variational problem

The problem of finding the plane supporting curve (brachistochrone) along which a heavy particle released from rest at the given starting point slides with dry (Coulomb) friction and reaches the given destination point in least time, is stated as a variational isoperimetrical problem. The finite parametric equations of the extremal curve are obtained. This curve is the sought-for brachistochrone if the solution of the problem exists. Several numerical examples are given.     

On the shear centre in Saint-Venant beam theory

The notion of shear centre in Saint-Venant beam theory was introduced by Robert Maillart who envisaged it to explain the results of experimental tests on beams with C-shaped sections. In literature, the location of the shear centre is provided in terms of flexure functions. The new result is a formula for the shear centre, based on the knowledge of the sole twist warping function of the cross-section.