A remedy for numerical oscillations in weakly compressible smoothed particle hydrodynamics

Weakly Compressible Smoothed Particle Hydrodynamics (WCSPH) can lead to non‐physical oscillations in the pressure and density fields when simulating incompressible flow problems. This in turn may result in tensile instability and sometimes divergence. In this paper, it is shown that this difficulty originates from the specific form of spatial discretization used for the pressure term when solving the mass conservation equation. After describing the pressure–velocity decoupling problem associated with the so‐called colocated grid methods, a modified approach is presented that overcomes this problem using a different discretization scheme for the second derivative of pressure. The modified scheme is employed for solving a number of benchmark problems including both single‐phase and two‐phase test cases. Copyright © 2010 John Wiley & Sons, Ltd.     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.     

Exact solutions of linearized equations of convection of a weakly compressible fluid

A mathematical model of fluid convection under microgravity conditions is considered. The equation of state is used in a form that allows considering the fluid as a weakly compressible medium. Based on the previously proposed mathematical model of convection of a weakly compressible fluid, unsteady convective motion in a vertical band, with a heat flux periodic in time set on the solid boundaries of this band, is considered. This model of convection allows one to study the problem with the boundary thermal model oscillating in an antiphase rather than in-phase mode, while the latter was required for the model of microconvection of an isothermally incompressible fluid. Exact solutions for velocity components and temperature are derived, and the trajectories of fluid particles are constructed. For comparison, the trajectories predicted by the classical Oberbeck-Boussinesq model of convection and by the microconvection model are presented.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 52–63, March–April, 2005.     

Differential model of turbulence: A numerical study of mixed convection in vertical pipes

A three-parameter model of turbulence with the transport equation for the transverse turbulent heat flux is complemented by terms taking thermogravity effects into account. The results of a numerical study carried out without using the Boussinesq approximation are compared with known experimental data on rising air flows in vertical heated pipes.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 73–86, March–April, 1996.     

Analytical and numerical study of buoyancy-driven convection in a vertical enclosure filled with nanofluids

This paper presents an analytical and numerical study of natural convection of nanofluids contained in a rectangular enclosure subject to uniform heat flux along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number Ra, the Prandtl number Pr, the aspect ratio of the cavity A and the solid volume fraction of nanoparticles, Φ. Three types of nanoparticles are taken into consideration: Cu, Al₂O₃ and TiO₂. Various models are used for calculating the effective viscosity and thermal conductivity of nanofluids. In the first part of the analytical study, a scale analysis is made for the boundary layer regime situation. In the second part, an analytical solution based on the parallel flow approximation is reported for tall enclosures (A ≫ 1). In the boundary layer regime a good agreement is obtained between the predictions of the scale analysis and those of the analytical solution. Solutions for the flow fields, temperature distributions and Nusselt numbers are obtained explicitly in terms of the governing parameters of the problem. A numerical study of the same phenomenon, obtained by solving the complete system of the governing equations, is also conducted. A good agreement is found between the analytical predictions and the numerical simulations.     

Natural convection in vertical cavities with internal heat generating porous medium

Steady natural convection heat transfer in a two dimensional cavity filled with a uniform heat generating, saturated porous medium has been studied. The boundary conditions were: Two isothermal walls at different temperatures, two horizontal adiabatic walls. The aspect ratio was varied from 0.1 to 10 and the Rayleigh number from 10⁰ to 10⁸. The results are presented in terms of the isotherms and stream functions, the temperature variation and maximum temperature in the cavity and heat transfer from the vertical walls. The study indicates that asymmetric vertical boundary conditions with _h >0 has an important effect on the temperature and flow fields as well as on the heat transfer characteristics of the cavity with highly asymmetric results. Various heat transfer modes are identified dependent on the Rayleigh number and the aspect ratio.Es wurde stetiger Wärmeübergang durch freie Konvektion in einen zweidimensionalen Hohlraum, gefüllt mit gleichmäßig wärmeerzeugendem porösen Medium, untersucht. Dabei wurden folgende Randbedingungen festgelegt: zwei isotherme Wände unterschiedlicher Temperatur, zwei horizontale adiabate Wände. Das Verhältnis vonH/L wurde zwischen 0,1 und 10, die Rayleigh-Zahl zwischen 10⁰ und 10⁸ variiert. Die Ergebnisse werden durch Terme der Temperaturverläufe, der maximalen Temperatur im Hohlraum, der Isothermen und der Strömungsfunktion, sowie der Wärmeübertragung der vertikalen Wände dargestellt. Die Studie zeigt, daß asymmetrische Randbedingungen mit _h >0 einen großen Einfluß auf Temperatur- und Strömungsfelder, als auch auf die Wärmeübertragungskennzahlen des Hohlraumes mit stark asymmetrischen Ergebnissen haben. Die verschiedenen Arten der Wärmeübertragung werden in Abhängigkeit von der Rayleigh-Zahl und dem Verhältnis vonH/L beschrieben.     

Experimental study on the magnetic convection in a vertical cylinder

The flow of paramagnetic fluid inside a cylinder placed in a bore of a superconducting magnet was studied experimentally. Single-phase closed thermosyphon configuration was employed. The lower side wall of the cylindrical enclosure was heated while the upper side wall was cooled with the thin adiabatic interface in between. The experiment was carried out with an aqueous solution of glycerol. The magnetic susceptibility of the working fluid was increased by adding Gd(NO₃)₃ × 6H₂O and was measured by a magnetic susceptibility balance. The encapsulated liquid crystal slurry (KWN-2025, Japan Capsular Product Inc.) dispersed in the working fluid was illuminated in the middle height horizontal cross-section of the enclosure to visualize the temperature field. The color images of flow mode were taken by a digital camera. The average heat transfer rates were also measured. Depending on the Rayleigh number, different spoke patterns were observed. The number of angular structures (spokes) increased with increase not only in the Rayleigh number but also in the strength of magnetic field. The heated fluid was repelled by the magnetic field, while the cooled fluid was attracted. The magnetic field enhanced the heat transfer rate.     

Performance of various RANS eddy-viscosity models for turbulent natural convection in tall vertical cavities

The present study is dedicated to the identification of turbulence models that are accurate and numerically economic for computing the natural air-flow and heat transfer by convection in tall cavities with differentially heated vertical walls. The eddy-viscosity models (EVM) are among the simplest to implement and the most economical to treat this problem. This study evaluated the dynamic, thermal and computational performances of twenty EVM turbulence models with one, two or three-equation closure. All the models were first implemented in several in-house codes using the finite volume method. The predictions of the retained models in terms of profiles of velocity, temperature and vertical velocity fluctuations in the cavity have been compared with those of experimental or numerical studies. The obtained results were used to identify the turbulence models that are accurate and numerically economic in predicting natural convection in vertical cavities with a high aspect ratio. The EVM models with three-equation (v2-f and ζ-f) provide the most accurate mean and fluctuating quantities, followed by the k-ε RNG (ReNormalization Group) and k-ω SST (Shear Stress Transport) models. The computing time of these four models is higher than that of the 2L (two-layer) and q-ω models, which provide fairly accurate results especially for the mean heat transfer between the vertical active walls. The other one-equation (Spalart and Allmaras model) and two-equation (k-ε, k-ω and hybrid models) turbulence models tested in this work, have a high computing time and/or predictions that are not sufficiently precise simultaneously for both velocity and temperature fields.     

Parameter-free numerical method for modeling thermal convection in square cavities in a wide range of Rayleigh numbers

Some numerical results for the two- and three-dimensional de Vahl Davis benchmark are presented. This benchmark describes thermal convection in a square (cubic) cavity with vertical heated walls in a wide range of Rayleigh numbers (10⁴ to 10¹⁴), which covers both laminar and highly turbulent f lows. Turbulent f lows are usually described using a turbulence model with parameters that depend on the Rayleigh number and require adjustment. An alternative is Direct Numerical Simulation (DNS) methods, but they demand extremely large computational grids. Recently, there has been an increasing interest in DNS methods with an incomplete resolution, which, in some cases, are able to provide acceptable results without resolving Kolmogorov scales. On the basis of this approach, the so-called parameter-free computational techniques have been developed. These methods cover a wide range of Rayleigh numbers and allow computing various integral properties of heat transport on relatively coarse computational grids. In this paper, a new numerical method based on the CABARET scheme is proposed for solving the Navier–Stokes equations in the Boussinesq approximation. This technique does not involve a turbulence model or any tuning parameters and has a second-order approximation scheme in time and space on uniform and nonuniform grids with a minimal computational stencil. Testing the technique on the de Vahl Davis benchmark and a sequence of refined grids shows that the method yields integral heat f luxes with a high degree of accuracy for both laminar and highly turbulent f lows. For Rayleigh numbers up to 10¹⁴, a several percent accuracy is achieved on an extremely coarse grid consisting of 20 × 20 cells refined toward the boundary. No definite or comprehensive explanation of this computational phenomenon has been given. Cautious optimism is expressed regarding the perspectives of using the new method for thermal convection computations at low Prandtl numbers typical of liquid metals.     

Experimental and numerical analysis of thermo-magnetic convection in a vertical annular enclosure

In the present paper an experimental and numerical analysis of a thermo-magnetic convective flow of paramagnetic fluid in an annular enclosure with a round rod core and a cylindrical outer wall is presented. It is complemented by an experimental analysis of natural convection depending on the inclination angle to show the stability of the present configuration. Convection in an annulus between two vertical co-axial cylinders resulting from gravitational and magnetic environments has been investigated. A strong magnetic field can be an alternative to heat transfer enhancement. The effect of the magnetic field on the convection of the paramagnetic fluid in the annular vessel in various positions was compared. The numerical analysis was done based on the continuity, momentum and energy equations. A term related to the magnetic buoyancy force was added to the momentum equation. The distributions of Nusselt number present minima in two positions of the enclosure, which depends on the reciprocal relationship between the gravitational and magnetising forces.     

Thin axisymmetric cavities in a subsonic compressible flow

The investigation of subsonic compressible flow past thin axisymmetric cavities carried out in [1–3] is continued by the method of asymptotic expansions. The dependence of the elongation of the cavity on the cavitation number and the Mach number is obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 174–177, September–October, 1987.The author is grateful to Yu. L. Yakimov for discussing his results.     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

A numerical study of three-dimensional laminar natural convection in a vented enclosure

A three-dimensional numerical investigation of steady laminar natural convection in vented enclosures is carried out. A discrete flush-type heat source mounted on the substrate is used to simulate an electronic component. Four different vent locations are investigated. Combined natural convection in the air and conduction in the heat source, the substrate, and the enclosure walls are solved. Solutions are obtained for Rayleigh numbers ranging from 10⁴ to 10⁶, different substrate thermal conductivity ratios, and varied vent sizes. The calculation domain is extended beyond the cubic enclosure in x-, y-, and z-directions. Appropriate boundary conditions are prescribed on the extended computational domain. The resulting flow and temperature patterns are discussed. Also, the local and overall heat transfer from the heat source and the substrate, in terms of Nusselt numbers and the surface temperatures, are presented to illustrate the vent effects.     

模拟微可压粘弹性流体的WCCBS_SU方法

针对微可压缩粘弹性流动问题,发展了微可压缩流的WCCBS方法,详细推导了基于Oldroyd-B本构模型的WCCBS_SU方法的求解过程。在流场微可压的条件下,分别对平面Poiseuille流和4：1粘弹性收缩流进行了数值模拟。Poiseuille流在不同We数下数值结果与解析解的比较,验证了本文方法具有较高的精度和较好...     

A numerical study for WENO scheme-based on different lattice Boltzmann flux solver for compressible flows

In this paper, the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows, including inviscid LBFS I, viscous LBFS II, hybrid LBFS III and hybrid LBFS IV. Hybrid LBFS can automatically realize the switch between inviscid LBFS I and viscous LBFS II through introducing a switch function. The resultant hybrid WENO–LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS. We investigate the performance of WENO scheme based on four kinds of LBFS systematically. Numerical results indicate that the devopled hybrid WENO–LBFS scheme has high accuracy, high resolution and no oscillations. It can not only accurately calculate smooth solutions, but also can effectively capture contact discontinuities and strong shock waves.     

Effect of Coriolis force on the thermosolutal convection threshold in a rotating annular Hele-Shaw cell

In this study, the linear stability analysis is used to determine the onset of thermosolutal convection in fluids confined in rotating annular Hele-Shaw cell. The fluid layer is submitted to radial gradients of temperature and concentration. The effects of both Coriolis force and curvature parameter on the stationary and oscillatory convection are investigated when the Prandtl number is of the order of unity or larger than unity.     

Numerical study of natural convection flow in a vertical slot

The stability of convective motion, generated by a lateral temperature difference across a vertical slot, is studied numerically over a range ofGr=5000 to 1.5 × 10⁵,Pr=0.01 to 10, andA=8,16 and 20. Various cellular flow structures and temperature patterns are illustrated. Several branches of solutions characterized by different numbers of the cells in the flow patterns as well as by both steady and unsteady multicellular patterns are found for low-Prandtl-number fluid in the vertical slot. Meanwhile, the behaviors of the temperature variation and heat transfer are also discussed. The project supported by the National Natural Science Foundation of China (59776011) and by the Returnee from Abroad Funding of Academia Sinica.     

A numerical technique for three-dimensional compressible boundary layers

The non-linear two-point boundary value problem for three-dimensional compressible boundary layers is solved through the application of a boundary value technique for a range of parameters characterizing the nature of stagnation point flows. The analytical boundary conditions, at infinity, are applied at the edge of the computational mesh with iterations on the size of the domain. The solutions obtained show excellent agreement with the established similarity solutions for three-dimensional flows. The present method has the potential advantage of yielding the wall values of f″_w, g″_w and θ′_w as a part of the solution, contrary to the previously used ‘shooting’ methods. The algorithm is computationally simple and numerically stable and extremely suitable for engineering design applications.