Zinc oxide is known to produce a wide variety of nanostructures that show promise for a number of applications. The use of electrochemical deposition techniques for growing ZnO nanostructures can allow tight control of the morphology of ZnO through the wide range of deposition parameters available. Here we model the growth of the rods under typical electrochemical conditions, using the Nernst–Planck equations in two dimensions to predict the growth rate and morphology of the nanostructures as a function of time. Generally good quantitative and qualitative agreement is found between the model predictions and recent experimental results.

Ignition or thermal explosion in an oxidizing porous body of material can be described by a dimensionless reaction–diffusion equation of the form ∂tu=∇2u+λe−1/u. Here such equations will be formulated in symmetrically shaped bounded regions Ω, effectively reducing the mathematical formulation to that of one dimension. This is critically re-examined from a modern perspective using numerical methods. A computer algorithm is constructed and used to carry out a broad-ranging evaluation of the watershed critical initial temperature conditions for thermal ignition of nonuniform assemblies. It is then shown how the resulting mathematical structure for the ignition threshold curves can be correlated by a hyperbolic conic section with a high degree of accuracy over the full range of positive ambient temperature values. However, this sometimes over-predicts (which is bad) and sometimes under-predicts (which is good) the critical initial condition. The definition of additional dimensionless parameters is found to generate further simplification, leading to a universal correlating form capable of collapsing the entire solution space onto a single line in the plane of the new variables. In addition, this study considers the physically intuitive conjecture that spatial moments of the initial temperature profile ought to possess a direct mathematical link to the critical ignition threshold. As such, the mth-order spatial moment of the critical total energy content integrals is defined, and an empirical result is derived stating that certain orders of this moment should be insensitive to changes in ambient temperature and initial shape profile and may be considered functionally dependent on the dimensionless eigenvalue only, within some quantifiable error band. Spatial moment integrals, based on computed critical threshold conditions, are found to support this conjecture, with the best accuracy obtained for the second-order moments.

This Special Issue of the ANZIAM Journal arose from a conference held at Industrial Research, Gracefield, Lower Hutt, New Zealand, on 29 August 2006, to celebrate the research career of Dr Stephen White. Stephen was an international expert in the development and application of numerical methods to quantify fluid flow in porous media, especially those involving heat and chemical transport. As a mathematical modeller, Stephen White’s work has been published extensively, and applied in many practical areas.

Mathematical modelling of Wairakei geothermal field is reviewed, both lumped-parameter and distributed-parameter models. In both cases it is found that reliable predictions require five to ten years of history for calibration. With such calibration distributed-parameter models are now used for field management. A prudent model of Wairakei, constructed without such historical data, would underestimate field capacity and provide only general projections of the type of changes in surface activity and subsidence.

An analytic solution is developed for the one-dimensional dissipational slip gradient equation first described by Gurtin [“On the plasticity of single crystals: free energy, microforces, plastic strain-gradients”, J. Mech. Phys. Solids48 (2000) 989–1036] and then investigated numerically by Anand et al. [“A one-dimensional theory of strain-gradient plasticity: formulation, analysis, numerical results”, J. Mech. Phys. Solids53 (2005) 1798–1826]. However we find that the analytic solution is incompatible with the zero-sliprate boundary condition (“clamped boundary condition”) postulated by these authors, and is in fact excluded by the theory. As a consequence the analytic solution agrees with the numerical results except near the boundary. The equation also admits a series of higher mode solutions where the numerical result corresponds to (a particular case of) the fundamental mode. Anand et al. also established that the one-dimensional dissipational gradients strengthen the material, but this proposition only holds if zero-sliprate boundary conditions can be imposed, which we have shown cannot be done. Hence the possibility remains open that dissipational gradient weakening may also occur.

A three-layer compartmental model of the geological structure in the Taupo Volcanic Zone of New Zealand is developed, based on the assumptions of isostasy (constant geostatic pressure at 25 km depth) and a constant rate of volcanism. The upper layer consists of volcanic infill to a depth of about 2500 m, then a middle layer of greywacke-like material to a depth of about 15 km, and a lower layer of andesitic-like material to a depth of 25 km. Our model assumptions predict that the area of each layer increases at a constant rate; that there is a constant ratio between the rate of energy production from volcanic activity and geothermal convection; and that there is the possibility of an abrupt change from rhyolitic to basaltic volcanism, if the middle layer becomes sufficiently thin. Two models are considered: a rifting and a spreading model. Both models predict the lower layer has an andesitic-like density. The spreading model has difficulty matching heat output with observed extension rates. The rifting model predicts the observed extension rates, but requires very deep circulation of groundwater to be consistent with observed chemical and isotopic properties of geothermal fluids.

We consider a model for the control of a linear network flow system with unknown but bounded demandand polytopic bounds on controlled flows. We are interested in the problem of finding a suitable objective functionthat makes robust optimal the policy represented by the so-called linear saturated feedback control. We regard the problem as a suitable differential game with switching cost and study it in the framework of the viscosity solutions theory for Bellman and Isaacs equations.

We investigate a global-in-time variational approach to abstract evolution by means of the weighted energy-dissipation functionals proposed by Mielke and Ortiz [ESAIM: COCV14 (2008) 494–516]. In particular, we focus on gradient flows in Hilbert spaces. The main result is the convergence of minimizers and approximate minimizers of these functionals to the unique solution of the gradient flow. Sharp convergence rates are provided and the convergence analysis is combined with time-discretization. Applications of the theory to various classes of parabolic PDE problems are presented. In particular, we focus on two examples of microstructure evolution from [S. Conti and M. Ortiz, J. Mech. Phys. Solids56 (2008) 1885–1904.].

We study the stabilization of global solutions of theKawahara (K) equation in a bounded interval, under the effect ofa localized damping mechanism. The Kawahara equation is a modelfor small amplitude long waves. Using multiplier techniques andcompactness arguments we prove the exponential decay of the solutions of the (K) model. The proofrequires of a unique continuation theorem and the smoothing effectof the (K) equation on the real line, which are proved in this work.

One shows that the linearized Navier-Stokes equation in ${\mathcal{O}}{\subset} R^d,\;d \ge 2$ , around an unstable equilibriumsolution is exponentially stabilizable in probability by aninternal noise controller $V(t,\xi)=\displaystyle\sum\limits_{i=1}^{N} V_i(t)\psi_i(\xi)\dot\beta_i(t)$ , $\xi\in{\mathcal{O}}$ , where $\{\beta_i\}^N_{i=1}$ areindependent Brownian motions in a probability space and $\{\psi_i\}^N_{i=1}$ is a system of functions on ${\mathcal{O}}$ withsupport in an arbitrary open subset ${\mathcal{O}}_0\subset {\mathcal{O}}$ . Thestochastic control input $\{V_i\}^N_{i=1}$ is found in feedbackform. One constructs also a tangential boundary noise controllerwhich exponentially stabilizes in probability the equilibriumsolution.

We investigate the evolution of an almost flat membrane driven by competition of the homogeneous, Frank, and bending energies as wellas the coupling of the local order of the constituent molecules of the membrane to its curvature.We propose an alternative to the model in [J.B. Fournier and P. Galatoa, J. Phys. II7 (1997) 1509–1520; N. Uchida, Phys. Rev. E66 (2002) 040902] which replacesa Ginzburg-Landau penalization for the length of theorder parameter by a rigid constraint.We introduce a fully discrete scheme, consisting of piecewise linearfinite elements, show that it is unconditionally stable for a large range of the elastic moduli in the model, and prove its convergence(up to subsequences) thereby proving the existence of a weak solutionto the continuous model. Numerical simulations are included that examine typical model situations, confirm our theory, and explore numerical predictions beyond that theory.

We provide an approximation of Mather variational problem by finite dimensional minimization problems in the framework of Γ-convergence. By a linear programming interpretation as done in [Evans and Gomes, ESAIM: COCV 8 (2002) 693–702] we state a duality theorem for the Mather problem, as well a finite dimensional approximation for the dual problem.

We provide a detailed analysis of the minimizers of the functional $u \mapsto \int_{\Bbb R^n} |\nabla u|^2 + D \int_{\Bbb R^n} |u|^{\gamma}$ , $\gamma \in (0, 2)$ , subject to the constraint $\|u\|_{L^2} = 1$ . This problem, e.g., describes the long-time behavior of the parabolic Anderson in probability theory or ground state solutions of a nonlinear Schrödinger equation. While existence can be proved with standard methods, we show that the usual uniqueness results obtained with PDE-methods can be considerably simplified by additional variational arguments. In addition, we investigate qualitative properties of the minimizers and also study their behavior near the critical exponent 2.

A new approach to irreversible quasistatic fracture growth is given, by means of Young measures.The study concerns a cohesive zone model with prescribed crack path, when the material gives different responses to loading and unloading phases.In the particular situation of constant unloading response,the result contained in [G. Dal Maso and C. Zanini,Proc. Roy. Soc. Edinburgh Sect. A137 (2007) 253–279] is recovered.In this case, the convergence of the discrete time approximationsis improved.

For a class of anisotropic integrodifferential operators ${\cal B}$ arising as semigroup generators of Markov processes, we present a sparse tensor product wavelet compression scheme for the Galerkin finite element discretization of the corresponding integrodifferential equations ${\cal B}$ u = f on [0,1]n with possibly large n. Under certain conditions on ${\cal B}$ , the scheme is of essentially optimal and dimension independent complexity $\mathcal{O}$ (h -1| log h |2(n-1)) without corrupting the convergence or smoothness requirements of the original sparse tensor finite element scheme. If the conditions on ${\cal B}$ are not satisfied, thecomplexity can be bounded by $\mathcal{O}$ (h -(1+ε)), whereε $\ll 1$ tends to zero with increasing number of the wavelets' vanishing moments. Here h denotes the width of the corresponding finite element mesh. The operators under consideration are assumed to be of non-negative (anisotropic) order and admit a non-standard kernel κ $(\cdot,\cdot)$ that can be singular on all secondary diagonals. Practical examples of such operators from Mathematical Finance are given and some numerical results are presented.

This paper is concerned with the numerical approximation of the solutions of a two-fluid two-pressure model used in the modelling of two-phase flows.We present a relaxation strategy for easily dealing with both the nonlinearities associated with the pressure laws and the nonconservative termsthat are inherently present in the set of convective equations and that couple the two phases.In particular, the proposed approximate Riemann solver is given by explicit formulas, preservesthe natural phase space, and exactly captures the coupling waves between the two phases.Numerical evidences are given to corroborate the validity of our approach.