The one-dimensional transient solidification of a binary alloy undergoing shrinkage is well-known as an invaluable benchmark for the testing of numerical codes that model macrosegregation. Here, recent work that considered the small-time behaviour of this problem is extended until complete solidification, thereby determining the solute profile across the entire solidified domain. The small-time solution is used as the initial condition for the numerical integration of a problem having three moving boundaries. Of particular significance is the so-called inverse segregation that is observed at the start of solidification, and the extreme segregation that is observed at the end; in the case of the example presented, which is for the often-cited Al–Cu system, the macrosegregation is found to be positive or negative, depending on whether Scheil’s equation or the lever rule is assumed at the microscale, respectively. The relevance of these results for the modelling of steady-state continuous casting processes – in particular, the phenomenon of centreline segregation – is also discussed.

In the one-dimensional solidification of a binary alloy undergoing shrinkage, there is a relative motion between solid and liquid phases in the mushy zone, leading to the possibility of macrosegregation; thus, the problem constitutes an invaluable benchmark for the testing of numerical codes that model these phenomena. Here, we revisit an earlier obtained solution for this problem, that was posed on a semi-infinite spatial domain and valid for the case of low superheat, with a view to extending it to the more general situation of a finite spatial domain, arbitrarily large superheat and both eutectic and non-eutectic solidification. We find that a similarity solution is available for short times which contains a boundary layer on the liquid side of the mush–liquid interface; this solution is believed to constitute the correct initial condition for the subsequent numerical solution of the full non-similar problem, which is deferred to future work.

The Brinkman model describes flow of fluid in complex porous media with a high-contrast permeability coefficient such that the flow is dominated by Darcy in some regions and by Stokes in others. A weak Galerkin (WG) finite element method for solving the Brinkman equations in two or three dimensional spaces by using polynomials is developed and analyzed. The WG method is designed by using the generalized functions and their weak derivatives which are defined as generalized distributions. The variational form we considered in this paper is based on two gradient operators which is different from the usual gradient-divergence operators for Brinkman equations. The WG method is highly flexible by allowing the use of discontinuous functions on arbitrary polygons or polyhedra with certain shape regularity. Optimal-order error estimates are established for the corresponding WG finite element solutions in various norms. Some computational results are presented to demonstrate the robustness, reliability, accuracy, and flexibility of the WG method for the Brinkman equations.

We critically compare the practicality and accuracy of numerical approximations of phase field models and sharp interface models of solidification. Here we focus on Stefan problems, and their quasi-static variants, with applications to crystal growth. New approaches with a high mesh quality for the parametric approximations of the resulting free boundary problems and new stable discretizations of the anisotropic phase field system are taken into account in a comparison involving benchmark problems based on exact solutions of the free boundary problem.

In this study, we present a new numerical model for crystal growth in a vertical solidification system. This model takes into account the buoyancy induced convective flow and its effect on the crystal growth process. The evolution of the crystal growth interface is simulated using the phase-field method. A semi-implicit lattice kinetics solver based on the Boltzmann equation is employed to model the unsteady incompressible flow. This model is used to investigate the effect of furnace operational conditions on crystal growth interface profiles and growth velocities. For a simple case of macroscopic radial growth, the phase-field model is validated against an analytical solution. The numerical simulations reveal that for a certain set of temperature boundary conditions, the heat transport in the melt near the phase interface is diffusion dominant and advection is suppressed.