An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.     

Mathematical Modeling of Czochralski Type Growth Processes for Semiconductor Bulk Single Crystals

This paper deals with the mathematical modeling and simulation of crystal growth processes by the so-called Czochralski method and related methods, which are important industrial processes to grow large bulk single crystals of semiconductor materials such as, e. g., silicon (Si) or gallium arsenide (GaAs) from the melt. In particular, we investigate a recently developed technology in which traveling magnetic fields are applied in order to control the behavior of the turbulent melt flow. Since numerous different physical effects like electromagnetic fields, turbulent melt flows, high temperatures, heat transfer via radiation, etc., play an important role in the process, the corresponding mathematical model leads to an extremely difficult system of initial-boundary value problems for nonlinearly coupled partial differential equations. In this paper, we describe a mathematical model that is under use for the simulation of real-life growth scenarios, and we give an overview of mathematical results and numerical simulations that have been obtained for it in recent years.     

A dynamic subgrid-scale model for the large eddy simulation of stratified flow

A new dynamic subgrid-scale (SGS) model, including subgrid turbulent stress and heat flux models for stratified shear flow is proposed by using Yoshizawa’s eddy viscosity model as a base model. Based on our calculated results, the dynamic subgrid-scale model developed here is effective for the large eddy simulation (LES) of stratified turbulent channel flows. The new SGS model is then applied to the large eddy simulation of stratified turbulent channel flow under gravity to investigate the coupled shear and buoyancy effects on the near-wall turbulent statistics and the turbulent heat transfer at different Richardson numbers. The critical Richardson number predicted by the present calculation is in good agreement with the value of theoretical analysis     

Flow rate measurements in turbulent liquid metal channel flow using Time-of-Flight Lorentz force velocimetry

Non-contact flow control and flow measurements in hot and aggressive metal melts are big challenges in metallurgical applications. Time-of-Flight Lorentz force velocimetry (ToF LFV) is an electromagnetic measurement technique to meet these challenges. Our experimental results demonstrate that this method is well suited to measure flow rate in turbulent liquid metal channel flow without knowledge of both melt and magnetic field properties. Moreover, the measured flow profiles are in very good agreement with predictions of numerical simulations using the commercial program Package FLUENT MHD. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

CFD of turbulent transport of particles behind a backward-facing step using a new model—k-ε-Sp

The k-ε-S_p model, describing two-dimensional gas–solid two-phase turbulent flow, has been developed. In this model, the diffusion flux and slip velocity of solid particles are introduced to represent the particle motion in two-phase flow. Based on this model, the gas–solid two-phase turbulent flow behind a vertical backward-facing step is simulated numerically and the turbulent transport velocities of solid particles with high density behind the step are predicted. The numerical simulation is validated by comparing the results of the numerical calculation with two other two-phase turbulent flow models (k-ε-A_p, k-ε-k_p) by Laslandes and the experimental measurements. This model, not only has the same virtues of predicting the longitudinal transport of the solid particles as the present practical two-phase flow models, but also can predict the lateral transport of the solid particles correctly.     

Particle-Laden Turbulent Channel Flow and Particle-Wall Interactions

The behavior of particle-laden gases in a turbulent channel flow is studied at moderate Reynolds number. Effects of the wallparticle interaction models on the velocity, statistics, and dispersion of the particles are presented. The results were obtained from a direct numerical simulation with particle feedback on the gas phase. The models were found to considerably influence the studied particle properties. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Mixed convection turbulent film condensation on a sphere

The present paper reports a research on condensation heat transfer of an isothermal sphere with an external flow of vapor. The high tangential velocity of the vapor flow is determined from potential flow theory. The transition criterion of the onset turbulence has been given in the local film Reynolds number (Re_Γ). An eddy diffusivity model along with an expression by [H. Kato, N.N. Shiwaki, M. Hirota, On the turbulent heat transfer by free convection from a vertical plate, Int. J. Heat Mass Transfer, 11(1968) 1117–1125] is used to model turbulence. And the local liquid–vapor interfacial shear which occurs for high velocity vapor flow across a sphere surface is defined by the Colburn analogy. The paper then presents analytical analysis for the local dimensionless film thickness and heat transfer characteristics for the film condensation. And a comparison with those generated by previous theoretical of laminar condensation is discussed. The comparison shows the heat transfer coefficient of turbulent film condensation is higher than laminar film condensation under the high vapor velocity.     

Analytical and numerical simulation investigation in effects of radiation and porosity on a non-orthogonal stagnation-point flow towards a stretching sheet

The steady two-dimensional non-orthogonal flow near the stagnation point on a stretching sheet embedded in a porous medium in the presence of radiation effects is studied. Using similarity variables, the nonlinear boundary-layer equations are solved analytically by homotopy perturbation method (HPM) employing Padé technique. Comparison between the results of HPM-Padé solution and numerical simulation as well as some other results which are available in the literature, demonstrates a very good agreement between them and the HPM-Padé solution provides a convenient way to control and adjust the convergence region of a system of nonlinear boundary-layer problems with high accurate. The effect of involved parameters such as striking angle, radiation parameter, porosity parameter and the Prandtl number on flow and heat transfer characteristics have been discussed with more details.     

Effects of CZSi furnace modification on density of grown-in defects

During large diameter Czochralski silicon growth, heat zone and argon flow influence the formation of defects in silicon crystal by changing the distribution of temperature. Different silicon crystals with various density of grown-in defects were grown by replacing the popular heater with the composite heater and changing the popular argon flow into a controlled flow. The experimental results have been explained well by the numeric simulation of argon flow.     

Evolution of weak discontinuities in a non-ideal radiating gas

A method of wavefront analysis is used to study the formation of shock waves in a two-dimensional steady supersonic flow of a non-ideal radiating gas past plane and axisymmetric bodies. The gas is taken to be sufficiently hot for the effect of thermal radiation to be significant, which is, of course, treated by the optically thin approximation to the radiative transfer equation. Transport equations, which lead to the determination of the shock formation distance and also to conditions which insure that no shock will ever evolve on the wavefront, is derived. The influence of the parameter of the non-idealness, upstream flow Mach number in the presence of thermal radiation on the behavior of the wavefront are examined.     

Direct Numerical Simulation of turbulent heat transfer behind a backward-facing step at low Prandtl number

The paper presents Direct Numerical Simulations of the turbulent flow of a low Prandtl number fluid over a backward-facing step with heat transfer. The backward-facing step flow is investigated as a generic configuration for sudden changes in cross section. Several simulations are reported: for isothermal conditions, for heat transfer with the Prandtl number of air, and for heat transfer with the Prandtl number of liquid sodium. The simulation for air is compared to results from literature. The differences induced by reduction of the Prandtl number are then assessed by comparison of the two cases. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Laminar film condensation in a vertical tube in the presence of noncondensable gas

A two-dimensional, steady state model of convective film-wise condensation of a vapor and noncondensable gas mixture flowing downward inside a vertical tube is developed. The mass transfer at the film and gas interface is treated as diffusion controlled process. The noncondensable effect on the condensation is taken into account through boundary layer analysis of species concentration and energy balance. Numerical predictions were obtained for the condensation heat transfer coefficient of turbulent vapor flow associated with laminar condensate. The predictions were compared with the experimental data in the literature to assess the model. Noncondensable mass fraction and vapor–noncondensable mixture temperature were presented in the form of radial and axial profiles.     

Einfluss eines Verbrennungsvorganges auf den Wärme- und Stoffaustausch in einer laminaren Grenzschicht

Summary The influence of combustion on heat and mass transfer is investigated on the following model. A mixture of an inert with a combustible gas (air) flows in steady, laminar flow over a flat plate. A mass flux of gaseous fuel away from the plate surface is produced by some means. Combustion is assumed to occur with very fast reaction rate so that the process is purely controlled by diffusion and the equilibrium is assumed as very close to complete combustion. It is studied under which conditions the combustion occurs at the surface or when the flame is displaced into the boundary layer. The influence of combustion on the heat transfer from a hot gas to the plate surface is calculated, for the condition that combustion occurs at the surface.      

Statistical modeling of compressible turbulent channel flow

Several questions that are relevant to turbulence modeling are addressed on the basis of recently obtained direct numerical simulation results of turbulent supersonic channel flow. In particular, this concerns the turbulence frequency production mechanism, wall damping effects on turbulence model parameters, and the relevance of compressibility effects. Limited support is found for usually applied models for the turbulence frequency production and wall damping effects. In contrast to that it is shown that turbulence frequency production mechanisms and wall damping effects may be explained very well on the basis of a frequency scaling that characterizes mean flow changes. The influence of compressibility is found to be relevant. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Prediction of Aerodynamic Noise Generated from an Aircraft in Low Mach Number Flight

The paper describes numerical prediction of aerodynamic noise generated from an Aircraft. Simulation of turbulent flow is done solving the incompressible Navier-Stokes equation, where turbulence is modeled applying the orthogonal subgrid scale (OSGS) method with dynamical subscales. Because of comparison, the same simulation is done using the LES (Large Eddy simulation). It is shown how simulation of turbulent flow affects the prediction of acoustic sources calculated using Lighthill's analogy. Translation from the time to frequency domain is done through DFT (Direct Fourier Transform), which gives smaller usage of memory. Acoustic sources are used in the inhomogeneous Helmholtz equation to simulate pressure wave propagation in the domain. It is shown that OSGS with dynamical subscales gives better representation of the spectrum. Overall, better prediction of energy transfer across large and small eddies will give better allocation and presentation of acoustics sources. These sources will change wave propagation of the pressure in the acoustic field. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of internal flow fields of swirl coaxial injector in a hot environment

Based on the fractional volume of fluid (VOF), a pure Eulerian model for defining and capturing the gas/liquid interface is developed in this paper. This model can describe gas/liquid interface in high refinement, which is better than the original VOF methodology. To validate the proposed model and the algorithm, the computational code is employed to predict the flow performance in a cylindrical swirl injector under cold-flow condition, and the predicted results agree well with experimental measurements. Furthermore, the proposed model is used to simulate gas-liquid reacting flows inside a gas/liquid coaxial swirl injector operating in a hot environment. The turbulent combustion process is simulated with the k−ε−f−g model. The numerical simulation is carried out under actual operating condition of the coaxial injector. The injector performances, such as liquid film thickness, liquid film injection velocity, spray angle, pressure drop, are obtained based on the detailed information of the internal flow field. The predicted results also show that droplets are shed from the liquid film in the recess cup of the coaxial injector because of the large velocity gradient between the gas and liquid streams, and a burning area, which is characterized by high temperature, is present inside the injector.     

Mathematical modeling of the processes occurring during high-temperature spray coating

We give a generalization of the mathematical model of the processes occurring in the preparation of gasthermal coatings in order to compute the radiation energy, thermal effects of turbulent flow of the gas jet of a plasmotron, spraying distance and required degree of preliminary heating of the base. The boundary conditions obtained describe the radiational/convective heat exchange of bodies with a medium through thin coatings taking account of the speed of flow of the gas jet.Translated fromMatematicheskie Metody i Fiziko-Mekhanicheskie Polya, Issue 33, 1991, pp. 13–18.     

Predicting gas–liquid flow in a mechanically stirred tank

Computational fluid dynamics (CFD) provides a method for investigating the highly complex fluid flow in mechanically stirred tanks. Although there are quite a number of papers in the literature describing CFD methods for modelling stirred tanks, most only consider single-phase flow. However, multiphase mixtures occur very frequently in the process industries, and these are more complex situations for which modelling is not as well developed. This paper reports on progress in developing CFD simulations of gas–liquid mixing in a baffled stirred tank. The model is three-dimensional and the impeller region is explicitly included using a Multiple Frames of Reference method to account for the relative movement between impeller and baffles. Fluid flow is calculated with a turbulent two-fluid model using a finite-volume method. Several alternative treatments of the multiphase equations are possible, including various expressions for drag and dispersion forces, and a number of these have been tested. Variation in bubble size due to coalescence and break-up is also modelled. The CFD simulation method has been used to model a gas-sparged tank equipped with a Rushton turbine, and simulation results are compared with experimental data. Results to date show the correct pattern of gas distribution and the correct trends in local bubble size in the tank. Further work is needed to improve the quantitative agreement with experimental data.     

Assessment of turbulence modeling for gas flow in two-dimensional convergent–divergent rocket nozzle

In the present study, the turbulent gas flow dynamics in a two-dimensional convergent–divergent rocket nozzle is numerically predicted and the associated physical phenomena are investigated for various operating conditions. The nozzle is assumed to have impermeable and adiabatic walls with a flow straightener in the upstream side and is connected to a plenum surrounding the nozzle geometry and extended in the downstream direction. In this integrated component model, the inlet flow is assumed a two-dimensional, steady, compressible, turbulent and subsonic. The physics based mathematical model of the considered flow consists of conservation of mass, momentum and energy equations subject to appropriate boundary conditions as defined by the physical problem stated above. The system of the governing equations with turbulent effects is solved numerically using different turbulence models to demonstrate their numerical accuracy in predicting the characteristics of turbulent gas flow in such complex geometry. The performance of the different turbulence models adopted has been assessed by comparing the obtained results of the static wall pressure and the shock position with the available experimental and numerical data. The dimensionless shear stress at the nozzle wall and the separation point are also computed and the flow field is illustrated. The various implemented turbulence models have shown different behavior of the turbulent characteristics. However, the shear-stress transport (SST) k–ω model exhibits the best overall agreement with the experimental measurements. In general, the proposed numerical procedure applied in the present paper shows good capability in predicting the physical phenomena and the flow characteristics encountered in such kinds of complex turbulent flow.