A method to reduce load imbalances in simulations of phase change processes with FS3D
|
|
- Jocelyn Williams
- 5 years ago
- Views:
Transcription
1 A method to reduce load imbalances in simulations of phase change processes with FS3D Johannes Müller Philipp Offenhäuser Martin Reitzle Workshop on Sustained Simulation Performance HLRS, Stuttgart, Germany October 9th and 10th, 2018
2 Introduction Institute of Aerospace Thermodynamics Droplet Dynamics Group Simulations of Multiphase flows with Free Surface 3D (FS3D) Phase changes Dynamic of Free Surface Droplet Splashing Droplet Dynamics Freezing Evaporation Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 2/18
3 Motivation Phase change process: supercooled water to ice Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 3/18
4 Motivation Numerical simulation of supercooled droplet Prozess 1 Prozess 2 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 4/18
5 Motivation Numerical simulation of supercooled droplet Prozess 1 Prozess 2 Load Imbalance: Prozess 1 Prozess 2 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 4/18
6 Outline Motivation Solidification simulations with Free Surface 3D Load-Balancing method Results Conclusions Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 5/18
7 Solidification simulation Volume of Fluid method (VOF) liquid Fractional volume: f( x, t) = V 0 in the liquid phase solid = 0 < f < 1 in the boundary cells V cell 1 in the solid phase solid Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 6/18
8 Solidification simulation Solidification model Calculation of normal vectors and curvature Evaluate surface temperatur Calculate temperature gradients at interface Solve heat conduction equation Advance interface Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 7/18
9 Solidification simulation Solidification model Calculation of normal vectors and curvature Evaluate surface temperatur Calculate temperature gradients at interface Solve heat conduction equation Advance interface Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 7/18
10 Solidification simulation Parallelization Discretized computational domain: Indexing of control volumina y n j (1, 4) (1, 3) (1, 2) (2, 4) (2, 3) (2, 2) (3, 4) (3, 3) (3, 2) (4, 4) (4, 3) (4, 2) Structured, equidistant grid Data structure is based on 3-dimensional arrays Decomposition in contiguous subarrays (1, 1) (2, 1) n i (3, 1) (4, 1) x Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 8/18
11 Solidification simulation Parallelization Decomposition of computational domain in subdomains cj max cj max cj min cj min ci min ci max ci min ci max Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 8/18
12 Solidification simulation Parallelization Evaluation of a stencil cj max cj max cj min cj min ci min ci max ci min ci max Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 8/18
13 Solidification simulation Parallelization Introduction of ghost cells cj max cj max cj min cj min ci min ci max ci min ci max Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 8/18
14 Solidification simulation Parallelization Data exchange between domains cj max cj max cj min cj min ci min ci max ci min ci max Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 8/18
15 Solidification simulation Load Imbalance Symmetrical domain decomposition #2 #1 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 9/18
16 Solidification simulation Load Imbalance Additional operations in domains with surface cells Symmetrical domain decomposition 1: 2: 3: 4: 5: 6: 7: 8: #2 9: 10: #1 11: 12: 13: 14: SUBROUTINE energy jump solid Communication tr_start = MPI_Wtime() loop k = 1,NK loop j = 1,NJ loop i = 1,NI if 0 < f (i, j, k) < 1 then Calculation end if end loop end loop end loop tr_end = MPI_Wtime() Communication Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 9/18
17 Solidification simulation Load Imbalance Additional operations in domains with surface cells Symmetrical domain decomposition 1: 2: 3: 4: 5: 6: 7: 8: #2 9: 10: #1 11: 12: 13: 14: SUBROUTINE energy jump solid Communication tr_start = MPI_Wtime() loop k = 1,NK loop j = 1,NJ loop i = 1,NI if 0 < f (i, j, k) < 1 then Calculation end if end loop end loop end loop tr_end = MPI_Wtime() Communication Work Load Imbalance: t2 > t1 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes (1) 9/18
18 Load-Balancing method Optimal Work Load Compute Blade of Cray XC40: 4 compute nodes with 2 processors each containing 12 cores. Optimal workload per core Pp 1 Lopt = Lmean = n=0 LT ptotal (2) Workload per Domain LT = XGN w (3) 1 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 10/18
19 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
20 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
21 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
22 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
23 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
24 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
25 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
26 Load-Balancing method Decomposition by Bisection 1: assign cells with computational weight w e(i, j, k) << w o(i, j, k) n j n i 2: while D temp < D total do 3: Find Domain with highest Work Load 4: Copy Bisection Coordinates 5: Calculate expansion of domain 6: Bisection of Domain 7: old subdomain c max halved 8: new subdomain c min corrected 9: end while Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 11/18
27 Load-Balancing method Load Balanced Domain Decomposition Equalized Work Load Implementation of process communication: Data exchange to at least 26 neighbors (27-point stencil) Construction of neighborhood (source and destination) Construction of coordinates for send and receive arrays Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 12/18
28 Results Testcases Name Total number of cells Initial diameter % Surface cells Test Case d = d Test Case d = 2d Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 13/18
29 Results Testcases Name Total number of cells Initial diameter % Surface cells Test Case d = d Test Case d = 2d Indicator for load imbalance: % Surface cells(pmax) % Surface cells(pmean) Imbalance = % Surface cells(p mean) (4) Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 13/18
30 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
31 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
32 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance LB 1 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
33 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance LB 1 LB 2 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
34 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance LB 1 LB 2 LB 3 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
35 Results Domain decompositions for Test Case 1 Load Balancing LB 1 LB 2 LB 3 LB 4 Weight per cell w e = w o w e < w o w e << w o w e << w o Number of processes with surface cells Imbalance increases LB 1 LB 2 LB 3 Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 14/18
36 Results Evaluation of domain decompositions t SY M : Time per timestep with symmetrical domain decomposition t LB : Time per timestep with load balanced domain decomposition Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 15/18
37 Results Evaluation of domain decompositions t SY M : Time per timestep with symmetrical domain decomposition t LB : Time per timestep with load balanced domain decomposition constant number of cores Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 15/18
38 Results Strong Scaling Test Case 2 Speedup = t SYM t P (5) Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 16/18
39 Results Strong Scaling Test Case 2 Speedup = t SYM t P (5) Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 16/18
40 Conclusions Load-Balanced Domain Decomposition Nearest neighbor process communication for solidification simulation Optimized performance if percentage of boundary cells is high Imbalance only in specific Subroutines Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 17/18
41 Conclusions Load-Balanced Domain Decomposition Nearest neighbor process communication for solidification simulation Optimized performance if percentage of boundary cells is high Imbalance only in specific Subroutines Further work: Optimization of domain decomposition Investigate optimal cell weight w opt = f(percentage of boundary cells) Optimization of the complex communication pattern Implementation for other phase change processes Johannes Müller, Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart: A method to reduce load imbalances in simulations of phase change processes 17/18
42 Thank you! Johannes Müller University of Stuttgart Institute of Aerospace Thermodynamics (ITLR) Pfaffenwaldring Stuttgart
Numerical modelling of phase change processes in clouds. Challenges and Approaches. Martin Reitzle Bernard Weigand
Institute of Aerospace Thermodynamics Numerical modelling of phase change processes in clouds Challenges and Approaches Martin Reitzle Bernard Weigand Introduction Institute of Aerospace Thermodynamics
More informationA framework for detailed multiphase cloud modeling on HPC systems
Center for Information Services and High Performance Computing (ZIH) A framework for detailed multiphase cloud modeling on HPC systems ParCo 2009, 3. September 2009, ENS Lyon, France Matthias Lieber a,
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationMolecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties
12 th HLRS Results and Review Workshop Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties High Performance Computing Center Stuttgart (HLRS), October 8,
More informationHeat Transfer Equations The starting point is the conservation of mass, momentum and energy:
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 On Computational Investigation of the Supercooled Stefan Problem
More informationEarth System Modeling Domain decomposition
Earth System Modeling Domain decomposition Graziano Giuliani International Centre for Theorethical Physics Earth System Physics Section Advanced School on Regional Climate Modeling over South America February
More informationNumerical simulation of polyurethane foaming processes on bubble scale
Numerical simulation of polyurethane foaming processes on bubble scale 7th OpenFOAM Workshop Darmstadt, 25-28 June 2012 Stephanie Geier and Manfred Piesche Institute of Mechanical Process Engineering University
More informationA Computation- and Communication-Optimal Parallel Direct 3-body Algorithm
A Computation- and Communication-Optimal Parallel Direct 3-body Algorithm Penporn Koanantakool and Katherine Yelick {penpornk, yelick}@cs.berkeley.edu Computer Science Division, University of California,
More informationScalable, performant, and resilient large-scale applications of molecular process engineering
Scalable, performant, and resilient large-scale applications of molecular process engineering M. Horsch,1 P. Gralka,2 C. Niethammer,3 N. Tchipev,4 J. Vrabec,5 H. Hasse1 1 University of Kaiserslautern,
More informationPerformance of WRF using UPC
Performance of WRF using UPC Hee-Sik Kim and Jong-Gwan Do * Cray Korea ABSTRACT: The Weather Research and Forecasting (WRF) model is a next-generation mesoscale numerical weather prediction system. We
More informationParallel Programming in C with MPI and OpenMP
Parallel Programming in C with MPI and OpenMP Michael J. Quinn Chapter 13 Finite Difference Methods Outline n Ordinary and partial differential equations n Finite difference methods n Vibrating string
More informationDetermination of the composition of multicomponent droplets by rainbow refractometry
Determination of the composition of multicomponent droplets by rainbow refractometry by J. Wilms*, N. Roth*, S. Arndt**, B. Weigand* *Institute of Aerospace Thermodynamics (ITLR), University of Stuttgart
More informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationCode Timings of a Bulk Liquid Simulation
1 Parallel Molecular Dynamic Code for Large Simulations using Truncated Octahedron Periodics M.M. Micci, L.N. Long and J.K. Little a a Aerospace Department, The Pennsylvania State University, 233 Hammond
More informationIntermolecular Forces and States of Matter AP Chemistry Lecture Outline
Intermolecular Forces and States of Matter AP Chemistry Lecture Outline Name: Chemical properties are related only to chemical composition; physical properties are related to chemical composition AND the
More informationZFS - Flow Solver AIA NVIDIA Workshop
ZFS - Flow Solver AIA NVIDIA Workshop Andreas Lintermann, Matthias Meinke, Wolfgang Schröder Institute of Aerodynamics and Chair of Fluid Mechanics RWTH Aachen University Outline Softwareproject ZFS (applications
More informationParallel programming using MPI. Analysis and optimization. Bhupender Thakur, Jim Lupo, Le Yan, Alex Pacheco
Parallel programming using MPI Analysis and optimization Bhupender Thakur, Jim Lupo, Le Yan, Alex Pacheco Outline l Parallel programming: Basic definitions l Choosing right algorithms: Optimal serial and
More informationVIRTUE / WP4 Delft Twisted Foil in steady flow Simulation with EOLE
6 th Framework Programme "Sustainable development, global change and ecosystems" Project No. 516201 / WP4 Delft Twisted Foil in steady flow Simulation with EOLE Workshop Oct. 17-18, 2007 R. Marcer Principia
More informationSupporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena
Supporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena Saurabh Nath and Jonathan B. Boreyko Department of Biomedical Engineering and Mechanics, Virginia
More informationThe coupled ocean atmosphere model at ECMWF: overview and technical challenges. Kristian S. Mogensen Marine Prediction Section ECMWF
The coupled ocean atmosphere model at ECMWF: overview and technical challenges Kristian S. Mogensen Marine Prediction Section ECMWF Slide 1 Overview of talk: Baseline: The focus of this talk is going to
More informationChapter 5 - Atmospheric Moisture
Chapter 5 - Atmospheric Moisture Understanding Weather and Climate Aguado and Burt Water Water Vapor - water in a gaseous form, not droplets. Water can also achieve solid and liquid phases on Earth Temperature
More informationWorksheet 1: Integrators
Simulation Methods in Physics I WS 2014/2015 Worksheet 1: Integrators Olaf Lenz Alexander Schlaich Stefan Kesselheim Florian Fahrenberger Peter Košovan October 22, 2014 Institute for Computational Physics,
More informationPhysics 127c: Statistical Mechanics. Application of Path Integrals to Superfluidity in He 4
Physics 17c: Statistical Mechanics Application of Path Integrals to Superfluidity in He 4 The path integral method, and its recent implementation using quantum Monte Carlo methods, provides both an intuitive
More informationRefrigeration. 05/04/2011 T.Al-Shemmeri 1
Refrigeration is a process of controlled removal of heat from a substance to keep it at a temperature below the ambient condition, often below the freezing point of water (0 O C) 05/04/0 T.Al-Shemmeri
More informationMassively parallel semi-lagrangian solution of the 6d Vlasov-Poisson problem
Massively parallel semi-lagrangian solution of the 6d Vlasov-Poisson problem Katharina Kormann 1 Klaus Reuter 2 Markus Rampp 2 Eric Sonnendrücker 1 1 Max Planck Institut für Plasmaphysik 2 Max Planck Computing
More informationHideyuki Usui 1,3, M. Nunami 2,3, Y. Yagi 1,3, T. Moritaka 1,3, and JST/CREST multi-scale PIC simulation team
Hideyuki Usui 1,3, M. Nunami 2,3, Y. Yagi 1,3, T. Moritaka 1,3, and JST/CREST multi-scale PIC simulation team 1 Kobe Univ., Japan, 2 NIFS,Japan, 3 JST/CREST, Outline Multi-scale interaction between weak
More informationPerformance and Quality Analysis of Interpolation Methods for Coupling
Performance and Quality Analysis of Interpolation Methods for Coupling Neda Ebrahimi Pour, Sabine Roller Outline Motivation Coupling Approaches Interpolation/Evaluation Results of Investigation Summary
More informationCRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel?
CRYSTAL in parallel: replicated and distributed (MPP) data Roberto Orlando Dipartimento di Chimica Università di Torino Via Pietro Giuria 5, 10125 Torino (Italy) roberto.orlando@unito.it 1 Why parallel?
More informationCoupled Lattice Boltzmann and Molecular Dynamics simulations on massively parallel computers
Coupled Lattice Boltzmann and Molecular Dynamics simulations on massively parallel computers Jens Harting 1,2, Stefan Frijters 1, Florian Janoschek 1,2, and Florian Günther 1 1 Department of Applied Physics,
More informationAn Experimental Investigation on Surface Water Transport and Ice Accreting Process Pertinent to Wind Turbine Icing Phenomena
An Experimental Investigation on Surface Water Transport and Ice Accreting Process Pertinent to Wind Turbine Icing Phenomena Dr. Hui HU Advanced Flow Diagnostics and Experimental Aerodynamics Laboratory
More informationAtomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec
Atomistic molecular simulations for engineering applications: methods, tools and results Jadran Vrabec Motivation Simulation methods vary in their level of detail The more detail, the more predictive power
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationParallel Programming & Cluster Computing Transport Codes and Shifting Henry Neeman, University of Oklahoma Paul Gray, University of Northern Iowa
Parallel Programming & Cluster Computing Transport Codes and Shifting Henry Neeman, University of Oklahoma Paul Gray, University of Northern Iowa SC08 Education Program s Workshop on Parallel & Cluster
More informationA Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries
A Framework for Hybrid Parallel Flow Simulations with a Trillion Cells in Complex Geometries SC13, November 21 st 2013 Christian Godenschwager, Florian Schornbaum, Martin Bauer, Harald Köstler, Ulrich
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationScalability Metrics for Parallel Molecular Dynamics
Scalability Metrics for arallel Molecular Dynamics arallel Efficiency We define the efficiency of a parallel program running on processors to solve a problem of size W Let T(W, ) be the execution time
More informationDependence of the surface tension on curvature
Dependence of the surface tension on curvature rigorously determined from the density profiles of nanodroplets Athens, st September M. T. Horsch,,, S. Eckelsbach, H. Hasse, G. Jackson, E. A. Müller, G.
More informationChapter 18 Temperature, Heat, and the First Law of Thermodynamics. Thermodynamics and Statistical Physics
Chapter 18 Temperature, Heat, and the First Law of Thermodynamics Thermodynamics and Statistical Physics Key contents: Temperature scales Thermal expansion Temperature and heat, specific heat Heat and
More informationAccelerating Three-Body Potentials using GPUs NVIDIA Tesla K20X
Using a Hybrid Cray Supercomputer to Model Non-Icing Surfaces for Cold- Climate Wind Turbines Accelerating Three-Body Potentials using GPUs NVIDIA Tesla K20X GE Global Research Masako Yamada Opportunity
More informationTOPS Contributions to PFLOTRAN
TOPS Contributions to PFLOTRAN Barry Smith Matthew Knepley Mathematics and Computer Science Division Argonne National Laboratory TOPS Meeting at SIAM PP 08 Atlanta, Georgia March 14, 2008 M. Knepley (ANL)
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
More informationCyclops Tensor Framework
Cyclops Tensor Framework Edgar Solomonik Department of EECS, Computer Science Division, UC Berkeley March 17, 2014 1 / 29 Edgar Solomonik Cyclops Tensor Framework 1/ 29 Definition of a tensor A rank r
More informationNumerical Study of Atmospheric Ice Accretion on Various Geometric Cross-sections
Numerical Study of Atmospheric Ice Accretion on Various Geometric Cross-sections by Muhammad S. Virk REPRINTED FROM WIND ENGINEERING VOLUME 35, NO. 5, 2011 MULTI-SCIENCE PUBLISHING COMPANY 5 WATES WAY
More informationQuiz Review Topical Questions
Quiz Review Topical Questions Kinetic Theory of Matter Expansion and Contraction Solids, Liquids, Gases States of Matter Phase Changes Distillation Water Properties Kinetic Theory 1. The kinetic theory
More information1. Thermal energy is transferred through the glass windows of a house mainly by. D. radiation and convection. (1)
1. Thermal energy is transferred through the glass windows of a house mainly by A. conduction. B. radiation. C. conduction and convection. D. radiation and convection. 2. The specific latent heat of vaporization
More informationPhase-Field Modeling of Technical Alloy Systems Problems and Solutions
Thermo-Calc User Meeting, Sept 8-9 2011, Aachen Phase-Field Modeling of Technical Alloy Systems Problems and Solutions B. Böttger ACCESS e.v., RWTH Aachen University, Germany Outline 1. Introduction 2.
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationAircraft Icing Icing Physics
Aircraft Icing Icing Physics Prof. Dr. Dept. Aerospace Engineering, METU Fall 2015 Outline Formation of ice in the atmosphere Supercooled water droplets Mechanism of aircraft icing Icing variations Ice
More informationConvective Vaporization and Burning of Fuel Droplet Arrays
Convective Vaporization and Burning of Fuel Droplet Arrays Guang Wu and William A. Sirignano Mechanical & Aerospace Engineering University of California, Irvine May 19, 2009 Background and Task In many
More informationImage Stitching II. Linda Shapiro CSE 455
Image Stitching II Linda Shapiro CSE 455 RANSAC for Homography Initial Matched Points RANSAC for Homography Final Matched Points RANSAC for Homography Image Blending What s wrong? Feathering + 1 0 1 0
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationPHASE FIELD MODELS, ADAPTIVE MESH REFINEMENT AND LEVEL SETS FOR SOLIDIFICATION PROBLEMS
PHASE FIELD MODELS, ADAPTIVE MESH REFINEMENT AND LEVEL SETS FOR SOLIDIFICATION PROBLEMS Department of Physics University of Illinois at Urbana-Champaign COWORKERS AND COLLABORATORS Phase-field calculations
More informationLarge-scale MD simulation of heterogeneous systems with ls1 mardyn
Large-scale MD simulation of heterogeneous systems with ls1 mardyn M. T. Horsch, R. Srivastava, S. J. Werth, C. Niethammer, C. W. Glass, W. Eckhardt, A. Heinecke, N. Tchipev, H.-J. Bungartz, S. Eckelsbach,
More informationImage Reconstruction And Poisson s equation
Chapter 1, p. 1/58 Image Reconstruction And Poisson s equation School of Engineering Sciences Parallel s for Large-Scale Problems I Chapter 1, p. 2/58 Outline 1 2 3 4 Chapter 1, p. 3/58 Question What have
More informationArches Part 1: Introduction to the Uintah Computational Framework. Charles Reid Scientific Computing Summer Workshop July 14, 2010
Arches Part 1: Introduction to the Uintah Computational Framework Charles Reid Scientific Computing Summer Workshop July 14, 2010 Arches Uintah Computational Framework Cluster Node Node Node Node Node
More informationCompact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems
Compact High Order Finite Difference Stencils for Elliptic Variable Coefficient and Interface Problems Daniel Ritter 1, Ulrich Rüde 1, Björn Gmeiner 1, Rochus Schmid 2 Copper Mountain, March 18th, 2013
More informationNumerical Modelling in Fortran: day 8. Paul Tackley, 2017
Numerical Modelling in Fortran: day 8 Paul Tackley, 2017 Today s Goals 1. Introduction to parallel computing (applicable to Fortran or C; examples are in Fortran) 2. Finite Prandtl number convection Motivation:
More informationImage Stitching II. Linda Shapiro EE/CSE 576
Image Stitching II Linda Shapiro EE/CSE 576 RANSAC for Homography Initial Matched Points RANSAC for Homography Final Matched Points RANSAC for Homography Image Blending What s wrong? Feathering + 1 0 ramp
More informationEfficiently solving large sparse linear systems on a distributed and heterogeneous grid by using the multisplitting-direct method
Efficiently solving large sparse linear systems on a distributed and heterogeneous grid by using the multisplitting-direct method S. Contassot-Vivier, R. Couturier, C. Denis and F. Jézéquel Université
More informationWhy Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important
Why Heteroepitaxy? Difficult to Model Multiple Species Dislocations and point defects important Species Flux important The Silicon Age is Ending New materials and growth required Useful in design of Devices
More informationPerformance Analysis of a List-Based Lattice-Boltzmann Kernel
Performance Analysis of a List-Based Lattice-Boltzmann Kernel First Talk MuCoSim, 29. June 2016 Michael Hußnätter RRZE HPC Group Friedrich-Alexander University of Erlangen-Nuremberg Outline Lattice Boltzmann
More informationOutline: 1 Motivation: Domain Decomposition Method 2 3 4
Multiscale Basis Functions for Iterative Domain Decomposition Procedures A. Francisco 1, V. Ginting 2, F. Pereira 3 and J. Rigelo 2 1 Department Mechanical Engineering Federal Fluminense University, Volta
More informationRobust Design Optimization of an Axial Compressor Johannes Einzinger ANSYS Germany GmbH
Robust Design Optimization of an Axial Compressor Johannes Einzinger ANSYS Germany GmbH 1 Motivation Turbo Machines show: Rotating and stationary Parts Transient Flow Field Choke, Stall Dynamic Blade Loading
More information1. Introduction. Student Submission for the 6 th ESI OpenFOAM User Conference 2018, Hamburg - Germany
Optimizing Load Balancing of Reacting Flow Solvers in OpenFOAM for High Performance Computing Thorsten Zirwes*, Feichi Zhang, Peter Habisreuther, Jordan A. Denev, Henning Bockhorn, Dimosthenis Trimis *Karlsruhe
More informationGRAPH PARTITIONING WITH MATRIX COEFFICIENTS FOR SYMMETRIC POSITIVE DEFINITE LINEAR SYSTEMS
GRAPH PARTITIONING WITH MATRIX COEFFICIENTS FOR SYMMETRIC POSITIVE DEFINITE LINEAR SYSTEMS EUGENE VECHARYNSKI, YOUSEF SAAD, AND MASHA SOSONKINA Abstract. Prior to the parallel solution of a large linear
More informationParticle in cell simulations
Particle in cell simulations Part III: Boundary conditions and parallelization Benoît Cerutti IPAG, CNRS, Université Grenoble Alpes, Grenoble, France. 1 Astrosim, Lyon, June 26 July 7, 2017. Plan of the
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationNumerical simulation of a droplet icing on a cold surface
Numerical simulation of a droplet icing on a cold surface Buccini Saverio Peccianti Francesco Supervisors: Prof. Alessandro Bottaro Prof. Dominique Legendre March 30, 2017 Università degli Studi di Genova
More informationNumerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach
Numerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach TFM Hybrid Interface Resolving Two-Fluid Model (HIRES-TFM) by Coupling of the Volume-of-Fluid (VOF)
More informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationAccelerating Three-Body Molecular Dynamics Potentials Using NVIDIA Tesla K20X GPUs. GE Global Research Masako Yamada
Accelerating Three-Body Molecular Dynamics Potentials Using NVIDIA Tesla K20X GPUs GE Global Research Masako Yamada Overview of MD Simulations Non-Icing Surfaces for Wind Turbines Large simulations ~ 1
More informationStatic-scheduling and hybrid-programming in SuperLU DIST on multicore cluster systems
Static-scheduling and hybrid-programming in SuperLU DIST on multicore cluster systems Ichitaro Yamazaki University of Tennessee, Knoxville Xiaoye Sherry Li Lawrence Berkeley National Laboratory MS49: Sparse
More informationInfluence of Low-g Conditions on Convection and Phase Transitions in a Vessel at High Hydrostatic Pressure
High Pressure Bioscience and Biotechnology 180 Proceedings of the 4 th International Conference on High Pressure Bioscience and Biotechnology, Vol. 1, 180 185, 2007 Influence of Low-g Conditions on Convection
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationNotation. Bounds on Speedup. Parallel Processing. CS575 Parallel Processing
Parallel Processing CS575 Parallel Processing Lecture five: Efficiency Wim Bohm, Colorado State University Some material from Speedup vs Efficiency in Parallel Systems - Eager, Zahorjan and Lazowska IEEE
More information1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING
1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING M. Fiocco, D. Borghesi- Mahindra Racing S.P.A. Outline Introduction
More informationJ.I. Aliaga 1 M. Bollhöfer 2 A.F. Martín 1 E.S. Quintana-Ortí 1. March, 2009
Parallel Preconditioning of Linear Systems based on ILUPACK for Multithreaded Architectures J.I. Aliaga M. Bollhöfer 2 A.F. Martín E.S. Quintana-Ortí Deparment of Computer Science and Engineering, Univ.
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationLecture 4. Writing parallel programs with MPI Measuring performance
Lecture 4 Writing parallel programs with MPI Measuring performance Announcements Wednesday s office hour moved to 1.30 A new version of Ring (Ring_new) that handles linear sequences of message lengths
More informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationTHERMODYNAMICS NOTES. These notes give a brief overview of engineering thermodynamics. They are based on the thermodynamics text by Black & Hartley.
THERMODYNAMICS NOTES These notes give a brief overview of engineering thermodynamics. They are based on the thermodynamics text by Black & Hartley. Topics covered include: concepts; properties; conservation
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationACCESS Beiratssitzung November Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany. FENET Workshop Copenhagen
Coupling FE Contact and Heat Transfer Analysis in Casting Simulations ACCESS Beiratssitzung November 2001 Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany FENET Workshop Copenhagen 27.02.2002 About
More informationNumerical analysis of natural convection in a latent heat thermal energy storage system containing rectangular enclosures
EUROTHERM99-01-074 Numerical analysis of natural convection in a latent heat thermal energy storage system containing rectangular enclosures Julian Vogel, Maike Johnson, Markus Eck, Dörte Laing German
More informationNumerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions
Numerical Simulation of a Complete Francis Turbine including unsteady rotor/stator interactions Ruprecht, A., Heitele, M., Helmrich, T. Institute for Fluid Mechanics and Hydraulic Machinery University
More information3-D Turbulent Multiphase Modeling of Molten Steel Flow and Heat Transfer in a Continuous Slab Caster. Abstract
Presented at FX User's onference, Wilmington, DE, Oct. 1, 1998. 3-D Turbulent Multiphase Modeling of Molten Steel Flow and Heat Transfer in a ontinuous Slab aster D.T. reech and B.G. Thomas, University
More informationRadial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance
1 / 27 Radial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance S. Milovanović, L. von Sydow Uppsala University Department of Information Technology Division
More informationFIELD TEST OF WATER-STEAM SEPARATORS FOR THE DSG PROCESS
FIELD TEST OF WATER-STEAM SEPARATORS FOR THE DSG PROCESS Markus Eck 1, Holger Schmidt 2, Martin Eickhoff 3, Tobias Hirsch 1 1 German Aerospace Center (DLR), Institute of Technical Thermodynamics, Pfaffenwaldring
More informationThin Film Behavior after Ink Transfer in Printing Processes N. Bornemann, H. M. Sauer, E. Dörsam
Thin Film Behavior after Ink Transfer in Printing Processes N. Bornemann, H. M. Sauer, E. Dörsam 15.04.2010 Institute of Printing Science and Technology Thin Film Behavior N. Bornemann Overview Thin Film
More informationSIMULATION OF THE FILM FORMATION AT A HIGH-SPEED ROTARY BELL ATOMIZER USED IN AUTOMOTIVE SPRAY PAINTING PROCESSES
Paper ID ILASS08-A009 ILASS08-2-14 ILASS 2008 Sep. 8-10, 2008, Como Lake, Italy SIMULATION OF THE FILM FORMATION AT A HIGH-SPEED ROTARY BELL ATOMIZER USED IN AUTOMOTIVE SPRAY PAINTING PROCESSES J. Domnick*,
More informationHIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS
HIGH-ORDER ACCURATE METHODS BASED ON DIFFERENCE POTENTIALS FOR 2D PARABOLIC INTERFACE MODELS JASON ALBRIGHT, YEKATERINA EPSHTEYN, AND QING XIA Abstract. Highly-accurate numerical methods that can efficiently
More informationDirect numerical simulation of compressible multi-phase ow with a pressure-based method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2902 Direct numerical simulation of compressible multi-phase ow with a pressure-based
More informationParallelization of Kinetic Theory Simulations
Parallelization of Kinetic Theory Simulations Jim Howell, Wolfgang Bauer, Dirk Colbry, Rodney Pickett, Alec Staber, Irina Sagert, and Terrance Strother Abstract Numerical studies of shock waves in large
More informationEfficient Molecular Dynamics on Heterogeneous Architectures in GROMACS
Efficient Molecular Dynamics on Heterogeneous Architectures in GROMACS Berk Hess, Szilárd Páll KTH Royal Institute of Technology GTC 2012 GROMACS: fast, scalable, free Classical molecular dynamics package
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
More informationLeveraging STAR-CCM+ for Aircraft Applications. Durrell Rittenberg, Ph.D.
Leveraging STAR-CCM+ for Aircraft Applications Durrell Rittenberg, Ph.D. Overview of Icing with STAR-CCM+ Icing in aerospace Common applications Impact of icing on Aircraft safety Common icing conditions
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationNWP Equations (Adapted from UCAR/COMET Online Modules)
NWP Equations (Adapted from UCAR/COMET Online Modules) Certain physical laws of motion and conservation of energy (for example, Newton's Second Law of Motion and the First Law of Thermodynamics) govern
More information