Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

A marker method for the solution of the damped Burgers' equation

A new method for the solution of the damped Burgers' equation is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details. The marker method is applicable to a general class of nonlinear dispersive partial differential equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

最大似然方法和Bayes方法在结构方程模型分析中的讨论(英文)

本文通过模拟研究,讨论了最大似然方法和Bayes方法在分析结构方程模型中的相似点和不同之处。     

Identification of nonstationary axisymmetric load distributed along a cylindrical shell

The paper outlines a procedure to identify the space-and time-dependent external nonstationary load acting on a closed circular cylindrical shell of medium thickness. Time-dependent deflections at several points of the shell are used as input data to solve the inverse problem. Examples of numerical identification of various nonstationary loads, including moving ones are presented. The relationship between the external load and the stress-strain state of the shell is described by the Volterra equation of the first kind. The identification problem is solved using Tikhonov's regularization method and Apartsin's h-regularization method __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 91–100, July 2008.     

Operator‐splitting method for the analysis of cavitation in liquid‐lubricated herringbone grooved journal bearings

This paper presents an operator‐splitting method (OSM) for the solution of the universal Reynolds equation. Jakobsson–Floberg–Olsson (JFO) pressure conditions are used to study cavitation in liquid‐lubricated journal bearings. The shear flow component of the oil film is first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component is computed by the Galerkin finite element method. Present OSM solutions for slider bearings are in good agreement with available analytical and experimental results. OSM is then applied to herringbone grooved journal bearings. The film pressure, cavitation areas, load capacity and attitude angle are obtained with JFO pressure conditions. The calculated load capacities are in agreement with available experimental data. However, a detailed comparison of the present results with those predicted using Reynolds pressure conditions shows some differences. The numerical results showed that the load capacity and the critical mass of the journal (linear stability indicator) are higher and the attitude angle is lower than those predicted by Reynolds pressure conditions for cases of high eccentricities. Copyright © 2004 John Wiley & Sons, Ltd.     

Linearized and rational approximation method for solving non‐linear Burgers' equation

A new numerical method called linearized and rational approximation method is presented to solve non‐linear evolution equations. The utility of the method is demonstrated for the case of differentiation of functions involving steep gradients. The solution of Burgers' equation is presented to illustrate the effectiveness of the technique for the solution of non‐linear evolution equations exhibiting nearly discontinuous solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Inviscid,laminar and turbulent opposed flows

This paper attempts to reproduce numerically previous experimental findings with opposed flows and extends their range to quantify the effects of upstream pipes and nozzles with inviscid, laminar and turbulent flows. The choice of conservation equations, boundary conditions, algorithms for their solution, the degree of grid dependence, numerical diffusion and the validity of numerical approximations are justified with supporting calculations where necessary. The results of all calculations on the stagnation plane show maximum strain rates close to the annular exit from the nozzles and pipes for lower separations and it can be expected that corresponding reacting flows will tend to extinguish in this region with the extinction moving towards the axis. With laminar flows, the maximum strain rate increased with Reynolds number and the maximum values were generally greater than with inviscid flows and smaller than with turbulent flows. With large separations, the strain rates varied less and this explains some results with reacting flows where the extinction appeared to begin on the axis. The turbulent‐flow calculations allowed comparison of three common variants of a two‐equation first‐moment closure. They provided reasonable and useful indications of strain rates but none correctly represented the rms of velocity fluctuations on the axis and close to the stagnation plane. As expected, those designed to deal with this problem produced results in better agreement with experiment but were still imperfect. Copyright © 2004 John Wiley & Sons, Ltd.