A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Influence of Azimuthal Structure of Surface Waves on Efficiency of their Excitation by Tubular Electron Beams

The excitation of eigen surface waves by tubular electron beams in cylindrical discharge devices is studied. The influence of the wave‐field azimuthal structure on the excitation efficiency and nonlinear stage of the plasmabeam instability is investigated both numerically and analytically. Analytical expressions for the saturation amplitude and excitation efficiency of the wave under study are derived. They are found to agree well with results obtained by numerical modelling of the plasma‐beam interaction presented in this paper. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Cochran's statistical theorem for outer inverses of matrices and matrix quadratic forms

We extend the matrix version of Cochran's statistical theorem to outer inverses of a matrix. As applications, we investigate the Wishartness and independence of matrix quadratic forms for Kronecker product covariance structures.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

一维粗糙表面散射特征的数值分析

使用Monte Carlo模拟的方法得出了随机粗糙表面，在Kirchhoff近似的基础上，利用数值分析的方法分析了一维随机粗糙表面的散射特征，得出了一维随机粗糙表面散射分布曲线，我们还考虑了遮蔽效应的影响，并且讨论了Kirchhoff近似的有效性。     

A posteriori error estimates for a time discrete scheme for a phase-field system of Penrose-Fife type

A time discrete scheme is used to approximate the solution toa phase field system of Penrose–Fife type with a non-conservedorder parameter. An a posteriori error estimate is presentedthat allows the estimation of the difference between continuousand semidiscrete solutions by quantities that can be calculatedfrom the approximation and given data.     

On the application of slip boundary condition on curved boundaries

Hydrodynamic simulations of sloshing phenomena often involve the application of slip boundary condition at the wetted surfaces. If these surfaces are curved, the ambiguous nature of the normal vector in the discretized problem can interfere with the application of such a boundary condition. Even the use of consistent normal vectors, preferred from the point of view of conservation, does not assure good approximation of the continuum slip condition in the discrete problem, and non‐physical recirculating flow fields may be observed. As a remedy, we consider the Navier slip condition, and more successfully, the so‐called BC‐free boundary condition. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Free boundary viscous flows at micro and mesolevel during liquid composites moulding process

Numerical simulation aspects, related to low Reynolds number free boundary viscous flows at micro and mesolevel during the resin impregnation stage of the liquid composite moulding process (LCM), are presented in this article. A free boundary program (FBP), developed by the authors, is used to track the movement of the resin front accurately by accounting for the surface tension effects at the boundary. Issues related to the global and local mass conservation (GMC and LMC) are identified and discussed. Unsuitable conditions for LMC and consequently GMC are uncovered at low capillary numbers, and hence a strategy for the numerical simulation of such flows is suggested. FBP encompasses a set of subroutines that are linked to modules in ANSYS. FBP can capture the void formation dynamics based on the analysis developed. We present resin impregnation dynamics in two dimensions. Extension to three dimensions is a subject for further research. Several examples are shown and efficiency of different stabilization techniques are compared. Copyright © 2004 John Wiley & Sons, Ltd.