Direct Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers.
|
|
- Mitchell Hodge
- 5 years ago
- Views:
Transcription
1 Direct Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers. G. Borrell, J.A. Sillero and J. Jiménez, Corresponding author: School of Aeronautics, Universidad Politécnica de Madrid, 284 Madrid, Spain. Center for Turbulence Research, Stanford University, Stanford, CA 9435, USA. Abstract: A new high-resolution code for the direct numerical simulation of zero pressure gradient turbulent boundary layers over a flat plate has been developed. Its purpose is to simulate a wide range of Reynolds numbers from Re θ =3 to 68 while showing a linear weak scaling up to cores in the BG/P architecture. Special attention has been payed to the generation of proper inflow boundary conditions. The results are in good agreement with existing numerical and experimental data sets. Keywords: DNS, Boundary Layer, Turbulence. Introduction. Turbulent boundary layers have an undeniable technological interest. Roughly half of the energy dissipation due to the movement of vehicles and transport through pipes is caused by or due to the presence of a turbulent boundary layer. This is the reason why turbulent boundary layers were among the first flows to be simulated[]. Our current research is focused on understanding the flow in the turbulent regions that are far from the wall where the Reynolds number plays a significant role. Our approach is to analyze data obtained from DNS simulations. While high-quality simulations exist for other wall bounded flows, mainly channels; similar resources are not available for boundary layers yet. We decided to generate them. We needed a high resolution code that is able to perform a DNS simulation of a boundary layer that has good performance and excellent scalability. The starting point is the code described in [4] but the result has seen severe modifications. Two simulations using this code are running at this moment in two BG/P supercomputers. A zero pressure gradient boundary layer over a flat plate with Re θ = 68 is running in Intrepid at Argonne National Laboratory and a forced boundary layer with artificial roughness with Re θ = 3 42 is running in Jugene at Jülich Forschungszentrum. The aggregated computational cost of these two simulations is roughly Million CPU hours.. Computational setup. Turbulent Boundary Layers are harder to compute than other wall parallel flows like channels because they have only one homogeneous direction and they are not universal. The non universality comes from
2 the fact that the boundary conditions in the inflow and wall-normal directions have to model a semi infinite domain. A schema of the computational domain can be seen in figure. Figure : Scheme of the computational domain and boundary conditions. The x, y, and z axis correspond to the streamwise, wall-normal and spanwise directions respectively. The integration method is a primitive-variable finite difference fractional-step. The temporal integration is a low-storage third order Runge-Kutta in which wall-normal second derivative terms use a Crank-Nicholson scheme to increase the time step. The spatial discretization is a fourth-order compact finite differences in the x and y directions and spectral Fourier in the z direction. The simulation is split in two concatenated domains with different boundary conditions. The planes π i and π i are given inflow boundary conditions, and outflow boundary conditions are assigned to π e and π e. The boundary conditions in π t and π t impose a zero pressure gradient on the domain. Finally, the spanwise direction is considered periodic. The mission of the first boundary layer (BL) is to provide accurate inflow boundary conditions to the second one (BL2). The inflow of BL is obtained from its own plane π that is rescaled using a method based on the one proposed by Lund, Wu and Squires[2]. The physical length of BL is chosen long enough to let the large scales recover from an unrealistic initial condition and, once this asymptotic state has been reached, the plane π 2 is used to give BL2 its inflow boundary condition. In consequence, a small portion of the BL simulation is thrown away. Given that the goal of BL is to allow the large scales to reach their asymptotic state and, given that the smaller scales take much shorter to reach a similar condition, BL is run at a coarser resolution than BL2. This setup allows to compute a single boundary layer with significantly less computational work..2 DNS and high order differentiation schemes. The goal of the present simulation is to achieve a Reynolds number based on the friction velocity up to Re τ = 2, a measure of the difference in size between the smallest and the largest scales present in wall bounded turbulent flows. To achieve a simulation up to Re τ = 2 a computational box of size is required with a total of 45 9 points per variable. In consequence it is crucial to optimize the memory usage. Being this a DNS simulation, the separation between adjacent collocation points is determined by the accuracy of the spatial discretization scheme and the local Kolmogorov scale. This scale changes depending on the distance to the wall, hence using a non uniform mesh in the wall normal direction is essential to save memory. Another important design decision is the order of the spatial discretization. Taking into account the modified wave number analysis, fourth order compact need the half of points than, for example, second order finite differences in each spatial coordinate. While a second order discretization scheme would probably get a higher flops count, it needs about 8 times more memory and twice the time steps.
3 2 Domain decomposition and communications. To take advantage of the distributed memory architectures, the computational domain must be partitioned. The only possible decomposition that guaranteed portability to the Blue Gene/P architecture was to use cross-stream planes schematized in the figure 2 as π i Figure 2: Elemental domains of the domain decomposition. To compute interpolations and derivatives over the x coordinate it is necessary to transpose the whole variable. This operation creates another elemental domain partition formed by lines in the streamwise direction, labeled in the figure 2 as ϖ i. Once these computations are finished the result is transposed back to planes π i. This global transpose is the bottleneck of the whole computational algorithm and has been carefully tuned. 2. The Blue Gene/P architecture. The most severe restriction of the Blue Gene/P architecture is the amount of memory installed per node being only 2 GiB. Taking the largest case to be run, only two π i planes fit in a single node. Given that every node has 4 processors, OpenMP is used to take advantage of all the resources of the node efficiently. Another important benefit from this approach is that the number of MPI processes is four times smaller than a MPI-only paradigm reducing the risk of performance degradation due to latency. The OpenMP runtime implementation available in the Blue Gene/P has, in our case, an efficiency of.92 running with all the available cores. It is important to emphasize that the use of a hybrid parallelism approach in the implementation is not a performance driven decision, but a mandatory feature to overcome the limitations of the architecture. 3 Performance and scalability. All the simulations run show a linear weak scaling up to cores and the same code without modifications is expected to scale further. The figure 3a shows that the communications time is typically 44% of the total run time, and that both computation and communications are scaling as expected. It is important to mention from figure 3b that the linear scaling is kept even when the MPI message size is close to the kb limit. This is an improvement from previous experiences with other network architectures where the latency is noticed below 3kB. Another interesting feature introduced in the code is parallel I/O operations using the parallel HDF5 library, now present in almost every supercomputer. Read and write operations on the Blue Gene s gpfs scale properly with the number of I/O nodes. The collective reading has achieved aggregated speeds of 25 GB/s, but collective writing is consistently two orders of magnitude slower.
4 -2 a) 2 b) Total Com m unications Time permessage [s] -3 Time [s] Size of the message [Bytes] Millions of points pernode Figure 3: (a) Scalability of total time and communication time for different test cases. The solid lines are linear regressions computed before taking logarithms of both axis. (b) Latency analysis. The solid line is a linear regression computed before taking logarithms of both axis. 4 Conclusions and future work. The code is giving results comparable to previous DNS simulations and experiments at similar Reynolds numbers. Figures 4, a-d are the mean and fluctuations of the velocity profiles of the present simulation compared with other available experimental [5],[6] and numerical [7] data sets. The contours of the last figure compare the obtained results with a channel simulation [3]. Some more work on tuning has been explored but a small speedup is expected. For example, overlapping computation and communications won t increase the performance more than 5%. This research is funded by the DoE s Incite program, the Prace initiative, CICYT TRA and Consolider CSD27-5. References [] P. R. Spalart. (988). Direct simulation of a turbulent boundary layer up to Re θ = 4. J. Fluid Mech. 87: 6-98 [2] T. S. Lund, X. Wu and K. D. Squires. (998) Generation or turbulent inflow data for spatiallydeveloping boundary layer simulations. J. Comput. Phys., 4: [3] S. Hoyas and J. Jiménez. (26) Scaling of the velocity fluctuations in turbulent channels up to Re τ = 23. Phys.Fluids, 8:72. [4] M. P. Simens, J. Jiménez, S. Hoyas and Y. Mizuno. (29) A high-resolution code for turbulent boundary layers. J. Comput. Phys., 228: [5] D. B. De Graaf and J. K. Eaton (2) Reynolds number scaling of the flat-plate turbulent boundary layer J. Fluid Mech. 422: [6] J. M Österlund, A. V Johansson, H. M. Nagib and M. Hites (2) A note on the overlap region in turbulent boundary layers Phys. Fluids 2: -4. [7] P. Schlatter and R. Örlü. (2) Assessment of direct numerical simulation data of turbulent boundary layers. J. Fluid Mech. 659: 6-26
5 U + v + (a) y + 3 (c).5.5 w u (b) 2 y + 3 (d) 2 3 y y + 3 λ z λ x + Figure 4: Experiments by [5],, Re θ = 526; and [6],, Re θ = 556. Simulations by [7], - -, Re θ = 46; present,, Re θ = 46,526. The law log(y + )/.4+5 is the discontinuous straight line in (a). (e) Twodimensional spectral densities from channels (solid) and boundary layers (dashed), present and [8]. [8] J. Jiménez, S. Hoyas, M. P. Simens and Y. Mizuno. (2) Turbulent boundary layers and channels at moderate Reynolds numbers. J. Fluid Mech. 657:
A code for direct numerical simulation of turbulent boundary layers at high Reynolds numbers in BG/P supercomputers
A code for direct numerical simulation of turbulent boundary layers at high Reynolds numbers in BG/P supercomputers Guillem Borrell *, Juan A. Sillero, Javier Jimenez School of Aeronautics, Universidad
More informationHybrid OpenMP-MPI Turbulent boundary Layer code over 32k cores
Hybrid OpenMP-MPI Turbulent boundary Layer code over 32k cores T/NT INTERFACE y/ x/ z/ 99 99 Juan A. Sillero, Guillem Borrell, Javier Jiménez (Universidad Politécnica de Madrid) and Robert D. Moser (U.
More informationNIA CFD Futures Conference Hampton, VA; August 2012
Petascale Computing and Similarity Scaling in Turbulence P. K. Yeung Schools of AE, CSE, ME Georgia Tech pk.yeung@ae.gatech.edu NIA CFD Futures Conference Hampton, VA; August 2012 10 2 10 1 10 4 10 5 Supported
More informationPARALLEL DNS USING A COMPRESSIBLE TURBULENT CHANNEL FLOW BENCHMARK
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 ECCOMAS PARALLEL DNS USING A COMPRESSIBLE
More informationLarge Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS)
The 14 th Asian Congress of Fluid Mechanics - 14ACFM October 15-19, 2013; Hanoi and Halong, Vietnam Large Eddy Simulation of Flow over a Backward Facing Step using Fire Dynamics Simulator (FDS) Md. Mahfuz
More informationDirect numerical simulation. in an annular pipe. of turbulent flow. Paolo Luchini & Maurizio Quadrio
P.Luchini & M.Quadrio SIMAI 2000 - Ischia - 9.6.2000 Direct numerical simulation of turbulent flow in an annular pipe Paolo Luchini & Maurizio Quadrio Dipartimento di Ingegneria Aerospaziale del Politecnico
More informationLES Analysis on Shock-Vortex Ring Interaction
LES Analysis on Shock-Vortex Ring Interaction Yong Yang Jie Tang Chaoqun Liu Technical Report 2015-08 http://www.uta.edu/math/preprint/ LES Analysis on Shock-Vortex Ring Interaction Yong Yang 1, Jie Tang
More informationLARGE-EDDY EDDY SIMULATION CODE FOR CITY SCALE ENVIRONMENTS
ARCHER ECSE 05-14 LARGE-EDDY EDDY SIMULATION CODE FOR CITY SCALE ENVIRONMENTS TECHNICAL REPORT Vladimír Fuka, Zheng-Tong Xie Abstract The atmospheric large eddy simulation code ELMM (Extended Large-eddy
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 informationTransition modeling using data driven approaches
Center for urbulence Research Proceedings of the Summer Program 2014 427 ransition modeling using data driven approaches By K. Duraisamy AND P.A. Durbin An intermittency transport-based model for bypass
More informationSolidWorks Flow Simulation 2014
An Introduction to SolidWorks Flow Simulation 2014 John E. Matsson, Ph.D. SDC PUBLICATIONS Better Textbooks. Lower Prices. www.sdcpublications.com Powered by TCPDF (www.tcpdf.org) Visit the following websites
More informationAn Introduction to SolidWorks Flow Simulation 2010
An Introduction to SolidWorks Flow Simulation 2010 John E. Matsson, Ph.D. SDC PUBLICATIONS www.sdcpublications.com Schroff Development Corporation Chapter 2 Flat Plate Boundary Layer Objectives Creating
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationFlat-Plate As stated earlier, the purpose of the flat-plate study was to create an opportunity for side-by-side comparison of ÒfastÓ RNG and
Chapter Six: Comparison of Turbulence Models Performance Comparisons There are a number of criteria by which to judge turbulence models. One criterion sometimes important to mathematicallyminded model
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationTurbulent boundary layer structure identification via POD
Center for Turbulence Research Proceedings of the Summer Program 21 55 Turbulent boundary layer structure identification via POD By J. R. Baltzer, R. J. Adrian AND X. Wu Proper orthogonal decomposition
More informationDirect Numerical Simulation of a Low Pressure Turbine Cascade. Christoph Müller
Low Pressure NOFUN 2015, Braunschweig, Overview PostProcessing Experimental test facility Grid generation Inflow turbulence Conclusion and slide 2 / 16 Project Scale resolving Simulations give insight
More informationTera-scalable Fourier Spectral Element Code for DNS of Channel Turbulent Flow at High Reynolds Number
Tera-scalable Fourier Spectral Element Code for DNS of Channel Turbulent Flow at High Reynolds Number Jin Xu Physics Division Argonne National Laboratory Argonne, IL 60439 USA jin xu@anl.gov Abstract Due
More informationChapter 4: Immersed Body Flow [pp (8e), or (9e)]
Chapter 4: Immersed Body Flow [pp. 445-459 (8e), or 374-386 (9e)] Dr. Bing-Chen Wang Dept. of Mechanical Engineering Univ. of Manitoba, Winnipeg, MB, R3T 5V6 When a viscous fluid flow passes a solid body
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
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 informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More 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 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 informationOptimizing TELEMAC-2D for Large-scale Flood Simulations
Available on-line at www.prace-ri.eu Partnership for Advanced Computing in Europe Optimizing TELEMAC-2D for Large-scale Flood Simulations Charles Moulinec a,, Yoann Audouin a, Andrew Sunderland a a STFC
More informationGeneration of Turbulent Inflow Data for Spatially-Developing Boundary Layer Simulations
JOURNAL OF COMPUTATIONAL PHYSICS 140, 233 258 (1998) ARTICLE NO. CP975882 Generation of Turbulent Inflow Data for Spatially-Developing Boundary Layer Simulations Thomas S. Lund Center for Turbulence Research,
More informationCOMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE
More informationA B C D E. Settings Choose height, H, free stream velocity, U, and fluid (dynamic viscosity and density ) so that: Reynolds number
Individual task Objective To derive the drag coefficient for a 2D object, defined as where D (N/m) is the aerodynamic drag force (per unit length in the third direction) acting on the object. The object
More informationThe Level Set Method THE LEVEL SET METHOD THE LEVEL SET METHOD 203
The Level Set Method Fluid flow with moving interfaces or boundaries occur in a number of different applications, such as fluid-structure interaction, multiphase flows, and flexible membranes moving in
More informationHighly Scalable Parallel Sorting. Edgar Solomonik and Laxmikant Kale University of Illinois at Urbana-Champaign April 20, 2010
Highly Scalable Parallel Sorting Edgar Solomonik and Laxmikant Kale University of Illinois at Urbana-Champaign April 20, 2010 1 Outline Parallel sorting background Histogram Sort overview Histogram Sort
More informationTurbulent Wall-Pressure Fluctuations: New Model for Off-Axis Cross-Spectral Density
NASA Contractor Report 9897 (Corrected Copy Dec. 99) Turbulent Wall-Pressure Fluctuations: New Model for Off-Axis Cross-Spectral Density Bart A. Singer High Technology Corporation, Hampton, Virginia Contract
More informationInternational Supercomputing Conference 2009
International Supercomputing Conference 2009 Implementation of a Lattice-Boltzmann-Method for Numerical Fluid Mechanics Using the nvidia CUDA Technology E. Riegel, T. Indinger, N.A. Adams Technische Universität
More informationSimulation of Laminar Pipe Flows
Simulation of Laminar Pipe Flows 57:020 Mechanics of Fluids and Transport Processes CFD PRELAB 1 By Timur Dogan, Michael Conger, Maysam Mousaviraad, Tao Xing and Fred Stern IIHR-Hydroscience & Engineering
More informationGFD 2012 Lecture 4 Part I: Structures in 3D quasi-geostrophic flows
GFD 2012 Lecture 4 Part I: Structures in 3D quasi-geostrophic flows Jeffrey B. Weiss; notes by Alban Sauret and Vamsi Krishna Chalamalla June 21, 2012 1 Introduction Geostrophic fluid dynamics describes
More informationINTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS
INTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS Yoshinobu Yamamoto Department of Nuclear Engineering, Kyoto University Yoshida Sakyo Kyoto, 66-85,
More informationAnalyzing the Performance of IWAVE on a Cluster using HPCToolkit
Analyzing the Performance of IWAVE on a Cluster using HPCToolkit John Mellor-Crummey and Laksono Adhianto Department of Computer Science Rice University {johnmc,laksono}@rice.edu TRIP Meeting March 30,
More informationDislocation Image Stresses at Free Surfaces by the Finite Element Method
Mat. Res. Soc. Symp. Proc. Vol. 795 2004 Materials Research Society U2.4.1 Dislocation Image Stresses at Free Surfaces by the Finite Element Method Meijie Tang, Guanshui Xu*, Wei Cai, Vasily Bulatov, Lawrence
More informationResearch Directions in. Unsteady Aerodynamics
APS DFD 2004 Seattle Research Directions in Unsteady Aerodynamics P. Spalart* Boeing Commercial Airplanes *with Strelets and Squires Outline Theory: )-: Computational Fluid Dynamics Resolved and modeled
More information3D Modeling of Urban Areas for Built Environment CFD Applications
3D Modeling of Urban Areas for Built Environment CFD Applications using C A.W.M. (Jos) van Schijndel Eindhoven University of Technology P.O. Box 513; 5600 MB Eindhoven; Netherlands, A.W.M.v.Schijndel@tue.nl
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 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 informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationTurbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics
Turbulent Premixed Combustion with Flamelet Generated Manifolds in COMSOL Multiphysics Rob J.M Bastiaans* Eindhoven University of Technology *Corresponding author: PO box 512, 5600 MB, Eindhoven, r.j.m.bastiaans@tue.nl
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationMultigrid 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 informationAnalytical Modeling of Parallel Systems. To accompany the text ``Introduction to Parallel Computing'', Addison Wesley, 2003.
Analytical Modeling of Parallel Systems To accompany the text ``Introduction to Parallel Computing'', Addison Wesley, 2003. Topic Overview Sources of Overhead in Parallel Programs Performance Metrics for
More informationsimulation framework for piecewise regular grids
WALBERLA, an ultra-scalable multiphysics simulation framework for piecewise regular grids ParCo 2015, Edinburgh September 3rd, 2015 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler
More informationModeling a Nozzle in a Borehole
Modeling a Nozzle in a Borehole E. Holzbecher, F. Sun Georg-August Universität Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY; E-mail: eholzbe@gwdg.de Abstract: A nozzle, installed in an injecting
More informationNumerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas
Numerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas Avi Levy, Dept. Mech. Eng., Ben Gurion University, Beer Sheva, Israel. Mohamed Sayed, CHC, National Research
More informationP. K. Yeung Georgia Tech,
Progress in Petascale Computations of Turbulence at high Reynolds number P. K. Yeung Georgia Tech, pk.yeung@ae.gatech.edu Blue Waters Symposium NCSA, Urbana-Champaign, May 2014 Thanks... NSF: PRAC and
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 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 informationParallel Algorithm Design. Parallel Algorithm Design p. 1
Parallel Algorithm Design Parallel Algorithm Design p. 1 Overview Chapter 3 from Michael J. Quinn, Parallel Programming in C with MPI and OpenMP Another resource: http://www.mcs.anl.gov/ itf/dbpp/text/node14.html
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationRyian Hunter MAE 598
Setup: The initial geometry was produced using the engineering schematics provided in the project assignment document using the ANSYS DesignModeler application taking advantage of system symmetry. Fig.
More informationHigh Performance Computing. Introduction to Parallel Computing
High Performance Computing Introduction to Parallel Computing Acknowledgements Content of the following presentation is borrowed from The Lawrence Livermore National Laboratory https://hpc.llnl.gov/training/tutorials
More informationτ-extrapolation on 3D semi-structured finite element meshes
τ-extrapolation on 3D semi-structured finite element meshes European Multi-Grid Conference EMG 2010 Björn Gmeiner Joint work with: Tobias Gradl, Ulrich Rüde September, 2010 Contents The HHG Framework τ-extrapolation
More informationComputational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy. Jordan A.
Computational Fluid Dynamics as an advanced module of ESP-r Part 1: The numerical grid - defining resources and accuracy Jordan A. Denev Abstract: The present paper is a first one from a series of papers
More informationEffects of Buildings Layout on the Flow and Pollutant Dispersion in Non-uniform Street Canyons ZHANG Yunwei, PhD candidate
Effects of Buildings Layout on the Flow and Pollutant Dispersion in Non-uniform Street Canyons ZHANG Yunwei, PhD candidate May 11, 2010, Xi an, Shannxi Province, China E_mail: zhangyunwe@gmail.com 1 Contents:
More informationMotion Estimation for Video Coding Standards
Motion Estimation for Video Coding Standards Prof. Ja-Ling Wu Department of Computer Science and Information Engineering National Taiwan University Introduction of Motion Estimation The goal of video compression
More informationValidation of the Pressure Code using JHU DNS Database
Validation of the Pressure Code using JHU DNS Database JHU turbulence DNS database (http://turbulence.pha.jhu.edu ) The data is from a direct numerical simulation of forced isotropic turbulence on a 1024
More informationTutorial: Modeling Liquid Reactions in CIJR Using the Eulerian PDF transport (DQMOM-IEM) Model
Tutorial: Modeling Liquid Reactions in CIJR Using the Eulerian PDF transport (DQMOM-IEM) Model Introduction The purpose of this tutorial is to demonstrate setup and solution procedure of liquid chemical
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 informationHierarchical Hybrid Grids
Hierarchical Hybrid Grids IDK Summer School 2012 Björn Gmeiner, Ulrich Rüde July, 2012 Contents Mantle convection Hierarchical Hybrid Grids Smoothers Geometric approximation Performance modeling 2 Mantle
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More informationBoundary layer transition in high-speed flows due to roughness
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Eposition 9-2 January 22, Nashville, Tennessee AIAA 22-6 Boundary layer transition in high-speed flows due to roughness
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 informationAn Object-Oriented Serial and Parallel DSMC Simulation Package
An Object-Oriented Serial and Parallel DSMC Simulation Package Hongli Liu and Chunpei Cai Department of Mechanical and Aerospace Engineering, New Mexico State University, Las Cruces, New Mexico, 88, USA
More informationAdvanced ANSYS FLUENT Acoustics
Workshop Modeling Flow-Induced (Aeroacoustic) Noise 14.5 Release Advanced ANSYS FLUENT Acoustics 2011 ANSYS, Inc. November 7, 2012 1 Introduction This tutorial demonstrates how to model 2D turbulent flow
More informationSan Diego Supercomputer Center. Georgia Institute of Technology 3 IBM UNIVERSITY OF CALIFORNIA, SAN DIEGO
Scalability of a pseudospectral DNS turbulence code with 2D domain decomposition on Power4+/Federation and Blue Gene systems D. Pekurovsky 1, P.K.Yeung 2,D.Donzis 2, S.Kumar 3, W. Pfeiffer 1, G. Chukkapalli
More informationPerformance Metrics of a Parallel Three Dimensional Two-Phase DSMC Method for Particle-Laden Flows
Performance Metrics of a Parallel Three Dimensional Two-Phase DSMC Method for Particle-Laden Flows Benzi John* and M. Damodaran** Division of Thermal and Fluids Engineering, School of Mechanical and Aerospace
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 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 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 informationPressure Correction Scheme for Incompressible Fluid Flow
AALTO UNIVERSITY School of Chemical Technology CHEM-E7160 Fluid Flow in Process Units Pressure Correction Scheme for Incompressible Fluid Flow Ong Chin Kai 620503 Lee De Ming Benedict 620448 Page 1 Abstract
More informationSimulation and Validation of Turbulent Pipe Flows
Simulation and Validation of Turbulent Pipe Flows ENGR:2510 Mechanics of Fluids and Transport Processes CFD LAB 1 (ANSYS 17.1; Last Updated: Oct. 10, 2016) By Timur Dogan, Michael Conger, Dong-Hwan Kim,
More informationAteles performance assessment report
Ateles performance assessment report Document Information Reference Number Author Contributor(s) Date Application Service Level Keywords AR-4, Version 0.1 Jose Gracia (USTUTT-HLRS) Christoph Niethammer,
More informationSupplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations
Supplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations 1. Extended Results 1.1. 2-D Smoke Plume Additional results for the 2-D smoke plume example are shown in Figures
More informationEVALUATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE URBAN FLOW
EVALATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE RBAN FLOW Jarkko Saloranta and Antti Hellsten Helsinki niversity of Technology, Laboratory of Aerodynamics, Finland INTRODCTION In this work we study
More informationParallel Fourier Transform A Practical Guide
Parallel Fourier Transform A Practical Guide Dhrubaditya Mitra Indian Institute of Science, Bangalore, 560012 Parallel Fourier Transform p.1/29 Outline Motivation Serial FFT Serial FFT : Basic Algorithm
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationDevelopment of Hybrid Fluid Jet / Float Polishing Process
COMSOL Conference - Tokyo 2013 Development of Hybrid Fluid Jet / Float Polishing Process A. Beaucamp, Y. Namba Dept. of Mechanical Engineering, Chubu University, Japan Zeeko LTD, United Kingdom Research
More informationHPC Usage for Aerodynamic Flow Computation with Different Levels of Detail
DLR.de Folie 1 HPCN-Workshop 14./15. Mai 2018 HPC Usage for Aerodynamic Flow Computation with Different Levels of Detail Cornelia Grabe, Marco Burnazzi, Axel Probst, Silvia Probst DLR, Institute of Aerodynamics
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationITTC Recommended Procedures and Guidelines
Page 1 of 7 Table of Contents 1. PURPOSE OF GUIDELINE... 2 2. PARAMETERS... 2 2.1 Basic Measurement Quantities... 2 2.2 Derived Parameters... 2 3. DIRECT SIMULATION OF A TARGET WAKE FIELD... 3 4. EXPERIMENTAL
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationLarge scale Imaging on Current Many- Core Platforms
Large scale Imaging on Current Many- Core Platforms SIAM Conf. on Imaging Science 2012 May 20, 2012 Dr. Harald Köstler Chair for System Simulation Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen,
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationCloud Cavitating Flow around an Axisymmetric Projectile in the shallow water
Cloud Cavitating Flow around an Axisymmetric Projectile in the shallow water 1,2 Chang Xu; 1,2 Yiwei Wang*; 1,2 Jian Huang; 1,2 Chenguang Huang 1 Key Laboratory for Mechanics in Fluid Solid Coupling Systems,
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 informationReproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software
Reports of Research Institute for Applied Mechanics, Kyushu University No.150 (71 83) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 3: For the Case
More informationSecond Symposium on Hybrid RANS-LES Methods, 17/18 June 2007
1 Zonal-Detached Eddy Simulation of Transonic Buffet on a Civil Aircraft Type Configuration V.BRUNET and S.DECK Applied Aerodynamics Department The Buffet Phenomenon Aircraft in transonic conditions Self-sustained
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More informationAssessment of dual plane PIV measurements in wall turbulence using DNS data
Experiments in Fluids (26) 41: 265 278 DOI 1.17/s348-6-168-z RESEARCH ARTICLE Neelakantan Saikrishnan Æ Ivan Marusic Ellen K. Longmire Assessment of dual plane PIV measurements in wall turbulence using
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering
More information3D Finite Element Software for Cracks. Version 3.2. Benchmarks and Validation
3D Finite Element Software for Cracks Version 3.2 Benchmarks and Validation October 217 1965 57 th Court North, Suite 1 Boulder, CO 831 Main: (33) 415-1475 www.questintegrity.com http://www.questintegrity.com/software-products/feacrack
More informationComputational Simulation of the Wind-force on Metal Meshes
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia 2-7 December 2007 Computational Simulation of the Wind-force on Metal Meshes Ahmad Sharifian & David R. Buttsworth Faculty
More informationTopology basics. Constraints and measures. Butterfly networks.
EE48: Advanced Computer Organization Lecture # Interconnection Networks Architecture and Design Stanford University Topology basics. Constraints and measures. Butterfly networks. Lecture #: Monday, 7 April
More information