We use a framework that takes into account the effects of deformation of both the solid and fluid in the solidification process, to study the solidification of a semi-infinite layer of fluid. It is shown that the time required for solidification, and the final location of the interface are significantly different form the predictions of the classical Stefan problem. A detailed numerical solution of the initial-boundary value problem is provided for a variety of values for non-dimensional parameters relevant for freezing water.

We develop a numerical analysis of the buoyancy driven natural convection of a fluid in a three dimensional shallow cavity (4.1.1) with a horizontal gradient of temperature along the larger dimension. The fluid is a liquid metal (Prandtl number equal to 0.015) while the Grashof number (Gr) varies in the range 100,000-300,000. The Navier-Stokes equations in vorticity-velocity formulation have been integrated by means of a linearized fully implicit scheme. The evaluation of fractal dimension of the attractors in the phase space has allowed the detection of the chaotic regime.

We deal with the emergence of the horizontal three-dimensional convection flow from an asymptotic mechanical equilibrium in a parallelepipedic box with rigid walls and a very small horizontal temperature gradient. The non-linear stability bound is associated with a variational problem. It is proved that this problem is equivalent to the eigenvalue problem governing the linear stability pf the asymptotic basic conduction state and so the two bounds, the linear one and the non-linear one, coincide.

A vast literature exists on the Benard flow, the vertical thermal convection flow, but almost no result is known on the horizontal counterpart. On account of the wide range of applications in geophysics, astrophysics, metereology, and material science; we think that the horizontal thermal convection flow deserves as much consideration as the Benard problem. The present study is the first step towards the description of the bifurcation pattern of the horizontal thermal convection flow.

We propose a time-advancing scheme for the discretization of the unsteady incompressible Navier-Stokes equations. At any time step, we are able to decouple velocity and pressure by solving some suitable elliptic problems. In particular, the problem related with the determination of the pressure does not require boundary conditions. The divergence free condition is imposed as a penalty term, according to an appropriate restatement of the original equations. Some experiments are carried out by approximating the space variables with the spectral Legendre collocation method.

A generalization of mesoscopic Lattice-Boltzmann models aimed at describing flows with solid/liquid phase transitions is presented. It exhibits lower computational costs with respect to the numerical schemes resulting from differential models, Moreover it is suitable to describe chaotic motions in the mushy zone.

A recent approach to generate a zero divergence velocity field by operating directly on the discretized Navier-Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the methodology suggested by Amit, Hall, and Porsching the feasibility of treating three dimensional flows and multiply connected domains is analyzed. The present model keeps the main features of the classical vector potential method in that it generates a divergence-free velocity field through an algebraic manipulation of the discrete equations.

The non-inertial flow of a shear thinning fluid between intersecting planes is studied using a multi-parameter continuation technique. Unlike the classical linearly viscous fluid, it is found that boundary layers develop even in the case of non-inertial flows in both converging and diverging flow. The boundary layers develop due to the non-linearities in the equation which reflect the fact that the fluid can shear thin. © 1991.