Numerical simulation of the horizontal Bridgman growth. Part II: Three-dimensional flow

We study the steady-state three-dimensional flow which occurs in a horizontal crucible of molten metal under the action of a horizontal temperature gradient. The geometry and the boundary conditions are similar to those encountered in the Bridgman growth process of semiconductor crystals. We find that three-dimensional effects can have a dramatic influence upon the flow, which, before the onset of periodic disturbances, differs appreciably from its two-dimensional counterpart. We also investigate the sensitivity of the flow to non-symmetric disturbances.     

Numerical simulation of the horizontal Bridgman growth. Part III: Calculation of the interface

We study the transient motion of the solidification front during the growth of semiconductor crystals in the horizontal Bridgman geometry. The calculation is based on a two-dimensional flow. We use finite elements which deform with the motion of the interface. The energy equation is coupled with the isothermal constraint of the interface in an implicit transient algorithm. Several examples show the oscillatory motion of the interface caused by the periodic flow of the melt, and they reveal the importance of the growth rate on the shape of the interface.     

用LDA测量U形水槽内的振荡流

本文介绍LDA用于测量U形水槽内的振荡流。水槽直管段为正方形,弯管段内侧弯面及外侧弯面是专门计算设计的。对于水平直管段、竖管段及弯管段共七个截面,测量了边界层外振荡流速的振幅分布。对于水平直管段上底面及下底面各三个截面,测量了边界层內振荡流速的振幅分布及相位分布。测量结果和无限长槽道的理论计算结果作了比较,二者基本符合。     

A transonic quasi-3D analysis for gas turbine engines including split-flow capability for turbofans

A numerical approximation is taken to the solution of the complex flows existing in gas turbine engines with transonic blading. The quasi-3D approach decouples the problem into through-flow and blade-to-blade solutions. An industrially practical finite element through-flow solution is developed and for blade-to-blade solutions a transonic finite areas method is utilized. The finite element code developed is capable of operating in an analysis or a design mode. In both modes a dynamic relaxation factor is employed and considerable reduction in solution time can be achieved. Comparisons to streamline curvature methods are carried out for simple analytical and complex industrial problems.     

The addition of quasi-three-dimensional terms into a finite element method for transonic turbomachinery blade-to-blade flows

This paper describes the extension of a purely two-dimensional finite element method for the calculation of transonic turbomachinery blade-to-blade flows to include the quasi-three-dimensional terms. These terms account for the effect of variations in streamline radius, stream-tube height and blade rotation. By approximating the stream surface as a piecewise linear function, then using a local developed cone transformation on an element basis, the finite element equations are shown to remain of the same form as the two-dimensional equations. The numerical results presented demonstrate that the stream-tube height, streamline radius and blade rotation terms must be included if the prediction of the Mach number distribution around a gas turbine blade is to be calculated correctly.     

A finite element procedure for viscoelastic flows

A finite element method for the simulation of viscoelastic flows has been developed. It uses a weak formulation of the method of characteristics to treat the viscoelastic constitutive law. Numerical results in a 4:1 contraction are presented and are discussed with respect to previous computations. New phenomena are put in evidence and new questions are opened in this already controversial problem.     

Finite element method for unsteady MHD flow through pipes with arbitrary wall conductivity

A finite element method is given to obtain the numerical solution of the coupled equations in velocity and magnetic field for unsteady MHD flow through a pipe having arbitrarily conducting walls. Pipes of rectangular, circular and triangular sections have been taken for illustration. Computations have been carried out for different Hartmann numbers and wall conductivity at various time levels. It is found that if the wall conductivity increases, the flux through a section decreases. The same is the effect of increasing the Hartmann number. It is also observed that the steady state is approached at a faster rate for larger Hartmann numbers or larger wall conductivity. Selected graphs are given showing the behaviour of velocity, induced magnetic field and flux across a section.     

Coating flow theory by finite element and asymptotic analysis of the navier-stokes system

Coating flows are laminar free surface flows, preferably steady and two-dimensional, by which a liquid film is deposited on a substrate. Their theory rests on mass and momentum accounting for which Galerkin's weighted residual method, finite element basis functions, isoparametric mappings, and a new free surface parametrization prove particularly well-suited, especially in coping with the highly deformed free boundaries, irregular flow domains, and the singular nature of static and dynamic contact lines where fluid interfaces intersect solid surfaces. Typically, short forming zones of rapidly rearranging two-dimensional flow merge with simpler asymptotic regimes of developing or developed flow upstream and downstream. The two-dimensional computational domain can be shrunk in size by imposing boundary conditions from asymptotic analysis of those regimes or by matching to one-dimensional finite element solutions of asymptotic equations. The theory is laid out with special attention to conditions at free surfaces, contact lines, and open inflow and outflow boundaries. Efficient computation of predictions is described with emphasis on a grand Newton iteration that converges rapidly and brings other benefits. Sample results for curtain coating and roll coating flows of Newtonian liquids illustrate the power and effectiveness of the theory.     

Fundamental solutions and numerical methods for flow problems

An approach for the numerical solution of flow problems based on the concept of fundamental solutions of differential equations is described. This approach uses the finite element methodology but does not rely on the concept of variational principle or that of residuals. The approach is shown to be well-suited for many types of flow problems. Various applications of this approach are discussed in this paper, with particular emphasis placed on the solution of potential flows and viscous flows containing appreciable regions of separation.     

Simulation of forming processes by FEM with a Bingham fluid model

We model the forming process as a fluid flow. A finite element program, FIDAP, which analyses flow problems, was used to calculate velocity and strain rates at points throughout the material during the deformation process. This allows predictions to be made on the shape and quality of the resulting part. The stress-strain relation we used models the plastic flow of metals (Bingham fluids). The FEM approximation of such a fluid is tested by comparing results for a simple analytical example. In forming processes provision must be made for friction between dye and workpiece, and the program was modified accordingly. Two classical ring forming simulations are compared to published results.     

Numerical Analysis of Thermally Coupled Flow Problems with Interfaces and Phase-change Effects

In this work a fixed mesh finite element approach is presented to solve thermally coupled flow problems including moving interfaces between immiscible fluids and phase-change effects. The weak form of the full incompressible Navier-Stokes equations is obtained using a generalized streamline operator (GSO) technique that enables the use of equal order interpolation of the primitive variables of the problem: velocity, pressure and temperature. The interfaces are defined with a mesh of marker points whose motion is obtained applying a Lagrangian scheme. Moreover, a temperature-based formulation is considered to describe the phase-change phenomena. The proposed methodology is used in the analysis of a filling of a step mould and a gravity-driven flow of an aluminium alloy in an obstructed vertical channel.     

Magnetohydrodynamic flow in an infinite channel

The magnetohydrodynamic (MHD) flow of an incompressible, viscous, electrically conducting fluid in an infinite channel, under an applied magnetic field has been investigated. The MHD flow between two parallel walls is of considerable practical importance because of the utility of induction flowmeters. The walls of the channel are taken perpendicular to the magnetic field and one of them is insulated, the other is partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 200. All the infinite integrals obtained are transformed to finite integrals which contain modified Bessel functions of the second kind. So, the difficulties associated with the computation of infinite integrals with oscillating integrands which arise for large M have been avoided. It is found that, as M increases, boundary layers are formed near the nonconducting boundaries and in the interface region for both velocity and magnetic fields, and a stagnant region in front of the conducting boundary is developed for the velocity field. Selected graphs are given showing these behaviours.     

Horizontal convection in a rectangular enclosure driven by a linear temperature profile

The horizontal convection in a square enclosure driven by a linear temperature profile along the bottom boundary is investigated numerically by using a finite difference method. The Prandtl number is fixed at 4.38, and the Rayleigh number Ra ranges from107 to 1011. The convective flow is steady at a relatively low Rayleigh number, and no thermal plume is observed, whereas it transits to be unsteady when the Rayleigh number increases beyond the critical value. The scaling law for the Nusselt number Nu changes from Rossby's scaling Nu ～ Ra~(1/5) in a steady regime to Nu ～ Ra~(1/4) in an unsteady regime, which agrees well with the theoretically predicted results. Accordingly,the Reynolds number Re scaling varies from Re ～ Ra~(3/11) to Re ～ Ra~(2/5). The investigation on the mean flows shows that the thermal and kinetic boundary layer thickness and the mean temperature in the bulk zone decrease with the increasing Ra. The intensity of fluctuating velocity increases with the increasing Ra.     

Magnetohydrodynamic flow in a rectangular duct

The magnetohydrodynamic flow of an incompressible, viscous, electrically conducting fluid in a rectangular duct, with an external magnetic field applied transverse to the flow, has been investigated. One of the duct's boundaries which is perpendicular to the magnetic field is taken partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 100. All the infinite series obtained are transformed to infinite integrals first and then to finite integrals which contain modified Bessel functions of the second kind. In this way, the difficulties associated with the computation of infinite integrals with oscillating integrands and slowly converging infinite series, the convergence of which is further affected for large values of M, have been avoided. It is found that, as M increases, boundary layers are formed near the non-conducting boundaries and in the interface region, and a stagnant region is developed in front of the conducting boundary for velocity field. The maximm value of magnetic field takes place on the conducting part. These behaviours are shown on some graphs.     

A numerical method for the calculation of incompressible,steady, separated flows around aerofoils

The present work deals with the numerical calculation of the incompressible turbulent flow around aerofoils. An orthogonal curvilinear grid of ‘C’ type is used for the solution of the time averaged equations and Reynolds stresses are modelled according to the κ-ε turbulence model. PISO and SIMPLE algorithms are used to solve the strongly coupled system of the derived finite volume equations and convergence is improved by applying the method of variable local underrelaxation factors. Comparisons between the calculated and measured pressure distributions are presented for NACA 0012 and NACA 4412 wing sections. The formation of separation bubbles according to calculations is also shown.     

Simulation of polymeric flows in the injection moulding process

Recent progress in the simulation of polymeric flows of two key problems in the injection moulding process, carried out by a team at Cornell University, is briefly described. For the filling of cooled thin cavities, the fluid is characterized by a power-law viscosity with exponential temperature dependence, and interaction between the transient thermal boundary-layer and the core flow in a domain with moving boundary is essential. The earlier procedure of Hieber and Shen is modified in two aspects: a boundary-integral formulation replaces the finite-element treatment of the pressure, and an ‘energy integral’ approach is used for the transient temperature. The second problem is the steady visco-elastic flow in the juncture region where sudden changes of the geometry and large strain rates occur. The constitutive equation is postulated according to the Leonov model. The main features in the numerical implementation are: integration along a streamline to determine the elastic deformation tensors for a given velocity field, and finite-element treatment (in time-dependent form) of the pressure and fields for given stresses. In an example where the contraction ratio is 7:1, results for nominal Deborah number exceeding 100 show no numerical instability. (However, for this problem, the true Weissenberg number, i.e. the ratio of local first-normal-stress difference to shear stress turns out to be generally O(10).) The predictions also correlate very well with experimental birefringence measurements.     

Horizontal impact of a sphere half immersed in a liquid of finite depth

We consider the problem of bottom influence during horizontal hydrodynamic impact of a spherical solid body of diameter 2a which is half submerged in a liquid layer of finite depth. The sphere is subjected to the action of a shock pulse, as a result of which it acquires an initial translational velocity u directed along the x-axis.The influence of the spherical bottom on the sphere impact phenomenon was first studied by Zhukovskii [1]. Vertical impact in a layer of finite thickness was considered in [2]; horizontal impact in a halfspace was examined in [3].     

Laminar mixed convection on a horizontal plate of finite length in a channel of finite width

The steady two-dimensional laminar mixed-convection flow past a horizontal plate of finite length is analysed for large Péclet numbers, small Prandtl numbers and weak buoyancy effects. The plate is placed in a channel of finite width, with the plane walls of the channel being parallel to the plate. The temperature of the plate is assumed to be constant. The hydrostatic pressure difference across the wake behind the plate is compensated by a perturbation of the inviscid channel flow. This outer flow perturbation affects the temperature distribution in the thermal boundary layer at the plate and the heat transfer rate, respectively. Solutions in closed form are given. The forces acting on the plate due to the potential flow perturbation are also determined.     

Oscillatory and steady convection in a dielectric viscoelastic layer subjected to a temperature gradient in the presence of an electric field

The linear stability, both stationary and oscillatory, of a horizontal layer of a dielectric Oldroyd fluid under the simultaneous action of a vertical a.c. field and a vertical temperature gradient is analysed. The stability regions in the space of the characteristic parameters are determined by means of a finite difference method. The physical mechanisms relevant to the problem are expounded using a heuristic argument.