Numerical Methods for PDEs Homework 2
|
|
- Rose Butler
- 5 years ago
- Views:
Transcription
1 Numerical Methods for PDEs Homework 2 Instructor : Prof. D. Willis david willis@uml.edu In this homework problem you will investigate finite difference methods on rectangular and nonrectangular 2D domains using MATLAB (or any other programming language of your choosing). You will approximate the porous media flow inside of a coffee pod, similar to those used by singlecup coffee brewers: Figure 1: (left) An example coffee pod and (right) A depiction of how it works. ( ( add.jpg) The diagram below illustrates the coffee pod. porous coffee media: Higher pressure water is driven through the Figure 2: (left) The flow inside a single brew coffee pod will be modeled as a porous-media flow. 1
2 The 2D physical domain will approximate one of these pods (in reality this should be modeled as a three-dimensional problem here, we will assume that the 2D model is adequate to model the basic physical principles). Since a pod is symmetric about the vertical centerline (x = 0), we shall only consider half of the domain in our analysis. The domain will be defined by the pod height, H, the bottom width W B, the top width, W T, and the hole size W H, as shown in the image below: Figure 3: (left) The physical domain. The physics of the flow through the coffee in the domain will be considered as pressure driven porous-media flow. To derive the equation, we start with the differential conservation of mass equation: q = q x x + q y y = 0, (1) where q is the velocity at a given point, q x is the x component of the velocity and q y is the y component of the velocity. Introducing Darcy s law to relate the flow velocity field to the pressure gradient field through the porous media: q = κ p (2) µ where, p is the pressure at a given location in the domain, µ is the viscosity of the fluid and κ is the permeability of the porous media. Substituting equation 2 into equation 1, the result is: q = ( κµ ) p = κ µ 2 p = 2 p = 0. (3) The equation governing the pressure in the porous media is the Laplace equation, which should be familiar to you at this point: 2 p = 0 (4) 2
3 The solution which you are finding is the pressure as a function of position in the coffee pod, p(x, y). In order to solve this equation numerically, boundary conditions are required. The pressure, p, at the inlet and exit holes are known. For the sake of simplicity, we will perform the analysis in gauge pressure, so: p = P on Γ T op p = 0 on Γ Hole The boundary conditions on the symmetry plane must reflect the symmetry assumption, hence: p n = 0 The boundary conditions on all solid boundaries must prevent flow from passing through the surface, hence: q ˆn = 0 ( κµ ) n p = 0 p n = 0 The following image illustrates the domain and boundary: Figure 4: (left) The physical domain with governing partial differential equation and boundary conditions. 3
4 Please answer the following questions. Submit written answers on paper. Plots and calculated results should be printed on paper. Computer code should be delivered in your dropbox.com folder. 1. Write the bilinear mapping functions, x(ξ, η) and y(ξ, η) that you will use. Present the linearsystem (matrix and vectors) uses H, W T, and W B to define the mapping coefficients for the bilinear mapping functions. 2. Write comprehensive pseudo code to illustrate how you will determine numerical values of the bilinear mapping coefficients. 3. Implement your mapping (only the bilinear mapping) using Matlab or other coding language. Demonstrate that your mapping works by mapping a collection of points in the ξ η space to the physical space for a collection of different H, W T and W B values (present 2-3 plots to confirm your mapping process works). 4. Show using a drawing and a pseudo code how you will go between a point with indices i, j in your (computational or physical) domain to a node number (ie. how you will take a point from i, j index definition to a node number so that you can access the appropriate row & column of a matrix). 5. List ALL of the second order finite difference equations that will be used to approximate the partial derivatives in the equation governing the computational domain solution. 6. Write the finite difference equation represents the governing equation for each internal node in the domain. Draw the associated 9-point finite difference stencil. Create a 3 3 table containing the coefficients (AA, BB, CC, etc) for each point used in the finite difference stencil for node i, j. 7. Write a detailed pseudo-code that describes the Matlab/computer code process for setting up and solving this finite difference problem. 8. Membrane deflection FD-code to start: Write your own finite difference solver (you may closely follow the example from class, but for your sake, try to write a code that differs from that in class) to solve the deflection of a membrane on a non-rectangular domain (set H = 2, W B = 0.5 and W T = 1.0). You will solve the following equation: 2 u = 1 (6) with boundary conditions, u = 0 on all domain boundaries. represents the deflection of the membrane. For this question, u 9. Illustrate that the deflection calculated using this finite difference solver has second order error convergence as a function of grid size (convergence is O( x 2 )). Please plot this result using an appropriate plot (eg. log, semilog, plot, etc.) 4
5 10. Simple(r) coffee pod problem: You will next try to solve the coffee pod problem for a rectangular domain. The coffee pod problem is different than the membrane problem due to boundary conditions. First describe on paper how the boundary conditions will be applied in the A-matrix for the rectangular domain finite difference problem (you should ensure a second order accurate solution is maintained so, second order FD expressions should be used for all derivatives). 11. Still using the rectangular domain (H = 2, W T = W B = 1 and W H = 0.25). Introduce in your code the appropriate boundary conditions for the porous-media flow in the relevant locations of the A-matrix. Solve the problem using Matlab (or other programming language of choice) twice. Once for a grid with 65 nodes and second for 129 nodes respectively in the x-direction (I only prescribe x-direction nodes to allow you flexibility to use whatever you want for NNodes in y). Assume P = 1 at the top of the domain. Plot the pressure distribution solution. Does the solution look reasonable & converged? 12. Now, describe on paper how you would introduce the boundary conditions when a nonrectangular domain is used (ie. the pod domain). You don t need to implement the nonrectangular domain boundary conditions to receive full credit for this assignment. 13. Bonus Points 1: Indicate how you would calculate the flow rate of coffee through the pod using your finite difference code. 14. Bonus Points 2: Implement the boundary conditions for a non-rectangular domain. Show that the flow rate calculation converges with refinement for a non-rectangular domain defined by: Set the pressure at the top to P = 1. Set the width of the bottom to W B = 1 Set the width of the hole to W H = 0.25 Set the width of the top to W T = 1.5 Set the height to H = Bonus Points 3: Implement the boundary conditions for a non-rectangular domain. Perform the pod analysis for a range of geometries defined by: Set the pressure at the top to P = 1. Set the width of the bottom to W B = 0.5 Set the width of the hole to W H = 0.25 Vary the length of the top boundary and height to maintain a constant pod area A pod = 2. This ensures a constant amount of coffee in the pod. Calculate and plot the coffee flow rate as a function of the pod height. What do you notice? 5
DUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
More informationMATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall
MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall James Cherry, Mehmet Sözen Grand Valley State University, cherryj1@gmail.com,
More informationREDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION
REDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION Lindaura Maria Steffens Dara Liandra Lanznaster lindaura.steffens@udesc.br daraliandra@gmail.com Santa Catarina State University CESFI, Av. Central,
More informationTable 17 1 Some general field equation terms. Heat Power. Current Source. 0 0 Boundary Current Porous Media Flow. Flow Source
17 Related Analogies 17.1 Basic Concepts The differential equation used in a finite element study in one discipline often appears in a different discipline, but with a different physical meaning for the
More informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
More informationWorksheet Week 7 Section
Worksheet Week 7 Section 8.. 8.4. This worksheet is for improvement of your mathematical writing skill. Writing using correct mathematical epression and steps is really important part of doing math. Please
More informationEfficient numerical solution of the Biot poroelasticity system in multilayered domains
Efficient numerical solution of the Biot poroelasticity system Anna Naumovich Oleg Iliev Francisco Gaspar Fraunhofer Institute for Industrial Mathematics Kaiserslautern, Germany, Spain Workshop on Model
More informationPage 1. Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.)
Page 1 Neatly print your name: Signature: (Note that unsigned exams will be given a score of zero.) Circle your lecture section (-1 point if not circled, or circled incorrectly): Prof. Vlachos Prof. Ardekani
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationAdditive Manufacturing Module 8
Additive Manufacturing Module 8 Spring 2015 Wenchao Zhou zhouw@uark.edu (479) 575-7250 The Department of Mechanical Engineering University of Arkansas, Fayetteville 1 Evaluating design https://www.youtube.com/watch?v=p
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationModeling as a tool for understanding the MEA. Henrik Ekström Utö Summer School, June 22 nd 2010
Modeling as a tool for understanding the MEA Henrik Ekström Utö Summer School, June 22 nd 2010 COMSOL Multiphysics and Electrochemistry Modeling The software is based on the finite element method A number
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationCS 4700: Foundations of Artificial Intelligence Homework 2 Solutions. Graph. You can get to the same single state through multiple paths.
CS 4700: Foundations of Artificial Intelligence Homework Solutions 1. (30 points) Imagine you have a state space where a state is represented by a tuple of three positive integers (i,j,k), and you have
More informationCourse Information Course Overview Study Skills Background Material. Introduction. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics
Introduction CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James CS 205A: Mathematical Methods Introduction 1 / 16 Instructor Prof. Doug James Office: Gates 363 Telephone: (650)
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationSound Propagation in Porous Media
Final Project Report for ENGN34 Sound Propagation in Porous Media ---Numerical simulation based on MATLAB Name: Siyuan Song Department: Engineering Date: Dec.15 17 1 Name: Siyuan Song Department: Engineering
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 informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationTurbomachinery Aerodynamics Prof. Bhaskar Roy Prof. A M Pradeep Department of Aerospace Engineering Indian Institute of Technology, Bombay
Turbomachinery Aerodynamics Prof. Bhaskar Roy Prof. A M Pradeep Department of Aerospace Engineering Indian Institute of Technology, Bombay Lecture No. # 40 CFD for Turbomachinery: 2D and 3D Blade Generation
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationLesson 13: More Factoring Strategies for Quadratic Equations & Expressions
: More Factoring Strategies for Quadratic Equations & Expressions Opening Exploration Looking for Signs In the last lesson, we focused on quadratic equations where all the terms were positive. Juan s examples
More information2 Regularized Image Reconstruction for Compressive Imaging and Beyond
EE 367 / CS 448I Computational Imaging and Display Notes: Compressive Imaging and Regularized Image Reconstruction (lecture ) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement
More informationMath Partial Differential Equations
Math 531 - Partial Differential Equations to Partial Differential Equations Joseph M. Mahaffy, jmahaffy@sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences
More informationLet s Do Algebra Tiles
Let s Do Algebra Tiles Algebra Tiles Algebra tiles can be used to model operations involving integers. Let the small green square represent +1 and the small pink square represent -1. The green and pink
More informationDESIGN OF MICRO-FLUIDIC BIO-REACTORS USING TOPOLOGY OPTIMIZATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 DESIGN OF MICRO-FLUIDIC BIO-REACTORS USING TOPOLOGY OPTIMIZATION
More information1 Binary Classification
CS6501: Advanced Machine Learning 017 Problem Set 1 Handed Out: Feb 11, 017 Due: Feb 18, 017 Feel free to talk to other members of the class in doing the homework. You should, however, write down your
More informationB005 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION
1 B5 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION LUDOVIC RICARD, MICAËLE LE RAVALEC-DUPIN, BENOÎT NOETINGER AND YVES GUÉGUEN Institut Français du Pétrole, 1& 4 avenue Bois Préau, 92852
More informationMachine Learning: Assignment 1
10-701 Machine Learning: Assignment 1 Due on Februrary 0, 014 at 1 noon Barnabas Poczos, Aarti Singh Instructions: Failure to follow these directions may result in loss of points. Your solutions for this
More informationUse the Rational Zero Theorem to list all the possible rational zeros of the following polynomials. (1-2) 4 3 2
Name: Math 114 Activity 1(Due by EOC Apr. 17) Dear Instructor or Tutor, These problems are designed to let my students show me what they have learned and what they are capable of doing on their own. Please
More informationBring a printed copy of this lab to your lab section. We will answer the lettered questions in class.
Lab 2 GEO 302C Week of January 30, 2006. Bring a printed copy of this lab to your lab section. We will answer the lettered questions in class. Goal for today: Be able to understand the greenhouse effect
More informationLecture 29: Forced Convection II
Lecture 29: Forced Convection II Notes by MIT Student (and MZB) As discussed in the previous lecture, the magnitude of limiting current can be increased by imposing convective transport of reactant in
More informationAP Calculus AB Information and Summer Assignment
AP Calculus AB Information and Summer Assignment General Information: Competency in Algebra and Trigonometry is absolutely essential. The calculator will not always be available for you to use. Knowing
More informationHomework #1. Denote the sum we are interested in as To find we subtract the sum to find that
Homework #1 CMSC351 - Spring 2013 PRINT Name : Due: Feb 12 th at the start of class o Grades depend on neatness and clarity. o Write your answers with enough detail about your approach and concepts used,
More information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More informationAn a posteriori error estimator for the weak Galerkin least-squares finite-element method
An a posteriori error estimator for the weak Galerkin least-squares finite-element method James H. Adler a, Xiaozhe Hu a, Lin Mu b, Xiu Ye c a Department of Mathematics, Tufts University, Medford, MA 02155
More informationAPPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell
APPENDIX Tidally induced groundwater circulation in an unconfined coastal aquifer modeled with a Hele-Shaw cell AaronJ.Mango* Mark W. Schmeeckle* David Jon Furbish* Department of Geological Sciences, Florida
More information18.303: Introduction to Green s functions and operator inverses
8.33: Introduction to Green s functions and operator inverses S. G. Johnson October 9, 2 Abstract In analogy with the inverse A of a matri A, we try to construct an analogous inverse  of differential
More informationHomework 3 solutions Math 136 Gyu Eun Lee 2016 April 15. R = b a
Homework 3 solutions Math 136 Gyu Eun Lee 2016 April 15 A problem may have more than one valid method of solution. Here we present just one. Arbitrary functions are assumed to have whatever regularity
More information*** ***! " " ) * % )!( & ' % # $. 0 1 %./ +, - 7 : %8% 9 ) 7 / ( * 7 : %8% 9 < ;14. " > /' ;-,=. / ١
١ ******!" #$ % & '!( ) % * ") +,-./ % 01. 3 ( 4 56 7/4 ) 8%9 % : 7 ;14 < 8%9 % : *7./ = ;-, >/'." Soil Permeability & Seepage ٢ Soil Permeability- Definition ٣ What is Permeability? Permeability is the
More information8 A pseudo-spectral solution to the Stokes Problem
8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,
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 informationFluid Flow Fluid Flow and Permeability
and Permeability 215 Viscosity describes the shear stresses that develop in a flowing fluid. V z Stationary Fluid Velocity Profile x Shear stress in the fluid is proportional to the fluid velocity gradient.
More informationNo prep assignment to do, but here are four questions anyway.
Preparation Assignments for Homework #3 Due at the start of class. Reading Assignments Please see the handouts for each lesson for the reading assignments. 3,4 February Lesson 2.5 No prep assignment to
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit VII Sparse Matrix
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 1: Direct Methods Dianne P. O Leary c 2008
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition
More informationPetroleum reservoir simulation using EbFVM: the negative transmissibility issue
Computational Methods in Multhase Flow IV 133 etroleum reservoir simulation using EbFVM: the negative transmissibility issue C. R. Maliska, J. Cordazzo & A. F. C. Silva Computational Fluid Dynamics Laboratory,
More informationA Lagrangian approach to the study of the kinematic dynamo
1 A Lagrangian approach to the study of the kinematic dynamo Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ October
More informationMath 1132 Practice Exam 1 Spring 2016
University of Connecticut Department of Mathematics Math 32 Practice Exam Spring 206 Name: Instructor Name: TA Name: Section: Discussion Section: Read This First! Please read each question carefully. Show
More informationThe effect of Entry Region on Thermal Field
The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving
More informationFormal Definition of a Finite Automaton. August 26, 2013
August 26, 2013 Why a formal definition? A formal definition is precise: - It resolves any uncertainties about what is allowed in a finite automaton such as the number of accept states and number of transitions
More informationAnalysis of oil displacement by water in oil reservoirs with horizontal wells
Analysis of oil displacement by water in oil reservoirs with horizontal wells Paulo Dore Fernandes, Thiago Judson L. de Oliveira and Rodrigo A. C. Dias Problem Description This work involves near-well
More informationThe Convergence of Mimetic Discretization
The Convergence of Mimetic Discretization for Rough Grids James M. Hyman Los Alamos National Laboratory T-7, MS-B84 Los Alamos NM 87545 and Stanly Steinberg Department of Mathematics and Statistics University
More informationREVIEW: Waves on a String
Lecture 14: Solution to the Wave Equation (Chapter 6) and Random Walks (Chapter 7) 1 Description of Wave Motion REVIEW: Waves on a String We are all familiar with the motion of a transverse wave pulse
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Solution to Final Examination Date: May 4, 2011 1:00 3:00 PM Instructor: J. Murthy Open Book, Open Notes Total: 100 points Use the finite volume
More informationMath 334 Midterm III KEY Fall 2006 sections 001 and 004 Instructor: Scott Glasgow
Math 334 Midterm III KEY Fall 6 sections 1 and 4 Instructor: Scott Glasgow Please do NO write on this exam No credit will be given for such work Rather write in a blue book, or on your own paper, preferably
More informationLab Week 6. Quiz #3 Voltage Divider Homework P11, P12 Kirchhoff's Voltage Law (KVL) Kirchhoff's Current Law (KCL) KCL + KVL Module Report tips
Lab Week 6 Quiz #3 Voltage Divider Homework P11, P12 Kirchhoff's Voltage Law (KVL) Kirchhoff's Current Law (KCL) KCL + KVL Module Report tips Quiz 3 Voltage Divider (20 pts.) Please clear desks and turn
More informationCS446: Machine Learning Spring Problem Set 4
CS446: Machine Learning Spring 2017 Problem Set 4 Handed Out: February 27 th, 2017 Due: March 11 th, 2017 Feel free to talk to other members of the class in doing the homework. I am more concerned that
More informationAltered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems
Altered Jacobian Newton Iterative Method for Nonlinear Elliptic Problems Sanjay K Khattri Abstract We present an Altered Jacobian Newton Iterative Method for solving nonlinear elliptic problems Effectiveness
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 informationImpact of a Jet. Experiment 4. Purpose. Apparatus. Theory. Symmetric Jet
Experiment 4 Impact of a Jet Purpose The purpose of this experiment is to demonstrate and verify the integral momentum equation. The force generated by a jet of water deflected by an impact surface is
More informationRevision of the thermal chip simulation, J.E. Akin, Rice University
Revision of the thermal chip simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationMachine Learning Applied to 3-D Reservoir Simulation
Machine Learning Applied to 3-D Reservoir Simulation Marco A. Cardoso 1 Introduction The optimization of subsurface flow processes is important for many applications including oil field operations and
More informationPhilosophy 220. Truth-Functional Equivalence and Consistency
Philosophy 220 Truth-Functional Equivalence and Consistency Review Logical equivalency: The members of a pair of sentences [of natural language] are logically equivalent if and only if it is not [logically]
More informationCapacitance Matrix Method
Capacitance Matrix Method Marcus Sarkis New England Numerical Analysis Day at WPI, 2019 Thanks to Maksymilian Dryja, University of Warsaw 1 Outline 1 Timeline 2 Capacitance Matrix Method-CMM 3 Applications
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More information. =. a i1 x 1 + a i2 x 2 + a in x n = b i. a 11 a 12 a 1n a 21 a 22 a 1n. i1 a i2 a in
Vectors and Matrices Continued Remember that our goal is to write a system of algebraic equations as a matrix equation. Suppose we have the n linear algebraic equations a x + a 2 x 2 + a n x n = b a 2
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 informationCOMPARISON OF CPU AND GPU IMPLEMENTATIONS OF THE LATTICE BOLTZMANN METHOD
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) c CIMNE, Barcelona, 2010 COMPARISON OF CPU AND GPU IMPLEMENTATIONS OF THE LATTICE BOLTZMANN METHOD James.E. McClure, Jan F. Prins
More informationAn investigation of mechanisms involved in the web-to-drum heat transfer coefficient when using gas injection between the web & drum.
An investigation of mechanisms involved in the web-to-drum heat transfer coefficient when using gas injection between the web & drum. Mike McCann, McCann Science Charles Bishop, Bishop Consulting Ltd 1
More informationElements of Matlab and Simulink Lecture 7
Elements of Matlab and Simulink Lecture 7 Emanuele Ruffaldi 12th May 2009 Copyright 2009,Emanuele Ruffaldi. This work is licensed under the Creative Commons Attribution-ShareAlike License. PARTIAL DIFFERENTIAL
More information10-701/15-781, Machine Learning: Homework 4
10-701/15-781, Machine Learning: Homewor 4 Aarti Singh Carnegie Mellon University ˆ The assignment is due at 10:30 am beginning of class on Mon, Nov 15, 2010. ˆ Separate you answers into five parts, one
More informationCSCI 150 Discrete Mathematics Homework 5 Solution
CSCI 150 Discrete Mathematics Homework 5 Solution Saad Mneimneh Computer Science Hunter College of CUNY Problem 1: Happy Birthday (if it applies to you)! Based on the size of the class, there is approximately
More informationRepresentations of Motion in One Dimension: Speeding up and slowing down with constant acceleration
Representations of Motion in One Dimension: Speeding up and slowing down with constant acceleration Name: Group Members: Date: TA s Name: Apparatus: Aluminum track and supports, PASCO Smart Cart, two cart
More informationOptimization of DPF Structures with a 3D-Unit Cell Model
Optimization of DPF Structures with a 3D-Unit Cell Model Wieland Beckert, Marcel Dannowski, Lisabeth Wagner, Jörg Adler, Lars Mammitzsch Fraunhofer IKTS, Dresden, Germany *Corresponding author: FhG IKTS,
More informationEmden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability Porous Layers
4 Int. J. Open Problems Compt. Math., Vol. 0, No., June 07 ISSN 998-66; Copyright ICSRS Publication, 07 www.i-csrs.org Emden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability
More informationTurbomachinery Aerodynamics Prof. Bhaskar Roy Prof. A M Pradeep Department of Aerospace Engineering Indian Institute of Technology, Bombay
Turbomachinery Aerodynamics Prof. Bhaskar Roy Prof. A M Pradeep Department of Aerospace Engineering Indian Institute of Technology, Bombay Lecture No. # 26 Tutorial 4: 3D Flows in Axial Flow Turbines We
More informationHOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, Maximum score 6.0 pts.
HOMEWORK 4: Numerical solution of PDEs for Mathmatical Models, Analysis and Simulation, Fall 2011 Report due Mon Nov 12, 2012. Maximum score 6.0 pts. Read Strang s (old) book, Sections 3.1-3.4 and 6.4-6.5.
More informationCFD analysis of the transient flow in a low-oil concentration hydrocyclone
CFD analysis of the transient flow in a low-oil concentration hydrocyclone Paladino, E. E. (1), Nunes, G. C. () and Schwenk, L. (1) (1) ESSS Engineering Simulation and Scientific Software CELTA - Rod SC-41,
More informationChapter 2 CONTINUUM MECHANICS PROBLEMS
Chapter 2 CONTINUUM MECHANICS PROBLEMS The concept of treating solids and fluids as though they are continuous media, rather thancomposedofdiscretemolecules, is one that is widely used in most branches
More informationCOMPUTATIONAL FLUID DYNAMICS (CFD) FOR THE OPTIMIZATION OF PRODUCTS AND PROCESSES
THE INTERNATIONAL CONFERENCE OF THE CARPATHIAN EURO-REGION SPECIALISTS IN INDUSTRIAL SYSTEMS 7 th EDITION COMPUTATIONAL FLUID DYNAMICS (CFD) FOR THE OPTIMIZATION OF PRODUCTS AND PROCESSES Franz, Haas DI
More informationNAPC Numerical investigation of axisymmetric underexpanded supersonic jets. Pratikkumar Raje. Bijaylakshmi Saikia. Krishnendu Sinha 1
Proceedings of the 1 st National Aerospace Propulsion Conference NAPC-2017 March 15-17, 2017, IIT Kanpur, Kanpur NAPC-2017-139 Numerical investigation of axisymmetric underexpanded supersonic jets Pratikkumar
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationFINITE DIFFERENCE METHODS (II): 1D EXAMPLES IN MATLAB
1.723 - COMPUTATIONAL METHODS FOR FLOW IN POROUS MEDIA Spring 2009 FINITE DIFFERENCE METHODS (II): 1D EXAMPLES IN MATLAB Luis Cueto-Felgueroso 1. COMPUTING FINITE DIFFERENCE WEIGHTS The function fdcoefs
More informationCHEN 7100 FA16 Final Exam
CHEN 7100 FA16 Final Exam Show all work and state all assumptions for full credit. The exam is closed book, notes, and homework. Only the course reader and your mind should be open. No electronic devices
More informationSignature. Printed Name. Math 312 Hour Exam 1 Jerry L. Kazdan March 5, :00 1:20
Signature Printed Name Math 312 Hour Exam 1 Jerry L. Kazdan March 5, 1998 12:00 1:20 Directions: This exam has three parts. Part A has 4 True-False questions, Part B has 3 short answer questions, and Part
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationPerformance Analysis of a Two phase Non-isothermal PEM Fuel Cell
Performance Analysis of a Two phase Non-isothermal PEM Fuel Cell A. H. Sadoughi 1 and A. Asnaghi 2 and M. J. Kermani 3 1, 2 Ms Student of Mechanical Engineering, Sharif University of Technology Tehran,
More informationBoundary Conditions - Inlet
Boundary Conditions - Inlet Boundary Conditions Need? What is physical and mathematical significance of a boundary condition? The boundary conditions of any problem is used to define the upper and lower
More informationAgain we will consider the following one dimensional slab of porous material:
page 1 of 7 REVIEW OF BASIC STEPS IN DERIVATION OF FLOW EQUATIONS Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations,
More informationPart III (15 credits)
PHYSICS continued Part III (15 credits) Please refer to the Department publication Regents Examination in Physics: Rating Guide f Part III. Teachers should become familiar with this guide befe rating students
More informationThe dependence of the cross-sectional shape on the hydraulic resistance of microchannels
3-weeks course report, s0973 The dependence of the cross-sectional shape on the hydraulic resistance of microchannels Hatim Azzouz a Supervisor: Niels Asger Mortensen and Henrik Bruus MIC Department of
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More information