The classical model for studying one-phase Hele-Shaw flows is based on a highly nonlinear moving boundary problem with the fluid velocity related to pressure gradients via a Darcy-type law. In a standard configuration with the Hele-Shaw cell made up of two flat stationary plates, the pressure is harmonic. Therefore, conformal mapping techniques and boundary integral methods can be readily applied to study the key interfacial dynamics, including the Saffman–Taylor instability and viscous fingering patterns. As well as providing a brief review of these key issues, we present a flexible numerical scheme for studying both the standard and nonstandard Hele-Shaw flows. Our method consists of using a modified finite-difference stencil in conjunction with the level-set method to solve the governing equation for pressure on complicated domains and track the location of the moving boundary. Simulations show that our method is capable of reproducing the distinctive morphological features of the Saffman–Taylor instability on a uniform computational grid. By making straightforward adjustments, we show how our scheme can easily be adapted to solve for a wide variety of nonstandard configurations, including cases where the gap between the plates is linearly tapered, the plates are separated in time, and the entire Hele-Shaw cell is rotated at a given angular velocity.

We consider the optimal portfolio and consumption problem for a jump-diffusion process with regime switching. Under the criterion of maximizing the expected discounted total utility of consumption, two methods, namely, the dynamic programming principle and the stochastic maximum principle, are used to obtain the optimal result for the general objective function, which is the solution to a system of partial differential equations. Furthermore, we investigate the power utility as a specific example and analyse the existence and uniqueness of the optimal solution. Under the constraints of no-short-selling and nonnegative consumption, closed-form expressions for the optimal strategy and the value function are derived. Besides, some comparisons between the optimal results for the jump-diffusion model and the pure diffusion model are carried out. Finally, we discuss our optimal results in some special cases.

We consider the explicit solution to the axisymmetric diffusion equation. We recast the solution in the form of a Mellin inversion formula, and outline a method to compute a formula for $u(r,t)$ as a series using the Cauchy residue theorem. As a consequence, we are able to represent the solution to the axisymmetric diffusion equation as a rapidly converging series.

We investigate the interaction between a singular surface and a strong shock in the self-gravitating interstellar gas clouds with the assumption of spherical symmetry. Using the method of the Lie group of transformations, a particular solution of the flow variables and the cooling–heating function for an infinitely strong shock is obtained. This paper explores an application of the singular surface theory in the evolution of an acceleration wave front propagating through an unperturbed medium. We discuss the formation of an acceleration, considering the cases of compression and expansion waves. The influence of the cooling–heating function on a shock formation is explained. The results of a collision between a strong shock and an acceleration wave are discussed using the Lax evolutionary conditions.

We study some properties of integro splines. Using these properties, we design an algorithm to construct splines $S_{m+1}(x)$ of neighbouring degrees to the given spline $S_m(x)$ with degree m. A local integro-sextic spline is constructed with the proposed algorithm. The local integro splines work efficiently, that is, they have low computational complexity, and they are effective for use in real time. The construction of nonlocal integro splines usually leads to solving a system of linear equations with band matrices, which yields high computational costs.