An axisymmetric model of the Mush-Chimney system

We develop a model of the mush-chimney system based on the assumption that the order of magnitude of the vertical velocity component in the chimney is much greater than the speed of upward advance of the mush, and making use of the fact that the radius of chimney is much less than the (dimensionless) thickness of the mush. We determine the fluid flow structure in the chimney by utilizing the knowledge of the mass fraction of light constituent, the vertical velocity component, and the pressure that are obtained from the mush. We find a relation between a parameter measuring the ratio of viscous and buoyancy forces in the chimney and the vertical velocity component on the top of the mush, and estimate numerically the value of this velocity.     

Homogenization of compressible two-phase two-component flow in porous media

The paper is devoted to the homogenization of immiscible compressible two-phase two-component flow in heterogeneous porous media. We consider liquid and gas phases, two-component (water and hydrogen) flow in a porous reservoir with periodic microstructure, modeling the hydrogen migration through engineered and geological barriers for a deep repository for radioactive waste. Phase exchange, capillary effects included by the Darcy–Muskat law and Fickian diffusion are taken into account. The hydrogen in the gas phase is supposed compressible and could be dissolved into the water obeying the Henry law. The flow is then described by the conservation of the mass for each component. The microscopic model is written in terms of the phase formulation, i.e. the liquid saturation phase and the gas pressure phase are primary unknowns. This formulation leads to a coupled system consisting of a nonlinear parabolic equation for the gas pressure and a nonlinear degenerate parabolic diffusion–convection equation for the liquid saturation, subject to appropriate boundary and initial conditions. The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system. Under some realistic assumptions on the data, we obtain a nonlinear homogenized problem with effective coefficients which are computed via a cell problem. We rigorously justify this homogenization process for the problem by using the two-scale convergence.     

Three‐dimensional solidification with two possible crystallization states: Existence of solutions with flow in the melt

In this paper we discuss the existence of solutions of a system of nonlinear and singular partial differential equations constituting a phase field model with convection for non‐isothermal solidification/melting of certain metallic alloys in the case where two different kinds of crystallization are possible. Each one of these crystallization states is described by its own phase field, while the liquid phase is described by another one. The model also allows the occurrence of fluid flow in non‐solid regions, which are a priori unknown, and then we have a free‐boundary value problem. Thus, the model relates the evolutions of these three phase fields, the temperature of the solidification/melting process and the fluid flow in non‐solid regions. Copyright © 2009 John Wiley & Sons, Ltd.     

A cellular automaton technique for modelling of a binary dendritic growth with convection

A two-dimensional model for the simulation of a binary dendritic growth with convection has been developed in order to investigate the effects of convection on dendritic morphologies. The model is based on a cellular automaton (CA) technique for the calculation of the evolution of solid/liquid (s/l) interface. The dynamics of the interface controlled by temperature, solute diffusion and Gibbs–Thomson effects, is coupled with the continuum model for energy, solute and momentum transfer with liquid convection. The solid fraction is calculated by a governing equation, instead of some approximate methods such as lever rule method [A. Jacot, M. Rappaz, Acta Mater. 50 (2002) 1909–1926.] or interface velocity method [L. Nastac, Acta Mater. 47 (1999) 4253; L. Beltran-Sanchez, D.M. Stefanescu, Mat. and Mat. Trans. A 26 (2003) 367.]. For the dendritic growth without convection, mesh independency of simulation results is achieved. The simulated steady-state tip velocity are compared with the predicted values of LGK theory [Lipton, M.E. Glicksmanm, W. Kurz, Metall. Trans. 18(A) (1987) 341.] as a function of melt undercooling, which shows good agreement. The growth of dendrite arms in a forced convection has been investigated. It was found that the dendritic growth in the upstream direction was amplified, due to larger solute gradient in the liquid ahead of the s/l interface caused by melt convection. In the isothermal environment, the calculated results under very fine mesh are in good agreement with the Oseen–Ivanstov solution for the concentration-driven growth in a forced flow.     

泥石流固液分相流速计算方法研究

泥石流固液分相流速是泥石流对岸坡、防治结构冲击、磨损机理的核心问题．将泥石流体简化为具有相同粒径的固相和具有相同力学性质的液相,基于泥石流体为沿流动方向的一维两相流体,运用两相流理论建立了泥石流固液分相流速控制方程．构建了泥石流平均压力、彻体力及平均表面力的计算方法,尤其通过浆体的Binhanm体流变方程、Bagnold颗粒相互作用试验成果建立了控制体平均表面力计算方法;建立了固液两相流速比例系数,以及理论固相流速与实际流速的比例系数．据此求解控制方程得到了固液分相流速计算方法,该方法既可同时适用于粘性泥石流和稀性泥石流,也可在泥石流爆发以后通过现场采集沉积物分析反求泥石流爆发期间的分相流速．工程实例分析显示,该方法计算结果与实测结果吻合较好．     

A steady-state Boussinesq-Stefan problem with continuous extraction

Summary One establishes an existence result for the weak solution to a steady-state strongly coupled system between a nonlinear two phases heat equation with convection and the Navier-Stokes equation in the liquid phase. The two phases Rayleigh-Bénard problem is included as the particular case corresponding to a zero extraction velocity.     

沉降过程的扩散机理及其数学表达式

本文分析了固液二相流体中的扩散机理，并从水流的紊动出发，指导出固体颗粒在二维水流中分布的理论公式。     

A phase field approach to solidification and solute separation in water solutions

We propose a phase field model for the solid–liquid phase transition in a water-salt (sodium chloride) solution in the absence of macroscopic motion, under possibly non-isothermal conditions. A thermodynamic approach based on a free energy functional is assumed. The model consists of three evolution equations: a time-dependent Ginzburg–Landau equation for the solid–liquid phase change, a diffusion equation of the Cahn–Hilliard kind for the solute dynamics and the heat equation for the temperature change. The proposed system is aimed to contribute to the modelling of the brine channels formation in the ice of the polar seas.     

An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary

A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l² norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.     

Coupled lattice Boltzmann method for simulating electrokinetic flows: A localized scheme for the Nernst–Plank model

We present a coupled lattice Boltzmann method (LBM) to solve a set of model equations for electrokinetic flows in micro-/nano-channels. The model consists of the Poisson equation for the electrical potential, the Nernst–Planck equation for the ion concentration, and the Navier–Stokes equation for the flows of the electrolyte solution. In the proposed LBM, the electrochemical migration and the convection of the electrolyte solution contributing to the ion flux are incorporated into the collision operator, which maintains the locality of the algorithm inherent to the original LBM. Furthermore, the Neumann-type boundary condition at the solid/liquid interface is then correctly imposed. In order to validate the present LBM, we consider an electro-osmotic flow in a slit between two charged infinite parallel plates, and the results of LBM computation are compared to the analytical solutions. Good agreement is obtained in the parameter range considered herein, including the case in which the nonlinearity of the Poisson equation due to the large potential variation manifests itself. We also apply the method to a two-dimensional problem of a finite-length microchannel with an entry and an exit. The steady state, as well as the transient behavior, of the electro-osmotic flow induced in the microchannel is investigated. It is shown that, although no external pressure difference is imposed, the presence of the entry and exit results in the occurrence of the local pressure gradient that causes a flow resistance reducing the magnitude of the electro-osmotic flow.     

On weakly nonlinear evolution of convective flow in a passive mushy layer

The problem of weakly nonlinear convective flow in a mushy layer, with a permeable mush–liquid interface and constant permeability, is studied under operating conditions for an experiment. A Landau type nonlinear evolution equation for the amplitude of the secondary solutions, which is based on the Landau theory and formulation for the Rayleigh, $R$, number close to its critical value, $R_c$, is developed. Using numerical and analytical methods, the solutions to the evolution equation are calculated for both supercritical and subcritical conditions. We found, in particular, that for $R<R_c$, the amplitude of the secondary solutions decays with time. For $R>R_c$, the tendency for chimney formation in the mushy layer increases with time. In addition, in such a supercritical regime, the basic flow is linearly unstable and we see the presence of steady flow for large values of time. These results suggest a possible slow transition to turbulence in such a flow system.     

Numerical Simulation of Two-Phase Flow through Heterogeneous Porous Media

A mixed finite element method is combined to finite volume schemes on structured and unstructured grids for the approximation of the solution of incompressible flow in heterogeneous porous media. A series of numerical examples demonstrates the effectiveness of the methodology for a coupled system which includes an elliptic equation and a nonlinear degenerate diffusion–convection equation arising in modeling of flow and transport in porous media.     

Taylor dispersion in a packed tube

Presented in this paper is a theoretical analysis for longitudinal scalar spread of mean concentration under a fully developed flow in a tube packed with porous media. A general form of momentum equation for superficial flow in porous media is introduced as a combination of the Navier–Stokes equation and Darcy’s law plus a superficial dispersion term due to phase discontinuity between the fluid flow and solid frame. The analytical solution presented for the fully developed superficial flow includes that for the Poiseulle flow in an evacuated tube as a limiting case. As an extension of Taylor’s classical work on dispersion of soluble matter in solvent flowing slowly through an evacuated tube, a one-dimensional dispersion equation valid for overall environmental assessment of contaminant is rigorously derived by cross-sectionally averaging the superficial mass equation and introducing a closure relation for a new unknown out of the averaging procedure, and corresponding Taylor dispersivity determined is shown to be a generalization of Taylor’s well-known result for the Poiseulle flow.     

Integrated simulation of thermo-solutal convection and pattern formation in directional solidification

Numerical analysis is carried out to examine the effects of thermo-solutal convection on the formation of complex patterns in directionally solidified binary alloys. A finite-difference analysis is used for dynamic modeling of a two-dimensional prototype of the vertical Bridgman system that takes into account heat transfer in the melt, crystal, and the ampoule, as well as the melt flow and solute transport. Actual temperature data from experimental measurements are used for accurately describing the thermal boundary conditions. A range of complex dynamical behavior is predicted in the melt flow due to flow transitions and this is found to be directly related to the spatial patterns observed experimentally in the solidified alloys. The model is applied to single phase solidification in the Al–Cu and Pb–Sn systems to characterize the effect of convection on the macroscopic shape of the interface. The application of the model to hyper-peritectic alloys in the Sn–Cd system shows that the presence of oscillating flow can give rise to a novel convection induced microstructure in which a tree-like primary phase in the center of the sample is embedded in the surrounding peritectic matrix.     

Numerical simulation of immersion quenching process of an engine cylinder head

In this article, we present the numerical simulations of a real cylinder head quench cooling process employing a newly developed boiling phase change model using the commercial CFD code AVL-FIRE v8.5. Separate computational domains constructed for the solid and liquid regions are numerically coupled at the interface of the solid–liquid boundaries using the AVL-Code-Coupling-Interface (ACCI) feature. The boiling phase change process triggered by the dipping hot metal and the ensuing two-phase flow is handled using an Eulerian two-fluid method. Multitude of flow features such as vapor pocket generation, bubble clustering and their disposition, are captured very effectively during the computation, in addition to the variation of the temperature pattern within the solid region. A comparison of the registered temperature readings at different monitoring locations with the numerical results generates an overall very good agreement and indicates the presence of intense non-uniformity in the temperature distribution within the solid. Overall, the predictive capability of the new boiling model is well demonstrated for real-time quenching applications.     

Towards Riccati-Feedback Control of Complex Flows with Moving Interfaces

Problems featuring moving interfaces appear in many applications. They can model solidification and melting of pure materials, crystal growth and other multi-phase problems. The control of the moving interface enables to, for example, influence production processes and, thus, the product material quality. We consider the two-phase Stefan problem that models a solid and a liquid phase separated by the moving interface. In the liquid phase, the heat distribution is characterized by a convection-diffusion equation. The fluid flow in the liquid phase is described by the Navier–Stokes equations which introduces a differential algebraic structure to the system. The interface movement is coupled with the temperature through the Stefan condition, which adds additional algebraic constraints. Our formulation uses a sharp interface representation and we define a quadratic tracking-type cost functional as a target of a control input. We compute an open loop optimal control for the Stefan problem using an adjoint system. For a feedback representation, we linearize the system about the trajectory defined by the open loop control. This results in a linear-quadratic regulator problem, for which we formulate the differential Riccati equation with time varying coefficients. This Riccati equation defines the corresponding feedback gain. Further, we present the feedback formulation that takes into account the structure and the differential algebraic components of the problem. Also, we discuss how the complexities that come, for example, with mesh movements, can be handled in a feedback setting. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exact and approximate expressions for resistance coefficients of a colloidal sphere approaching a solid surface at intermediate Reynolds numbers

A mathematical model has been developed to describe the force of liquid flow acting on a colloidal spherical particle as it approaches a solid surface at intermediate-Reynolds-number-flow regime. The model has incorporated bispherical coordinates to determine a stream function for the flow disturbed by the sphere. The stream function was then used to derive the flow force on the particle as a function of the inter-surface separation distance. The force equation was related to the modified Stokes equation to obtain an exact analytical expression for the correction factor to the Stokes law. Finally, a rational approximation is presented, which is in good agreement with the exact numerical result, and can be readily applied to more general particle–surface interactions involving short-range hydrodynamics associated with colloidal particles in the near vicinity of a large solid collector surface at intermediate Reynolds number of the supporting flow.