Transcript Visualizing Large Eddy Simulations

April 13, 2015
DNS of Spatially Developing Turbulent
Boundary Layers
PI: Antonino Ferrante, University of Washington
ECSS: Darren Adams, David Bock, NCSA
Science Overview
•
•
•
•
Computational fluid flow at boundary layers
Solution of 3D, unsteady Navier-Stokes (NS)
Generate initial inflow conditions to DNS
Velocity components must satisfy NS
Reference: http://dx.doi.org/10.1016/j.jcp.2004.01.016
2
XSEDE project
•
•
•
•
•
SU hrs and TB on Kraken, Ranger, tape storage
Simulation using largest Reynolds number ever reported
Share generated DNS database with scientific community
Perform first DNS of fully-resolved droplet-laden turbulence
ECSS consultants: Darren Adams, David Bock (NCSA)
3
ECSS goals
• Petascale scaling and deployment on XSEDE resources
• High-level HDF5 parallel IO library layer (H5DNS)
• H5DNS provides
– Data creation via petascale simulation runs
– Data analysis
– Simplified I/O data conversion (e.g. visualization)
• Visualization tests per simulation run
• Production visualization imagery (e.g. journals, conferences)
4
H5DNS software
• IO library written in C on top of HDF5 library.
• Encapsulates HDF5 File structure and metadata particular to
the DNS simulation.
• For simulation code, H5DNS provides:
–
–
–
–
Large data dumps of field variables
Restart capability
Particle IO
Metadata (run vars., input parameters etc.)
5
H5DNS software
• Analysis and visualization
– VisIt Reader (C++)
– Other direct HDF5 access (Matlab,TecPlot,…)
– Swig interfaces ??? (Perl , Python,…)
6
Visualization
Data input
•
•
•
•
•
Before H5DNS, data generation cumbersome, error-prone
Required editing and executing Fortran post-processor
Simple H5DNS reader now used to read and extract data
Converted to format used by custom visualization system
Interested in scalar Lambda2 (vortices using velocity tensors)
7
Visualization
Initial representation
• Tests verifying data conversion using slice planes
• Lambda2 (left) and velocity magnitude (right)
• Note small structures in lambda2
8
Visualization
Isolating structure
• Transfer function to isolate narrow Lambda2 range
• Evident are very small tubular structures
• Colormap with gaussian hue and alpha functions
9
Visualization
Volume rendering tests
• Ray-traced, volume renderings using same colormap
• Looking down chamber (left) and close-up side view (right)
• Note small, wispy tubular structures
10
Visualization
Volume rendering tests
• Testing different visual representations
• Top-down view (left), three-quarter side view (right)
• Added data boundary grid for reference
11
Visualization
Production sequences
• Lighting and shadows added
• Increased absorption coefficients
• Multiple viewing angles
12
Visualization
Production sequences
• Lighting and shadows added
• Increased absorption coefficients
• Multiple viewing angles
13
Visualization
Future work
• Integration of particle flow visualizing velocity
• Example of previous particle flow using custom system
14