A new family of mathematical functions to fit longitudinal growth data is described. All members derive from the differential equation dh/dt = s(t). (h1-h) where h1 is adult size and s(t) is a function of time. The form of s(t) is given by one of many functions, all solutions of differential equations, thus generating a family of different models. Three versions were compared. All were superior to previously described models. Model 1, in which s(t) was defined by ds/dt = (s1 - s)(s - s0) was especially accurate and robust, containing only five parameters to describe growth in stature from age two to maturity. Derived "biological" parameters such as Peak Height Velocity were very consistent between these three members of the family but, in some cases, differed signficantly from previous estimates.

No abstract

No abstract

The migration of liquids in porous media, such as sand, has been commonly considered at high saturation levels with liquid pathways at pore dimensions. In this Letter, we reveal a low saturation regime observed in our experiments with droplets of extremely low volatility liquids deposited on sand. In this regime, the liquid is mostly found within the grain surface roughness and in the capillary bridges formed at the contacts between the grains. The bridges act as variable-volume reservoirs and the flow is driven by the capillary pressure arising at the wetting front according to the roughness length scales. We propose that this migration (spreading) is the result of interplay between the bridge volume adjustment to this pressure distribution and viscous losses of a creeping flow within the roughness. The net macroscopic result is a special case of nonlinear diffusion described by a superfast diffusion equation for saturation with distinctive mathematical character. We obtain solutions to a moving boundary problem defined by superfast diffusion equation that robustly convey a time power law of spreading as seen in our experiments.

Abstract. This article describes a number of velocity-based moving mesh numerical methods for multidimensional nonlinear time-dependent partial differential equations (PDEs). It consists of a short historical review followed by a detailed description of a recently developed multidimensional moving mesh finite element method based on conservation. Finite element algorithms are derived for both mass-conserving and non mass-conserving problems, and results shown for a number of multidimensional nonlinear test problems, including the second order porous medium equation and the fourth order thin film equation as well as a two-phase problem. Further applications and extensions are referenced.

No abstract

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold.

Abstract. Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving-point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly. Our approach is also well suited to capture waiting-time behaviour efficiently. A finite-difference moving-point scheme is derived and applied in a simplified context (continental radially symmetrical shallow ice approximation). The scheme, which is inexpensive, is verified by comparing the results with steady states obtained from an analytic solution and with exact moving-margin transient solutions. In both cases the scheme is able to track the position of the ice sheet margin with high accuracy.