A stable and efficient hybrid method for aeroacoustic sound generation and propagation

We discuss how to combine the node based unstructured finite volume method widely used to handle complex geometries and nonlinear phenomena with very efficient high order finite difference methods suitable for wave propagation dominated problems. This fully coupled numerical procedure reflects the coupled character of the sound generation and propagation problem. The coupling procedure is based on energy estimates and stability can be guaranteed. Numerical experiments using finite difference methods that shed light on the theoretical results are performed. To cite this article: J. Nordström, J. Gong, C. R. Mecanique 333 (2005).     

Simulation du front d'écoulement dans les procédés de moulage des composites liquides

This work proposes a procedure of finite difference discretized in a system of curvilinear coordinates adapted to the shape of the saturated zone, to simulate a flow in the LCM (Liquid Composites Molding) process. Formulation and the numerical application of the procedure are described. We describe two configurations of injections. A good agreement is found between numerical, analytical and literature experimental results. To cite this article: M. Hattabi et al., C. R. Mecanique 333 (2005).     

A coupled solution of the vorticity–velocity formulation of the incompressible Navier–Stokes equations

A relatively novel formulation of the Navier-Stokes equations is used for obtaining solutions of two dimensional incompressible fluid flow and convective heat transfer problems. A vorticity transport equation along with two Poisson equations for the velocity components and the energy equation are solved by a finite difference scheme. A coupled solution procedure is used for solving simultaneously the dependent variables along a line, using a block tridiagonal matrix algorithm. The formulation is found to be stable and has features that may be desirable for solving a wide variety of flow and heat transfer problems.     

Study of stabilization methods for computational aeroacoustics

The use of high-order centered finite difference to solve the Euler equations commonly requires a stabilization procedure. The present work is a theoretical analysis of these stabilization methods that make the whole algorithm (i) still consistent with the continuous problem and (ii) able to run long time simulations. In the present study, a theoretical analysis of the three commonly used methods resorting to the application of high-order filters is performed. An extension to non-periodic boundary conditions is studied to avoid numerical reflection and numerical instabilities due to the use of specific boundary schemes. To cite this article: R. Guénanff, M. Terracol, C. R. Mecanique 333 (2005).     

A novel non‐upwind,interconnected, multi‐grid,overlapping numerical procedure for problems involving fluid flow

A novel numerical procedure for heat, mass and momentum transfer in fluid flow is presented. The new scheme is passed on a non‐upwind, interconnected, multi‐grid, overlapping (NIMO) finite‐difference algorithm. In 2D flows, the NIMO algorithm solves finite‐difference equations for each dependent variable on four overlapping grids. The finite‐difference equations are formulated using the control‐volume approach, such that no interpolations are needed for computing the convective fluxes. For a particular dependent variable, four fields of values are produced. The NIMO numerical procedure is tested against the exact solution of two test problems. The first test problem is an oblique laminar 2D flow with a double step abrupt change in a passive scalar variable for infinite Peclet number. The second test problem is a rotating radial flow in an annular sector with a single step abrupt change in a passive scalar variable for infinite Peclet number. The NIMO scheme produced essentially the exact solution using different uniform and non‐uniform square and rectangular grids for 45 and 30° angle of inclination. All other schemes were unable to capture the exact solution, especially for the rectangular and non‐uniform grids. The NIMO scheme was also successful in predicting the exact solution for the rotating radial flow, using a uniform cylindrical‐polar coordinate grid. Copyright © 2008 John Wiley & Sons, Ltd.     

A numerical model for thermal systems with helically-coiled tubes and its experimental verification

 A conjugate numerical model proposed by Nakayama et al. for the steady problem of cooling a fluid flowing through a coiled tube, has been successfully extended to investigate two distinctive thermal problems, namely, the transient cooling processes associated with a beer dispenser, and the transient processes of heat storage and recovery associated with a packed bed saturated with a molten salt. An axisymmetric numerical procedure is adopted for determining the velocity and temperature fields within the chilled water bath of the beer dispenser. A simplified one-dimensional heat transfer model is introduced for coupling the tube flow with the recirculating flow in the bath. A similar axisymmetric finite difference procedure is applied for the heat transfer analysis of the packed bed saturated with a molten salt. For the heat recovery process, a one-dimensional heat balance equation for the two-phase flow with a helically-coiled tube is introduced to update the wall surface temperatures, which are needed to calculate the temperature field in the saturated packed bed. The numerical results for both thermal systems associated with coiled tubes agree very well with the corresponding velocity and temperature data obtained from the experiments. Received on 28 August 2000 / Published online: 29 November 2001     

A computational method for geometric optimization of enhanced heat transfer devices based upon entropy generation minimization

A computational fluid dynamics‐based optimization methodology is developed, appropriate for the geometric optimization of enhanced heat transfer devices based upon the principle of entropy generation minimization, in which the objective function is evaluated from a flow field obtained by computational simulation. A quasi‐Newton optimization procedure is employed, with computation of the objective function gradients based upon a finite difference approach. The optimization procedure is developed to be general with regard to the choice of objective function, the details of the problem under consideration, and the computational methodology employed in solving the fluid flow and heat transfer problems. A novel implementation of a Taylor series‐based procedure for the fast solution of nearby problems is presented, which is found to greatly benefit the efficiency of the present methodology. Finally, a numerical experiment is presented, illustrating the use of the present method in the geometric optimization of a practical enhanced heat transfer device on the basis of the criterion of entropy generation minimization. The optimization of the fin spacing of a simple plate fin heat sink is considered, and a comparison of the computational results with results obtained by analytical optimization based upon empirical friction factor and Nusselt number correlations is given. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional finite element model for metal displacement and heat transfer in squeeze casting processes

A three-dimensional finite element model for the numerical simulation of metal displacement and heat transfer in the squeeze casting process has been developed. In the model, a numerical approach, termed as ‘Quasi-static Eulerian’, is proposed, in which the dynamic metal displacement process is divided into a certain number of sub-cycles. In each of the sub-cycles, the dieset configuration is assumed to be static and a fixed finite element mesh is created, thus making the Eulerian approach applicable to the solution of metal flow and heat transfer. Mesh-to-mesh data mapping is carried out for any two adjacent sub-cycles to ensure that the physical continuity of the real metal displacement process is represented. A numerical example is presented, which shows the application of the present model to geometrically complex three-dimensional squeeze casting problems. To cite this article: R.W. Lewis et al., C. R. Mecanique 335 (2007).     

A model of a subcritical Joule–Thomson cryocooler with condensation inside the recuperator

To develop a tool for predicting of heat and mass transfer in Joule–Thomson cryocoolers working at subcritical pressures, we study a counter flow heat exchanger with condensation by employing the integral method. The effects of inlet pressure and working fluid are predicted. We also show that there is an optimal value of the enthalpy difference along the heat exchanger for which its length is minimal.

---

On the stability of a Bingham fluid flow in an annular channelSur la stabilité de l'écoulement d'un fluide de Bingham dans une conduite annulaire

Linear stability of a fully developed Bingham fluid flow between two coaxial cylinders subject to infinitesimal axisymetric perturbations is investigated. The analysis leads to two uncoupled Orr–Sommerfeld equations with appropriate boundary conditions. The numerical solution is obtained using fourth order finite difference scheme. The computations were performed for various plug flow dimensions and radii ratios. Within the range of the parameters considered in this paper, the Poiseuille flow of Bingham fluid is found to be linearly stable. To cite this article: N. Kabouya, C. Nouar, C. R. Mecanique 331 (2003).     

Entropy generation during fluid flow in a channel under the effect of transverse magnetic field

Entropy generation due to fluid flow and heat transfer inside a horizontal channel made of two parallel plates under the effect of transverse magnetic field is numerically investigated. The flow is assumed to be steady, laminar, hydro-dynamically and thermally fully developed of electrically conducting fluid. Both horizontal walls are maintained at constant temperatures higher than that of the fluid. The governing equations in Cartesian coordinate are solved by an implicit finite difference technique. After the flow field and the temperature distributions are obtained, the entropy generation profiles are computed and presented graphically. The factors, which were found to affect the problem under consideration are the magnetic parameter, Eckert number, Prandtl number, and the temperature parameter (θ_∞). It was found that, entropy generation increased as all parameters involved in the present problem increased.     

Velocity–pressure coupling in finite difference formulations for the Navier–Stokes equations

A new numerical procedure for solving the two‐dimensional, steady, incompressible, viscous flow equations on a staggered Cartesian grid is presented in this paper. The proposed methodology is finite difference based, but essentially takes advantage of the best features of two well‐established numerical formulations, the finite difference and finite volume methods. Some weaknesses of the finite difference approach are removed by exploiting the strengths of the finite volume method. In particular, the issue of velocity–pressure coupling is dealt with in the proposed finite difference formulation by developing a pressure correction equation using the SIMPLE approach commonly used in finite volume formulations. However, since this is purely a finite difference formulation, numerical approximation of fluxes is not required. Results presented in this paper are based on first‐ and second‐order upwind schemes for the convective terms. This new formulation is validated against experimental and other numerical data for well‐known benchmark problems, namely developing laminar flow in a straight duct, flow over a backward‐facing step, and lid‐driven cavity flow. Copyright © 2010 John Wiley & Sons, Ltd.     

Application of a 3-D FEM/FDM Overlapping Scheme to Heat and Mass Transfer and Moving-Boundary Problems

A 3-D FEM/FDM overlapping scheme for viscous, incompressible flow problems is presented that combines the finite element method, which is best suited to analyze flow in any arbitrarily shaped flow geometry, with the finite difference method, which is advantageous in both computing time and computer storage. The combination of both methods enables large-scale viscous flow to be analyzed, which is crucial both for detailed analysis of 3-D flows and for solving flow problems around moving bodies, A modified ABM AC method is used as the basic algorithm, to which a sophisticated time integration scheme, proposed by the present authors, has been applied. In this paper, some numerical results including 3-D heat and mass transfer problem and moving-boundary problems are presented.     

A comparative study of two fast nonlinear free‐surface water wave models

This paper presents a comparison in terms of accuracy and efficiency between two fully nonlinear potential flow solvers for the solution of gravity wave propagation. One model is based on the high‐order spectral (HOS) method, whereas the second model is the high‐order finite difference model OceanWave3D. Although both models solve the nonlinear potential flow problem, they make use of two different approaches. The HOS model uses a modal expansion in the vertical direction to collapse the numerical solution to the two‐dimensional horizontal plane. On the other hand, the finite difference model simply directly solves the three‐dimensional problem. Both models have been well validated on standard test cases and shown to exhibit attractive convergence properties and an optimal scaling of the computational effort with increasing problem size. These two models are compared for solution of a typical problem: propagation of highly nonlinear periodic waves on a finite constant‐depth domain. The HOS model is found to be more efficient than OceanWave3D with a difference dependent on the level of accuracy needed as well as the wave steepness. Also, the higher the order of the finite difference schemes used in OceanWave3D, the closer the results come to the HOS model. Copyright © 2011 John Wiley & Sons, Ltd.     

A heat transfer measurement of jet impingement with high injection temperature

Local heat transfer from an impinging high temperature jet is studied using a method based on the heat thin foil technique and on the infrared thermography. Heat thin foil technique is used to impose several heat fluxes. For each flux, the temperature distribution is recorded using infrared imaging. Local heat transfer coefficients and adiabatic wall temperatures are determined by means of a linear regression method. This procedure is validated for a single round jet impinging on a flat plate for a range of injection temperatures. To cite this article: M. Fénot et al., C. R. Mecanique 333 (2005).     

Influence of longitudinal conduction in the matrix on effectiveness of rotary heat regenerator used in air-conditioning

The paper presents an efficient analytical method of solving the problem of a rotary heat regenerator taking into account longitudinal heat conduction in the matrix. The small parameter method, Laplace transform as well as one of the spline functions have been applied for approximation of an initial condition in the reversion time. In the application part, the solution for a model in analysis of an influence of longitudinal conduction in the matrix on effectiveness of rotary heat regenerator in a wide range of dimensionless parameters as well as for the particular matrix applied in air-conditioning was used.

---

An Immersed Boundary Method for Complex Flow and Heat Transfer

The need to predict flow and heat transfer problems requires a flexible and fast tool able to simulate complex geometries without increasing the complexity of the flow solver architecture. Here we use a finite volume code that uses a direct solver with pressure correction. A new immersed boundary method (IBM) is used for a geometry consisting of a square body in a flow. The method is applied to flow cases with and without heat transfer. The obstacle simulated in the domain is implemented by local forcing of the flow with a procedure that adjusts locally the shear stress at the position of the object in conjunction with a non-penetration condition on the body walls. This approach has already been successfully applied by Breugem and Boersma (Phys. Fluids 17:15, 2005). We extend it for the case of heat transfer between body and flow. Comparison with other methods has been carried out as well. However, the proposed method can not be simply extended to immersed boundaries not aligned with the grid.     

On the effectiveness and limitations of local criteria for the identification of a vortex

Three proposed local criteria for the identification of vortices are analyzed and discussed; they are based on the analysis of invariants of the velocity gradient tensor ▿u or invariants of the tensor S² +Ω², where S and Ω are the symmetric and antisymmetric parts of ▿u.Moreover, a tentative non-local procedure is proposed, which takes advantage of the observation that vortices tend to be made up of the same fluid particles; this leads to the definition of a Galilean invariant quantity, which can be computed and used to identify vortical structures.Three analytical flow fields are used for a comparative evaluation of both local and non-local criteria, which allows a deeper understanding of the physical meaning of the considered techniques.     

Modelling of the free surface with heat and mass transfer considerations on a liquid filling a porous slab heated from its sides

The position of the free surface is calculated numerically for a porous slab which is partly filled with a liquid and differentially heated from its sides. A coordinate transformation is used to transform the original problem from a physical coordinate system to a non‐orthogonal system where the free surface becomes a fixed straightline. The transformed problem is then solved using a finite difference method. Results are obtained for Rayleigh numbers up to 1000. The Nusselt numbers increase slightly with medium Rayleigh numbers (convection‐dominated region) as expected. Since at low Ra the conduction is dominant and at high Ra radiation is dominant. Hadizadeh and Tien (Int. J. Heat Mass Transfer 2004; 17 (6):799–804) studied the forced convection on the surface of porous layer. In that paper they dealt with in detail the boundary regime of liquid in the channel and modelled the flow and heat transfer. Copyright © 2007 John Wiley & Sons, Ltd.     

Viscous dissipative fluid flow past a semi‐infinite vertical plate with variable surface temperature

An analysis of the effect of viscous dissipative heat on two‐dimensional viscous incompressible fluid flow past a semi‐infinite vertical plate with variable surface temperature is carried out. The dimensionless governing equations are unsteady, two‐dimensional, coupled, and non‐linear governing equations. A most accurate, unconditionally stable and fast converging implicit finite‐difference scheme is used to solve the non‐dimensional governing equations. Velocity and temperature of the flow have been presented graphically for various parameters occurring in the problem. The local and average skin friction and Nusselt number are also shown graphically. It is observed that greater viscous dissipative heat causes a rise in the temperature. Copyright © 2007 John Wiley & Sons, Ltd.