Numerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
|
|
- Deborah Kerry Shepherd
- 5 years ago
- Views:
Transcription
1 Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1
2 Overview These notes are meant to be an overview of the various memory access patterns in typical Partial Differential Equation (PDE) solvers in the physical sciences. Although our experience has been mostly in fluids and structures, the comments should be applicable to other systems governed by PDEs. The comments in these slides are based on our experience with the following types of flow solvers: 1. Single-block structured. 2. Multiblock structured. 3. Tetrahedral unstructured. SSC Working Group Meetings, JJA 2
3 4. Cartesian unstructured. Only a few comments about Cartesian unstructured methods will be made (for DoE Adaptive Mesh Refinement - AMR applications.) SSC Working Group Meetings, JJA 3
4 ABird seyeviewofapde Regardless of the physical system you want to solve, a large majority of the PDEs in the physical sciences can be written in the following form w t + R(w) =0, (1) where w is the vector of dependent variables and R is a typically non-linear vector function called the residual. The equation is solved in a domain of finite size subject to some boundary conditions. The domain is discretized using a mesh of some type and the flow variables are stored at either the nodes or the cell centers. R(w) results from the process of semidiscretization by which all of the spatial derivatives are discretized, while the time derivatives are not. All of the first, second, or higher order spatial derivatives of the dependent variables are discretized using some stencil and lumped into R(w). SSC Working Group Meetings, JJA 4
5 How Do You Write a PDE solver? 2 Min Crash Course Just follow the recipe: 1. Choose a suitable spatial discretization. 2. For all nodes in the mesh, calculate the spatial residual according to the discretization chosen. 3. Update our current guess of the solution w using either a time-stepping scheme or some sort of matrix factorization. 4. Repeat until converged. SSC Working Group Meetings, JJA 5
6 Euler Equations in 2D For example, the 2D compressible, Euler (inviscid) equations can be used as an approximation to the behavior of a flow. Let p, ρ, u, v, H, and E denote the pressure, density, cartesian velocity components, total enthalpy, and total energy respectively. Consider a control volume Ω with boundary Ω. The equations of motion of the fluid can then be written in integral form as d dt Ω w dx dy + Ω (f dy g dx) =0, (2) SSC Working Group Meetings, JJA 6
7 where w is the vector of flow variables ρ ρu w = ρv ρe, and f, g are the Euler flux vectors f = ρ(u x t ) ρu(u x t )+p ρv(u x t ) ρe(u x t )+pu, g = ρ(v y t ) ρu(v y t ) ρv(v y t )+p ρe(v y t )+pv. Also, for an ideal gas, the equation of state may be written as p =(γ 1) ρ [ E 1 ] 2 (u2 + v 2 ). SSC Working Group Meetings, JJA 7
8 Applying Equation 2 independently to each cell in the mesh we obtain a set of ordinary differential equations of the form d dt (w ij V ij )+R(w ij )=0, (3) where V ij is the volume of the i, j cell and the residual R(w ij ) is obtained by evaluating the flux integral in Equation 2. Note that if you used the Reynolds Averaged Navier-Stokes (RANS) equations, the only thing you need to do is to add another portion of the residual R(w). Once the residual has been calculated, it can be used to update the current guess for the solution w using Equation 1. This process is repeated iteratively until convergence is reached. SSC Working Group Meetings, JJA 8
9 Applicability Notice that this decription of semi-discrete PDEs is applicable to 1. Euler and RANS solvers. 2. LES (Large Eddy Simulation) procedures. 3. Finite Element methods. and largely any PDE solution method that is matrix-free and that doesn t use numerical techniques such as FFTs, etc. For purposes of SSC, the most important thing is to ask ourselves the following questions: SSC Working Group Meetings, JJA 9
10 1. How is the memory access pattern affected by the stencil of the semidiscretization? 2. How does the choice of mesh topology affect these memory accesses? 3. Once a node in the mesh (discrete memory location) is touched, what are the operations that typically need to be performed on it? It turns out that the stencil of the spatial discretization and the mesh choices are fundamental to understanding how memory is accessed. SSC Working Group Meetings, JJA 10
11 Single-Block Structured Discretization Figure 1: Structured Mesh Usually for domains that are not geometrically complex (although flow may be). SSC Working Group Meetings, JJA 11
12 Typical codes: FLO87, FLO107, DNS, development codes for testing various ideas, etc. Mesh ordering follows logical structure of an n-dimensional array: neighbors (important in computing the residual) can be easily accessed by unit shifts on array indices. Connectivity of the mesh is implicit in the data structure. For cache optimization purposes, varying stride in each coordinate direction is important. SSC Working Group Meetings, JJA 12
13 Multi-Block Structured Discretization Figure 2: Multiblock Structured Mesh Around an Airplane Geometrically complex domains. SSC Working Group Meetings, JJA 13
14 Typical codes: FLO107-MB, TFLO, etc. Mesh ordering within each block follows logical structure of an n- dimensional array: neighbors (important in computing the residual) can be easily accessed by unit shifts on array indices. Connectivity of each block is implicit in the data structure. Inter-block connectivity must be supplied by user and defines the way in which all blocks are attached to each other. Halo communication typically used to isolate blocks. SSC Working Group Meetings, JJA 14
15 Tetrahedral Unstructured Discretization Figure 3: Unstructured Tetrahedral Mesh Around A320 Used for domains with high geometric complexity (full aircraft, combustors, etc.) SSC Working Group Meetings, JJA 15
16 Typical codes: AIRPLANE, combustor LES code, etc. Mesh locations follow no particular order. In fact, mesh can be reordered to improve cache performance. Connectivity of the mesh must be supplied by user. Mapping of mesh to data structure fundamental to achieve high performance in a variety of computers. SSC Working Group Meetings, JJA 16
17 AMR Cartesian Discretization Figure 4: AMR Cartesian Mesh Used for domains with high geometric complexity (full aircraft, combustors, etc.) but typically only for either Euler flows or domains with complex fluid phenomena but Cartesian boundaries. SSC Working Group Meetings, JJA 17
18 Typical codes: DoE codes, ASCI Chicago, CART3D. Mesh locations are structured but the underlying data structure is usually of the tree type. Connectivity of the mesh is inherent in data structure but not trivial. Mapping of mesh to data structure fundamental to achieve high performance in a variety of computers. SSC Working Group Meetings, JJA 18
19 Stencils in Various Types of Meshes Single-block structured and multi-block structured codes usually have a one- or two layer stencil structure. That means that in 2D all nodes/cells that adjoin the node/cell in question are likely to belong to the stencil. Examples later on. The situation in unstructured meshes is similar, although, due to the nature of the mesh, a little more random SSC Working Group Meetings, JJA 19
20 Note that in a node-based unstructured mesh solver, if a double halo is required, one may need to reach 2 edges over. Depending on the edge connectivity (and its incidence order) one may need a large amount of data. SSC Working Group Meetings, JJA 20
21 Operations/Work Required For a typical solver like TFLO, the following operations need to be performed, every iteration, for every cell in the mesh (TFLO is a cellcentered code.) Inviscid fluxes Artificial dissipation fluxes SSC Working Group Meetings, JJA 21
22 Viscous fluxes Time-step computation Residual smoothing SSC Working Group Meetings, JJA 22
23 Solution update Other operations need to be considered such as 1. Multigrid 2. Interpolation between meshes for multiphysics computations (DoE) 3. Matrix factorization steps (ADI and similar). Similar to Residual Smoothing. 4. Local and global preconditioning (fall in the classes described above). SSC Working Group Meetings, JJA 23
24 5. Advanced artificial dissipation algorithms. SSC Working Group Meetings, JJA 24
25 Code Excerpts Subroutine EFLUX from FLO82 (2D Euler solver for airfoils) C C C FLUX IN THE J DIRECTION DO 30 J=2,JL DO 30 I=2,IL XX = X(I,J,1) -X(I-1,J,1) YX = X(I,J,2) -X(I-1,J,2) PA = P(I,J+1) +P(I,J) QSP = (XX*W(I,J+1,3) -YX*W(I,J+1,2))/W(I,J+1,1) QSM = (XX*W(I,J,3) -YX*W(I,J,2))/W(I,J,1) FS(I,J,1) = QSP*W(I,J+1,1) +QSM*W(I,J,1) FS(I,J,2) = QSP*W(I,J+1,2) +QSM*W(I,J,2) -YX*PA FS(I,J,3) = QSP*W(I,J+1,3) +QSM*W(I,J,3) +XX*PA FS(I,J,4) = QSP*(W(I,J+1,4) +P(I,J+1)) +QSM*(W(I,J,4) +P(I,J)) SSC Working Group Meetings, JJA 25
26 C C C 30 CONTINUE ACCUMULATE THE FLUX IN THE J DIRECTION DO 40 N=1,4 DO 40 J=2,JL DO 40 I=2,IL DW(I,J,N) = DW(I,J,N) +FS(I,J,N) -FS(I,J-1,N) 40 CONTINUE SSC Working Group Meetings, JJA 26
27 Key Issues/Questions In order to determine the features of a language that allows for all of these types of operations, we must be able to handle 1. All types (or a generalization) of the local stencils described earlier. 2. Ability to handle not just scalars, but also vector fields. 3. Rewrite code in such a way that all operations are performed one point at a time (outer loop on mesh nodes/cells, not operations!) 4. Must conceptualize the ability to do non-local operations (multigrid, residual smoothing, ADI). SSC Working Group Meetings, JJA 27
28 5. Must efficiently implement shifts in all coordinate directions so that data accesses with non-unity stride are fast. 6. For unstructured/cartesian AMR meshes, only real option is to copy small chunks of data obtained by domain decomposition (graph partitioning/graph coloring) and make sure that operations are performed efficiently within each chunk. 7. Others? SSC Working Group Meetings, JJA 28
Multigrid Algorithms for Three-Dimensional RANS Calculations - The SUmb Solver
Multigrid Algorithms for Three-Dimensional RANS Calculations - The SUmb Solver Juan J. Alonso Department of Aeronautics & Astronautics Stanford University CME342 Lecture 14 May 26, 2014 Outline Non-linear
More informationAirfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition
Airfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition.5.5.5.5.5.5.5..5.95.9.85.8.75.7 Patrick A. LeGresley and Juan J. Alonso Dept. of Aeronautics & Astronautics
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationRealistic Animation of Fluids
Realistic Animation of Fluids p. 1/2 Realistic Animation of Fluids Nick Foster and Dimitri Metaxas Realistic Animation of Fluids p. 2/2 Overview Problem Statement Previous Work Navier-Stokes Equations
More informationIntroduction to Multigrid and its Parallelization
Introduction to Multigrid and its Parallelization! Thomas D. Economon Lecture 14a May 28, 2014 Announcements 2 HW 1 & 2 have been returned. Any questions? Final projects are due June 11, 5 pm. If you are
More informationImplementation of an integrated efficient parallel multiblock Flow solver
Implementation of an integrated efficient parallel multiblock Flow solver Thomas Bönisch, Panagiotis Adamidis and Roland Rühle adamidis@hlrs.de Outline Introduction to URANUS Why using Multiblock meshes
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationAn Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid
An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid J. V. Lassaline Ryerson University 35 Victoria St, Toronto, ON, M5B 2K3, Canada D. W. Zingg University
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationA Hybrid Cartesian Grid and Gridless Method for Compressible Flows
rd AIAA Aerospace Sciences Meeting and Exhibit,, January 5, Reno, Nevada A Hybrid Cartesian Grid and Gridless Method for Compressible Flows Hong Luo and Joseph D. Baum Science Applications International
More informationComputation of Fictitious Gas Flow with Euler Equations
1 Computation of Fictitious Gas Flow with Euler Equations Pei Li and Helmut Sobieczky DLR Göttingen, Germany Abstract The Fictitious Gas Concept supports some computational design methods to construct
More informationThe gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last
Parallel Implementation of Gas-Kinetic BGK Scheme on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer Middle East Technical
More informationComputational Fluid Dynamics for Engineers
Tuncer Cebeci Jian P. Shao Fassi Kafyeke Eric Laurendeau Computational Fluid Dynamics for Engineers From Panel to Navier-Stokes Methods with Computer Programs With 152 Figures, 19 Tables, 84 Problems and
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More information3D ADI Method for Fluid Simulation on Multiple GPUs. Nikolai Sakharnykh, NVIDIA Nikolay Markovskiy, NVIDIA
3D ADI Method for Fluid Simulation on Multiple GPUs Nikolai Sakharnykh, NVIDIA Nikolay Markovskiy, NVIDIA Introduction Fluid simulation using direct numerical methods Gives the most accurate result Requires
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationPROGRAMMING OF MULTIGRID METHODS
PROGRAMMING OF MULTIGRID METHODS LONG CHEN In this note, we explain the implementation detail of multigrid methods. We will use the approach by space decomposition and subspace correction method; see Chapter:
More informationA Two-Dimensional Multigrid-Driven Navier-Stokes Solver for Multiprocessor Architectures
A wo-dimensional Multigrid-Driven Navier-Stokes Solver for Multiprocessor Architectures Juan J. Alonso, odd J. Mitty, Luigi Martinelli, and Antony Jameson Department of Mechanical and Aerospace Engineering
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationA MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES
A MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES Khalil Bensassi, and Herman Deconinck Von Karman Institute for Fluid Dynamics Aeronautics & Aerospace
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 informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationAcceleration of a 2D Euler Flow Solver Using Commodity Graphics Hardware
Acceleration of a 2D Euler Flow Solver Using Commodity Graphics Hardware T. Brandvik and G. Pullan Whittle Laboratory, Department of Engineering, University of Cambridge 1 JJ Thomson Avenue, Cambridge,
More informationA fully implicit Navier-Stokes algorithm for unstructured grids incorporating a two-equation turbulence model
Copyright 1996, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1996 A9618376, AIAA Paper 96-0414 A fully implicit Navier-Stokes algorithm for unstructured
More informationSolution of 2D Euler Equations and Application to Airfoil Design
WDS'6 Proceedings of Contributed Papers, Part I, 47 52, 26. ISBN 8-86732-84-3 MATFYZPRESS Solution of 2D Euler Equations and Application to Airfoil Design J. Šimák Charles University, Faculty of Mathematics
More informationMultigrid Pattern. I. Problem. II. Driving Forces. III. Solution
Multigrid Pattern I. Problem Problem domain is decomposed into a set of geometric grids, where each element participates in a local computation followed by data exchanges with adjacent neighbors. The grids
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 informationRobustness improvement of polyhedral mesh method for airbag deployment simulations. TU Delft
Robustness improvement of polyhedral mesh method for airbag deployment simulations. TU Delft Santiago Alagon Carrillo, Numerical Analysis Daily Supervisor: Prof. dr. ir. C. Vuik, Company Supervisor: Ir.
More informationWe present a high-order accurate space-time discontinuous Galerkin method for solving two-dimensional
A High-Order Discontinuous Galerkin Method with Unstructured Space-Time Meshes for Two-Dimensional Compressible Flows on Domains with Large Deformations Luming Wang a, Per-Olof Persson a, a Department
More informationAnalysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization
Analysis of the Adjoint Euler Equations as used for Gradient-based Aerodynamic Shape Optimization Final Presentation Dylan Jude Graduate Research Assistant University of Maryland AMSC 663/664 May 4, 2017
More informationNumerical Methods in Aerodynamics. Fluid Structure Interaction. Lecture 4: Fluid Structure Interaction
Fluid Structure Interaction Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark
More informationDiscontinuous Galerkin Spectral Element Approximations for CFD
Discontinuous Galerkin Spectral Element Approimations for CFD D.A. Kopriva Florida State Universit Tallahassee, FL 3236 G.B. Jacobs San Diego State Universit San Diego, CA 92182 September 3, 211 1 Code
More informationRecent progress in the area of meshless methods for CFD computations has shown great promise in terms
46th AIAA Aerospace Sciences Meeting and Exhibit 7 - January 28, Reno, Nevada AIAA 28-699 Edge-based Meshless Methods for ompressible Flow Simulations Aaron Katz Antony Jameson Department of Aeronautics
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationNumerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods
Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,
More informationHybrid Simulation of Wake Vortices during Landing HPCN-Workshop 2014
Hybrid Simulation of Wake Vortices during Landing HPCN-Workshop 2014 A. Stephan 1, F. Holzäpfel 1, T. Heel 1 1 Institut für Physik der Atmosphäre, DLR, Oberpfaffenhofen, Germany Aircraft wake vortices
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 informationInterface and Boundary Schemes for High-Order Methods
19th AIAA Computational Fluid Dynamics 22-25 June 29, San Antonio, Texas AIAA 29-3658 Interface and Boundary Schemes for High-Order Methods Xun Huan, Jason E. Hicken, and David W. Zingg Institute for Aerospace
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 informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
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 informationMulti-Element High-Lift Configuration Design Optimization Using Viscous Continuous Adjoint Method
JOURNAL OF AIRCRAFT Vol. 41, No. 5, September October 2004 Multi-Element High-Lift Configuration Design Optimization Using Viscous Continuous Adjoint Method Sangho Kim, Juan J. Alonso, and Antony Jameson
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationImmersed Boundary Method and Chimera Method applied to Fluid-
The numericsacademy Fixed Colloquium IBM on Moving Immersed IBM Boundary Applications Methods : Conclusion Current Status and Future Research Directions 15-17 June 2009, Academy Building, Amsterdam, the
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 informationAlgorithmic Developments in TAU
Algorithmic Developments in TAU Ralf Heinrich, Richard Dwight, Markus Widhalm, and Axel Raichle DLR Institute of Aerodynamics and Flow Technology, Lilienthalplatz 7, 38108, Germany ralf.heinrich@dlr.de,
More informationA Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids
A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids Ruo Li, Xin Wang and Weibo Zhao January 26, 27 Abstract We proposed an efficient and robust algorithm to solve the steady Euler
More informationThe purpose of this tutorial is to illustrate how to set up and solve a problem using the. Moving Deforming Mesh (MDM) using the layering algorithm.
Tutorial: Introduction The purpose of this tutorial is to illustrate how to set up and solve a problem using the following two features in FLUENT. Moving Deforming Mesh (MDM) using the layering algorithm.
More informationTwo-dimensional laminar shock wave / boundary layer interaction
Two-dimensional laminar shock wave / boundary layer interaction J.-Ch. Robinet (), V. Daru (,) and Ch. Tenaud () () SINUMEF Laboratory, ENSAM-PARIS () LIMSI-CNRS 5, Bd. de l Hôpital, PARIS 753, France
More informationApplication of coupled solver technology to CFD modeling of multiphase flows with CFX
Application of coupled solver technology to CFD modeling of multiphase flows with CFX A. Bums, A. Splawski, S. Lo & C. Guetari AEA Technology, 2000 Oxford Drive, Suite 61 0, Bethel Park, PA 151 02, USA
More informationStrö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 informationMesh-Free Applications for Static and Dynamically Changing Node Configurations
Mesh-Free Applications for Static and Dynamically Changing Node Configurations Natasha Flyer Computational Information Systems Lab National Center for Atmospheric Research Boulder, CO Meshes vs. Mesh-free
More informationExploring unstructured Poisson solvers for FDS
Exploring unstructured Poisson solvers for FDS Dr. Susanne Kilian hhpberlin - Ingenieure für Brandschutz 10245 Berlin - Germany Agenda 1 Discretization of Poisson- Löser 2 Solvers for 3 Numerical Tests
More informationOptimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers
Optimization with Gradient and Hessian Information Calculated Using Hyper-Dual Numbers Jeffrey A. Fike and Juan J. Alonso Department of Aeronautics and Astronautics, Stanford University, Stanford, CA 94305,
More informationITU/FAA Faculty of Aeronautics and Astronautics
S. Banu YILMAZ, Mehmet SAHIN, M. Fevzi UNAL, Istanbul Technical University, 34469, Maslak/Istanbul, TURKEY 65th Annual Meeting of the APS Division of Fluid Dynamics November 18-20, 2012, San Diego, CA
More informationVerification of Moving Mesh Discretizations
Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions
More informationImplementing third order compressible flow solver for hexahedral meshes in OpenFoam
Tutorial/Report in OpenFoam Course 8 Implementing third order compressible flow solver for hexahedral meshes in OpenFoam Martin Olausson, Chalmers University of Technology, SE-1 9 Gothenburg, Sweden Abstract
More informationAn Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm
An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
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 informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationA high order moving boundary treatment for compressible inviscid flows 1. Abstract
A high order moving boundary treatment for compressible inviscid flows Sirui Tan and Chi-Wang Shu Abstract We develop a high order numerical boundary condition for compressible inviscid flows involving
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationcomputational Fluid Dynamics - Prof. V. Esfahanian
Three boards categories: Experimental Theoretical Computational Crucial to know all three: Each has their advantages and disadvantages. Require validation and verification. School of Mechanical Engineering
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationSensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method
2th AIAA Computational Fluid Dynamics Conference 27-3 June 2, Honolulu, Hawaii AIAA 2-348 Sensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method Li Wang
More informationParallel Algorithms: Adaptive Mesh Refinement (AMR) method and its implementation
Parallel Algorithms: Adaptive Mesh Refinement (AMR) method and its implementation Massimiliano Guarrasi m.guarrasi@cineca.it Super Computing Applications and Innovation Department AMR - Introduction Solving
More informationAnalysis, extensions and applications of the Finite-Volume Particle Method (FVPM) PN-II-RU-TE Synthesis of the technical report -
Analysis, extensions and applications of the Finite-Volume Particle Method (FVPM) PN-II-RU-TE-2011-3-0256 - Synthesis of the technical report - Phase 1: Preparation phase Authors: Delia Teleaga, Eliza
More informationAn adaptive discretization of compressible flow using a multitude of moving Cartesian grids
An adaptive discretization of compressible flow using a multitude of moving Cartesian grids Linhai Qiu, Wenlong Lu, Ronald Fedkiw Stanford University, 353 Serra Mall Room 27, Stanford, CA 9435 Abstract
More informationADJOINT methods for sensitivity analysis involving partial
AIAA JOURNAL Vol. 46, No. 4, April 2008 ADjoint: An Approach for the Rapid Development of Discrete Adjoint Solvers Charles A. Mader and Joaquim R. R. A. Martins University of Toronto, Toronto, Ontario
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More informationFOURTH ORDER COMPACT FORMULATION OF STEADY NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS
International Journal of Mechanical Engineering and Technology (IJMET Volume 9 Issue 10 October 2018 pp. 179 189 Article ID: IJMET_09_10_11 Available online at http://www.iaeme.com/ijmet/issues.asp?jtypeijmet&vtype9&itype10
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 informationDevelopment of mesh refinement methods at CFD codes for computational fluid mechanics problems DIPLOMA THESIS
NATIONAL TECHNICAL UNIVERSITY OF ATHENS DEPARTMENT OF MECHANICAL ENGINEERING SECTION OF FLUID MECHANICS Development of mesh refinement methods at CFD codes for computational fluid mechanics problems DIPLOMA
More informationAn Upwind Multigrid Method for Solving Viscous Flows on Unstructured Triangular Meshes
An Upwind Multigrid Method for Solving Viscous Flows on Unstructured Triangular Meshes by Daryl Lawrence Bonhaus B.S. June 1990, University of Cincinnati A Thesis submitted to The Faculty of The School
More informationDevelopment of a Computational Framework for Block-Based AMR Simulations
Procedia Computer Science Volume 29, 2014, Pages 2351 2359 ICCS 2014. 14th International Conference on Computational Science Development of a Computational Framework for Block-Based AMR Simulations Hideyuki
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationACCELERATION OF A COMPUTATIONAL FLUID DYNAMICS CODE WITH GPU USING OPENACC
Nonlinear Computational Aeroelasticity Lab ACCELERATION OF A COMPUTATIONAL FLUID DYNAMICS CODE WITH GPU USING OPENACC N I C H O L S O N K. KO U K PA I Z A N P H D. C A N D I D AT E GPU Technology Conference
More informationAn Embedded Cartesian Grid Euler Solver with Radial Basis Function for Boundary Condition Implementation
46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada AIAA 2008-532 An Embedded Cartesian Grid Euler Solver with Radial Basis Function for Boundary Condition Implementation L.
More informationUnstructured Grid Numbering Schemes for GPU Coalescing Requirements
Unstructured Grid Numbering Schemes for GPU Coalescing Requirements Andrew Corrigan 1 and Johann Dahm 2 Laboratories for Computational Physics and Fluid Dynamics Naval Research Laboratory 1 Department
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationAdjoint Solver Workshop
Adjoint Solver Workshop Why is an Adjoint Solver useful? Design and manufacture for better performance: e.g. airfoil, combustor, rotor blade, ducts, body shape, etc. by optimising a certain characteristic
More informationA linear solver based on algebraic multigrid and defect correction for the solution of adjoint RANS equations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2014; 74:846 855 Published online 24 January 2014 in Wiley Online Library (wileyonlinelibrary.com)..3878 A linear solver
More informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
More informationNumerical Simulations of Fluid-Structure Interaction Problems using MpCCI
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI François Thirifay and Philippe Geuzaine CENAERO, Avenue Jean Mermoz 30, B-6041 Gosselies, Belgium Abstract. This paper reports
More informationSHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY
SHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY A. N. Ryabinin Department of Hydroaeromechanics, Faculty of Mathematics and Mechanics, Saint-Petersburg State University, St. Petersburg, Russia E-Mail: a.ryabinin@spbu.ru
More informationSTABILIZED FINITE ELEMENT METHOD WITH AN ALE STRATEGY TO SOLVE MOVING BOUNDARIES PROBLEMS
STABILIZED FINITE ELEMENT METHOD WITH AN ALE STRATEGY TO SOLVE MOVING BOUNDARIES PROBLEMS M. Viale and N. Nigro Universidad Nacional de Rosario Pelegrini 5, () Rosario, Argentina nnigro@intec.unl.edu.ar
More informationLarge-scale Gas Turbine Simulations on GPU clusters
Large-scale Gas Turbine Simulations on GPU clusters Tobias Brandvik and Graham Pullan Whittle Laboratory University of Cambridge A large-scale simulation Overview PART I: Turbomachinery PART II: Stencil-based
More information39th AIAA Aerospace Sciences Meeting and Exhibit January 8 11, 2001/Reno, NV
AIAA 1 717 Static Aero-elastic Computation with a Coupled CFD and CSD Method J. Cai, F. Liu Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697-3975 H.M. Tsai,
More informationFlux Vector Splitting Methods for the Euler Equations on 3D Unstructured Meshes for CPU/GPU Clusters
Flux Vector Splitting Methods for the Euler Equations on 3D Unstructured Meshes for CPU/GPU Clusters Manfred Liebmann Technische Universität München Chair of Optimal Control Center for Mathematical Sciences,
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 informationCUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata
CUDA Fluid simulation Lattice Boltzmann Models Cellular Automata Please excuse my layout of slides for the remaining part of the talk! Fluid Simulation Navier Stokes equations for incompressible fluids
More information