In this paper, we investigate the thermal evolution in a one-dimensional bagasse stockpile. The mathematical model involves four unknowns: the temperature, oxygen content, liquid water content and water vapour content. We first nondimensionalize the model to identify dominant terms and so simplify the system. We then calculate solutions for the approximate and full system. It is shown that under certain conditions spontaneous combustion will occur. Most importantly, we show that spontaneous combustion can be avoided by sequential building. To be specific, in a situation where, say, a $4.7\,$m stockpile can spontaneously combust, we could construct a $3\,$m pile and then some days later add another $1.7\,$m to produce a stable $4.7\,$m pile.

The classical half-line Robin problem for the heat equation may be solved via a spatial Fourier transform method. In this work, we study the problem in which the static Robin condition $$bq(0,t) + {q_x}(0,t) = 0$$ is replaced with a dynamic Robin condition; $$b = b(t)$$ is allowed to vary in time. Applications include convective heating by a corrosive liquid. We present a solution representation and justify its validity, via an extension of the Fokas transform method. We show how to reduce the problem to a variable coefficient fractional linear ordinary differential equation for the Dirichlet boundary value. We implement the fractional Frobenius method to solve this equation and justify that the error in the approximate solution of the original problem converges appropriately. We also demonstrate an argument for existence and unicity of solutions to the original dynamic Robin problem for the heat equation. Finally, we extend these results to linear evolution equations of arbitrary spatial order on the half-line, with arbitrary linear dynamic boundary conditions.

We consider an extension of the classical Fisher–Kolmogorov equation, called the “Fisher–Stefan” model, which is a moving boundary problem on $0<x<L(t)$. A key property of the Fisher–Stefan model is the “spreading–vanishing dichotomy”, where solutions with $L(t)>L_{\text{c}}$ will eventually spread as $t\rightarrow \infty$, whereas solutions where $L(t)\ngtr L_{\text{c}}$ will vanish as $t\rightarrow \infty$. In one dimension it is well known that the critical length is $L_{\text{c}}=\unicode[STIX]{x1D70B}/2$. In this work, we re-formulate the Fisher–Stefan model in higher dimensions and calculate $L_{\text{c}}$ as a function of spatial dimensions in a radially symmetric coordinate system. Our results show how $L_{\text{c}}$ depends upon the dimension of the problem, and numerical solutions of the governing partial differential equation are consistent with our calculations.

Pressure-correction projection finite element methods (FEMs) are proposed to solve nonstationary natural convection problems in this paper. The first-order and second-order backward difference formulas are applied for time derivative, the stability analysis and error estimates of the semi-discrete schemes are presented using energy method. Compared with characteristic variational multiscale FEM, pressure-correction projection FEMs are more efficient and unconditionally energy stable. Ample numerical results are presented to demonstrate the effectiveness of the pressure-correction projection FEMs for solving these problems.

Weighted interior penalty discontinuous Galerkin method is developed to solve the two-dimensional non-equilibrium radiation diffusion equation on unstructured mesh. There are three weights including the arithmetic, the harmonic, and the geometric weight in the weighted discontinuous Galerkin scheme. For the time discretization, we treat the nonlinear diffusion coefficients explicitly, and apply the semi-implicit integration factor method to the nonlinear ordinary differential equations arising from discontinuous Galerkin spatial discretization. The semi-implicit integration factor method can not only avoid severe timestep limits, but also takes advantage of the local property of DG methods by which small sized nonlinear algebraic systems are solved element by element with the exact Newton iteration method. Numerical results are presented to demonstrate the validity of discontinuous Galerkin method for high nonlinear and tightly coupled radiation diffusion equation.