Verification of Moving Mesh Discretizations
|
|
- Marshall McDaniel
- 5 years ago
- Views:
Transcription
1 Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018
2 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions on moving domains and meshes Considerations: Inviscid versus viscous term discretization Solution field norm versus scalar outputs Integration quadrature requirements Mapping singularities Geometric conservation Boundary condition consistency Mesh Motion Verification Introduction 2/19
3 The arbitrary Lagrangian-Eulerian (ALE) method ALE Idea: solve transformed PDE on a static reference domain Reference domain: X, u X, H X ux t NdA + X H X X(u X, q X) = 0 Mapping X, t x( X, t) G = x X g = det(g) u X = gu q X = gg T q v G = x t H X = gg 1 H u X G 1 v G nda = gg T NdA NdA = g 1 G T nda Physical domain: x, u, H u t nda x + H(u, q) = 0 Key definitions X = reference-domain coordinates x = physical-domain coordinates g = determinant of Jacobian matrix v G = grid velocity, x/ t u = physical state u X = reference state H = physical flux vector H X = reference flux vector Mesh Motion Verification Introduction 3/19
4 The transformed equations Integrate the evolution PDE over a time-varying volume v(t), apply the divergence theorem, transform terms to reference space, and apply the divergence theorem again, u X t + X H X (u X, X u X ) = 0, X where u X = gu H X = gg 1 H ux G 1 v G X is the gradient with respect to the reference coordinates The transformed flux, H X, separates into inviscid and viscous contributions, H X = F X + G X, FX = gg 1 F + ux G 1 v G, GX = gg 1 G Mesh Motion Verification Introduction 4/19
5 ALE implementation To minimize code intrusion, express the reference-space fluxes in terms of the physical fluxes FX = gg 1 F ux G 1 v G = gg 1 ( F u vg ) G X = gg 1 G = gg 1 K q = } G 1 {{ KG T } q X K X Linearizations must be performed w.r.t reference states Viscous stabilization does not change (same as physical)! Boundary conditions (e.g. no slip) must incorporate the boundary velocity, v G. Mesh Motion Verification Introduction 5/19
6 Outputs General form of boundary flux integral outputs: J = o T Ĥ n da, o R s is a weight function that defines the output For 2D Navier-Stokes: o = [0; cos α; sin α; 0] Drag o = [0; y; x; 0] Moment about origin o = [0; v Gx ; v Gy ; 0] Power Ω For example, if v G = v 0 + ω r, where r is a position vector relative to some origin of rotation, J = ( v 0 + ω r) f surf da = v 0 f surf da + ω ( r f surf ) da Ω Ω } {{ } Ω } {{ } v 0 F net ω T net Mesh Motion Verification Introduction 6/19
7 Inviscid field norm check: Euler vortex Analytical vortex solution to the Euler equations Sinusoidal interior mesh deformation for testing x(t) = X + 2 sin(2πx/20) sin(2πy/15) sin(2πt) y(t) = Y sin(2πx/20) sin(2πy/15) sin(4πt) Pressure at final time Deformed mesh Mesh Motion Verification Introduction 7/19
8 Euler vortex results Look at L 2 error norm at end of the simulation 10 0 DIRK4 with small t 2 L 2 error norm h = sqrt(1/nelem) DG, p=1 DG, p=2 DG, HDG, p=1 HDG, p=2 HDG, Mesh Motion Verification Introduction 8/19
9 Inviscid boundary output: bump M=0.2, structured mesh, sinusoidal mesh motion Mach contours mesh at t = 0.0 mesh at t = 0.5 mesh at t = 1.5 Mesh Motion Verification Introduction 9/19
10 Inviscid boundary output: bump Monitor instantaneous drag output on the bump Compare moving mesh results to a static mesh Coarse (p = 3, N t = 40) and fine (p = 4, N t = 60) solutions Verification: error decreases with increasing resolution drag e-05,N_t=40, motion,n_t=40,n_t=60, motion,n_t=60 6e-05 4e time Mesh Motion Verification Introduction 10/19
11 Viscous boundary output: flat plate M=0.2, Re = 1000, adapted mesh, sinusoidal mesh motion Mach contours mesh at t = 0.0 mesh at t = 0.5 mesh at t = 1.5 Mesh Motion Verification Introduction 11/19
12 Viscous boundary output: flat plate Monitor instantaneous drag output on the flat plate Compare moving mesh results to a static mesh Coarse (p = 2, N t = 20) and fine (p = 3, N t = 40) solutions Verification: error decreases with increasing resolution drag p=2,n_t=20, motion p=2,n_t=20,n_t=40, motion,n_t= time Mesh Motion Verification Introduction 12/19
13 NACA 0012 airfoil in pitch/plunge motion Smooth plunge, h(t), and pitch, θ(t), motion, 0 t 1 h(t) = 1 cos πt, θ(t) = π cos 2πt. 2 NACA 0012, Re c = 5000, M = 0.2, steady-state IC t = 0 t = 0.5 t = 1.0 Mesh Motion Verification Introduction 13/19
14 Outputs Impulse = Work = F y dt = 1 0 airfoil ŷ f surf ds dt ] [ v 0 F net + ω T net dt = Work done on airfoil appears as total energy deficit in flow 1 0 airfoil v G f surf ds dt Power/Work Instantaneous power Time integrated power Delta total energy in flow Time Mesh Motion Verification Introduction 14/19
15 Outputs Impulse = Work = F y dt = 1 0 airfoil ŷ f surf ds dt ] [ v 0 F net + ω T net dt = Work done on airfoil appears as total energy deficit in flow 1 0 airfoil v G f surf ds dt Power/Work Instantaneous power Time integrated power Delta total energy in flow Time Mesh Motion Verification Introduction 14/19
16 Worskhop runs meshes: quasi-uniform refinement Mesh 0: 664 elements (farfield at 100c) Mesh 2: 7389 elements Mesh 1: 2137 elements Mesh 3: elements Mesh Motion Verification Introduction 15/19
17 Worskhop runs: details Considered spatial orders p = 1, 2, 3, 4, 5 ESDIRK5 (5th order) time stepping Number of time steps depends on mesh and order: p = 1 p = 2 p = 3 p = 4 p = 5 Mesh Mesh Mesh (these values were determined empirically) Error computed relative to a truth solution from a p = 5, N t = 800, ESDIRK5 run on Mesh 3 Mesh Motion Verification Introduction 16/19
18 Worskhop result: case 1, pure heaving h(t) = t2 (3 t), θ(t) = LiftInt error p= p= /(dof) 1/2 LiftInt error p= p= work units Truth lift integral = Mesh Motion Verification Introduction 17/19
19 Worskhop result: case 1, pure heaving h(t) = t2 (3 t), θ(t) = 0 4 PowerInt error 10 3 p=1 p= /(dof) 1/2 PowerInt error 10 3 p=1 p= work units Truth power integral = Mesh Motion Verification Introduction 17/19
20 Worskhop result: case 2, flow aligning h(t) = t2 (3 t), θ(t) = π 4 3 t2 (t 2 4t + 4) LiftInt error p=1 p=2 LiftInt error p=1 p= /(dof) 1/ work units Truth lift integral = Mesh Motion Verification Introduction 18/19
21 Worskhop result: case 2, flow aligning h(t) = t2 (3 t), θ(t) = π 4 3 t2 (t 2 4t + 4) PowerInt error p=1 p= /(dof) 1/2 PowerInt error 10 3 p=1 p= work units Truth power integral = Mesh Motion Verification Introduction 18/19
22 Worskhop result: case 3, energy extracting h(t) = t2 16 ( 8t3 + 51t 2 111t + 84), θ(t) = 4π 9 t2 (t 2 4t + 4) LiftInt error 10 3 p=1 p= /(dof) 1/2 LiftInt error 10 3 p=1 p= work units Truth lift integral = Mesh Motion Verification Introduction 19/19
23 Worskhop result: case 3, energy extracting h(t) = t2 16 ( 8t3 + 51t 2 111t + 84), θ(t) = 4π 9 t2 (t 2 4t + 4) PowerInt error 10 3 p=1 p= /(dof) 1/2 PowerInt error 10 3 p=1 p= work units Truth power integral = Mesh Motion Verification Introduction 19/19
Studies 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 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 informationHigh-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA )
High-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA9550-07-0195) Sachin Premasuthan, Kui Ou, Patrice Castonguay, Lala Li, Yves Allaneau,
More informationOutput-Adaptive Solution Strategies for Unsteady Aerodynamics on Deformable Domains
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, awaii, July 9-13, 2012 ICCFD7-3802 Output-Adaptive Solution Strategies for Unsteady Aerodynamics on Deformable Domains
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 informationA High Order Discontinuous Galerkin Method for Fluid-Structure Interaction
A High Order Discontinuous Galerkin Method for Fluid-Structure Interaction Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A. Javier Bonet University
More informationCase 3.1: Turbulent Flow over a 2D Multi-Element Airfoil. Summary of Results. Marco Ceze
Case 3.1: Turbulent Flow over a 2D Multi-Element Airfoil Summary of Results Marco Ceze (mceze@umich.edu) 2 nd International Workshop on High-Order CFD Methods, May 27-28, Cologne, Germany C3.1 1/14 Case
More informationDevelopment of a Consistent Discrete Adjoint Solver for the SU 2 Framework
Development of a Consistent Discrete Adjoint Solver for the SU 2 Framework Tim Albring, Max Sagebaum, Nicolas Gauger Chair for Scientific Computing TU Kaiserslautern 16th Euro-AD Workshop, Jena December
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 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 informationThree-dimensional numerical simulations of flapping wings at low Reynolds numbers
Three-dimensional numerical simulations of flapping wings at low Reynolds numbers OpenFOAM Workshop, Zagreb Frank Bos, Bas van Oudheusden, Hester Bijl 7-9 June 2007 1/22 Delft University of Technology
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 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 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 informationSummary of the C2.3 test case results
1st International Workshop on High-Order CFD Methods, Nashville, January 7-8, 2012 Institute of Aerodynamics and Flow Technology German Aerospace Center 8. Jan. 2012 1 / 19 Test case C2.3 Compressible
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
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 informationAn Embedded Boundary Method with Adaptive Mesh Refinements
An Embedded Boundary Method with Adaptive Mesh Refinements Marcos Vanella and Elias Balaras 8 th World Congress on Computational Mechanics, WCCM8 5 th European Congress on Computational Methods in Applied
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 informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
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 informationHigh-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson
More informationComparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils
Comparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils Moritz Kompenhans 1, Esteban Ferrer 2, Gonzalo Rubio, Eusebio Valero E.T.S.I.A. (School of Aeronautics)
More informationOverview. Applications of DEC: Fluid Mechanics and Meshing. Fluid Models (I) Part I. Computational Fluids with DEC. Fluid Models (II) Fluid Models (I)
Applications of DEC: Fluid Mechanics and Meshing Mathieu Desbrun Applied Geometry Lab Overview Putting DEC to good use Fluids, fluids, fluids geometric interpretation of classical models discrete geometric
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 informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
More 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 informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationIntegral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach
Integral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach J. Bevan, UIUC ICERM/HKUST Fast Integral Equation Methods January 5, 2016 Vorticity and Circulation Γ =
More informationShape optimisation using breakthrough technologies
Shape optimisation using breakthrough technologies Compiled by Mike Slack Ansys Technical Services 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Introduction Shape optimisation technologies
More informationLecture 1.1 Introduction to Fluid Dynamics
Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More 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 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 informationUnsteady CFD Simulations of Flow-Control Valves
Unsteady CFD Simulations of Flow-Control Valves by an Unstructured Overset-Grid Method Ph.D. dissertation presentation Jingxin Liu Advisory Committee: Dr. Akin Ecer (co-chair) Dr. Steve Frankel (co-chair)
More informationFlow Structures Extracted from Visualization Images: Vector Fields and Topology
Flow Structures Extracted from Visualization Images: Vector Fields and Topology Tianshu Liu Department of Mechanical & Aerospace Engineering Western Michigan University, Kalamazoo, MI 49008, USA We live
More informationAil implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations
Ail implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations Vimala Dutta Computational and Theoretical Fluid Dynamics Division National Aerospace
More information김태희, 권형일 1, 최성임 2* T.H. Kim, H.I. Kwon, and S.I. Choi
김태희, 권형일 1, 최성임 2* A STUDY ON INVERSE DESIGN OF AIRFOIL USING e-science BASED AERODYNAMICS DESIGN OPTIMIZATION FRAMEWORK T.H. Kim, H.I. Kwon, and S.I. Choi Recently, with advanced computational performance,
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 informationCurves: We always parameterize a curve with a single variable, for example r(t) =
Final Exam Topics hapters 16 and 17 In a very broad sense, the two major topics of this exam will be line and surface integrals. Both of these have versions for scalar functions and vector fields, and
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 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 informationThe Role of the Kutta-Joukowski Condition in the Numerical Solution of Euler Equations for a Symmetrical Airfoil.
The Role of the Kutta-Joukowski Condition in the Numerical Solution of Euler Equations for a Symmetrical Airfoil. M.Z. Dauhoo Dept. of Math, University of Mauritius, Rep. of Mauritius. E-mail : m.dauhoo@uom.ac.mu
More informationCase C1.3: Flow Over the NACA 0012 Airfoil: Subsonic Inviscid, Transonic Inviscid, and Subsonic Laminar Flows
Case C1.3: Flow Over the NACA 0012 Airfoil: Subsonic Inviscid, Transonic Inviscid, and Subsonic Laminar Flows Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts
More informationRunge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh 1
Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh Jun Zhu, inghui Zhong 3, Chi-Wang Shu 4 and Jianxian Qiu 5 Abstract In this paper we generalize a new type
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 informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationAirfoil shape optimization using adjoint method and automatic differentiation. Praveen. C
'th Annual AeSI CFD Symposium, -2 August 2009, Bangalore Airfoil shape optimization using adjoint method and automatic differentiation Praveen. C TIFR Center for Applicable Mathematics Post Bag No. 6503,
More informationThe use of gas-kinetic schemes for the simulation of compressible flows become widespread in the two last
A Gas-Kinetic BGK Scheme for Parallel Solution of 3-D Viscous Flows on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer
More informationValidation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons Stores Separation Introduction Flight Test Expensive, high-risk, sometimes catastrophic
More informationA MESH ADAPTATION METHOD FOR SIMULATION OF UNSTEADY FLOWS
A MESH ADAPTATION METHOD FOR SIMULATION OF UNSTEAD FLOWS C. H. Zhou* * Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing, 6, China Keywords: mesh adaptation, unsteady
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
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 informationSolution for Euler Equations Lagrangian and Eulerian Descriptions
Solution for Euler Equations Lagrangian and Eulerian Descriptions Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com Abstract We find an exact solution for the system of Euler equations, supposing
More informationFLOWING FLUIDS AND PRESSURE VARIATION
Chapter 4 Pressure differences are (often) the forces that move fluids FLOWING FLUIDS AND PRESSURE VARIATION Fluid Mechanics, Spring Term 2011 e.g., pressure is low at the center of a hurricane. For your
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
More informationStructure and Synthesis of Robot Motion
Structure and Synthesis of Robot Motion Dynamics: Constraints, Continua, etc. Subramanian Ramamoorthy School of Informatics 5 February, 2009 Recap Last time, we discussed two major approaches to describing
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 informationSU 2 : Advanced Analysis Topics
SU 2 : Advanced Analysis Topics OpenMDAO-SU 2 Joint Workshop Stanford University Tuesday, October 1, 2013 Thomas D. Economon Department of Aeronautics & Astronautics Stanford University 1 Three Main Topics
More informationfunctional prediction
An anisotropic hp-adaptation framework for functional prediction Marco Ceze 1 and Krzysztof J. Fidkowski 2 Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 This paper presents
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 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 informationEFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationNumerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns
More informationSimulation of Turbulent Flow in an Asymmetric Diffuser
Simulation of Turbulent Flow in an Asymmetric Diffuser 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 3 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University of Iowa C.
More informationAdaptive Cartesian Immersed Boundary Method for Simulation of Flow over Flexible Lifting Surfaces
20th AIAA Computational Fluid Dynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3384 Adaptive Cartesian Immersed Boundary Method for Simulation of Flow over Flexible Lifting Surfaces Robert
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationMesh deformation and shock capturing techniques for high-order simulation of unsteady compressible flows on dynamic meshes
AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1741 Mesh deformation and shock capturing techniques for high-order simulation of unsteady compressible
More informationSolution of final examination
of final examination Math 20, pring 201 December 9, 201 Problem 1 Let v(t) (2t e t ) i j + π cos(πt) k be the velocity of a particle with initial position r(0) ( 1, 0, 2). Find the accelaration at the
More informationIntroduction to Aerodynamic Shape Optimization
Introduction to Aerodynamic Shape Optimization 1. Aircraft Process 2. Aircraft Methods a. Inverse Surface Methods b. Inverse Field Methods c. Numerical Optimization Methods Aircraft Process Conceptual
More informationThe Divergence Theorem
The Divergence Theorem MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Summer 2011 Green s Theorem Revisited Green s Theorem: M(x, y) dx + N(x, y) dy = C R ( N x M ) da y y x Green
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 informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
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 informationComparison of h- and p- Adaptations for Spectral Difference Methods
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4435 Comparison of h- and p- Adaptations for Spectral Difference Methods Yi Li, Sachin Premasuthan and Antony
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 informationModelling the Unsteady Loads of Plunging Airfoils in Attached, Light and Deep Stall Conditions
Modelling the Unsteady Loads of Plunging Airfoils in Attached, Light and Deep Stall Conditions N. Chiereghin, D.J.Cleaver, I. Gursul, S.Bull DiPart 2017 Partnership with Contents Introduction Aims and
More informationA Local Sampling Approach to Anisotropic Metric-Based Mesh Optimization
A Local Sampling Approach to Anisotropic Metric-Based Mesh Optimization Krzysztof J. Fidkowski Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 4819, USA This work presents an
More informationExpress Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 -
More informationA Higher-Order Accurate Unstructured Finite Volume Newton-Krylov Algorithm for Inviscid Compressible Flows
A Higher-Order Accurate Unstructured Finite Volume Newton-Krylov Algorithm for Inviscid Compressible Flows by AMIR NEJAT B.Sc. (Aerospace Engineering), AmirKabir University of Technology, 1996 M.Sc. (Aerospace
More informationState of the art at DLR in solving aerodynamic shape optimization problems using the discrete viscous adjoint method
DLR - German Aerospace Center State of the art at DLR in solving aerodynamic shape optimization problems using the discrete viscous adjoint method J. Brezillon, C. Ilic, M. Abu-Zurayk, F. Ma, M. Widhalm
More informationSemi-automatic transition from simulation to one-shot optimization with equality constraints
Semi-automatic transition from simulation to one-shot optimization with equality constraints Lisa Kusch, Tim Albring, Andrea Walther, Nicolas Gauger Chair for Scientific Computing, TU Kaiserslautern, www.scicomp.uni-kl.de
More informationCase C3.1: Turbulent Flow over a Multi-Element MDA Airfoil
Case C3.1: Turbulent Flow over a Multi-Element MDA Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
More information10. Cartesian Trajectory Planning for Robot Manipulators
V. Kumar 0. Cartesian rajectory Planning for obot Manipulators 0.. Introduction Given a starting end effector position and orientation and a goal position and orientation we want to generate a smooth trajectory
More informationCase C2.2: Turbulent, Transonic Flow over an RAE 2822 Airfoil
Case C2.2: Turbulent, Transonic Flow over an RAE 2822 Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
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 information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationRunge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes
Runge Kutta discontinuous Galerkin method using WENO limiters II: unstructured meshes Jun Zhu, Jianxian Qiu,Chi-WangShu 3 and Michael Dumbser 4 Abstract In [], Qiu and Shu investigated using weighted essentially
More informationImproving High-Order Finite Element Approximation Through Geometrical Warping
Improving High-Order Finite Element Approximation Through Geometrical Warping Devina P. Sanjaya and Krzysztof J. Fidkowski Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 489,
More informationBASICS OF FLUID MECHANICS AND INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS
BASICS OF FLUID MECHANICS AND INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS Numerical Methods and Algorithms Volume 3 Series Editor: Claude Brezinski Université des Sciences et Technologies de Lille, France
More informationHigh-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder
High-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder Aerospace Application Areas Aerodynamics Subsonic through Hypersonic Aeroacoustics Store release & weapons bay analysis High lift devices
More informationAPPROACH FOR NUMERICAL MODELING OF AIRFOIL DYNAMIC STALL
APPROACH FOR NUMERICAL MODELING OF AIRFOIL DYNAMIC STALL Cvetelina Velkova, Ivan Dobrev, Michael Todorov, Fawaz Massouh To cite this version: Cvetelina Velkova, Ivan Dobrev, Michael Todorov, Fawaz Massouh.
More informationTheoretical and Computational Modeling of Dynamic Stall for Rotorcraft Applications
Theoretical and Computational Modeling of Dynamic Stall for Rotorcraft Applications Ashok Gopalarathnam and Jack Edwards Department of Mechanical and Aerospace Engineering North Carolina State 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 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 information