scholarly journals On the thermodynamic boundary conditions of a solidifying mushy layer with outflow

2014 ◽  
Vol 762 ◽  
Author(s):  
David W. Rees Jones ◽  
M. Grae Worster
Keyword(s):  
Transition Region ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Darcy Number ◽  
Forced Flow ◽  
Mushy Layer ◽  
Liquid Region ◽  
Precise Knowledge ◽  
Steady Solutions ◽  
Velocity Changes

AbstractThe free-boundary problem between a liquid region and a mushy layer (a reactive porous medium) must respect both thermodynamic and fluid dynamical considerations. We develop a steady two-dimensional forced-flow configuration to investigate the thermodynamic condition of marginal equilibrium that applies to a solidifying mushy layer with outflow and requires that streamlines are tangent to isotherms at the interface. We show that a ‘two-domain’ approach in which the mushy layer and liquid region are distinct domains is consistent with marginal equilibrium by extending the Stokes equations in a narrow transition region within the mushy layer. We show that the tangential fluid velocity changes rapidly in the transition region to satisfy marginal equilibrium. In convecting mushy layers with liquid channels, a buoyancy gradient can drive this tangential flow. We use asymptotic analysis in the limit of small Darcy number to derive a regime diagram for the existence of steady solutions. Thus we show that marginal equilibrium is a robust boundary condition and can be used without precise knowledge of the fluid flow near the interface.

Numerical Simulation of the Interplay of Electrical Double Layers, Electrode Reactions, and Pressure-Driven Flows in Microchannels

2008 ◽  
Author(s):  
Bettina Wa¨lter ◽  
Peter Ehrhard
Keyword(s):  
Flow Field ◽  
Stokes Equations ◽  
Electrical Potential ◽  
Planck Equation ◽  
Forced Flow ◽  
Double Layers ◽  
Rectangular Microchannel ◽  
Small Width ◽  
Electrode Reactions ◽  
Precise Knowledge

We investigate the influence of flow field and electrode reactions onto an electrical double layer (EDL), which is located in the immediate vicinity of the walls of a rectangular microchannel. The precise knowledge of the EDL appears to be important for many technical applications in microchannels of small width, since the electrokinetic effects, as electroosmosis or electrophoresis, in such cases depend on the detailed charge distribution. The mathematical model for the numerical treatment relies on a first–principle description of the EDL and the electrical forces caused by the electrical field between internal electrodes. Hence, the so–called Debye–Hu¨ckel approximation is avoided. The governing system of equations consists of a Poisson equation for the electrical potential, the Navier–Stokes equations for the flow field, species transport equations, based on the Nernst–Planck equation, and a model for the electrode reactions, based on the Butler–Volmer equation. The simulations are time–dependent and two–dimensional (plane) in nature and employ a finite–volume method. It is discussed, e.g., how the thickness of the EDL expands at the stagnation point of a forced flow, as the velocity (or Reynolds number) is increased. Furthermore, the effect of electrode reactions on the ionic strength and, hence, on the EDL and the electroosmotic flow, are discussed.

Weak convection, liquid inclusions and the formation of chimneys in mushy layers

1999 ◽  
Vol 388 ◽  
pp. 197-215 ◽  
Author(s):  
T. P. SCHULZE ◽  
M. GRAE WORSTER
Keyword(s):  
Solid Fraction ◽  
Stokes Equations ◽  
Numerical Study ◽  
Solidification Rate ◽  
Mushy Region ◽  
Mushy Layer ◽  
Liquid Region ◽  
Interfacial Conditions ◽  
Buoyancy Driven Convection ◽  
Steady Convection

We present a numerical study of steady convection in a two-dimensional mushy layer during solidification of a binary mixture at a constant speed V. The mushy layer is modelled as a reactive porous medium whose permeability is a function of the local solid fraction. The flow in the liquid region above the mushy layer is modelled using the Stokes equations (i.e. the Prandtl number is taken to be infinite). The calculations follow the development of buoyancy-driven convection as the flow amplitude is increased to the level where the solid fraction is driven to zero at some point within the mushy region. We show that this event cannot occur before the local buoyancy-driven volume flux exceeds the solidification rate V. Further increases in the flow amplitude lead to the formation of a region with negative solid fraction, indicating the need to switch from the Darcy approximation to the Stokes flow approximation. These regions ultimately become what are known as chimneys. We exhibit solutions which give the detailed structure of the temperature, solute, flow and solid fraction fields within the mushy layer. A key finding of the numerics is that these fledgling chimneys emerge from the interior of the mushy layer, rather than eating their way down from the top of the layer, as the amplitude of the steady convection is increased. We discuss some qualitative features of the resulting liquid inclusions and, in the light of these, reassess the interfacial conditions between mushy and liquid regions.

Permeability in the Mushy Zone

2010 ◽  
Vol 649 ◽  
pp. 399-408 ◽  
Author(s):  
R.G. Erdmann ◽  
D.R. Poirier ◽  
A.G. Hendrick
Keyword(s):  
Mushy Zone ◽  
Volume Fraction ◽  
Primary Dendrite ◽  
Fluid Velocity ◽  
Interdendritic Liquid ◽  
Creeping Flow ◽  
Liquid Region ◽  
Local Volume ◽  
Local Volume Fraction ◽  
Deterministic Function

When modeled at macroscopic length scales, the complex dendritic network in the solid-plus-liquid region of a solidifying alloy (the “mushy zone”) has been modeled as a continuum based on the theory of porous media. The most important property of a porous medium is its permeability, which relates the macroscopic pressure gradient to the throughput of fluid flow. Knowledge of the permeability of the mushy zone as a function of the local volume-fraction of liquid and other morphological parameters is thus essential to successfully modeling the flow of interdendritic liquid during alloy solidification. In current continuum models, the permeability of the mushy zone is given as a deterministic function of (1) the local volume fraction of liquid and (2) a characteristic length scale such as the primary dendrite arm spacing or the reciprocal of the specific surface area of the solid-liquid interface. Here we first provide a broad overview of the experimental data, mesoscale numerical flow simulations, and resulting correlations for the deterministic permeability of both equiaxed and columnar mushy zones. A extended view of permeability in mushy zones which includes the stochastic nature of permeability is discussed. This viewpoint is the result of performing extensive numerical simulations of creeping flow through random microstructures. The permeabilities obtained from these simulations are random functions with spatial autocorrelation structures, and variations in the local permeability are shown to have dramatic effects on the flow patterns observed in such microstructures. Specifically, it is found that “lightning-like” patterns emerge in the fluid velocity and that the flows in such geometries are strongly sensitive to small variations in the solid structure. We conclude with a comparison of deterministic and stochastic permeabilities which suggests the importance of incorporating stochastic descriptions of the permeability of the mushy zone in solidification modeling.

Modeling of the Vortex-Structure in a Particle-Laden Mixing-Layer

2005 ◽  
Author(s):  
Jens A. Melheim ◽  
Stefan Horender ◽  
Martin Sommerfeld
Keyword(s):  
Vortex Structure ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Mixing Layer ◽  
Equation Model ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Calculated Concentration ◽  
Langevin Model ◽  
Particle Velocities

Numerical calculations of a particle-laden turbulent horizontal mixing-layer based on the Eulerian-Lagrangian approach are presented. Emphasis is given to the determination of the stochastic fluctuating fluid velocity seen by the particles in anisotropic turbulence. The stochastic process for the fluctuating velocity is a “Particle Langevin equation Model”, based on the Simplified Langevin Model. The Reynolds averaged Navier-Stokes equations are closed by the standard k-epsilon turbulence model. The calculated concentration profile and the mean, the root-mean-square (rms) and the cross-correlation terms of the particle velocities are compared with particle image velocimetry (PIV) measurements. The numerical results agree reasonably well with the PIV data for all of the mentioned quantities. The importance of the modeled vortex structure “seen” by the particles is discussed.

Effects of Disturbance of Current Field on Power Characteristics of a Floating Type Pitch-Controlled VAMT in a Real Sea

2016 ◽  
Author(s):  
Tomoki Ikoma ◽  
Koichi Masuda ◽  
Hiroaki Eto ◽  
Chang-Kyu Rheem ◽  
Osamu Enomoto
Keyword(s):  
Fluid Velocity ◽  
Vertical Axis ◽  
Ocean Waves ◽  
Floating Body ◽  
Maximum Output ◽  
Power Characteristics ◽  
Turbine Performance ◽  
Velocity Changes ◽  
Floating Type ◽  
Japanese Islands

While a type of marine turbine for tidal current generation can be chosen from several types, a vertical axis marine turbine (VAMT) should be better in Japan because sea areas around Japanese islands where current velocity is sufficient are limited. This study conducted a sea test of a VAMT of a floating type installed with six straight pitch-controllable blades. The cycloidal mechanism was adapted for the pitch control. The purpose of the study is to understand effects of ocean waves and motion of a floating body on turbine performance and behaviours of the VAMT in unideal current conditions. Besides, the data taken should be effective to consider that effects in order to design VAMTs. The setup with the setting angle of −30 degrees suggested highest performance from the sea tests, then 15% in maximum turbine power and maximum output was 40W. Ocean waves strongly affected on the turbine performance because fluid velocity changes due to ocean waves and it was unable to neglect the variation of the velocity in spite of small. The characteristics of the turbine sensitively varied because of ocean waves. The results suggested that during accelerating and decelerating incoming fluid speed, characteristics of the turbine were different in each case.

scholarly journals Modeling Free Surface Flows Using Stabilized Finite Element Method

2018 ◽  
Vol 2018 ◽  
pp. 1-9 ◽  
Author(s):  
Deepak Garg ◽  
Antonella Longo ◽  
Paolo Papale
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Free Surface ◽  
Numerical Scheme ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Computer Code ◽  
Free Surface Flows ◽  
Surface Flows ◽  
Element Method

This work aims to develop a numerical wave tank for viscous and inviscid flows. The Navier-Stokes equations are solved by time-discontinuous stabilized space-time finite element method. The numerical scheme tracks the free surface location using fluid velocity. A segregated algorithm is proposed to iteratively couple the fluid flow and mesh deformation problems. The numerical scheme and the developed computer code are validated over three free surface problems: solitary wave propagation, the collision between two counter moving waves, and wave damping in a viscous fluid. The benchmark tests demonstrate that the numerical approach is effective and an attractive tool for simulating viscous and inviscid free surface flows.

On the drag-out problem in liquid film theory

2008 ◽  
Vol 617 ◽  
pp. 283-299 ◽  
Author(s):  
E. S. BENILOV ◽  
V. S. ZUBKOV
Keyword(s):  
Surface Tension ◽  
Liquid Film ◽  
Constant Velocity ◽  
Stokes Equations ◽  
Film Theory ◽  
Infinite Plate ◽  
Self Similar Solution ◽  
Self Similar ◽  
The Mean ◽  
Steady Solutions

We consider an infinite plate being withdrawn (at an angle α to the horizontal, with a constant velocity U) from an infinite pool of viscous liquid. Assuming that the effects of inertia and surface tension are weak, Derjaguin (C. R. Dokl. Acad. Sci. URSS, vol. 39, 1943, p. 13.) conjectured that the ‘load’ l, i.e. the thickness of the liquid film clinging to the plate, is l=(μU/ρgsinα)1/2, where ρ and μ are the liquid's density and viscosity, and g is the acceleration due to gravity.In the present work, the above formula is derived from the Stokes equations in the limit of small slopes of the plate (without this assumption, the formula is invalid). It is shown that the problem has infinitely many steady solutions, all of which are stable – but only one of these corresponds to Derjaguin's formula. This particular steady solution can only be singled out by matching it to a self-similar solution describing the non-steady part of the film between the pool and the film's ‘tip’.Even though the near-pool region where the steady state has been established expands with time, the upper, non-steady part of the film (with its thickness decreasing towards the tip) expands faster and, thus, occupies a larger portion of the plate. As a result, the mean thickness of the film is 1.5 times smaller than the load.

Two-dimensional flow of a viscous fluid in a channel with porous walls

1991 ◽  
Vol 227 ◽  
pp. 1-33 ◽  
Author(s):  
Stephen M. Cox
Keyword(s):  
Finite Time ◽  
Stokes Equations ◽  
Pitchfork Bifurcation ◽  
Time Dependent ◽  
High Rate ◽  
Navier Stokes ◽  
Recent Theory ◽  
Navier Stokes Equations ◽  
Two Dimensional Flow ◽  
Steady Solutions

We consider the flow of a viscous incompressible fluid in a parallel-walled channel, driven by steady uniform suction through the porous channel walls. A similarity transformation reduces the Navier-Stokes equations to a single partial differential equation (PDE) for the stream function, with two-point boundary conditions. We discuss the bifurcations of the steady solutions first, and show how a pitchfork bifurcation is unfolded when a symmetry of the problem is broken.Then we describe time-dependent solutions of the governing PDE, which we calculate numerically. We analyse these unsteady solutions when there is a high rate of suction through one wall, and the other wall is impermeable: there is a limit cycle composed of an explosive phase of inviscid growth, and a slow viscous decay. The inviscid phase ‘almost’ has a finite-time singularity. We discuss whether solutions of the governing PDE, which are exact solutions of the Navier-Stokes equations, may develop mathematical singularities in a finite time.When the rates of suction at the two walls are equal so that the problem is symmetrical, there is an abrupt transition to chaos, a ‘homoclinic explosion’, in the time-dependent solutions as the Reynolds number is increased. We unfold this transition by perturbing the symmetry, and compare direct numerical integrations of the governing PDE with a recent theory for ‘Lorenz-like’ dynamical systems. The chaos is found to be very sensitive to symmetry breaking.

scholarly journals Mathematical Model for the Electrokinetic Instability of Electrolyte Flow

2019 ◽  
Vol 224 ◽  
pp. 02003
Author(s):  
Andrey Shobukhov
Keyword(s):  
Electric Potential ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Fluid Velocity ◽  
Steady State Solution ◽  
Stability Loss ◽  
Navier Stokes ◽  
Ion Distribution ◽  
Navier Stokes Equations ◽  
Dilute Aqueous Solution

We study a one-dimensional model of the dilute aqueous solution of KCl in the electric field. Our model is based on a set of Nernst-Planck-Poisson equations and includes the incompressible fluid velocity as a parameter. We demonstrate instability of the linear electric potential variation for the uniform ion distribution and compare analytical results with numerical solutions. The developed model successfully describes the stability loss of the steady state solution and demonstrates the emerging of spatially non-uniform distribution of the electric potential. However, this model should be generalized by accounting for the convective movement via the addition of the Navier-Stokes equations in order to substantially extend its application field.

Instability and Transition of a Plane Bubble Plume

2002 ◽  
Author(s):  
Ophe´lie Caballina ◽  
Eric Climent ◽  
Jan Dusˇek
Keyword(s):  
Stokes Equations ◽  
Fluid Layer ◽  
Fluid Velocity ◽  
Continuous Phase ◽  
Collective Effects ◽  
Recent Experimental Data ◽  
Fluid Viscosity ◽  
Bubble Plume ◽  
Carrier Fluid ◽  
Bubble Cloud

When bubbles are continuously released from a located source at the bottom of a fluid layer initially at rest, a plume is produced. The motion of the carrier fluid is initiated and driven by buoyancy of the bubble cloud. In the present study, a detailed analysis of the bubble plume transition is investigated. The continuous phase flow is obtained by direct numerical resolution of Navier-Stokes equations forced by the presence of bubbles. Collective effects induced by the presence of bubbles are modelled by a spatio-temporal distribution of momentum. Time evolution of the dispersed phase is solved by lagrangian tracking of all the bubbles. Focused on the description of plume transition, several configurations (plume widths, fluid viscosity, injection rate) are investigated. During the laminar ascension of the plume, fluid velocity profiles can be non-dimensionalised on a single auto-similar evolution. Dimensional analysis provides a prediction of the limit rising velocity of the plume top. This prediction has been confirmed by our numerical simulations. Furthermore, our first results point out the symmetry breaking induced by plume instability which appears beyond a critical transition height. Various data show that the Grashof number based on injection conditions is the key parameter to predict the transition of the plume. Our results agree very well with recent experimental data. Comparison with experiments on thermal plumes in air shows that the bubble plume is more unstable. This feature should be related to the lack of diffusion in the lagrangian transport of density gradient by the bubble cloud and to the slip velocity between the two phases.

Sign in / Sign up

Export Citation Format

Share Document