FEMLAB Exercise 1 for ChE366
|
|
- Candice Mason
- 5 years ago
- Views:
Transcription
1 FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let η and ρ denote the viscosity and density of the fluid respectively. The center of the sphere always lies on the axis of the cylinder. Compute the velocity and pressure fields by using FEMLAB for the following parameters: r s = 0.2 mm, U = 1 mm/s, and η = 10 cp and ρ = 1 g/cm 3 for R values (in mm) of 2, 1, 0.6, 0.4, 0.3 and Compare the results obtained for different values of κ r s /R. Show that for κ >>1, friction coefficient ζ F/(6πηUr s ), where F denotes the drag force, approaches 1. Construct a plot of ζ vs. κ. Key Steps Step 1: Defining the geometry (Simulating an infinitely long cylinder ) In practice, we need to consider a finite cylinder with length L >> r s such that the end effects do not influence our final answer. This requires a trial and error procedure as given below. 1. Choose a length L >>r s, say L = 10 r s. 2. Solve the problem and check whether the solution (e.g. velocity and pressure) changes appreciably at the ends of the cylinder. 3. If the solution changes appreciably at the cylinder ends, increase L by 10 percent and go back to step 2. If not, you have succeeded in simulating an infinite cylinder. 1 Prepared by Kartik Arora based on guidance from Prof. R. Sureshkumar Privileged material. The use may be subject to software licensing agreements. Do not circulate or modify without written permission from Sureshkumar. 1
2 You would require boundary conditions at the cylinder ends. These boundary conditions can be assigned in the form of pressure, i.e., maintaining the same pressure at the two ends of the cylinder: see the attached handout and figure 1 in this document. Step 2: Changing the frame of reference (Simulating a translating sphere) It is not possible to solve a problem with moving boundaries in FEMLAB. Therefore, it makes sense to solve the problem in a new of frame of reference in which the sphere is stationary and the cylindrical wall moves with a constant velocity U. Step 3: Exploiting symmetries to reduce problem size Cylindrical symmetry can be used to reduce the problem geometry to an axisymmetric one, see figure 4 of the attached handout and figure below. p = p a Line of symmetry r s No slip z r L u z = U u r = 0 p = p a Figure 1: Schematic of the flow geometry along with the boundary conditions R Step 4: Non-dimensionalization of the Navier-Stokes equations The fluid motion inside the cylinder can be described by Navier-Stokes equations which in its dimensional form can be written as follows: 2
3 u ρ + u u = p + η t 2 r r u + ρg + F. (1) In the above equation, u is the velocity vector, p is the hydrodynamic pressure, g r is the acceleration due to gravity, and F r is the external force applied per unit volume on the fluid. For this problem, we neglect gravity and assume no external force, i.e. F r = 0. The equation is solved along with the incompressibility constraint (i.e., conservation of mass or the continuity equation) on the velocity vector, which is, u = 0. (2) Equations (1) and (2) can be non-dimensionalized by using quantities, r s, U, r s and U the scales for length, velocity, time and pressure respectively. Let u, t and p denote the dimensionless velocity vector, time and pressure respectively. Then equation (1) can be written in the dimensionless form as (verify this): ρu r s 2 u ηu ηu + u u = p + t r r 2 s 2 s 2 u. Note that in the above equation, the gradient and the Laplacian operators are in dimensionless forms. Dividing equation (3) by ηu 2, we obtain r s ρur s η u + u u = p + 2 u. t The quantity ρur s appearing as the coefficient of the left hand side term, is referred to as the η Reynolds number and is denoted by Re, which simplifies equation (4) as, u 2 Re + u u = p + u. t (5) Recall that Re represents the ratio of inertial to viscous forces in the flow. Equation (5) looks similar to equation (1), except that it is in dimensionless form, the coefficient of the Laplacian of u is 1 instead of η, and the coefficient of the convection term in the l.h.s. is Re instead of ρ. Note that we have ignored the gravitational and the external force terms in equation (1). ηu (3) (4) r s as 3
4 FEMLAB only allows equations to be written in dimensional form. However, by cleverly renaming the variables, a non-dimensional equation (e.g. Eq. 5) can be represented in FEMLAB. Specifically, this can be achieved in our example if we assign Re to ρ, keep η as 1, and define the geometry in non-dimensional coordinates. Note that in non-dimensional coordinates, the radius of the sphere will always be 1, length of the cylinder will be L/r s and radius of the cylinder will be R/r s. The dimensionless geometry and the dimensionless boundary conditions are shown in figure 2 below. Note that the non-dimensional form allows us to determine the flow regime, e.g. Re << 1 implies creeping flow. Step 5: Calculating Drag The drag on the sphere can be calculated by integrating the z-component of the force along the surface of the sphere. You can compute this quantity by using the boundary integration tool in FEMLAB (see attached handout). p = 0 Line of symmetry 1 No slip z r L/r s u z = 1 u r = 0 p = 0 Figure 2: Schematic of the dimensionless flow geometry with the boundary conditions (dimensionless coordinates). R/r s 4
5 Notice that when R/r s >> 1, the particle will not feel the presence of the cylinder wall. In this case, the drag on the sphere can be compared with the Stokes law valid in the limit Re 0. According to Stokes law, the force on a spherical particle moving with a velocity U in an infinite medium, is given by, F = 6πr s Uη. (6) This force can be expressed in its dimensionless form as, F = 6π. (7) You can compare the force evaluated using FEMLAB to the equation (7). You can also study the influence of R/r s and Re on the force (drag). Theoretically, when Re << 1 and R/r s >> 1, the drag on the particle has to be equal to 6π. 5
6 Attachment: Handout on the use of FEMLAB FEMLAB utilizes the Finite Element Method in the MATLAB interface to solve physical and mathematical problems. It has a user-friendly interface that allows to prescribe process geometry and built-in models for processes from the library of chemical/mechanical/electrical processes. You can also create and customize applications in FEMLAB. The steps involved in solving a problem in FEMLAB. 1. Choose the physics: The moment you open the FEMLAB software, it will ask you to choose a process from its built-in library. Take some time to browse through the library; see what all processes you can identify. You also need to choose the dimensionality of the problem, i.e., whether it is 1-dimensional, 2-dimensional, 3-dimensional, 1-d axisymmetric, or 2-d axisymmetric. You can always go back and change the physics of the problem, but changing the dimensionality will not be easy. In some cases you might have to start over the whole problem. For our example, choose 2-d axisymmetric geometry + Incompressible Navier-Stokes equation as shown in figure 1. Figure 1. Choosing the physics (process). Note that in the Space dimension field, it shows Axial symmetry (2D). 6
7 2. Specify the details of the process geometry: Once you have chosen the physics, depending on the dimensionality, FEMLAB will bring up a 1-d, 2-d, or 3-d interface as shown in figure 2. At this point you are required to define the geometry, i.e., the domain of the process. There are various preset geometries in FEMLAB that you can use to draw the domain, or you can draw an arbitrary shape by using the 2 nd and 3 rd degree brazier curve object in menu Draw Draw Objects. Also, the toolbar on the left has all the short-cut keys to these objects. Figure 2. The axisymmetric interface. Note that the vertical and horizontal axes represent the z and the r directions respectively. For our example, draw a rectangle with its vertices at (0, -20), (0, 20), (20, 20), (20, -20), and a circle with its center at (0,0) and radius of 1 2. After drawing these objects, select them by keeping the shift key pressed and clicking the objects individually. Once the objects are selected, press the difference button on the left panel/toolbar (one that has two intersecting circles with the non-intersected region of one of the circles colored) as shown in figure 3. This should result in the desired axisymmetric geometry for our problem, as shown in figure 4. 2 You can view and change the properties of an object by double clicking it. 7
8 Figure 3. Taking the difference of two objects to construct the geometry. Figure 4. Desired Geometry: Rotating along the vertical axis gives the 3-D geometry for sphere in a cylindrical vessel. 3. Assign the parameters of the problem: The Navier-Stokes equation, in its dimensional form, contains two parameters, η and ρ, denoting the viscosity and the density of the fluid respectively. However, we will only work with non-dimensional quantities, as explained in the 8
9 problem definition. In this case, η will always be 1 and ρ will represent the Reynolds number, as FEMLAB doesn t allow equations in their non-dimensional forms. Refer to the attached document (pgs 2-3) for details. The parameter, µ and ρ, can be assigned by opening menu physics Subdomain Setting As a test problem run the following case, enter 0 in the field for ρ, and 1 in the field for η, as shown in figure 5. Figure 5. Assigning the parameters. 4. Specify the boundary conditions: Once you are sure that you have assigned the right values to all the parameters, it is time to specify the boundary conditions. It is important to note that N-S equation is a second order PDE; hence it requires either the flux (derivative) specification, or the variable(s) specification all around the boundary. The boundary conditions can be specified by opening menu physics Boundary Settings 9
10 For our geometry, FEMLAB divides the boundary into 7 different segments, which consist of 5 lines (black) and two quarter-circles (red), as shown in figure 6. You can change the conditions on each of the boundary segments by highlighting its number on the left side, i.e., Boundary selection field. For segments 1 and 3, in the boundary condition field, choose slip/symmetry. For segments 2 and 4, in the boundary condition field, choose normal flow/pressure and let p 0 = 0. For segment 5, in the boundary condition field, choose inflow/outflow velocity and set z-velocity to 1 and r-velocity to 0. For segments 6 and 7, choose the no slip field as shown in figure 6. Figure 6. Boundary segments: There are seven different segments; condition on each one of them can be specified individually. 5. Construct a mesh: Finite element method divides the domain into a finite number of elements, which can in principle have arbitrary shapes. However, in most of the problems, it is convenient to have them as polygons, i.e., triangles, rectangles, pentagons etc. FEMLAB employs only triangular elements. You can draw the mesh (a mesh consists of several elements) by choosing the following option menu Mesh Initialize mesh. 10
11 After initializing the mesh, the FEMLAB will draw all the elements as shown in figure 7. You can further refine the mesh in the whole domain or you can focus on a desired region, and refine that part only. Do not refine the mesh unnecessarily because your system can run out of memory. Besides, it takes more time for finer meshes. However, you will need to ensure that your results are mesh-independent or mesh-converged. Figure 7. Mesh generation. If N is the number of nodes and there are 3 dependent variables (u, v, and p), then FEMLAB calculates the number of degree of freedom as 3N. 6. Solving the problem: First go to menu Solver Parameter, and make sure that in the General tab, your solver is set as Stationary nonlinear. Do not change any other fields at this point. Congratulations! You have completed the problem definition; go ahead and choose menu solve Solve Problem. For the initialized mesh, this should take 2-5 seconds to run on CELL computers, unless you are running some other programs in background. 11
12 7. Understanding and presenting the solution: By default, FEMLAB will display the surface plot (color-coded) of the absolute velocity. You can plot streamlines and vector plots of the flow, or you can plot contour plots of individual velocities (u r or u z ), by going to menu Post processing Plot Parameters as shown in figure 8. Figure 8. Displaying the solution. The figure illustrates the surface plot (color-coded) of the absolute velocity. You can also plot the solution (variables or derived quantities) along one of domain s boundaries by going to menu Post processing Domain Plot Parameters: Tab = Line/Extrusion, and selecting the desired boundary and variable. Or you can draw a line plot by clicking on the lineplot short cut in the left toolbar, as shown in figure 9. 12
13 Short-cut for line plots Figure 9. Constructing the line plot of the absolute velocity. 8. Calculating the drag on the sphere: The total drag, or in other words the force exerted by the fluid, on the sphere will be in the positive z-direction. The net force in the z-direction can be evaluated by performing a boundary integral of Total force per area, z-component, on the sphere s perimeter, i.e., segments 6 and 7. You can do this by going to menu Post Processing Boundary Integration. The two segments should be selected simultaneously in order to calculate the integral. 13
14 Remember to check this box Result for the surface integral Figure 10. Calculating the drag on the sphere. For our problem, remember to check the box below the expression fields in the Boundary Integration window, that reads Compute surface integral (for axisymmetric modes), as shown in figure 10. The results appear only after clicking the OK button, which shows up in the log generated below the figure. Compare the drag with equation (7) of the handout FEMLAB Exercise for ChE
Strömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationTerminal Falling Velocity of a Sand Grain
Terminal Falling Velocity of a Sand Grain Introduction The first stop for polluted water entering a water work is normally a large tank, where large particles are left to settle. More generally, gravity
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationSteady Flow: Lid-Driven Cavity Flow
STAR-CCM+ User Guide Steady Flow: Lid-Driven Cavity Flow 2 Steady Flow: Lid-Driven Cavity Flow This tutorial demonstrates the performance of STAR-CCM+ in solving a traditional square lid-driven cavity
More informationAppendix: To be performed during the lab session
Appendix: To be performed during the lab session Flow over a Cylinder Two Dimensional Case Using ANSYS Workbench Simple Mesh Latest revision: September 18, 2014 The primary objective of this Tutorial is
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationVerification of Laminar and Validation of Turbulent Pipe Flows
1 Verification of Laminar and Validation of Turbulent Pipe Flows 1. Purpose ME:5160 Intermediate Mechanics of Fluids CFD LAB 1 (ANSYS 18.1; Last Updated: Aug. 1, 2017) By Timur Dogan, Michael Conger, Dong-Hwan
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationSolved with COMSOL Multiphysics 4.0a. COPYRIGHT 2010 COMSOL AB.
Journal Bearing Introduction Journal bearings are used to carry radial loads, for example, to support a rotating shaft. A simple journal bearing consists of two rigid cylinders. The outer cylinder (bearing)
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationThe viscous forces on the cylinder are proportional to the gradient of the velocity field at the
Fluid Dynamics Models : Flow Past a Cylinder Flow Past a Cylinder Introduction The flow of fluid behind a blunt body such as an automobile is difficult to compute due to the unsteady flows. The wake behind
More informationModule D: Laminar Flow over a Flat Plate
Module D: Laminar Flow over a Flat Plate Summary... Problem Statement Geometry and Mesh Creation Problem Setup Solution. Results Validation......... Mesh Refinement.. Summary This ANSYS FLUENT tutorial
More informationMicrowell Mixing with Surface Tension
Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical
More informationCHAPTER 3. Elementary Fluid Dynamics
CHAPTER 3. Elementary Fluid Dynamics - Understanding the physics of fluid in motion - Derivation of the Bernoulli equation from Newton s second law Basic Assumptions of fluid stream, unless a specific
More informationFLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016
FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
More informationSimulation of Turbulent Flow over the Ahmed Body
1 Simulation of Turbulent Flow over the Ahmed Body ME:5160 Intermediate Mechanics of Fluids CFD LAB 4 (ANSYS 18.1; Last Updated: Aug. 18, 2016) By Timur Dogan, Michael Conger, Dong-Hwan Kim, Maysam Mousaviraad,
More informationSupersonic Flow Over a Wedge
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 2 Supersonic Flow Over a Wedge Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification A uniform supersonic stream encounters a wedge
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
1 Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose ME:5160 Intermediate Mechanics of Fluids CFD LAB 2 (ANSYS 19.1; Last Updated: Aug. 7, 2018) By Timur Dogan, Michael Conger,
More informationTUTORIAL#3. Marek Jaszczur. Boundary Layer on a Flat Plate W1-1 AGH 2018/2019
TUTORIAL#3 Boundary Layer on a Flat Plate Marek Jaszczur AGH 2018/2019 W1-1 Problem specification TUTORIAL#3 Boundary Layer - on a flat plate Goal: Solution for boudary layer 1. Creating 2D simple geometry
More informationVelocity and Concentration Properties of Porous Medium in a Microfluidic Device
Velocity and Concentration Properties of Porous Medium in a Microfluidic Device Rachel Freeman Department of Chemical Engineering University of Washington ChemE 499 Undergraduate Research December 14,
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationSimulation of Turbulent Flow in an Asymmetric Diffuser
Simulation of Turbulent Flow in an Asymmetric Diffuser 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 3 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University of Iowa C.
More informationSTAR-CCM+: Wind loading on buildings SPRING 2018
STAR-CCM+: Wind loading on buildings SPRING 2018 1. Notes on the software 2. Assigned exercise (submission via Blackboard; deadline: Thursday Week 3, 11 pm) 1. NOTES ON THE SOFTWARE STAR-CCM+ generates
More informationSimulation of Laminar Pipe Flows
Simulation of Laminar Pipe Flows 57:020 Mechanics of Fluids and Transport Processes CFD PRELAB 1 By Timur Dogan, Michael Conger, Maysam Mousaviraad, Tao Xing and Fred Stern IIHR-Hydroscience & Engineering
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationAnalysis of Flow Dynamics of an Incompressible Viscous Fluid in a Channel
Analysis of Flow Dynamics of an Incompressible Viscous Fluid in a Channel Deepak Kumar Assistant Professor, Department of Mechanical Engineering, Amity University Gurgaon, India E-mail: deepak209476@gmail.com
More informationIsotropic Porous Media Tutorial
STAR-CCM+ User Guide 3927 Isotropic Porous Media Tutorial This tutorial models flow through the catalyst geometry described in the introductory section. In the porous region, the theoretical pressure drop
More informationAn Introduction to SolidWorks Flow Simulation 2010
An Introduction to SolidWorks Flow Simulation 2010 John E. Matsson, Ph.D. SDC PUBLICATIONS www.sdcpublications.com Schroff Development Corporation Chapter 2 Flat Plate Boundary Layer Objectives Creating
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationMiddle East Technical University Mechanical Engineering Department ME 485 CFD with Finite Volume Method Fall 2017 (Dr. Sert)
Middle East Technical University Mechanical Engineering Department ME 485 CFD with Finite Volume Method Fall 2017 (Dr. Sert) ANSYS Fluent Tutorial Developing Laminar Flow in a 2D Channel 1 How to use This
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More informationENERGY-224 Reservoir Simulation Project Report. Ala Alzayer
ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................
More informationUsing Periodic Boundary Conditions
1 of 6 2004 11 08 15:20 Use periodic boundary conditions, periodic edge conditions, and periodic point conditions to define a constraint that makes two quantities equal on two different (but usually equally
More informationSolidWorks Flow Simulation 2014
An Introduction to SolidWorks Flow Simulation 2014 John E. Matsson, Ph.D. SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following websites
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY. Analyzing wind flow around the square plate using ADINA Project. Ankur Bajoria
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Analyzing wind flow around the square plate using ADINA 2.094 - Project Ankur Bajoria May 1, 2008 Acknowledgement I would like to thank ADINA R & D, Inc for the full
More informationSolved with COMSOL Multiphysics 4.2
Peristaltic Pump Solved with COMSOL Multiphysics 4.2 Introduction In a peristaltic pump, rotating rollers squeeze a flexible tube. As the pushed-down rollers move along the tube, fluids in the tube follow
More informationNumerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers
Numerical studies for Flow Around a Sphere regarding different flow regimes caused by various Reynolds numbers R. Jendrny, H. Damanik, O. Mierka, S. Turek Institute of Applied Mathematics (LS III), TU
More informationANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step
ANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Governing Equation Start-Up Geometry
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
More informationModeling and Simulation of Single Phase Fluid Flow and Heat Transfer in Packed Beds
Modeling and Simulation of Single Phase Fluid Flow and Heat Transfer in Packed Beds by:- Balaaji Mahadevan Shaurya Sachdev Subhanshu Pareek Amol Deshpande Birla Institute of Technology and Science, Pilani
More informationSolution Recording and Playback: Vortex Shedding
STAR-CCM+ User Guide 6663 Solution Recording and Playback: Vortex Shedding This tutorial demonstrates how to use the solution recording and playback module for capturing the results of transient phenomena.
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationand to the following students who assisted in the creation of the Fluid Dynamics tutorials:
Fluid Dynamics CAx Tutorial: Pressure Along a Streamline Basic Tutorial #3 Deryl O. Snyder C. Greg Jensen Brigham Young University Provo, UT 84602 Special thanks to: PACE, Fluent, UGS Solutions, Altair
More informationTechnical Report TR
Technical Report TR-2015-09 Boundary condition enforcing methods for smoothed particle hydrodynamics Arman Pazouki 1, Baofang Song 2, Dan Negrut 1 1 University of Wisconsin-Madison, Madison, WI, 53706-1572,
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationArbitrary Lagrangian-Eulerian method for coupled Navier-Stokes and convection-diffusion equations with moving boundaries
1th WSEAS Int. Conf. on APPLIED MATHEMATICS, Cairo, Egypt, December 9-31, 007 31 Arbitrary Lagrangian-Eulerian method for coupled Navier-Stokes and convection-diffusion equations with moving boundaries
More informationMEAM 550 Modeling and Design of MEMS Spring Solution to homework #3. In our notation and values, k = = =
MEAM 550 Modeling and Design of MEMS Spring 004 Solution to homework # Problem 1 A fixed-guided beam (length = l, width = b, depth = h ) with a transverse tip load of F has the following formulas for maximum
More informationSimulation of Turbulent Flow around an Airfoil
Simulation of Turbulent Flow around an Airfoil ENGR:2510 Mechanics of Fluids and Transfer Processes CFD Pre-Lab 2 (ANSYS 17.1; Last Updated: Nov. 7, 2016) By Timur Dogan, Michael Conger, Andrew Opyd, Dong-Hwan
More informationSTAR-CCM+ User Guide 6922
STAR-CCM+ User Guide 6922 Introduction Welcome to the STAR-CCM+ introductory tutorial. In this tutorial, you explore the important concepts and workflow. Complete this tutorial before attempting any others.
More information3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic
ISCS 2013 Selected Papers Dam-break effect on a Solid Wall 1 3D Simulation of Dam-break effect on a Solid Wall using Smoothed Particle Hydrodynamic Suprijadi a,b, F. Faizal b, C.F. Naa a and A.Trisnawan
More informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
More informationAIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN ABSTRACT
AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN Adeel Khalid *, Daniel P. Schrage + School of Aerospace Engineering, Georgia Institute of Technology
More informationFinite Volume Discretization on Irregular Voronoi Grids
Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling
More informationSolved with COMSOL Multiphysics 4.2
Backstep Introduction This tutorial model solves the incompressible Navier-Stokes equations in a backstep geometry. A characteristic feature of fluid flow in geometries of this kind is the recirculation
More informationA Comparative CFD Analysis of a Journal Bearing with a Microgroove on the Shaft & Journal
Proceedings of International Conference on Innovation & Research in Technology for Sustainable Development (ICIRT 2012), 01-03 November 2012 182 A Comparative CFD Analysis of a Journal Bearing with a Microgroove
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE
More informationNumerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert
Numerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert A. Kadyyrov 1 1 Research center for power engineering problems Federal government budgetary institution of science Kazan
More informationand to the following students who assisted in the creation of the Fluid Dynamics tutorials:
Fluid Dynamics CAx Tutorial: Channel Flow Basic Tutorial # 4 Deryl O. Snyder C. Greg Jensen Brigham Young University Provo, UT 84602 Special thanks to: PACE, Fluent, UGS Solutions, Altair Engineering;
More informationImplementation in COMSOL
Implementation in COMSOL The transient Navier-Stoke equation will be solved in COMSOL. A text (.txt) file needs to be created that contains the velocity at the inlet of the common carotid (calculated as
More informationSolved with COMSOL Multiphysics 4.2
Laminar Static Mixer Introduction In static mixers, also called motionless or in-line mixers, a fluid is pumped through a pipe containing stationary blades. This mixing technique is particularly well suited
More informationCFD Post-Processing of Rampressor Rotor Compressor
Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements
More informationA B C D E. Settings Choose height, H, free stream velocity, U, and fluid (dynamic viscosity and density ) so that: Reynolds number
Individual task Objective To derive the drag coefficient for a 2D object, defined as where D (N/m) is the aerodynamic drag force (per unit length in the third direction) acting on the object. The object
More informationTutorial 2: Particles convected with the flow along a curved pipe.
Tutorial 2: Particles convected with the flow along a curved pipe. Part 1: Creating an elbow In part 1 of this tutorial, you will create a model of a 90 elbow featuring a long horizontal inlet and a short
More informationANSYS FLUENT. Airfoil Analysis and Tutorial
ANSYS FLUENT Airfoil Analysis and Tutorial ENGR083: Fluid Mechanics II Terry Yu 5/11/2017 Abstract The NACA 0012 airfoil was one of the earliest airfoils created. Its mathematically simple shape and age
More informationSimulation of Turbulent Flow around an Airfoil
1. Purpose Simulation of Turbulent Flow around an Airfoil ENGR:2510 Mechanics of Fluids and Transfer Processes CFD Lab 2 (ANSYS 17.1; Last Updated: Nov. 7, 2016) By Timur Dogan, Michael Conger, Andrew
More informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationANSYS AIM Tutorial Steady Flow Past a Cylinder
ANSYS AIM Tutorial Steady Flow Past a Cylinder Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Solution Domain Boundary Conditions Start-Up
More informationANSYS AIM Tutorial Compressible Flow in a Nozzle
ANSYS AIM Tutorial Compressible Flow in a Nozzle Author(s): Sebastian Vecchi Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Pre-Analysis Start-Up Geometry Import Geometry Mesh
More informationDemo problem: Steady finite-reynolds-number flow through an iliac bifurcation
Chapter 1 Demo problem: Steady finite-reynolds-number flow through an iliac bifurcation The purpose of this tutorial is to demonstrate the simulation of cardiovascular fluid mechanics problems with oomph-lib.
More informationFree Convection Cookbook for StarCCM+
ME 448/548 February 28, 2012 Free Convection Cookbook for StarCCM+ Gerald Recktenwald gerry@me.pdx.edu 1 Overview Figure 1 depicts a two-dimensional fluid domain bounded by a cylinder of diameter D. Inside
More informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
More informationSPC 307 Aerodynamics. Lecture 1. February 10, 2018
SPC 307 Aerodynamics Lecture 1 February 10, 2018 Sep. 18, 2016 1 Course Materials drahmednagib.com 2 COURSE OUTLINE Introduction to Aerodynamics Review on the Fundamentals of Fluid Mechanics Euler and
More informationSETTLEMENT OF A CIRCULAR FOOTING ON SAND
1 SETTLEMENT OF A CIRCULAR FOOTING ON SAND In this chapter a first application is considered, namely the settlement of a circular foundation footing on sand. This is the first step in becoming familiar
More informationFirst Steps - Ball Valve Design
COSMOSFloWorks 2004 Tutorial 1 First Steps - Ball Valve Design This First Steps tutorial covers the flow of water through a ball valve assembly before and after some design changes. The objective is to
More informationµ = Pa s m 3 The Reynolds number based on hydraulic diameter, D h = 2W h/(w + h) = 3.2 mm for the main inlet duct is = 359
Laminar Mixer Tutorial for STAR-CCM+ ME 448/548 March 30, 2014 Gerald Recktenwald gerry@pdx.edu 1 Overview Imagine that you are part of a team developing a medical diagnostic device. The device has a millimeter
More informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationComputational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent
MEGR 7090-003, Computational Fluid Dynamics :1 7 Spring 2015 Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent Rahul R Upadhyay Master of Science, Dept of Mechanical
More informationThe Transient Modeling of Bubble Pinch-Off Using an ALE Moving Mesh
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston The Transient Modeling of Bubble Pinch-Off Using an ALE Moving Mesh Christopher J. Forster, Marc K. Smith * School of Mechanical Engineering,
More informationA 3D VOF model in cylindrical coordinates
A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved
More informationExcerpt from the Proceedings of the COMSOL Conference 2008 Boston
Excerpt from the Proceedings of the COMSOL Conference 008 Boston CFD-based Evaluation of Drag Force on a Sphere Unsteadily Moving Perpendicularly toward a Solid Surface: a Simple Model of a Biological
More informationSimLab 14.3 Release Notes
SimLab 14.3 Release Notes Highlights SimLab 14.0 introduced new graphical user interface and since then this has evolved continuously in subsequent versions. In addition, many new core features have been
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 Outline
More informationChapter 24. Creating Surfaces for Displaying and Reporting Data
Chapter 24. Creating Surfaces for Displaying and Reporting Data FLUENT allows you to select portions of the domain to be used for visualizing the flow field. The domain portions are called surfaces, and
More informationTUTORIAL#4. Marek Jaszczur. Turbulent Thermal Boundary Layer on a Flat Plate W1-1 AGH 2018/2019
TUTORIAL#4 Turbulent Thermal Boundary Layer on a Flat Plate Marek Jaszczur AGH 2018/2019 W1-1 Problem specification TUTORIAL#4 Turbulent Thermal Boundary Layer - on a flat plate Goal: Solution for Non-isothermal
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationStep 1: Create Geometry in GAMBIT
Step 1: Create Geometry in GAMBIT If you would prefer to skip the mesh generation steps, you can create a working directory (see below), download the mesh from here (right click and save as pipe.msh) into
More informationSimulation of Turbulent Flow over the Ahmed Body
Simulation of Turbulent Flow over the Ahmed Body 58:160 Intermediate Mechanics of Fluids CFD LAB 4 By Timur K. Dogan, Michael Conger, Maysam Mousaviraad, and Fred Stern IIHR-Hydroscience & Engineering
More informationAnswer Key: Three-Dimensional Cross Sections
Geometry A Unit Answer Key: Three-Dimensional Cross Sections Name Date Objectives In this lesson, you will: visualize three-dimensional objects from different perspectives be able to create a projection
More information0 Graphical Analysis Use of Excel
Lab 0 Graphical Analysis Use of Excel What You Need To Know: This lab is to familiarize you with the graphing ability of excels. You will be plotting data set, curve fitting and using error bars on the
More information