Implicit Finite Volume Method to Simulate Reacting Flow
|
|
- Annabel Dixon
- 5 years ago
- Views:
Transcription
1 43rd I erospace Sciences Meeting and Ehibit, -3 Jan. 5, Reno, Nevada Implicit Finite Volume Method to Simulate Reacting Flow Masoud Darbandi and raz Banaeizadeh Sharif University of Technology, Tehran, P.O. Bo , Iran Gerry E. Schneider University of Waterloo, Waterloo, Ontario, NL 3G, Canada In this work, an efficient bi-implicit strategy is suitably developed within the contet of a finite volume element approach in order to solve turbulent reactive flow governing equations. Based on the essence of control-volume-based finite-element methods, the formulation retains the geometrical fleibility of the pure finite element methods while derives the discrete algebraic governing equations through using the conservation balance applied to discrete control volumes distributed all over the solution domain. The physical influence upwinding scheme is used to approimate the advection flues at all cell faces. While respecting the physics of flow, this scheme also provides the necessary coupling of velocity and pressure fields. The two-equation k ɛ turbulence model and one step miture fraction chemistry equation are simultaneously solved in a semi-coupled manner in order to achieve a better prediction of both the transport of turbulent species and the transport of mass fraction species. The validation of the current numerical results is fulfilled by comparing them with eperimental data and other available numerical results. I. Nomenclature B body force C pi specific heat of species i at constant pressure D diffusion coefficient f miture fraction h enthalpy h F heat reaction N finite element shape function n total number of species P pressure R gas constant T temperature u, v velocity components ū friction velocity V velocity vector ssociate Professor, Department of erospace Engineering. Graduate Student, Department of erospace Engineering. Professor and Chair, Department of Mechanical Engineering, I Fellow. Copyright c 5 by M. Darbandi. Published by the merican Institute of eronautics and stronautics, Inc. with permission. of
2 ū friction velocity, y Cartesian coordinates Y i mass fraction of species i z a reative of the etensive quantities Υ dissipation rate of turbulent kinetic energy at nodes ɛ dissipation rate of turbulent kinetic energy κ turbulent kinetic energy µ laminar viscosity µ e effective viscosity µ t turbulent viscosity ρ miture density σ viscous shear tensor. Subscripts i p st chemical species the grid point adjacent to the solid wall stoichiometric condition B. Superscripts + non-dimensional magnitudes II. Introduction The reacting flow in combustion chambers mostly occurs in comple geometries which in turn need robust numerical tools to treat the resulting comple turbulent reactive flow fields. Finite-difference-based methods do not straightforwardly provide efficient means to solve the flow fields occurred in comple geometries. On the contrary, the finite element methods provide easy treatment of complicated domains. The past eperience has shown that a robust approach toward handling the comple geometries is the finitevolume-based finite-element one. This method guarantees the great benefits of both finite volume and finite element methods. Finite volume element method incorporated with a suitable physical upwinding influence scheme effectively enhances the capabilities of the current dual-based method. In this work, the primitive implicit strategy introduced by Darbandi and Schneider 3 is suitably etended in order to solve the turbulent reactive flow governing equations. Back to the primary problem with the continuity equation, Darbandi and Bostandoost 4 have developed a fully implicit finite volume method which eliminates the need for the pressure Poisson equation and takes the pressure role directly into the continuity equation. The approach is fully implicit and produces a 7-diagonal matri in treating the D Navier-Stokes equations. Darbandi, et al. 5 etended the formulation for solving the diffusive flame. In this work, the original approach is applied to discretize turbulence and fast chemistry equations in a way that continuity, -momentum and y-momentum are solved implicitly in one matri and κ, ɛ, miture fraction and energy equations in another matri. This approach leads to a bi-implicit treatment of the governing equations. III. The Governing Equations The conservative equations for a reacting flow can be categorized into fluid flow and species transport equations. The radiation heat transfer effect is ignored in this study. The fluid flow governing equations consist of the conservation statements for mass, momentums, and energy. The incompressible steady governing of
3 equations are written as (ρ V ) = () (ρ V u) = p + (µ e u) () (ρ V v) = p y + (µ e v) (3) where µ e is the effective viscosity defined by µ e = µ + µ t. The eddy turbulent viscosity µ t is calculated from the turbulent kinetic energy κ equation and its dissipation rate ɛ equation. The transport equations for these two components are written as (ρ V κ) = ( µ e σ κ κ) + (Gκ ρɛ) (4) (ρ V ɛ) = ( µ e σ ɛ ɛ) + ɛ κ (c G κ c ρɛ) (5) In these equations, G κ and c G κ are the production terms and ρɛ and c ρɛ are the destruction terms. The epression for G κ is given by G κ = µ t {[( u ) + ( v ) ] + [ u y + v ] } (6) The derived κ and ɛ magnitudes from Eqs. (4)-(5) are used to calculate the eddy viscosity definition, i.e., µ t = c d ρ κ ɛ The constants in κ and ɛ equations are σ κ =., σ ɛ =.3, c =.44, c =.9, and c d =.9. s is known, the calculated turbulent energy and its dissipation rate are valid sufficiently far from the solid walls and mainly in the interior region. The viscous effects dominate near the solid walls which subsequently require etra considerations. Depending on the flow Reynolds number, there are different choices such as low Reynolds number model and wall function approimation to cure the problem. In the work, the wall function choice is used to describe the flow close to the wall. The elliptic nature of κ and ɛ transport equations directs the users to specify boundary conditions at all domain boundaries. s was mentioned, we apply wall function at solid walls. Indeed, the region close to the solid walls can be divided into two layers ) a laminar sublayer or viscous sublayer where viscous effects are dominant and ) turbulent sublayer. In our model, it is supposed that the grid node p adjacent to the solid wall is far enough from the wall to be in turbulent sublayer region where the velocity is parallel to the wall at this node. The condition near the wall is calculated from the known logarithmic law. The turbulent kinetic energy and its dissipation rate at node p are then calculated from κ p = ɛ p = (7) ū cd (8) ū3 k y p (9) where c d and k constants are.9 and.4, respectively. y p is the actual distance between node p and wall and ū is the friction velocity. To specify the friction velocity at node p, the non-dimensional distance from the solid wall is defined as y p + = ρūy p /µ. Using this definition, the non-dimensional velocity u + is determined from { y p +.63 u + p = k ln y+ p () y p + >.63 u + p = y p + 3 of
4 Then, the friction velocity is derived from ū p = u p u + p For the combustion model, we consider a diffusion flame and a fast chemistry with one step irreversible chemical reaction. It is assumed that the chemistry is sufficiently fast and that the intermediate species do not play a significant role. Considering a mass flow rate of M = (Kg/s) for the air/fuel miture, it consists of a fuel rate of f(kg/s) and an air rate of ( f)(kg/s). Considering the above mass flow rates, any etensive property, such as z, of the miture resulting from the mi of these two streams can be written as () f = z M z z F z () If F indicates the fuel and indicates the oidizer, the miture fraction yields f = [Y F (F/O) st Y O ] M + (F/O) st Y O, + (F/O) st Y O, (3) From the chemical equilibrium assumption, no oidant s if a miture fraction richer than stoichiometric condition is employed. dditionally, no fuel s if the miture fraction is lower than the stoichiometric condition. 6 For any adiabatic operation under the assumption of unit Lewis number, the enthalpy equation is similar to the miture fraction equation and the enthalpy can be calculated directly from the miture fraction and the inlet enthalpy value. Thus, for all combustion related variables such as enthalpy, mass fraction of fuel, oygen, and the products, it is sufficient to solve the transport equation for miture fraction f. 7 The transport equation for f is then given by (ρ V f) + (D f) = (4) ssuming a unit Lewis number, the reactive energy equation is reduced to where the enthalpy is defined by h = Y F h F + (ρ V h) = (µ h) (5) n i= T T ref Y i C pi (T )dt + V (6) Eventually, the density is calculated from the equation of state which is written as P = ρrt n Y i W i (7) IV. The Domain Discretization The solution domain is broken into a huge number of quadrilateral elements. The elements fully cover the solution domain with no overlapping. Figure shows a small part of the solution domain. Nodes are located at the corners of elements and are shown by circles. The nodes are the locations of the unknown variables. Each node belongs to four neighboring elements. There are four quadrilaterals which enclose node P in figure. To utilize the benefits of cell-centered schemes, each element is divided into four quadrilaterals by the help of its medians. The median is demonstrated by dashed-line in this figure. The cells are then constructed from the proper assemblage of these sub-quadrilaterals. s is seen, irrespective of the shape and distribution of the elements, each node is surrounded by a number of sub-quadrilaterals. The proper assemblage of neighboring sub-quadrilaterals around any non-boundary node creates a complete cell. 4 of
5 Integration Point Control Volume Element Node t s y Figure. part of the solution domain illustrating four elements, one complete finite volumes, siteen sub-volumes, and eight cell faces. V. The Computational Modelling To utilize the advantages of finite element volume methods, the governing equations are initially integrated over an arbitrary volume, e.g., the shaded area or the cell face shown in Fig.. The employment of Gauss divergence theorem to the governing equations leads to ρv d = (8) u(ρv ) d = v(ρv ) d = p d + p d y + µ e ( u d) (9) µ e ( v d) () The above integrals are evaluated over the surface which encloses each cell. The surface area is indicated by. The above equations are suitably discretized using finite difference scheme and finite element interpolations. In the above epressions, d=d i d y j is a normal vector to the edges of cell. Using this definition, the above integrals can be evaluated by summation over the faces that enclose the cell center, i.e., ns ns ns ns ns [ρû u d + ρˆv u d y )] i = (p d ) i + [ρ(u d + v d y )] i = () ns ns [ρû v d + ρˆv v d y )] i = (p d y ) i + [ ( u µ e d + u )] y d y i [ ( v µ e d + v )] y d y where i counts the number of cell faces from to ns. There are 8 cell faces around non-boundary cells. To linearize the governing equations, the hat over û and ˆv indicates that these velocity components are approimated from the known magnitudes of the preceding iteration. Such approimation is essential to linearize the nonlinear momentum convection terms. The rest of procedure is to relate the cell face magnitudes i () (3) 5 of
6 (identified by lower case letters such as u, v and p variables) directly to the nodal magnitudes (identified by upper case letters such as U, V and P variables) which re the locations for the unknown variables in the current algorithm. simple idea for treating the right-hand-side terms is to use the finite element shape functions N j=...4. The treatment results in p i = φ ξ = i 4 N P j (4) j= 4 j= N ξ Φ j (5) where p i identifies the magnitude of p at the mid-point of ith edge of the cell face. The j notation counts the node numbers of an element where the ith cell face is located inside it. dditionally, the variable ξ res either or y coordinates and φ (and Φ) res either u (and U) or v (and V ) velocity components. s was mentioned, lower and upper case letters are used to re cell face and nodal magnitudes, respectively. The above approimations end the pressure and diffusion term treatments at all cell faces. However, more sophisticated epressions are required to treat the convection terms. In fact, the treatment should not disregard the convection-diffusion physics and concept. To respect the correct physics of the convection, Reference 5, 8 employ an upwind-based scheme (known as a physical influence scheme) within quadrilateral and triangular elements, respectively. The references show that a physical-based treatment of the -momentum governing equation can result in 4 4 φ i = α Φ j + β P j + γ i (6) j= where α, β, and γ re matri, matri, and vector coefficients, respectively. The above statement indicates that φ ( u, v) at cell face can be approimated by the proper assemblage of Φ ( U, V ) and P influences. In fact, this approimation can be regarded as a pressure-weighted upwind scheme. s is observed in Eq.(6), contrary to the problems raised in the Introduction section, the pressure field is not obtained through solving the pressure Poisson equation in this study. In fact, the continuity preserves its original identification and the direct contribution of the pressure terms are forced through the mass flu statements at the cell faces. The substitution of Eqs. (4)-(6) in Eqs. ()-(3) provides a set of algebraic equations for each cell. It is given by c pp c up c vp c pu c uu c vu c pv c uv c vv P j U j V j j= = where i and j count the global node numbers, i.e., i, j =... Nnode. It should be mentioned that the matri is a diagonal one which is normally encountered in implicit-based finite-element methods. Therefore, it is strongly sparse and need sparse solution strategies. The coefficients in the global assembled matri is identified by c. The first letter in each superscript depicts the type of equation, i.e., p, u and v indicate continuity, -momentum and y-momentum equations, respectively. The second letter in the superscripts indicates which unknown the coefficient belongs to. The right-hand-side vector is shown by d. Since all U, V, and P unknowns eplicitly appear in all the three governing equations, there is strong coupling between the pressure and velocity fields and the pressure checkerboard problem may not occur. In a similar manner, the procedure etended for the fluid flow governing equations can be repeated for the turbulent kinetic energy, Eq.(4), the dissipation rate of turbulent kinetic energy, Eq.(5), enthalpy of d p i d u i d v i (7) 6 of
7 chemical species, Eq.(5), and mass fraction of chemical species, Eq.(4). Similar to the fluid flow governing equations, the aforementioned four governing equations are solved simultaneously. This implicit procedure enhances the stability of the etended method. The consideration of this latter point is so crucial in solving turbulent reacting flow using the current pressure-weighted approach. If we follow the preceding procedure, another set of algebraic equations is derived. The new derived matri is given by c κκ c ɛκ c eκ c fκ c κɛ c ɛɛ c eɛ c fɛ c κe c ɛe c ee c fe c κf c ɛf c ef c ff K j Υ j H j F j = s was mentioned in the Introduction section, the algorithm is a bi-implicit one and solves two sets of algebraic equations in each iteration, one for the fluid flow governing equations and the other one for energy, turbulence, combustion, and mass fraction of species. d κ i d ɛ i d e i d f i (8) VI. The Results In this section, the new developed procedure is tested in two stages. t the first stage, the chemical reaction switch is off and only the transport of turbulence species is investigated using our new convectiondiffusion model. To eamine the etended model, turbulent flow over a backward-facing step is chosen as the test case. The chosen test case is known as a benchmark case which is mainly used to test newly developed numerical algorithms. There is a wide range of either eperimental 9 or numerical solutions for this test case. Figure shows the geometry (not scaled) of the calculation domain and the specified boundary conditions. Following, et al. 9 and, a uniform unit velocity is considered at the inlet. The flow Reynolds number based on the step height and inlet velocity is 696. t the eit, the flow is assumed to be fully developed. The specified wall functions are applied at all solid walls. Different fine grid resolutions have been used to ensure the grid independency of the solution..5. Wall Function Wall Function 4. Figure. Turbulent flow in backward-facing step. Figure 3 shows the turbulent kinetic energy distributions at two different locations downstream of the step, i.e., =. and 4.. The current results are compared with those of, et al. 9 and. Despite using general κ ɛ model, there are good agreement among the ed results. Figure 4 illustrates the mean aial velocity profiles at =.33 and There are ecellent agreements between the current solutions and those of the references. The mean velocity profiles at =.33 and 5.33 indicate that there is a recirculation zone behind the step., et al. 9 report that the length of this recirculation is about 6 to 8 following their eperimental study. The numerical solution of predicts a length of 5.59 using 7 of
8 .5.5 Y (X=) Y (X=4.) K Figure 3. Turbulent kinetic energy distributions at =. and 4. and comparing them with those of, et al. 9 and. K.5 Y (X=.33) Y (X=5.33).75 - U Figure 4. The velocity profiles at =.33 and 5.33 and comparing them with those of, et al. 9 and. U standard κ ɛ model. The method calculates a length of Finally, Fig. 5 shows both mean aial velocity profile and turbulent kinetic energy distribution far downstream of the step at =7.67 and 8., respectively. In this figure, the behavior and conclusions are very similar to those discussed over Figs. 3 and 4. t the second stage, the turbulent reacting flow and the transport of chemical species are investigated. The chosen test case is a confined diffusion flame. The solution domain geometry is observed in Fig. 6. This test case has been widely investigated either eperimentally or numerically in cylindrical coordinates. Unfortunately, the authors did not find any investigation in Cartesian coordinates for the non-premied turbulent flame. Therefore, we compare our results with those obtained in cylindrical coordinates. Of course, this can cause some uncertainty in the evaluation of the current results in the validation part. Because of the symmetry, only one half of the physical domain is studied. The turbulent-combustion interaction and the 8 of
9 .5.5 Y (X=7.67) Y (X=8) K.75 U Figure 5. Turbulent kinetic energy distributions at =7.67 and the velocity profiles at =8. and comparing them with those of, et al. 9 and. y IR D d=8 mm d=. mm d3=8.8 mm R=.6 mm L=54 mm FUEL d3 d d L Figure 6. Confined turbulent diffusion flame. fluctuation of species are not taken into account. Considering the dimensions given in Fig. 6, the boundary conditions for the chosen geometry are defined as ) y < d3 is the fuel (which is pure methane) inlet with a uniform velocity at this port equal to.3 m/s and a uniform temperature of 3 K and ) d < y < d is the air inlet with Y O =., Y N =.79, a uniform velocity of 34.3 m/s, and an average temperature of 589 K. t this stage of our investigations, we the results obtained after solving the above confined turbulent methane diffusion flame. Figures 7-9 the miture fraction distributions and molar concentrations at three different aial locations of =.95,.75, and.46 m. The agreement between the current solutions and those of references is good. The current solutions are almost over-specified in all cases. This might be because of either considering a rectangular burner geometry instead of an aisymmetry one (which has been investigated by the other workers) or ignoring the turbulent-combustion interaction in our algorithm. Figures 7-9(left) show the molar concentration of species which are determined from the miture fraction quantities. Generally speaking, s moving forward in the burner channel, the miture fraction decreases; thus, resulting in a decrease in fuel concentration. 9 of
10 =.95 m Miture Fraction =.95 m Present Elkaim et al. Eperimental Molar Concentration CH4 O CO HO N Figure 7. The miture fraction distributions and the distributions of molar concentration of species and a comparison with those reported in Smoot and Lewis and Elkaim, et al., =.95 m =.75 m Miture Fraction =.75 m Present Elkaim et al. Eperimental Molar Concentration CH4 O CO HO N Figure 8. The miture fraction distributions and the distributions of molar concentration of species and a comparison with those reported in Smoot and Lewis and Elkaim, et al., =.75 m. VII. Conclusion The Navier-Stokes equations were implicitly solved and the results were used to simulate both turbulent and combustion behavior in a reacting turbulent flow using a bi-implicit algorithm. The implementation of a new pressure-weighted upwinding scheme in a finite element volume contet improves the accuracy and performance of the etended algorithm. The fast chemistry assumption is used in our combustion model. Two test cases including the turbulent flow in a backward-facing step and the confined turbulent diffusion flame were used to validate the etended algorithm. The results perform very good agreement with the eperimental data and other numerical solutions. The results of chemical simulation indicates that the physical model has to be improved in order to enhance the overall accuracy of the algorithm. of
11 .9.9 Miture Fraction =.46 Present Elkaim et al. Eperimental Molar Concentration =.46 CH4 O CO HO N Figure 9. The miture fraction distributions and the distributions of molar concentration of species and a comparison with those reported in Smoot and Lewis and Elkaim, et al., =.46 m. References Schneider, G.E., and Raw, M.J., Control Volume Finite Element Method for Heat Transfer and Fluid Flow Using Colocated Variables -.Computational Procedure, Numerical Heat Transfer, Vol., No.4, 987, pp Darbandi, M. and Banaeizadeh,., Parallel Computation of the Navier-Stokes Equations Using Implicit Finite Volume Method, Proceedings of the 4th European Congress on Computational Methods in pplied Sciences and Engineering, ECCOMS 4, Jyvaskyla, Finland, July 4-8, 4. 3 Darbandi, M., and Schneider, G.E., nalogy-based Method for Solving Compressible and Incompressible Flows, Journal of Thermophysics and Heat Transfer, Vol., No., 998, pp Darbandi, M., and Bostandoost S.M., New Formulation Toward Unifying the Velocity Role in Collocated Variable rrangement, Numerical Heat Transfer, Part B Vol.47, 5, in press. 5 Darbandi, M., Banaeizadeh,., and Schneider, G.E., Parallel Computation of a Mied Convection Problem Using Fully- Coupled and Segregated lgorithms, Paper no. SME-IMECE , SME International Mechanical Engineering Congress, naheim, California, US, Nov. 4-9, 4. 6 Kuo, K.K., Principles of Combustion John Wiley & Sons, New York, Elkaim, D., Reggio, M., and Camarero, R., Numerical Solution of Reactive Laminar Flow by a Control-Volume Based Finite-Element Method and the Vorticity-Streamfunction Formulation, Numerical Heat Transfer, Part B, Vol., 99, pp Darbandi, M., Mazaheri-Body, K., and Vakilipour, S., Pressure-Weighted Upwind Scheme in Unstructured Finite- Element Grids, in Numerical Mathematics nd dvanced pplications, Editted by M. Feistauer, V. Dolejsi, P. Knobloch, and K. Najzar, Springer-Verlag, 4, pp , J., Kline, S.J., and Johnston, J.P., Investigation of a Reattaching Turbulent Shear Layer: Flow over a Backward- Facing Step, Journal of Fluids Engineering, SME Trans., Vol., 984, pp.7-74., J.L., Evaluation of FIDP on Some Classical Laminar and Turbulent Benchmarks, International Journal for Numerical Methods in Fluids, Vol.8, 988, pp Smoot, J.L., and Lewis, H.M., Turbulent Gaseous Combustion: Part I, Local Species Concentration Measurements, Combustion and Flame, Vol. 4, 98, pp Elkaim, D., Reggio, M., and Camarero, R., Control Volume Finite-Element Solution of Confined Turbulent Diffusion Flame, Numerical Heat Transfer, Vol. 3, 993, pp of
AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
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 informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationSecond Law Analysis of Forced Convective Cooling in a Channel with a Heated Wall Mounted Obstacle
Journal of Electronics Cooling and Thermal Control, 3, 3, - http://d.doi.org/.436/jectc.3.33 Published Online September 3 (http://www.scirp.org/journal/jectc) Second Law Analysis of Forced Convective Cooling
More informationA NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
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 informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
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 informationFall Exam II. Wed. Nov. 9, 2005
Fall 005 10.34 Eam II. Wed. Nov. 9 005 (SOLUTION) Read through the entire eam before beginning work and budget your time. You are researching drug eluting stents and wish to understand better how the drug
More informationChapter 6 Laminar External Flow
Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than
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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav
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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
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 informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationChapter 8: Flow in Pipes
8-1 Introduction 8-2 Laminar and Turbulent Flows 8-3 The Entrance Region 8-4 Laminar Flow in Pipes 8-5 Turbulent Flow in Pipes 8-6 Fully Developed Pipe Flow 8-7 Minor Losses 8-8 Piping Networks and Pump
More informationApplication of Computational Fluid Dynamics in Discontinuous Unsteady Flow with Large Amplitude Changes; the Shock Tube Problem
Application of Computational Fluid Dynamics in Discontinuous Unsteady Flow with Large Amplitude Changes; the Shock Tube Problem KHALED ALHUSSAN Research Assistant Professor, Space Research Institute, King
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationAdvanced Numerical Methods for non-premixed Flames
0 Advanced Numerical Methods for non-premixed Flames Annarita Viggiano Department of Environmental Engineering and Physics University of Basilicata Italy 1. Introduction Engine designers are under increasing
More informationfluctuations based on the resolved mean flow
Temperature Fluctuation Scaling in Reacting Boundary Layers M. Pino Martín CTR/NASA Ames, Moffett Field, CA 94035 Graham V. Candler Aerospace Engineering and Mechanics University of Minnesota, Minneapolis,
More informationNumerical modeling of a cutting torch
Journal of Physics: Conference Series OPEN ACCESS Numerical modeling of a cutting torch To cite this article: B R Mancinelli et al 2014 J. Phys.: Conf. Ser. 511 012071 View the article online for updates
More information4.2 Concepts of the Boundary Layer Theory
Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
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 informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationA Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM
206 9th International Conference on Developments in esystems Engineering A Computational Investigation of a Turbulent Flow Over a Backward Facing Step with OpenFOAM Hayder Al-Jelawy, Stefan Kaczmarczyk
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
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 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 informationABSTRACT I. INTRODUCTION
2016 IJSRSET Volume 2 Issue 4 Print ISSN : 2395-1990 Online ISSN : 2394-4099 Themed Section: Engineering and Technology Analysis of Compressible Effect in the Flow Metering By Orifice Plate Using Prasanna
More informationADVANCES in NATURAL and APPLIED SCIENCES
ADVANCES in NATURAL and APPLIED SCIENCES ISSN: 1995-0772 Published BYAENSI Publication EISSN: 1998-1090 http://www.aensiweb.com/anas 2017 September 11(11): 1-10 Open Access Journal Two-Dimensional Numerical
More informationPreliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data
35th AIAA Fluid Dynamics Conference, June 6 9, 2005/Toronto,Canada Preliminary Study of the Turbulence Structure in Supersonic Boundary Layers using DNS Data Ellen M. Taylor, M. Pino Martín and Alexander
More informationProcess Chemistry Toolbox - Mixing
Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationLaminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationB.1 NAVIER STOKES EQUATION AND REYNOLDS NUMBER. = UL ν. Re = U ρ f L μ
APPENDIX B FLUID DYNAMICS This section is a brief introduction to fluid dynamics. Historically, a simplified concept of the boundary layer, the unstirred water layer, has been operationally used in the
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 7
Numerical Fluid Mechanics Fall 2011 Lecture 7 REVIEW of Lecture 6 Material covered in class: Differential forms of conservation laws Material Derivative (substantial/total derivative) Conservation of Mass
More informationPredictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical and CFD Method
International Journal of Engineering and Technical Research (IJETR) ISSN: 2321-0869, Volume-3, Issue-5, May 2015 Predictionof discharge coefficient of Venturimeter at low Reynolds numbers by analytical
More informationNUMERICAL INVESTIGATION OF THE TURBULENCE- COMBUSTION INTERACTION IN NON-PREMIXED CH 4 /AIR FLAMES
UMERICAL IVESTIGATIO OF THE TURBULECE- COMBUSTIO ITERACTIO I O-PREMIXED CH 4 /AIR FLAMES K. M. Saqr, M. M. Sies, and M. A. Wahid High-Speed Reacting Flow Laboratory, Faculty of Mechanical Engineering,
More informationHelsinki University of Technology Laboratory of Applied Thermodynamics
Helsinki University of Technology Laboratory of pplied Thermodynamics umerical Methods for Coupling the Reynolds-veraged avier Stokes Equations with the Reynolds-Stress Turbulence Model Patrik Rautaheimo
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
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 informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
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 informationChapter 3. Formulation of FEM for Two-Dimensional Problems
Chapter Formulation of FEM for Two-Dimensional Problems.1 Two-Dimensional FEM Formulation Many details of 1D and 2D formulations are the same. To demonstrate how a 2D formulation works we ll use the following
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationA Comparative Analysis of Turbulent Pipe Flow Using k And k Models
A Comparative Analysis of Turbulent Pipe Flow Using k And k Models 1 B. K. Menge, 2 M. Kinyanjui, 3 J. K. Sigey 1 Department of Mathematics and Physics. Technical University of Mombasa, P.O BOX 90420-80100,Mombasa,
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationFINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS
IJRET: International Journal of Research in Engineering and Technology eissn: 39-63 pissn: 3-738 FINITE-VOLUME SOLUTION OF DIFFUSION EQUATION AND APPLICATION TO MODEL PROBLEMS Asish Mitra Reviewer: Heat
More informationCranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom. Cranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom
MultiScience - XXX. microcad International Multidisciplinary Scientific Conference University of Miskolc, Hungary, 21-22 April 2016, ISBN 978-963-358-113-1 NUMERICAL INVESTIGATION OF AN INCOMPRESSIBLE
More information2Dimensional CFD Simulation of Gas Fired Furnace during Heat Treatment Process
2Dimensional CFD Simulation of Gas Fired Furnace during Heat Treatment Process Vivekanand D. Shikalgar 1, Umesh C.Rajmane 2, Sanjay S. Kumbhar 3, Rupesh R.bidkar 4, Sagar M. Pattar 5, Ashwini P. Patil
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 informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER
ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER Ing. Vojtech Betak Ph.D. Aerospace Research and Test Establishment Department of Engines Prague, Czech Republic Abstract
More informationCoupled calculations in OpenFOAM -
Coupled calculations in OpenFOAM - Multiphysics handling, structures and solvers, Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg Chalmers University of Technology November 14, 2012 Outline
More informationME 144: Heat Transfer Introduction to Convection. J. M. Meyers
ME 144: Heat Transfer Introduction to Convection Introductory Remarks Convection heat transfer differs from diffusion heat transfer in that a bulk fluid motion is present which augments the overall heat
More informationMETHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION
TASKQUARTERLY9No1,17 35 METHOD OF LINES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN THE STREAM-FUNCTION FORMULATION ZBIGNIEW KOSMA Institute of Applied Mechanics, Radom University of Technology,
More informationAdvanced near-wall heat transfer modeling for in-cylinder flows
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress April 20, 2015 Detroit, MI S. Šarić, B. Basara AVL List GmbH Advanced near-wall heat transfer modeling for in-cylinder
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationEntropy 2011, 13, ; doi: /e OPEN ACCESS. Entropy Generation at Natural Convection in an Inclined Rectangular Cavity
Entropy 011, 13, 100-1033; doi:10.3390/e1305100 OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Article Entropy Generation at Natural Convection in an Inclined Rectangular Cavity Mounir
More informationSolution Methods. Steady convection-diffusion equation. Lecture 05
Solution Methods Steady convection-diffusion equation Lecture 05 1 Navier-Stokes equation Suggested reading: Gauss divergence theorem Integral form The key step of the finite volume method is to integrate
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 informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
More informationValidation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
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 informationChapter 6 An introduction of turbulent boundary layer
Chapter 6 An introduction of turbulent boundary layer T-S Leu May. 23, 2018 Chapter 6: An introduction of turbulent boundary layer Reading assignments: 1. White, F. M., Viscous fluid flow. McGraw-Hill,
More informationAnalysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data
45th AIAA Aerospace Science Meeting and Exhibit, January 8 11, 25/Reno, Nevada Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data M. Wu and M.P. Martin Mechanical
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationNumerical Modeling of Anode Baking Process with COMSOL Multiphysics
Numerical Modeling of Anode Baking Process with COMSOL Multiphysics P. Nakate 1,2, W. Voogd de 3, D.J.P. Lahaye 1, C. Vuik 1, M. Talice 4 1. Department of applied mathematics, Delft University of technology,
More informationSYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL
ISSN 1726-4529 Int j simul model 7 (2008) 4, 186-197 Original scientific paper SYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL Ternik, P. University
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER
THERMAL SCIENCE: Year 2018, Vol. 22, No. 2, pp. 963-972 963 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF A V-RIB WITH GAP ROUGHENED SOLAR AIR HEATER by Jitesh RANA, Anshuman SILORI, Rajesh MAITHANI *, and
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationChapter 8: Flow in Pipes
Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate the major and minor losses associated with pipe flow in piping networks
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationEnhancement of the momentum interpolation method on non-staggered grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 33: 1 22 Enhancement of the momentum interpolation method on non-staggered grids J. Papageorgakopoulos, G. Arampatzis,
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationStabilized Formulations and Smagorinsky Turbulence Model for Incompressible Flows
WCCM V Fifth World Congress on Computational Mechanics July 7 2, 22, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Stabilized Formulations and Smagorinsky Turbulence Model for
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationViscous Flow in Ducts
Dr. M. Siavashi Iran University of Science and Technology Spring 2014 Objectives 1. Have a deeper understanding of laminar and turbulent flow in pipes and the analysis of fully developed flow 2. Calculate
More information