The Control-Volume Finite-Difference Approximation to the Diffusion Equation
|
|
- Patricia Fox
- 5 years ago
- Views:
Transcription
1 The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering ME 448/548: D Convection-Diffusion Equation Overview Develop a numerical model of the two-dimensional Poisson equation Demonstrate application of the finite volume method for a cartesian mesh Allow for non-uniform mesh, diffusion coefficient, and source term Implement in Matlab Demonstrate truncation error for a model problem with a simple solution Appl to full-developed flow in a rectangle ME 448/548: D Convection-Diffusion Equation page
2 Model Problem The two dimensional diffusion equation in Cartesian coordinates is Γ φ + Γ φ + S = () where φ is the scalar field, Γ is the diffusion coefficient, and S is the source term ME 448/548: D Convection-Diffusion Equation page D Cartesian Finite Volume Mesh 3 Interior node + + Boundar node Ambiguous corner node j = i = 3 ME 448/548: D Convection-Diffusion Equation page 3
3 Finite Volume Mesh Notation for a tpical control volume w e Regardless of whether the control volumes sizes are uniform, the node is alwas located in the geometric center of each control volume Thus, δ n W N P E n P w = e P = P s = n P = δ s s δ w S δ e ME 448/548: D Convection-Diffusion Equation page 4 Convert the Differential Equation to a Discrete Equation Integrate over the control volume w e n e s w = Γ φ d d n s Γ φ Γ φ d e w Γ φ Γ φ e w δ n δ s W δ w N P S δ e E n s Γ e φ E φ P δ e φ P φ W Γ w δ w ME 448/548: D Convection-Diffusion Equation page 5
4 Convert the Differential Equation to a Discrete Equation Integrate the source term e n w s Sdd S P () Note: The Control Volume Finite Difference (CVFD) method treats the source term and diffusion coefficients as piecewise constants This is a rather crude approimation, sa, compared to allowing the source term and diffusion coefficient to var linearl within the control volume However, piecewise constant profiles of S and Γ allow the method to be conservative, ie, conserving mass or energ, automaticall The conservative nature of the CVFD method is one of its primar strengths ME 448/548: D Convection-Diffusion Equation page 6 Convert the Differential Equation to a Discrete Equation Putting pieces back into the model equation gives the discrete sstem of equations where a E = a S φ S a W φ W + a P φ P a E φ E a N φ N = b (3) Γ e δ e, a W = Γ w δ w, a N = Γ n δ n, a S = Γ s δ s a P = a E + a W + a N + a S (4) b = S P (5) This is not a tridiagonal sstem of equations ME 448/548: D Convection-Diffusion Equation page 7
5 Non-uniform Γ Continuit of flues at the interface requires φ Γ P φ = Γ E φ = Γ e e e+ Use central difference approimations e material material P E Γ e φ E φ P δ e = Γ P φ e φ P δ e (6) δ e- δ e+ Γ e φ E φ P δ e = Γ E φ E φ e δ e+ (7) δ e ME 448/548: D Convection-Diffusion Equation page 8 Non-uniform Γ Equations 6 and 7 can be rearranged as material material φ e φ P = δ e Γ P φ E φ e = δ e+ Γ E Γ e δ e (φ E φ P ) (8) Γ e δ e (φ E φ P ) (9) P E Add Equation 8 and Equation 9 φ E φ P = Γ e δe (φ E φ P ) + δ e+ δ e Γ P Γ E Cancel the factor of (φ E φ P ) and solve for Γ e /δ e to get δ e- δ e δ e+ Γ e δe = + δ e+ = δ e Γ P Γ E Γ E Γ P δ e Γ E + δ e+ Γ P ME 448/548: D Convection-Diffusion Equation page 9
6 Non-uniform Γ Thus, the diffusion coefficient at the interface that results in flu continuit is Γ e = Γ E Γ P βγ E + ( β)γ P () where β δ e δ e = e P E P () An analogous derivation gives formulas for Γ w, Γ n,andγ s ME 448/548: D Convection-Diffusion Equation page Solving the Sstem of Equations Regardless of uniform or variable Γ, the discrete equation has a five-point stencil, and the discrete equation for an interior node can be written a S φ S a W φ W + a P φ P a E φ E a N φ N = b () To set up the matri for this sstem of equations, we need to re-number the unknowns ME 448/548: D Convection-Diffusion Equation page
7 Solving the Sstem of Equations Order the nodes: n = i +(j )n With natural ordering the neighbors in the compass point notation have these indices np = i + (j-)*n ne = np + nw = np - nn = np + n ns = np - n j = i = 3 Interior node Boundar node Ambiguous corner node ME 448/548: D Convection-Diffusion Equation page Solving the Sstem of Equations This leads to a vector of unknowns φ, φ, φ n, φ, φ, φ i,j φ n,n φ φ φ n φ n+ φ n+ φ n φ N (3) ME 448/548: D Convection-Diffusion Equation page 3
8 Solving the Sstem of Equations A = a S a W ap a E a N ME 448/548: D Convection-Diffusion Equation page 4 Model Problems Choose the source term so that the eact solution is simple to compute Uniform source term Analtical solution eists, but requires an infinite series 3 Nonuniform diffusivit and concentrated source term This problem is meant to stress the numerical solution code ME 448/548: D Convection-Diffusion Equation page 5
9 Model Problem Choose a source term that ma be phsicall unrealistic, but one that gives an eact solution that is eas to evaluate π π π π S = + sin sin L L L L The eact solution is π π φ = sin sin L L ME 448/548: D Convection-Diffusion Equation page 6 Model Problem The eact solution is π π φ = sin sin L L ME 448/548: D Convection-Diffusion Equation page 7
10 Solutions to Model Problem Show Correct Truncation Error The local truncation error at each node is e i O( ) Since the eact solution is known we can compute e N = e i N Nē N = ē N where N = n n is the total number of interior nodes in the domain, and ē is the average truncation error per node Measured error Measured Theoretical error ~ ( Since e i O( ), N n,and = L /(n + ), wecanestimate e N ē = O( ) = O L /(n + ) 3 O = O( 3 ) N n n n ME 448/548: D Convection-Diffusion Equation page 8 Model Problem Uniform source term: S = ME 448/548: D Convection-Diffusion Equation page 9
11 Model Problem 3 Solution with α = Γ = αγ Γ = L 5L S = S = 5L L 5L L ME 448/548: D Convection-Diffusion Equation page Model Problem 4: Full Developed Flow in a Rectangular Duct For simple full-developed flow the governing equation for the aial velocit w is µ w + w dp dz = This corresponds to the generic model equation with φ = w, Γ = µ (= constant), S = dp dz ME 448/548: D Convection-Diffusion Equation page
12 Model Problem 4: Full Developed Flow in a Rectangular Duct The smmetr in the problem allows alternative was of defining the numerical model Full Duct Quarter Duct L L L L For the full duct simulation depicted on the left hand side, the boundar conditions are no slip conditions on all four walls w(, ) = w(, L )=w(,)=w(l,)= (full duct) ME 448/548: D Convection-Diffusion Equation page Model Problem 4: Full Developed Flow in a Rectangular Duct The smmetr in the problem allows alternative was of defining the numerical model Full Duct Quarter Duct L L L L For the quarter duct simulation depicted on the right hand side, the boundar conditions are no slip conditions on the solid walls ( = L and = L ) w(l,)=w(, L )= (quarter duct) and smmetr conditions on the other two planes u = u = = = ME 448/548: D Convection-Diffusion Equation page 3
MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
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 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 informationThe Plane Stress Problem
14 The Plane Stress Problem IFEM Ch 14 Slide 1 Plate in Plane Stress Thickness dimension or transverse dimension z Top surface Inplane dimensions: in, plane IFEM Ch 14 Slide 2 Mathematical Idealization
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More information1.1 The Equations of Motion
1.1 The Equations of Motion In Book I, balance of forces and moments acting on an component was enforced in order to ensure that the component was in equilibrium. Here, allowance is made for stresses which
More informationCONSERVATION OF ENERGY FOR ACONTINUUM
Chapter 6 CONSERVATION OF ENERGY FOR ACONTINUUM Figure 6.1: 6.1 Conservation of Energ In order to define conservation of energ, we will follow a derivation similar to those in previous chapters, using
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Models for CFD ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: 1D Convection-Diffusion
More informationChapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More informationwhere is the Laplace operator and is a scalar function.
Elliptic PDEs A brief discussion of two important elliptic PDEs. In mathematics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace who first studied its
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationExercise solutions: concepts from chapter 5
1) Stud the oöids depicted in Figure 1a and 1b. a) Assume that the thin sections of Figure 1 lie in a principal plane of the deformation. Measure and record the lengths and orientations of the principal
More informationApply mass and momentum conservation to a differential control volume. Simple classical solutions of NS equations
Module 5: Navier-Stokes Equations: Appl mass and momentum conservation to a differential control volume Derive the Navier-Stokes equations Simple classical solutions of NS equations Use dimensional analsis
More informationTriple Integrals in Cartesian Coordinates. Triple Integrals in Cylindrical Coordinates. Triple Integrals in Spherical Coordinates
Chapter 3 Multiple Integral 3. Double Integrals 3. Iterated Integrals 3.3 Double Integrals in Polar Coordinates 3.4 Triple Integrals Triple Integrals in Cartesian Coordinates Triple Integrals in Clindrical
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationAlternative Boundary Condition Implementations for Crank Nicolson Solution to the Heat Equation
Alternative Boundary Condition Implementations for Crank Nicolson Solution to the Heat Equation ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu
More informationNumerical Solution of Non-Darcian Effects on Natural Convection in a Rectangular Porous Enclosure with Heated Walls
International Journal of Advanced Mechanical Engineering. ISSN 2250-3234 Volume 8, Number 1 (18), pp. 71-86 Research India Publications http://www.ripublication.com Numerical Solution of Non-Darcian Effects
More informationEstimating h Boundary Layer Equations
Estimating h Boundar Laer Equations ChE 0B Before, we just assumed a heat transfer coefficient, but can we estimate them from first rinciles? Look at stead laminar flow ast a flat late, again: Clearl,
More informationNST1A: Mathematics II (Course A) End of Course Summary, Lent 2011
General notes Proofs NT1A: Mathematics II (Course A) End of Course ummar, Lent 011 tuart Dalziel (011) s.dalziel@damtp.cam.ac.uk This course is not about proofs, but rather about using different techniques.
More informationFIRST- AND SECOND-ORDER IVPS The problem given in (1) is also called an nth-order initial-value problem. For example, Solve: Solve:
.2 INITIAL-VALUE PROBLEMS 3.2 INITIAL-VALUE PROBLEMS REVIEW MATERIAL Normal form of a DE Solution of a DE Famil of solutions INTRODUCTION We are often interested in problems in which we seek a solution
More informationRancho Bernardo High School/Math Department Honors Pre-Calculus Exit Exam
Rancho Bernardo High School/Math Department Honors Pre-Calculus Eit Eam You are about to take an eam that will test our knowledge of the RBHS Honors Pre-Calculus curriculum. You must demonstrate genuine
More informationNUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR
NUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR Jenő SUDA, István BERTA, Gergel KRISTÓF +, Tamás LAJOS # Ph.D. Student, + Assistant, # Professor, Department
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationContents. Dynamics and control of mechanical systems. Focuses on
Dnamics and control of mechanical sstems Date Da (/8) Da (3/8) Da 3 (5/8) Da 4 (7/8) Da 5 (9/8) Da 6 (/8) Content Review of the basics of mechanics. Kinematics of rigid bodies - coordinate transformation,
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationTENSOR TRANSFORMATION OF STRESSES
GG303 Lecture 18 9/4/01 1 TENSOR TRANSFORMATION OF STRESSES Transformation of stresses between planes of arbitrar orientation In the 2-D eample of lecture 16, the normal and shear stresses (tractions)
More informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationSeparation of Variables in Cartesian Coordinates
Lecture 9 Separation of Variables in Cartesian Coordinates Phs 3750 Overview and Motivation: Toda we begin a more in-depth loo at the 3D wave euation. We introduce a techniue for finding solutions to partial
More informationCONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM
Chapter 4 CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Figure 4.1: 4.1 Conservation of Angular Momentum Angular momentum is defined as the moment of the linear momentum about some spatial reference
More informationTowards Higher-Order Schemes for Compressible Flow
WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN 8-867-84- MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague,
More informationThe Plane Stress Problem
. 4 The Plane Stress Problem 4 Chapter 4: THE PLANE STRESS PROBLEM 4 TABLE OF CONTENTS Page 4.. INTRODUCTION 4 3 4... Plate in Plane Stress............... 4 3 4... Mathematical Model.............. 4 4
More informationTwo-Dimensional Formulation.
Two-Dimensional Formulation mi@seu.edu.cn Outline Introduction( 概论 ) Two vs. Three-Dimensional Problems Plane Strain( 平面应变 ) Plane Stress( 平面应力 ) Boundar Conditions( 边界条件 ) Correspondence between Plane
More informationOne of the most common applications of Calculus involves determining maximum or minimum values.
8 LESSON 5- MAX/MIN APPLICATIONS (OPTIMIZATION) One of the most common applications of Calculus involves determining maimum or minimum values. Procedure:. Choose variables and/or draw a labeled figure..
More information6 = 1 2. The right endpoints of the subintervals are then 2 5, 3, 7 2, 4, 2 9, 5, while the left endpoints are 2, 5 2, 3, 7 2, 4, 9 2.
5 THE ITEGRAL 5. Approimating and Computing Area Preliminar Questions. What are the right and left endpoints if [, 5] is divided into si subintervals? If the interval [, 5] is divided into si subintervals,
More informationGreen s Theorem Jeremy Orloff
Green s Theorem Jerem Orloff Line integrals and Green s theorem. Vector Fields Vector notation. In 8.4 we will mostl use the notation (v) = (a, b) for vectors. The other common notation (v) = ai + bj runs
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationReview of the Fundamentals of Groundwater Flow
Review of te Fundamentals of Groundwater Flow Darc s Law 1 A 1 L A or L 1 datum L : volumetric flow rate [L 3 T -1 ] 1 : draulic ead upstream [L] A : draulic ead downstream [L] : draulic conductivit [L
More information17. Find the moments of inertia I x, I y, I 0 for the lamina of. 4. D x, y 0 x a, 0 y b ; CAS. 20. D is enclosed by the cardioid r 1 cos ; x, y 3
SCTION 2.5 TRIPL INTGRALS 69 2.4 XRCISS. lectric charge is distributed over the rectangle, 2 so that the charge densit at, is, 2 2 (measured in coulombs per square meter). Find the total charge on the
More informationA first contact with STAR-CCM+ Comparing analytical and finite volume solutions with STAR-CCM+ simulations
A first contact with STAR-CCM+ Comparing analtical and finite volume solutions with STAR-CCM+ simulations Michael Heer Analtical Finite volumes STAR-CCM+ What is ParisTech? ParisTech is a consortium of
More information7. Choice of approximating for the FE-method scalar problems
7. Choice of approimating for the FE-method scalar problems Finite Element Method Differential Equation Weak Formulation Approimating Functions Weighted Residuals FEM - Formulation Weak form of heat flow
More informationDynamics and control of mechanical systems
JU 18/HL Dnamics and control of mechanical sstems Date Da 1 (3/5) 5/5 Da (7/5) Da 3 (9/5) Da 4 (11/5) Da 5 (14/5) Da 6 (16/5) Content Revie of the basics of mechanics. Kinematics of rigid bodies coordinate
More informationCONTINUOUS SPATIAL DATA ANALYSIS
CONTINUOUS SPATIAL DATA ANALSIS 1. Overview of Spatial Stochastic Processes The ke difference between continuous spatial data and point patterns is that there is now assumed to be a meaningful value, s
More informationElectromagnetic Induction (EMI) Response from Conducting and Permeable Spheroidal Shells
Electromagnetic Induction EMI Response from Conducting and Permeable Spheroidal Shells Benjamin E. Barrowes, Kevin O Neill, Chi O. Ao, and J. A. Kong Department of Electrical Engineering and Computer Science
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationAircraft Structures Structural & Loading Discontinuities
Universit of Liège Aerospace & Mechanical Engineering Aircraft Structures Structural & Loading Discontinuities Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
More informationOverview of Convection Heat Transfer
Overview of Convection Heat Transfer Lecture Notes for ME 448/548 Gerald Recktenwald March 16, 2006 1 Introduction CFD is a powerful tool for engineering design involving fluid flow and heat transfer.
More informationChapter Adequacy of Solutions
Chapter 04.09 dequac of Solutions fter reading this chapter, ou should be able to: 1. know the difference between ill-conditioned and well-conditioned sstems of equations,. define the norm of a matri,
More informationSimplified Fast Multipole Methods for Micromagnetic Modeling
implified Fast Multipole Methods for Micromagnetic Modeling D. M. Apalkov and P. B.Visscher Department of Phsics and Astronom and MIT Center The Universit of Alabama This project was supported b F grants
More informationMultivariable Calculus Lecture #13 Notes. in each piece. Then the mass mk. 0 σ = σ = σ
Multivariable Calculus Lecture #1 Notes In this lecture, we ll loo at parameterization of surfaces in R and integration on a parameterized surface Applications will include surface area, mass of a surface
More informationGUIDED WAVES IN A RECTANGULAR WAVE GUIDE
GUIDED WAVES IN A RECTANGULAR WAVE GUIDE Consider waves propagating along Oz but with restrictions in the and/or directions. The wave is now no longer necessaril transverse. The wave equation can be written
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
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 informationLecture on: 3D Dislocation Dynamics: Numerical Treatment
Tri-La Short Course on Dislocations in in Materials Pleasanton, CA June 8-0, 998 Lecture on: 3D Dislocation Dnamics: Numerical Treatment H.M. Zi, M. hee & J.P. Hirth School of Mechanical and Materials
More informationAnti-diffusive Methods for Hyperbolic Heat Transfer
Anti-diffusive Methods for Hperbolic Heat Transfer Wensheng Shen and Leigh Little Department of Computational Science, SUNY Brockport Brockport, NY 44, USA Liangjian Hu Department of Applied Mathematics,
More informationPHY 5246: Theoretical Dynamics, Fall Assignment # 9, Solutions. y CM (θ = 0) = 2 ρ m
PHY 546: Theoretical Dnamics, Fall 5 Assignment # 9, Solutions Graded Problems Problem (.a) l l/ l/ CM θ x In order to find the equation of motion of the triangle, we need to write the Lagrangian, with
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationModeling of Failure along Predefined Planes in Fractured Reservoirs
PROCEEDINGS, Thirt-Ninth Workshop on Geothermal Reservoir Engineering Stanford Universit, Stanford, California, Februar 24-26, 2014 SGP-TR-202 Modeling of Failure along Predefined Planes in Fractured Reservoirs
More informationBoundary-Fitted Coordinates!
Computational Fluid Dnamics http:wwwndedu~gtrggvacfdcourse Computational Fluid Dnamics Computational Methods or Domains with Comple BoundariesI Grétar Trggvason Spring For most engineering problems it
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationOne Dimensional Convection: Interpolation Models for CFD
One Dimensional Convection: Interpolation Moels for CFD ME 448/548 Notes Geral Recktenwal Portlan State University Department of Mechanical Engineering gerry@p.eu ME 448/548: D Convection-Di usion Equation
More informationNumerical Integration (Quadrature) Another application for our interpolation tools!
Numerical Integration (Quadrature) Another application for our interpolation tools! Integration: Area under a curve Curve = data or function Integrating data Finite number of data points spacing specified
More informationBTCS Solution to the Heat Equation
BTCS Solution to the Heat Equation ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@mepdxedu ME 448/548: BTCS Solution to the Heat Equation Overview
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationOperations depend on pixel s Coordinates. Context free. Independent of pixel values. I(x,y) I (x,y )
Geometric Transformation Operations depend on piel s Coordinates. Contet free. Independent of piel values. f f (, ) = ' (, ) = ' I(, ) = I' ( f (, ), f ( ) ), (,) (, ) I(,) I (, ) Eample: Translation =
More information11.1 Double Riemann Sums and Double Integrals over Rectangles
Chapter 11 Multiple Integrals 11.1 ouble Riemann Sums and ouble Integrals over Rectangles Motivating Questions In this section, we strive to understand the ideas generated b the following important questions:
More informationReview Test 2. c ) is a local maximum. ) < 0, then the graph of f has a saddle point at ( c,, (, c ) = 0, no conclusion can be reached by this test.
eview Test I. Finding local maima and minima for a function = f, : a) Find the critical points of f b solving simultaneousl the equations f, = and f, =. b) Use the Second Derivative Test for determining
More informationMAE 323: Chapter 4. Plane Stress and Plane Strain. The Stress Equilibrium Equation
The Stress Equilibrium Equation As we mentioned in Chapter 2, using the Galerkin formulation and a choice of shape functions, we can derive a discretized form of most differential equations. In Structural
More informationMath 1431 Final Exam Review. 1. Find the following limits (if they exist): lim. lim. lim. lim. sin. lim. cos. lim. lim. lim. n n.
. Find the following its (if they eist: sin 7 a. 0 9 5 b. 0 tan( 8 c. 4 d. e. f. sin h0 h h cos h0 h h Math 4 Final Eam Review g. h. i. j. k. cos 0 n nn e 0 n arctan( 0 4 l. 0 sin(4 m. cot 0 = n. = o.
More informationRoe Scheme in Generalized Coordinates; Part I- Formulations
Roe Scheme in Generalized Coordinates; Part I- Formulations M.. Kermani. G. Plett Department Of Mechanical & Aerospace ngineering Carleton Universit, Ottawa, Canada Abstract. Most of the computations b
More informationChapter 3: BASIC ELEMENTS. solid mechanics)
Chapter 3: BASIC ELEMENTS Section 3.: Preliminaries (review of solid mechanics) Outline Most structural analsis FE codes are displacement based In this chapter we discuss interpolation methods and elements
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
More informationHandout for Adequacy of Solutions Chapter SET ONE The solution to Make a small change in the right hand side vector of the equations
Handout for dequac of Solutions Chapter 04.07 SET ONE The solution to 7.999 4 3.999 Make a small change in the right hand side vector of the equations 7.998 4.00 3.999 4.000 3.999 Make a small change in
More informationEntropy ISSN
Entrop 5, 7[4], 3 37 3 Entrop ISSN 99-43 www.mdpi.org/entrop/ Full Paper Second law analsis of laminar flow in a channel filled with saturated porous media: a numerical solution Kamel Hooman, Arash Ejlali
More informationChemical Reaction and Viscous Dissipation Effects on an Unsteady MHD Flow of Heat and Mass Transfer along a Porous Flat Plate
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 2 Issue 2 December 24 PP 977-988 ISSN 2347-37X (Print) & ISSN 2347-342 (Online) www.arcjournals.org Chemical Reaction
More information6.1 The Linear Elastic Model
Linear lasticit The simplest constitutive law or solid materials is the linear elastic law, which assumes a linear relationship between stress and engineering strain. This assumption turns out to be an
More informationPhysics 101. Vectors. Lecture 2. h0r33fy. EMU Physics Department. Assist. Prof. Dr. Ali ÖVGÜN
Phsics 101 Lecture 2 Vectors ssist. Prof. Dr. li ÖVGÜN EMU Phsics Department h0r33f www.aovgun.com Coordinate Sstems qcartesian coordinate sstem qpolar coordinate sstem qfrom Cartesian to Polar coordinate
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationVector Calculus. Vector Fields. Reading Trim Vector Fields. Assignment web page assignment #9. Chapter 14 will examine a vector field.
Vector alculus Vector Fields Reading Trim 14.1 Vector Fields Assignment web page assignment #9 hapter 14 will eamine a vector field. For eample, if we eamine the temperature conditions in a room, for ever
More informationFinite Difference Methods (FDMs) 1
Finite Difference Methods (FDMs) 1 1 st - order Approxima9on Recall Taylor series expansion: Forward difference: Backward difference: Central difference: 2 nd - order Approxima9on Forward difference: Backward
More informationINTERFACE CRACK IN ORTHOTROPIC KIRCHHOFF PLATES
Gépészet Budapest 4-5.Ma. G--Section-o ITERFACE CRACK I ORTHOTROPIC KIRCHHOFF PLATES András Szekrénes Budapest Universit of Technolog and Economics Department of Applied Mechanics Budapest Műegetem rkp.
More informationFirst Order Equations
10 1 Linear and Semilinear Equations Chapter First Order Equations Contents 1 Linear and Semilinear Equations 9 Quasilinear Equations 19 3 Wave Equation 6 4 Sstems of Equations 31 1 Linear and Semilinear
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at.sciencedirect.com ScienceDirect Procedia Engineering 9 (14 383 388 1th International Conference on Mechanical Engineering, ICME 13 Effects of volumetric heat source and temperature
More informationPhysics Gravitational force. 2. Strong or color force. 3. Electroweak force
Phsics 360 Notes on Griffths - pluses and minuses No tetbook is perfect, and Griffithsisnoeception. Themajorplusisthat it is prett readable. For minuses, see below. Much of what G sas about the del operator
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More informationIntroduction to Finite Di erence Methods
Introduction to Finite Di erence Methods ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@pdx.edu ME 448/548: Introduction to Finite Di erence Approximation
More information3.5 Finite Differences and Fast Poisson Solvers
3.5 Finite Differences and Fast Poisson Solvers It is extremely unusual to use eigenvectors to solve a linear system KU = F. You need to know all the eigenvectors of K, and (much more than that) the eigenvector
More informationBroken Extremals. Variational Methods & Optimal Control lecture 20. Broken extremals. Broken extremals
Broken etremals Variational Methods & Optimal Control lecture 2 Matthew Roughan Discipline of Applied Mathematics School of Mathematical Sciences Universit of Adelaide
More informationy=1/4 x x=4y y=x 3 x=y 1/3 Example: 3.1 (1/2, 1/8) (1/2, 1/8) Find the area in the positive quadrant bounded by y = 1 x and y = x3
Eample: 3.1 Find the area in the positive quadrant bounded b 1 and 3 4 First find the points of intersection of the two curves: clearl the curves intersect at (, ) and at 1 4 3 1, 1 8 Select a strip at
More informationMath 214 Spring problem set (a) Consider these two first order equations. (I) dy dx = x + 1 dy
Math 4 Spring 08 problem set. (a) Consider these two first order equations. (I) d d = + d (II) d = Below are four direction fields. Match the differential equations above to their direction fields. Provide
More informationMicroscopic Momentum Balance Equation (Navier-Stokes)
CM3110 Transport I Part I: Fluid Mechanics Microscopic Momentum Balance Equation (Navier-Stokes) Professor Faith Morrison Department of Chemical Engineering Michigan Technological University 1 Microscopic
More informationThe embedded finite difference method for the Poisson equation in a domain with an irregular boundary and Dirichlet boundary conditions
The embedded finite difference method for the Poisson equation in a domain with an irregular boundary and Dirichlet boundary conditions Z. Jomaa and C. Macaskill School of Mathematics & Statistics University
More informationDirectional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables.
Directional derivatives and gradient vectors (Sect. 14.5). Directional derivative of functions of two variables. Partial derivatives and directional derivatives. Directional derivative of functions of
More informationChapter 3: Three faces of the derivative. Overview
Overview We alread saw an algebraic wa of thinking about a derivative. Geometric: zooming into a line Analtic: continuit and rational functions Computational: approimations with computers 3. The geometric
More information