Taylor-Green Vortex LES-Smagorinsky (using Velocity Gradient)

In this example, we will investigate the creation and evolution of vortices in a simple pre-defined cube with periodic boundaries. For this case, the Reynolds number (Re) is set to 1600.

A detailed description of test case can be found in the parent directory Description of test case TGV.

A description about the LES cases can be found here: Description of test case TGV_LES.

The objectives of this example is to introduce how to: * Use a pre-defined geometry instead of creating a mesh with Seeder. * Use a turbulence model. * Simulate the Taylor-Green Vortex in the cube using Musubi. * Create 2D plots using the Gleaner tool. Gleaner is a Python tool that extracts data from Musubi ASCII output and uses the plotting library Matplotlib in Python to create a plot. * Post-process the results by calculating the volume average of the tracked kinetic energy. * Calculate the dissipation rate by means of the tracked quantities. * Compare the tracked and calculated dissipation rate with a reference solution.

Running simulation

For the Smagorinksy turbulence model, the strain rate can be computed from velocity gradient as in this example here or from the particle distribution function (PDF). For the setup of latter one see:

TGV-LES Smagorinksy from PDF

Define flow parametes

Define collision parameter

Post-processing

Here are the results from the simulation.

Kinetic energy over time compared:

Dissipation rate over time compared to reference solution from Brachet:

To create these plots, run python plot_track.py to create the plots. Before running the plot script, open 'plot_track.py' and update path to Gleaner script in 'glrPath'. Download Gleaner script using hg clone https://geb.sts.nt.uni-siegen.de/hg/gleaner

Documentation Bibliography Examples Weakly compressible flows Applications Nonreflecting Boundary Conditions Nozzle Nozzle flow inside an around Nozzle flow inside Benchmarks Channel 2D Channel 2D Boundary Conditions Channel 2D MfrBB PressExpol Channel 2D MfrEq PressEq Channel 2D PressExpol PressExpol Channel 2D VelBB PressExpol Channel 2D VelBFL PressExpol Channel 2D VelEq PressEq Channel 2D VelNonEqExpol PressNonEqExpol Poisueille flow in a channel 2D Poisueille flow in a channel 2D Poisueille flow in a channel 2D Poisueille flow in a channel 2D Poisueille flow in a channel 2D Poisueille flow in a channel 2D Poisueille flow in a channel 2D Channel 3D test cases Turbulent flow around the cylinder in a channel 3D at Re=3900 Turbulent flow around the sphere in a channel 3D Turbulent channel Turbulent channel wall model validation Turbulent channel wall model Musker Turbulent channel wall model Musker (Newton) Turbulent channel wall model Power-Law Turbulent channel wall model Reichardt Convected vortex in a domain with a coarser grid at the center Absorbing layers Acoustic Line Source Acoustic Line Source 2D Acoustic pulse 2D Acoustic pulse 2D with absorbing layer box 2d Acoustic pulse 2D with absorbing layer plane Acoustic pulse 2D with absorbing layer radial Gaussian Pulse in Pressure in a Cube Incompressible flows Applications Benchmarks Channel 2D Flow around the cylinder 2D single level Flow around the cylinder 2D single level Channel 3D Simple Flow in a 3D Channel Concentric Cylinders Concentric Cylinders Concentric Cylinders Lid driven cavity Lid driven cavity Lid driven cavity Pipe 3D Pipe Force Pipe LES Pipe 3D Pipe Simple Pipe Split Taylor Green Vortex Taylor-Green Vortex with LES Turbulence Models at Re = 1600 Taylor-Green Vortex with WALE Smagorinsky (GradU) model at Re = 1600 Taylor-Green Vortex with WALE Smagorinsky (PDF) model at Re = 1600 Taylor-Green Vortex with Vreman LES model at Re = 1600 Taylor-Green Vortex with WALE LES model at Re = 1600 Gaussian Pulse in Pressure in a Cube Tutorials Prerequisites Configuration Toolchain Tracking Boundary Conditions Restart Initial Conditions Abort Criteria Source Terms Multilevel Simulations Prerequisites Feature Overview Description of the Multi-Level Algorithm Data Structures Interpolation Methods Multispecies Musubi-Harvesting Non-Newtonian Fluid Simulations Vectorization Scheme Implementation Build and run Musubi