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.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

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.     

Pediatric Laryngotracheal Stenosis and Airway Reconstruction: A Review of Voice Outcomes, Assessment, and Treatment Issues

Laryngotracheal stenosis is defined as a congenital or acquired narrowing of the airway. Congenital causes may include subglottic membranous or cartilaginous narrowing. Acquired causes may include trauma due to prolonged endotracheal or tracheal intubation or laryngotracheal injury. Although advances have been made over the past 30 years in reconstructive surgeries to improve airway patency in these patients, long-term laryngeal function for voice production is not well defined in this population. This review examines causes, symptoms and signs, and methods for diagnosing laryngotracheal stenosis. Surgical management procedures are briefly summarized. The current literature on voice outcomes is summarized. The predominant voice characteristics in the population are presented, although results are challenged by the heterogeneity of voice presentation and paucity of data from instrumental measures. Considerations for subjective and instrumental assessment, measures of quality of life, instrumental methods, and treatment options specific to the needs of this population are discussed. Research strategies to identify long-term outcomes of surgical and behavioral treatments in this population are posed.     

A model study for the steady state error of numerical approximations of inviscid flow equations on Cartesian grid

The widely used locally adaptive Cartesian grid methods involve a series of abruptly refined interfaces. In this paper we consider the influence of the refined interfaces on the steady state errors for second‐order three‐point difference approximations of flow equations. Since the various characteristic components of the Euler equations should behave similarly on such grids with regard to refinement‐induced errors, it is sufficient enough to conduct the analysis on a scalar model problem. The error we consider is a global error, different to local truncation error, and reflects the interaction between multiple interfaces. The steady state error will be compared to the errors on smooth refinement grids and on uniform grids. The conclusion seems to support the numerical findings of Yamaleev and Carpenter (J. Comput. Phys. 2002; 181: 280–316) that refinement does not necessarily reduce the numerical error. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates

In this paper we consider the numerical approximation of steady and unsteady generalized Newtonian fluid flows using divergence free finite elements generated by the Powell–Sabin–Heindl elements. We derive a priori and a posteriori finite element error estimates and prove convergence of the method of successive approximations for the steady flow case. A priori error estimates of unsteady flows are also considered. These results provide a theoretical foundation and supporting numerical studies are to be provided in Part II. Copyright © 2003 John Wiley & Sons, Ltd.     

Predictability and model selection in the context of ARCH models

Most of the methods used in the ARCH literature for selecting the appropriate model are based on evaluating the ability of the models to describe the data. An alternative model selection approach is examined based on the evaluation of the predictability of the models in terms of standardized prediction errors. Copyright © 2005 John Wiley & Sons, Ltd.     

Krylov subspaces and the analytic grade

Typical behaviour of the solution of a linear system of equations obtained iteratively by Krylov methods can be characterized by three stages. Initially the residual diminishes steadily; this is followed by stagnation and finally rapid convergence near the algebraic grade. This study examines this behaviour in terms of the concepts of approximately invariant subspace and what we have called the analytic grade of a Krylov sequence. It is shown how the small Ritz values play a vital role in the convergence and how this knowledge helps in the construction of an effective preconditioner. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydroxylation-induced modifications of the Al2O3/NiAl(0 0 1) surface at low water vapour pressure

STM, STS, LEED and XPS data for crystalline θ-Al₂O₃ and non-crystalline Al₂O₃ ultra-thin films grown on NiAl(0 0 1) at 1025 K and exposed to water vapour at low pressure (1 × 10⁻⁷-1 × 10⁻⁵ mbar) and room temperature are reported. Water dissociation is observed at low pressure. This reactivity is assigned to the presence of a high density of coordinatively unsaturated cationic sites at the surface of the oxide film. The hydroxyl/hydroxide groups cannot be directly identify by their XPS binding energy, which is interpreted as resulting from the high BE positions of the oxide anions (O_1s signal at 532.5-532.8 eV). However the XPS intensities give evidence of an uptake of oxygen accompanied by an increase of the surface coverage by Al³⁺ cations, and a decrease of the concentration in metallic Al at the alloy interface. A value of ∼2 for the oxygen to aluminium ions surface concentration ratio indicates the formation of an oxy-hydroxide (AlO_xOH_y with x + y ∼ 2) hydroxylation product. STM and LEED show the amorphisation and roughening of the oxide film. At P(H₂O) = 1 × 10⁻⁷ mbar, only the surface of the oxide film is modified, with formation of nodules of ∼2 nm lateral size covering homogeneously the surface. STS shows that essentially the valence band is modified with an increase of the density of states at the band edge. With increasing pressure, hydroxylation is amplified, leading to an increased coverage of the alloy by oxy-hydroxide products and to the formation of larger nodules (∼7 nm) of amorphous oxy-hydroxide. Roughening and loss of the nanostructure indicate a propagation of the reaction that modifies the bulk structure of the oxide film. Amorphisation can be reverted to crystallization by annealing under UHV at 1025 K when the surface of the oxide film has been modified, but not when the bulk structure has been modified.