scholarly journals Direct Numerical Simulation of Liquid Flow in a Horizontal Microchannel

2005 ◽  
Author(s):  
Cenk Evren Ku¨krer ◽  
I˙lker Tarı

Numerical Simulations of liquid flow in a micro-channel between two horizontal plates are performed. The channel is infinite in streamwise and spanwise directions and its height is taken as 3.1×10−4 m which falls within the dimension ranges of micro-channels. The Navier-Stokes equations with the addition of Brinkman number (Br) to the energy equation are used as the governing equations and a spectral methods based approach is applied to obtain the required accuracy to handle liquid flow in the micro-channel. It is known for micro-channels that Br combines the effects of conduction and viscous dissipation in liquids and is also a way of comparing the importance of later relative to former. A laminar flow of a liquid in a micro-channel shows different characteristics compared to a similar flow in a macro-channel. To observe the differences, three different cases are run over each of a range of Reynolds numbers: one with no axial conduction assumption that correspond to a case similar to macro-channel flow, another case with axial conduction included in the energy equation to simulate one of the main differences and lastly a case with inclusion of Br number in the governing equations. The results are compared with each other to see the effects of axial conduction and Br inclusion. A qualitative comparison is made with the previous results in literature.

1997 ◽  
Vol 52 (4) ◽  
pp. 358-368 ◽  
Author(s):  
Michio Nishida ◽  
Masashi Matsumotob

Abstract • This paper describes a computational study of the thermal and chemical nonequilibrium occuring in a rapidly expanding flow of high-temperature air transported as a free jet from an orifice into low-density stationary air. Translational, rotational, vibrational and electron temperatures are treated separately, and in particular the vibrational temperatures are individually treated; a multi-vibrational temperature model is adopted. The governing equations are axisymmetric Navier-Stokes equations coupled with species vibrational energy, electron energy and species mass conservation equations. These equations are numerically solved, using the second order upwind TVD scheme of the Harten-Yee type. The calculations were carried out for two different orifice temperatures and also two different orifice diameters to investigate the effects of such parameters on the structure of a nonequilibrium free jet.


Author(s):  
Wolfgang Höhn

During the design of the compressor and turbine stages of today’s aeroengines, aerodynamically induced vibrations become increasingly important since higher blade load and better efficiency are desired. In this paper the development of a method based on the unsteady, compressible Navier-Stokes equations in two dimensions is described in order to study the physics of flutter for unsteady viscous flow around cascaded vibrating blades at stall. The governing equations are solved by a finite difference technique in boundary fitted coordinates. The numerical scheme uses the Advection Upstream Splitting Method to discretize the convective terms and central differences discretizing the viscous terms of the fully non-linear Navier-Stokes equations on a moving H-type mesh. The unsteady governing equations are explicitly and implicitly marched in time in a time-accurate way using a four stage Runge-Kutta scheme on a parallel computer or an implicit scheme of the Beam-Warming type on a single processor. Turbulence is modelled using the Baldwin-Lomax turbulence model. The blade flutter phenomenon is simulated by imposing a harmonic motion on the blade, which consists of harmonic body translation in two directions and a rotation, allowing an interblade phase angle between neighboring blades. Non-reflecting boundary conditions are used for the unsteady analysis at inlet and outlet of the computational domain. The computations are performed on multiple blade passages in order to account for nonlinear effects. A subsonic massively stalled unsteady flow case in a compressor cascade is studied. The results, compared with experiments and the predictions of other researchers, show reasonable agreement for inviscid and viscous flow cases for the investigated flow situations with respect to the Steady and unsteady pressure distribution on the blade in separated flow areas as well as the aeroelastic damping. The results show the applicability of the scheme for stalled flow around cascaded blades. As expected the viscous and inviscid computations show different results in regions where viscous effects are important, i.e. in separated flow areas. In particular, different predictions for inviscid and viscous flow for the aerodynamic damping for the investigated flow cases are found.


It is shown that the boundary layer approximation to the flow of a viscous fluid past a flat plate of length l , generally valid near the plate when the Reynolds number Re is large, fails within a distance O( lRe -3/4 ) of the trailing edge. The appropriate governing equations in this neighbourhood are the full Navier- Stokes equations. On the basis of Imai (1966) these equations are linearized with respect to a uniform shear and are then completely solved by means of a Wiener-Hopf integral equation. The solution so obtained joins smoothly on to that of the boundary layer for a flat plate upstream of the trailing edge and for a wake downstream of the trailing edge. The contribution to the drag coefficient is found to be O ( Re -3/4 ) and the multiplicative constant is explicitly worked out for the linearized equations.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Filipe S. Pereira ◽  
Luís Eça ◽  
Guilherme Vaz

The importance of the turbulence closure to the modeling accuracy of the partially-averaged Navier–Stokes equations (PANS) is investigated in prediction of the flow around a circular cylinder at Reynolds number of 3900. A series of PANS calculations at various degrees of physical resolution is conducted using three Reynolds-averaged Navier–Stokes equations (RANS)-based closures: the standard, shear-stress transport (SST), and turbulent/nonturbulent (TNT) k–ω models. The latter is proposed in this work. The results illustrate the dependence of PANS on the closure. At coarse physical resolutions, a narrower range of scales is resolved so that the influence of the closure on the simulations accuracy increases significantly. Among all closures, PANS–TNT achieves the lowest comparison errors. The reduced sensitivity of this closure to freestream turbulence quantities and the absence of auxiliary functions from its governing equations are certainly contributing to this result. It is demonstrated that the use of partial turbulence quantities in such auxiliary functions calibrated for total turbulent (RANS) quantities affects their behavior. On the other hand, the successive increase of physical resolution reduces the relevance of the closure, causing the convergence of the three models toward the same solution. This outcome is achieved once the physical resolution and closure guarantee the precise replication of the spatial development of the key coherent structures of the flow.


Symmetry ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 580 ◽  
Author(s):  
Zahra Abdelmalek ◽  
Mohammad Yaghoub Abdollahzadeh Jamalabadi

Micromixers are significant segments inside miniaturized scale biomedical frameworks. Numerical investigation of the effects of galloping cylinder characteristics inside a microchannel Newtonian, incompressible fluid in nonstationary condition is performed. Governing equations of the system include the continuity equation, and Navier–Stokes equations are solved within a moving mesh domain. The symmetry of laminar entering the channel is broken by the self-sustained motion of the cylinder. A parameter study on the amplitude and frequency of passive moving cylinder on the mixing of tiny particles in the fluid is performed. The results show a significant increase to the index of mixing uses of the galloping body in biomedical frameworks in the course of micro-electromechanical systems (MEMS) devices.


1998 ◽  
Vol 120 (1) ◽  
pp. 72-75 ◽  
Author(s):  
V. N. Kurdyumov ◽  
E. Ferna´ndez

A correlation formula, Nu = W0(Re)Pr1/3 + W1(Re), that is valid in a wide range of Reynolds and Prandtl numbers has been developed based on the asymptotic expansion for Pr → ∞ for the forced heat convection from a circular cylinder. For large Prandtl numbers, the boundary layer theory for the energy equation is applied and compared with the numerical solutions of the full Navier Stokes equations for the flow field and energy equation. It is shown that the two-terms asymptotic approximation can be used to calculate the Nusselt number even for Prandtl numbers of order unity to a high degree of accuracy. The formulas for coefficients W0 and W1, are provided.


1965 ◽  
Vol 87 (2) ◽  
pp. 252-258 ◽  
Author(s):  
Y. Peng ◽  
S. W. Yuan

The effect of foreign coolant injection at the wall on the temperature distribution of a laminar flow of a fluid with variable transport and thermodynamic properties in a porous-wall pipe has been investigated. The velocity components, mass concentration, and temperature distribution were obtained by the solution of the Navier-Stokes equations, the diffusion equation, and the energy equation. A perturbation method was used to solve the first equations for small flows through the porous wall, and the eigenvalues in the latter two equations were calculated with the aid of the CDC 1604 computer. The results from this investigation depict the significant differences in both velocity distribution and temperature distribution between the present case of hydrogen coolant and the case of air coolant [1]. The results also show that the heat transfer coefficient at the wall in the present case is considerably smaller than the case of air-coolant injection.


1989 ◽  
Vol 56 (1) ◽  
pp. 47-50 ◽  
Author(s):  
C. Y. Wang

Melting of a disk is facilitated by rotation. The problem is governed by a nondimensional parameter α which represents the relative importance of injection (melt) rate and rotation times viscosity. The nonlinear governing equations are solved by perturbations for small α and numerical integration for arbitrary α. Torque and heat transfer rates are found. The solution is one of the rare exact similarity solutions of the Navier-Stokes equations.


2009 ◽  
Vol 19 (09) ◽  
pp. 2939-2951 ◽  
Author(s):  
JACQUES HONA ◽  
ELISABETH NGO NYOBE ◽  
ELKANA PEMHA

In this paper, the axisymmetric flow of a viscous fluid through a porous annular tube with walls kept at different temperatures is studied theoretically. The physical properties of the fluid remain constant, notably its specific mass, its dynamic viscosity and its thermal diffusivity. The nondimensional parameters which the solutions of the problem depend on are defined. A numerical integration using the shooting method is applied for solving the Navier–Stokes equations and the energy equation. Bifurcation diagrams are presented and enable to highlight significant properties of the flow. Some thermal behaviors corresponding to specific values of the parameters are performed. Asymmetric solutions of the steady flow are described and some results about velocity components are also analyzed.


Sign in / Sign up

Export Citation Format

Share Document