This paper deals with the motion of a fluid in a closed loop under the effect of natural convection and a given external heat flux. More precisely, we show that the stationary solutions of a system describing the intermediate asymptotics of the previous problem are structurally linearly unstable.

Low permeability formations are often hydrofractured to increase the production rate of oil and gas. This process creates a thin, but highly permeable, fracture which provides an easy path for oil and gas to flow through the reservoir to the borehole. Here we examine the payoff of hydrofracturing by determining the increased production rate of a hydrofractured well. We find explicit formulas for the steady production rate in the three regimes of small, intermediate, and large (dimensionless) fracture conductivities. Previously, only the formula for the large fracture conductivity case was known.

We assume that Darcy flow pertains throughout the reservoir. Then, the steady fluid flow through the reservoir is governed by Laplace's equation with a second-order boundary condition along the fracture. We first analyze this boundary value problem for the case of small fracture conductivities. An explicit formula for the production rate is obtained for this case, essentially by combining singular perturbation methods with spectral methods in a function space which places the second-order boundary condition on the same footing as Laplace's equation. Next, we re-cast Laplace's equation as a variational principle which has the second-order boundary condition as its natural boundary condition. This allows us to use simple trial functions to derive accurate estimates of the production rate in the intermediate conductivity case. Then, an asymptotic analysis is used to find the production rate for the large fracture conductivity case. Finally, the asymptotic and variationally-derived production rate formulasare compared to exact values of the production rate, which have been obtained numerically.

It may be feasible to create more than a single fracture about a borehole. So we also develop similar asymptotic and variational formulas for the production rate of a well with multiple fractures.

Numerical methods for initial-value problems which develop singularities in finite time are analyzed. The objective is to determine simple strategies which produce the correct asymptotic behaviour and give an accurate approximation of the blow-up time. Fixed step methods for scalar ordinary differential equations are studied first and it is shown that there is a natural embedding of the discrete process in a continuous one. This shows clearly how and why the fixed-step strategy fails. A class of time-stepping strategies that correspond to a time- continuous re-scaling of the underlying differential equation is then proposed; this class is analyzed and criteria established to determine suitable choices for the re-scaling. Finally the ideas are applied to a partial differential equation arising from the study of a fluid with temperature-dependent viscosity. The numerical method involves re-formulating the equationas a moving boundary problem for the peak value and applying the ODE time-steppingstrategies based on this peak value.

This article is concerned with the structure and stability properties of a combustion front that propagates in the axial direction along the surface of a cylindrical solid fuel element. The fuel consists of a mixture of two finely ground metallic powders, which combine upon ignition in a one-step chemical reaction. The reaction is accompanied by a melting process, which in turn enhances the reaction rate. The combustion products are in the solid state. The reaction zone, inside which the melting occurs, is modelled as a front that propagates along the surface of the cylinder. The different modes of propagation that have been observed experimentally (such as single- and multiheaded spin combustion and multiple-point combustion) are explained as the results of bifurcations from a uniformly propagating plane circular front. The stability properties of the various modes are discussed.

This paper is concerned with the mathematical study of a nonlinear system modelling an irreversible phase change problem. Uniqueness of the solution is proved using the accretivity of the system in (L1)2. Expressing one of the two unknowns as an explicit functional of the other reduces the system to a single nonlinear evolution equation and ultimately leads to an existence theorem.

In this paper the existence and uniqueness of the solution of a nonlinear system modelling some irreversible phase changes is established.