An unstructured mesh finite volume method for modelling saltwater intrusion into coastal aquifers

In this paper, a two-dimensional finite volume unstructured mesh method (FVUM) based on a triangular background interpolation mesh is developed for analysing the evolution of the saltwater intrusion into single and multiple coastal aquifer systems. The model formulation consists of a ground-water flow equation and a salt transport equation. These coupled and non-linear partial differential equations are transformed by FVUM into a system of differential/algebraic equations, which is solved using backward differentiation formulas of order one through five. Simulation results are compared with previously published solutions where good agreement is observed.     

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.     

Multiple Wavelength Extinction Technique for Particle Characterization in Dense Particle Clouds

In utilizing the advantages of extinction measurements in micron and especially submicron particle characterization, the properties of a multiple wavelength extinction technique have been the subject of extended theoretical studies. Furthermore, an experimental set-up was designed which provides high flexibility owing to its modular design. The performance of the technique described is demonstrated by a large variety of applications in aerosol and combustion research and in large-scale industrial systems. It was found to be a reliable tool in characterizing dense particulate systems.     

Essentially optimal explicit Runge–Kutta methods with application to hyperbolic–parabolic equations

Optimal explicit Runge–Kutta methods consider more stages in order to include a particular spectrum in their stability domain and thus reduce time-step restrictions. This idea, so far used mostly for real-line spectra, is generalized to more general spectra in the form of a thin region. In thin regions the eigenvalues may extend away from the real axis into the imaginary plane. We give a direct characterization of optimal stability polynomials containing a maximal thin region and calculate these polynomials for various cases. Semi-discretizations of hyperbolic–parabolic equations are a relevant application which exhibit a thin region spectrum. As a model, linear, scalar advection–diffusion is investigated. The second-order-stabilized explicit Runge–Kutta methods derived from the stability polynomials are applied to advection–diffusion and compressible, viscous fluid dynamics in numerical experiments. Due to the stabilization the time step can be controlled solely from the hyperbolic CFL condition even in the presence of viscous fluxes.     

On a zero equation model of turbulence

In this paper we study a widely used zero equation model of turbulence. The governing equations are derived by applying to the incompressible Navier-Stokes equations the Reynolds time averaging procedure. We achieve closure by employing the eddy viscosity concept. Using the Implicit Function Theorem we obtain an existence and uniquencess result. We also discuss the existence of nonsingular solutions. Finally, we present an algorithm for solving the modeled equations.