Comparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel
|
|
- Samantha Cox
- 6 years ago
- Views:
Transcription
1 Comparison of some approximation schemes for convective terms for solving gas flow past a square in a micorchannel Kiril S. Shterev and Sofiya Ivanovska Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Block 4, Sofia 1113, Bulgaria Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, Acad. G. Bonchev St., Bl. 25A, Sofia 1113, Bulgaria Abstract. Rapidly emerging micro-electro-mechanical devices create new potential microfluidic applications. A simulation of an internal and external gas flows is important for their design. For small Knudsen number Kn < 0.1 (Kn = l 0 /L, where l 0 is the mean free path of the gas molecules and L is the characteristic length), a continuum approach based on modified Navier-Stokes-Fourier or extended hydrodynamic continuum models with corresponding velocity-slip and temperature-jump boundary conditions is still applicable and, respectively, preferable. We restrict ourself to the use of Navier-Stokes-Fourier continuum model. A development of the algorithm to solve a specific class of problems is closely related to numerical schemes used for approximation of equations terms. Higher-order approximation schemes can reduce the number of mesh nodes and respectively computational time, but it is possible to obtain physical unrealistic results. In this paper we study influence of some approximation schemes for convective terms over the spatial steps. It is compared upwind, central difference and total variation diminishing (TVD) schemes Min-Mod, QUICK and SUPERBEE. A test case is gas flow past a square in a microchannel at subsonic speed (Mach number M = 0.1) and supersonic speed (M = 2.43), available in a literature. Keywords: mesh convergence, TVD schemes, subsonic flow, supersonic flow, SIMPLE-TS, micro-gas flows, continuum model PACS: Mv, j, Df, x, n, Fg, Dx INTRODUCTION All devices with character dimensions between 1µm and 1mm are called micro-devices. Micro mechanical devices are rapidly emerging technologies, where new potential applications are continuously being developed. Microchannel fluid flow is currently modelling using either the continuum approach or the non-continuum approach most commonly called molecular approach. The Knudsen number (Kn), a nondimensional parameter, determines the degree of appropriateness of the continuum model. It is defined as the ratio of the mean free path l 0 to the macroscopic length scale of a physical system L. Kn = l 0 (1) L Micro-gas flow simulations should be fast and accurate. Higher order spatial accuracy scheme can reduce required computational resources. Correct calculation of supersonic microflows use first order upwind scheme for approximation of convective terms and density in the middle points requires a very fine mesh and enormous computational resources [1]. On the other hand, higher order schemes reduce the number of nodes in the mesh, but it is possible to be unstable and to obtain physical unrealistic results. TVD schemes are second order spatial accuracy schemes and are designed to prevent undesirable oscillatory behaviour of higher-order schemes. In TVD schemes, the tendency towards oscillation is counteracted by adding an artificial diffusion fragment or by adding a weighting towards upstream contribution. In the literature early schemes based on these ideas were called flux corrected transport (FCT) schemes: see Boris and Book [2, 3]. Further works by Van Leer [4, 5, 6, 7], Harten [8, 9], Sweby [10], Roe [11], Osher and Chakravarthy [12] and many others have contributed to the development of present-day TVD schemes. These schemes are widely used in Computational Fluid Dynamic (CFD) codes and software packages like OpenFOAM [13]. The main target of the present work is to compare mesh-dependent convergence of TVD schemes applied for calculation of gas microflows at different regimes. TVD schemes are implemented in the algorithm SIMPLE-TS [1]. There is comparison between first order upwind scheme and second order TVD schemes: Min-Mod [11], QUICK [14] and SUPERBEE [11]. Flow past a square is used as a test case in a microchannel at different regimes: subsonic (M = 0.1) and supersonic (M = 2.43) speed.
2 CONTINUUM MODEL EQUATIONS A two dimensional system of equations describing the unsteady flow of viscous, compressible, heat conductive fluid can be expressed in a general form as follows: where: (ρu) t (ρv) t (ρt ) t + (ρuu) x + (ρuv) x + (ρut ) x + (ρvu) + (ρvv) + (ρvt ) [ ( u ) 2 Φ = 2 + x ρ t + (ρu) + (ρv) = 0 x (2) =ρg x A p [ ( x + B Γ u ) Γ u )] x x { ( + B Γ u ) Γ v ) 2 [ ( u Γ x x x 3 x x + v )]} (3) =ρg y A p [ ( + B Γ v ) Γ v )] x x { ( + B Γ v ) Γ u ) 2 [ ( u Γ x 3 x + v )]} (4) [ ( =C T 1 Γ λ T ) Γ λ T )] +C T 2.Γ.Φ +C T 3 Dp x x Dt (5) p = ρt (6) ( ) ] v 2 ( v + x + u ) 2 2 ( u 3 x + v ) 2 (7) u is the horizontal component of velocity, v is the vertical component of velocity, p is pressure, T is temperature, ρ is density, t is time, x and y are coordinates of a Cartesian coordinate system. The reference quantities used to scale the system of equations (2) - (6) will be defined later. The parameters A, B, g x, g y, C T 1, C T 2, C T 3 and diffusion coefficients Γ and Γ λ, given in Eqs. (2)-(6), depend on the gas model and the equation non-dimensional form. The system of equations (2) - (6) is solved by using an algorithm SIMPLE-TS [1]. Second order central difference scheme is employed for the approximation of the diffusion terms. The main topic of this paper is comparison of mesh-dependent convergence of schemes for approximation of convective terms and density in middle points. In following sections are presented details of applications and results. The system of equations (2) - (6) is given in a general form of the Navier-Stokes-Fourier equations. For a gaseous microflow description we use the model of a compressible, viscous hard sphere gas with diffusion coefficients determined by the first approximation of the Chapman-Enskog theory for low Knudsen numbers [15]. For a hardsphere gas, the viscosity coefficient µ and the heat conduction coefficient λ read (first approximations are sufficient for our considerations) as: µ = µ h T, µh = (5/16)ρ 0 l 0 V th π (8) λ = λ h T, λh = (15/32)c p ρ 0 l 0 V th π (9) The Prandtl number is given by Pr = 2/3, γ = c p /c v = 5/3. The dimensionless system of equations (2) - (6) is scaled by the following reference quantities, as given in [15]: molecular thermal velocity V 0 = V th = 2RT 0 for velocity, for length - square size a (Fig. 3), for time - t 0 = a/v 0, the reference pressure (p 0 ) is the pressure at the inflow of the channel, the reference temperature (T 0 ) is equal to the channel walls, reference density (ρ 0 ) is calculated using equation of state (6), the calculated case neglects the influence of gravity field, therefore g x = g y = 0. The corresponding non-dimensional parameters in the equation system (2) - (6) read as follows: A = 0.5, B = 5 π 16 Kn, Γ = Γλ = T, C T 1 = 15 π π 32 Kn, CT 2 = 4 Kn, CT 3 = 2 5 (10)
3 TVD SCHEMES - ESSENCE AND IMPLEMENTATION Generalization of upwind-biased discretisation schemes is considered in many papers and books (see [16]). Here is presented briefly the application of TVD scheme in algorithm SIMPLE-TS. Consider the standard control volume discretisation of the one dimensional case, Fig. 1. The implementation for the two dimensional case is straightforward. Approximated variable φ is defined in nodes of a mesh, φ i 1, φ i and φ i+1, Fig. 1. TVD scheme is applied to FIGURE 1. 1D Structured mesh approximate variable in points, where the variable is not defined. The variable is not defined on control surfaces, φi CS and φi+1 CS, Fig. 1. The general form of the value φcs i+1 in a descretisation scheme for u > 0, Fig. 1, can be written as: φ CS i+1 = φ i Ψ(r)(φ i+1 φ i ) (11) where r is the ratio of upwind-side gradient to downwind-side gradient r = (φ i φ i 1 )/(φ i+1 φ i ) determines the value of function Ψ and the nature of the scheme. To prevent division by zero, in code is included condition that r = 0, if (φ i+1 φ i )/ < In the literature ([8, 9, 10]) the total variation has been considered for transient one-dimensional transport equation and are used explicit TVD schemes. TVD schemes are successfully applied also for calculation of steady problems (see [16]). Here are used implicit TVD schemes for calculation of steady microfluidic problems. The algorithm SIMPLE- TS, with the implemented TVD scheme, is also applicable for calculation of unsteady problems. Sweby [10] has given necessary and sufficient conditions for a scheme to be first- (Fig. 2. (a)) and second-order (a) FIGURE 2. Region in grey for a first-order TVD schemes (a) and second-order TVD schemes (Fig. 2. ) TVD in terms of r Ψ relationship. Fig. 2. (a) confirm that upwind is a first order TVD scheme, while the central difference is second-order TVD scheme for r 0.5, but is not TVD for r < 0.5. The quadratic upstream interpolation for convective kinetics (QUICK) scheme [17], Ψ(r) = (3 + r)/4, also leave the TVD region, Fig. 2. (a). The idea of designing a TVD scheme is to force the r Ψ relationship to remain within the shaded region for all values of r. QUICK schemes with limiters (13), [14], is in a TVD region. In this paper are implemented and used TVD schemes Min-Mod (12), [11], QUICK (13) and SUPERBEE (14), [11]. The Min-Mod limiter function exactly traces the lower limit of the TVD region, whereas SUPERBEE scheme follows the upper limit, (Fig. 2. ). Ψ(r) Min Mod = max(0,min(r,1)) (12) Ψ(r) QUICK = max(0,min(2r,(3 + r)/4,2)) (13) Ψ(r) SUPERBEE = max(0,min(2r,1),min(r,2)) (14) TVD schemes are implemented for approximation of convective terms of equations (3), (4), (5). When density in middle points is approximated, using central differences in algorithm SIMPLE-TS, [1], unphysical oscillations of fluid flow at supersonic speed are obtained. The approximation of upwind first-order scheme is used to prevent these oscillations. Here are applied the same approach: density approximation in middle points is the same as an approximation for convective terms.
4 PROBLEM FORMULATION As a example we consider a 2D steady-state laminar flow around a small square particle with size a confined in a FIGURE 3. Flow geometry for a square-shaped particle with size a confined in a channel with length L ch and height H ch. plane microchannel (height H ch ) as shown in Fig. 3. The blockage ratio B = a/h ch is equal to B = 10, the inflow length is L a. The problem is considered in a local Cartesian coordinate system, which is moving with the particle. Thus for an observer moving along with the particle the problem is transformed to a consideration of a gas flow past a stationary square confined in a microchannel with moving walls. Reference parameters for this problem are: p 0 = p in, T 0 = T in and L = a. Temperature of the square and channel walls are equal to T 0. On the walls of the channel and the square velocity-slip and temperature-jump boundary conditions are imposed. The velocity-slip BC is given as: v s v w = ζ v n, where v s is velocity of the gas at the solid wall surface, v w is velocity of the wall, s v ζ = Kn local = Kn/ρ local, Kn local is the local Knudsen number, ρ local is the local density, n is the s derivative of velocity normal to the wall surface. The temperature-jump boundary condition is: T s T w = τ T, where s T s is temperature of the gas at the wall surface, T w is temperature of the wall, τ = Kn local = Kn/ρ local, T n is the derivative of temperature normal to the wall surface. The Knudsen number is Kn = s Detailed problem formulation and validation of SIMPLE-TS using upwind first-order scheme to approximate convective terms and density in middle points are presented in [1]. Subsonic gas flow In the subsonic case the Mach number is 0.1, L ch = 40, L a = The channel walls are moving with a constant velocity equal to u = (according to Mach number 0.1). The inflow boundary conditions (BC in ) are: u in is calculated from the continuity equation for first control volume on OX, similarly to the pressure driven flow in a channel, v/ x = 0, p in = 1 and T in = 1. The outflow boundary conditions (BC out ) are: u out is calculated using the continuity equation defined for the last control volume on OX. Here are compared upwind first-order scheme, results are available in [1], second-order TVD schemes Min-Mod (12) and QUICK (13), and central difference scheme. Mesh-dependent convergence is investigated on spatial steps = base, base /2, base /4 (corresponding uniform mesh: 400x100, x100 and x100), where = x = y, x is step on OX, y is step on OY, base = 0.1 is based spatial step. Here are compared profiles of horizontal velocity and temperature at the centre line along the channel, Fig. 4, obtained at different spatial steps and schemes. The profiles for horizontal velocity, obtained using different spatial steps and approximations are very close to each other, Fig. 4. (a). The profiles of temperature obtained, using different approximations at a given spatial steps are very close to each other, Fig. 4.. The groups of profiles obtained at a different spatial steps are different. On the other hand, the temperature deviations of fluid flow is negligible therefore fluid flow can be considered isothermal. The results show no advantages of second order schemes for a considered case. The velocity of subsonic case is very low (Mach number is 0.1). When Mach number is sufficiently low (close to zero) the fluid flow became creeping (Stokes flow) and convective terms can be neglected in a mathematical model of fluid. The low flow velocity is the most likely reason of independence of the numerical solution of the order of approximation of convective members. n
5 FIGURE 4. Profiles of horizontal velocity (a) and temperature obtained by upwind first-order, central difference, QUICK and Min-Mod schemes at the centre line along the channel (y = H ch /2), M = 0.1. Supersonic gas flow In the supersonic calculations the Mach number is The channel length is L ch = 50, the distance between the square and the channel inlet is L a = 5.5. The channel walls are moving with a constant velocity u = (according to Mach number ). The inflow boundary conditions (BC in ) are: u in = , v = 0, p in = 1 and T in = 1. The outflow boundary conditions (BC out ) are: u/ x = 0, v/ x = 0, p/ x = 0 and T/ x = 0. Here are compared results, obtained using first-order and second order schemes for calculation of flow past square particle in a microchannel at supersonic speed. The supersonic speed leads to the existence of a bow shock wave, which reflects from the channel walls. As a result, past the particle, we obtain a complex picture of the interaction of reflected shock waves (see Fig. 5). The shock waves have significant gradients of velocities, pressure and temperature. The calculation of this flow using first-order upwind schemes accurately requires a very fine mesh (8000x1600) cells, [1], or adaptive one. Here are compared results obtained using first-order upwind scheme, available in [1] with results obtained using second-order TVD schemes Min-Mod (12), QUICK (13) and SUPERBEE (14). Mesh-dependent convergence is investigated for a range of spatial steps for first- and second order schemes. Range FIGURE 5. Fields of the horizontal component of velocity (left part) and temperature (right part), M = 2.43, obtained using TVD scheme SUPERBEE, spatial step is = base /2 = of spacial steps for upwind first-order schemes are = base /2, base /4, base /8, base /16 (corresponding uniform mesh: x100, x100, x100 and x100), where = x = y, base = 0.1 is based spatial step. Range of spacial steps for second-order TVS schemes are = base, base /2, base /4, base /8 (corresponding uniform mesh: 500x100, x100, x100 and x100). It is expected that second-order schemes will reduce computational time in comparison with first-order schemes, but it is possible to bring instability and to obtain physical unrealistic results. Here are compared profiles of horizontal velocity and temperature at the centre line along the channel obtained using first- and second-order schemes for the approximation of convective terms and density in middle points, Fig The second-order Min-Mod scheme mesh convergence is much better than the first-order upwind scheme, Fig. 6. We accepted that the numerical solution obtained using Min-Mod scheme is accurate enough for a spatial step base /4 = 0.025, while for upwind scheme is base /16. Therefore Min-Mod scheme needs 4 times less cells in each spatial direction, which is 4 2 = 16 less cells for the two dimensional case. Mesh convergence of QUICK scheme is also rapider compared to those of upwind scheme, Fig. 7. However mesh convergence of QUICK scheme is slightly better then Min-Mod scheme. The convergence
6 of SUPERBEE scheme is better than upwind scheme, Fig. 8 and is even better than convergences of Min-Mod and QUICK schemes. SUPERBEE schemes obtain an accurate solution on a spatial step base /2 = 0.05, therefore needs 8 2 = 64 times less cells in the computational domain, then the upwind scheme. The results show many times reduction of the number of cells in a computational domain of TVD schemes, Fig. 10, compared to first-order upwind scheme. On the other hand, investigated TVD schemes required smaller time step and larger number of iterations to be algorithm converge against upwind first-order scheme. TVD scheme with best convergence SUPERBEE obtain small over-shoot in a shock wave in front of the square for spatial steps = base = base /2, Fig. 9. Other investigated TVD schemes Min-Mod and QUICK do not obtain such kind of physical unrealistic results. Nevertheless the requirements for smaller time steps and increase the number of internal iterations compared to first-order upwind scheme TVD schemes reduce many times computational time. Investigated TVD schemes are applied using their standard format and do not obtain a convergent solution for spatial step = base /16 = A round-off error is the most probable reason TVD schemes to be unstable using small spatial steps. The implementation of TVD schemes in algorithm has to be changed so the round-off error to be reduced. Results show that mesh convergence of TVD schemes closer to the upper TVD region (Fig. 2, ) is faster than TVD schemes closer to lower TVD region limit, Fig. 6, Fig. 7 and Fig. 8. On the other hand TVD schemes closer to the upper TVD limit (SUPERBEE) needs smaller time step and more internal iteration to be algorithm SIMPLE-TS converged according to TVD schemes closer to lower TVD limit region (Min-Mod). TVD schemes at the lower limit of the TVD region (Min-Mod) required less time step and more internal iterations then upwind first-order scheme, which correspond to Ψ(r) = 0. This investigation shows that schemes, closer to the upper limit of TVD region have rapider mesh convergence, then schemes, closer to to lower limit of TVD region and first-order upwind scheme (Ψ(r) = 0). On the other hand, schemes, closer to the upper limit of TVD region obtain small over-shoot (SUPERBEE) and are less stable then schemes, closer to to lower limit of TVD region and first-order upwind scheme (Ψ(r) = 0). The results show that use of schemes, closer to the upper TVD region will give us rapider mesh convergence small reducing of stability of the algorithm and possible appearance of small over-shoots. (a) FIGURE 6. Profiles of horizontal velocity (a) and temperature obtained by upwind first-order and Min-Mod schemes at the centre line along the channel (y = H ch /2), M = (a) FIGURE 7. Profiles of horizontal velocity (a) and temperature obtained by upwind first-order and QUICK schemes at the centre line along the channel (y = H ch /2), M = 2.43.
7 (a) FIGURE 8. Profiles of horizontal velocity (a) and temperature obtained by upwind first-order and SUPERBEE schemes at the centre line along the channel (y = H ch /2), M = FIGURE 9. Profiles of temperature obtained by upwind first-order and SUPERBEE schemes in front of the square at the centre line along the channel (y = H ch /2), M = (a) FIGURE 10. Profiles of horizontal velocity (a) and temperature obtained by upwind first-order and Min-Mod, QUICK and SUPERBEE schemes at the centre line along the channel (y = H ch /2), M = CONCLUSIONS It this paper is investigated mesh convergence of first-order and second order schemes for the approximation of convective terms of Navier-Stokes-Fourier system of equations described gas microflow. A test case is flow past a square in a microchannel at subsonic (M = 0.1) and supersonic (M = 2.43) speeds. Results, obtained from subsonic case show no advantages of any of compared schemes (first-order upwind and second order central difference, TVD
8 Min-Mod and TVD QUICK schemes. Results, obtained from supersonic case is that TVD schemes obtained the same results as first-order upwind schemes with 4 2 = 16 (Min-Mod and QUICK) to 8 2 = 64 (SUPERBEE) times less cells. The disadvantages of TVD schemes then first-order upwind scheme are: the usage of smaller time step and the increase of iterations in internal loops. SUPERBEE scheme, which mesh convergence is rapider according to investigated schemes, obtains small over-shoot in a shock wave in front of the square. Investigated TVD schemes (Min-Mod, QUICK, SUPERBEE) are not converged for small spatial steps = base /16 = Nevertheless the requirements for smaller time steps and increase the number of internal iterations compared to first-order upwind scheme TVD schemes reduce many times resources to calculate supersonic fluid flow and are highly recommended for this flow regime. ACKNOWLEDGMENTS The author appreciates the financial support by the NSF of Bulgaria under Grant No DMU 03/37, since This work makes use of results produced by the High-Performance Computing Infrastructure for South East Europe s Research Communities (HP-SEE), a project co-funded by the European Commission (under contract number ) through the Seventh Framework Programme. HP-SEE involves and addresses specific needs of a number of new multidisciplinary international scientific communities (computational physics, computational chemistry, life sciences, etc.) and thus stimulates the use and expansion of the emerging new regional HPC infrastructure and its services. Full information is available at: REFERENCES 1. K. S. Shterev, and S. K. Stefanov, Journal of Computational Physics 229, (2010), ISSN , URL 2. J. P. Boris, and D. L. Book, Journal of Computational Physics 11, (1973), ISSN , URL http: // 3. J. P. Boris, and D. L. Book, Solution of continuity equations by the method of flux-corrected transport, 1976, pp B. van Leer, Journal of Computational Physics 14, (1974), ISSN , URL sciencedirect.com/science/article/pii/ B. V. Leer, Journal of Computational Physics 23, (1977), ISSN , URL sciencedirect.com/science/article/pii/ B. V. Leer, Journal of Computational Physics 23, (1977), ISSN , URL sciencedirect.com/science/article/pii/ x. 7. B. van Leer, Journal of Computational Physics 32, (1979), ISSN , URL sciencedirect.com/science/article/pii/ A. Harten, Journal of Computational Physics 49, (1983), ISSN , URL sciencedirect.com/science/article/pii/ A. Harten, SIAM Journal on Numerical Analysis 21, 1 23 (1984), URL / , P. K. Sweby, SIAM Journal on Numerical Analysis, vol. 21, no. 5, pp , (1984). 11. P.L.Roe, Lectures in Applied Mathematics 22, (1985). 12. S. Osher, and S. Chakravarthy, SIAM Journal on Numerical Analysis 21, (1984), URL org/doi/abs/ / , Openfoam (????), URL B. P. Leonard, International Journal for Numerical Methods in Fluids 8, (1988), ISSN , URL S. Stefanov, V. Roussinov, and C. Cercignani, Physics of Fluids 14, (2002). 16. H. K. Versteeg, and W. Malalasekra, An Introduction to Computational Fluid Dynamics: The Finite Volume Method, Prentice Hall, Pearson, 2007, 2nd edn. 17. B. Leonard, Computer Methods in Applied Mechanics and Engineering 19, (1979), ISSN , URL
Gaseous Slip Flow in Three-Dimensional Uniform Rectangular Microchannel
Gaseous Slip Flow in Three-Dimensional Uniform Rectangular Microchannel Khaleel Khasawneh, Hongli Liu and Chunpei Cai Department of Mechanical and Aerospace Engineering, New Mexico State University, Las
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
More informationBehaviour of microscale gas flows based on a power-law free path distribution function
Behaviour of microscale gas flows based on a power-law free path distribution function Nishanth Dongari, Yonghao Zhang and Jason M Reese Department of Mechanical Engineering, University of Strathclyde,
More informationA MICROBEARING GAS FLOW WITH DIFFERENT WALLS TEMPERATURES
Mili}ev, S. S., et al.: A Microbearing Gas Flow with Different Walls Temperatures THERMAL SCIENCE, Year 01, Vol. 16, No. 1, pp. 119-13 119 A MICROBEARING GAS FLOW WITH DIFFERENT WALLS TEMPERATURES by Snežana
More informationMicro-Scale Gas Transport Modeling
Micro-Scale Gas Transport Modeling Continuum & Slip Flow Regimes: Navier-Stokes Equations Slip Boundary Conditions U g U w = ( σ ) σ U Kn n Slip, Transitional & Free Molecular: Direct Simulation Monte
More informationChapter 5. The Finite Volume Method for Convection-Diffusion Problems. The steady convection-diffusion equation is
Chapter 5 The Finite Volume Method for Convection-Diffusion Problems Prepared by: Prof. Dr. I. Sezai Eastern Mediterranean University Mechanical Engineering Department Introduction The steady convection-diffusion
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationStefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran
Stefan Stefanov Bulgarian Academy of Science, Bulgaria Ali Amiri-Jaghargh Ehsan Roohi Hamid Niazmand Ferdowsi University of Mashhad, Iran Outlines: Introduction DSMC Collision Schemes Results Conclusion
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationRAREFACTION EFFECT ON FLUID FLOW THROUGH MICROCHANNEL
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationDevelpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow
Develpment of NSCBC for compressible Navier-Stokes equations in OpenFOAM : Subsonic Non-Reflecting Outflow F. Piscaglia, A. Montorfano Dipartimento di Energia, POLITECNICO DI MILANO Content Introduction
More informationWP8: Parallel algorithms and GRID applications for computer modeling of processes and phenomena in Micro-Electro-Mechanical Systems (MEMS). 1.
WP8: Parallel algorithms and GRID applications for computer modeling of processes and phenomena in Micro-Electro-Mechanical Systems (MEMS). 1. Major Activities and Results Task 8.1: Computer Simulation
More informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More informationTHE EFFECTS OF LONGITUDINAL RIBS ON ENTROPY GENERATION FOR LAMINAR FORCED CONVECTION IN A MICROCHANNEL
THE EFFECTS OF LONGITUDINAL RIBS ON ENTROPY GENERATION FOR LAMINAR FORCED CONVECTION IN A MICROCHANNEL Nader POURMAHMOUD, Hosseinali SOLTANIPOUR *1,, Iraj MIRZAEE Department of Mechanical Engineering,
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationDSMC simulation of the gas flow through a bend and a short microchannel
Journal of Physics: Conference Series DSMC simulation of the gas flow through a bend and a short microchannel To cite this article: N K Kulakarni and S K Stefanov 2012 J. Phys.: Conf. Ser. 362 012014 View
More informationResearch Article. Slip flow and heat transfer through a rarefied nitrogen gas between two coaxial cylinders
Available online wwwjocprcom Journal of Chemical and Pharmaceutical Research, 216, 8(8):495-51 Research Article ISSN : 975-7384 CODEN(USA) : JCPRC5 Slip flow and heat transfer through a rarefied nitrogen
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.
More informationHEAT TRANSFER COEFFICIENT CHARACTERIZATION AT THE SOLAR COLLECTOR WALL-FLUID INTERFACE
SASEC15 Third Southern African Solar Energy Conference 11 13 May 15 Kruger National Park, South Africa HEAT TRANSFER COEFFICIENT CHARACTERIZATION AT THE SOLAR COLLECTOR WALL-FLUID INTERFACE Mébarki Ghazali*
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationContribution to the gas flow and heat transfer modelling in microchannels
Applied and Computational Mechanics 3 (2009) 63 74 Contribution to the gas flow and heat transfer modelling in microchannels H. Klášterka a,, J. Vimmr b,m.hajžman b a Faculty of Mechanical Engineering,
More informationENHANCED HEAT TRANSFER IN LAMINAR PULSATING FLOW PAST A FLAT PLATE
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta ENHANCED HEAT TRANSFER IN LAMINAR PULSATING FLOW PAST A FLAT PLATE Cárdenas Miranda, A.*
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationThe Use of CFD Simulations in Learning Fluid Mechanics at the Undergraduate Level
Presented at the COMSOL Conference 2009 Boston The Use of CFD Simulations in Learning Fluid Mechanics at the Undergraduate Level Marc K. Smith Woodruff School of Mechanical Engineering Georgia Institute
More informationMicrochannel flow in the slip regime: gas-kinetic BGK Burnett solutions
J. Fluid Mech. (2004), vol. 513, pp. 87 110. c 2004 Cambridge University Press DOI: 10.1017/S0022112004009826 Printed in the United Kingdom 87 Microchannel flow in the slip regime: gas-kinetic BGK Burnett
More informationAmir-Mehran Mahdavi 1, Nam T.P. Le 1, Ehsan Roohi 1,*, Craig White 2
Thermal rarefied gas flow investigations through micro/nano backward-facing step: Comparison of DSMC and CFD subject to hybrid slip and jump boundary conditions Amir-Mehran Mahdavi 1, Nam T.P. Le 1, Ehsan
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationTHE CURVATURE EFFECT ON THERMOACOUSTIC WAVES IN A COUETTE RAREFIED GAS FLOW * Peter Gospodinov, Vladimir Roussinov, Dobri Dankov
МАТЕМАТИКА И МАТЕМАТИЧЕСКО ОБРАЗОВАНИЕ, 2015 MATHEMATICS AND EDUCATION IN MATHEMATICS, 2015 Proceedings of the Forty Fourth Spring Conference of the Union of Bulgarian Mathematicians SOK Kamchia, April
More informationNatural convection adjacent to a sidewall with three fins in a differentially heated cavity
ANZIAM J. 48 (CTAC2006) pp.c806 C819, 2007 C806 Natural convection adjacent to a sidewall with three fins in a differentially heated cavity F. Xu 1 J. C. Patterson 2 C. Lei 3 (Received 31 August 2006;
More informationNumerical Methods for Problems with Moving Fronts Orthogonal Collocation on Finite Elements
Electronic Text Provided with the Book Numerical Methods for Problems with Moving Fronts by Bruce A. Finlayson Ravenna Park Publishing, Inc., 635 22nd Ave. N. E., Seattle, WA 985-699 26-524-3375; ravenna@halcyon.com;www.halcyon.com/ravenna
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationOn the Effect of the Boundary Conditions and the Collision Model on Rarefied Gas Flows
On the Effect of the Boundar Conditions and the Collision Model on Rarefied Gas Flows M. Vargas a, S. Stefanov a and D. Valougeorgis b a Institute of Mechanics, Bulgarian Academ of Sciences, Acad. G. Bonchev
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationDSMC Modeling of Rarefied Flow through Micro/Nano Backward-Facing Steps
DSMC Modeling of Rarefied Flow through Micro/Nano Backward-Facing Steps Amir-Mehran Mahdavi 1, Ehsan Roohi 2 1,2- Department of Mechanical Engineering, Faculty of Engineering, Ferdowsi university of Mashhad,
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationRAREFIED GAS HEAT CONDUCTION BETWEEN TWO COAXIAL CYLINDERS UNDER PULSATING TEMPERATURE BOUNDARY CONDITION *
МАТЕМАТИКА И МАТЕМАТИЧЕСКО ОБРАЗОВАНИЕ, 2017 MATHEMATICS AND EDUCATION IN MATHEMATICS, 2017 Proceedings of the Forty-sixth Spring Conference of the Union of Bulgarian Mathematicians Borovets, April 9 13,
More informationDSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations
DSMC-Based Shear-Stress/Velocity-Slip Boundary Condition for Navier-Stokes Couette-Flow Simulations J. R. Torczynski and M. A. Gallis Engineering Sciences Center, Sandia National Laboratories, P. O. Box
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationA wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations
A wall-function approach to incorporating Knudsen-layer effects in gas micro flow simulations D. A. Lockerby 1, J. M. Reese 2 and M. A. Gallis 3 1 Department of Mechanical Engineering, King s College London,
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationThe effect of Entry Region on Thermal Field
The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving
More informationNon-Equilibrium Flow Modeling Using High- Order Schemes for the Boltzmann Model Equations
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 28 Non-Equilibrium Flow Modeling Using High- Order Schemes for the Boltzmann
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
More informationImpact of Magnetic Field Strength on Magnetic Fluid Flow through a Channel
ISSN: 2278-8 Vol. 2 Issue 7, July - 23 Impact of Magnetic Field Strength on Magnetic Fluid Flow through a Channel S. Saha, S. Chakrabarti 2 Dept. of Mechanical Engineering, Dr. Sudhir Chandra Sur Degree
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationSummary We develop an unconditionally stable explicit particle CFD scheme: Boltzmann Particle Hydrodynamics (BPH)
Summary Steps in BPH Space is divided into Cartesian cells Finite number of particles ~10 6-10 7 Particles fly freely between t n and t n+1 Mass, momentum and total energy are conserved Relax into a (LTE)
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationScientific Computing I
Scientific Computing I Module 10: Case Study Computational Fluid Dynamics Michael Bader Winter 2012/2013 Module 10: Case Study Computational Fluid Dynamics, Winter 2012/2013 1 Fluid mechanics as a Discipline
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationApplication of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows
Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows Timothy R. Deschenes and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor,
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationContents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1
Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationViscous dissipation and temperature dependent shear thinning rheology
Viscous dissipation and temperature dependent shear thinning rheology M. Filippucci 1, A. Tallarico 1, M. Dragoni 2 1 Dipartimento di Scienze della Terra e Geoambientali, Università degli Studi di Bari
More informationSIMULATION OF MIXED CONVECTIVE HEAT TRANSFER USING LATTICE BOLTZMANN METHOD
International Journal of Automotive and Mechanical Engineering (IJAME) ISSN: 2229-8649 (Print); ISSN: 2180-1606 (Online); Volume 2, pp. 130-143, July-December 2010 Universiti Malaysia Pahang DOI: http://dx.doi.org/10.15282/ijame.2.2010.3.0011
More informationNUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00394 NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP
More informationThe limits of Navier-Stokes theory and kinetic extensions for describing small-scale gaseous hydrodynamics
PHYSICS OF FLUIDS 18, 111301 2006 The limits of Navier-Stokes theory and kinetic extensions for describing small-scale gaseous hydrodynamics Nicolas G. Hadjiconstantinou a Department of Mechanical Engineering,
More informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationMIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM
THERMAL SCIENCE, Year 015, Vol. 19, No. 1, pp. 119-18 119 MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM by Gurminder SINGH *a and Oluwole Daniel MAKINDE
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationSimulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code
1 st. Annual (National) Conference on Industrial Ventilation-IVC2010 Feb 24-25, 2010, Sharif University of Technology, Tehran, Iran IVC2010 Simulation and improvement of the ventilation of a welding workshop
More informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationLattice Boltzmann Simulations of 2D Laminar Flows past Two Tandem Cylinders
Lattice Boltzmann Simulations of 2D Laminar Flows past Two Tandem Cylinders Alberto Mussa a, Pietro Asinari a,, and Li-Shi Luo b a Department of Energetics, Politecnico di Torino, Torino 10129, Italy Tel.:
More information