DNS 1-3: Difference between revisions

From KBwiki
Jump to navigation Jump to search
Line 23: Line 23:
<div id="figure1"></div>
<div id="figure1"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_fig1.png|740px]]
|[[Image:DNS_1_3_fig1.png|1000px]]
|-
|-
|'''Figure 1:''' Stanford double diffuser, Re=10000. Non-dimensional instantaneous velocity at the two middle planes using Alya with 250 Million DoF and P2 polynomials.
|'''Figure 1:''' Stanford double diffuser, Re=10000. Non-dimensional instantaneous velocity at the two middle planes using Alya with 250 Million DoF and P2 polynomials.
Line 31: Line 31:
<div id="figure2"></div>
<div id="figure2"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_sigma.png|740px]]
|[[Image:DNS_1_3_sigma.png|1000px]]
|-
|-
|'''Figure 2:''' Stanford double diffuser, Re=10000. Sigma (see EARSM variables) at the middle plane using Alya with 250 Million DoF and P2 polynomials.
|'''Figure 2:''' Stanford double diffuser, Re=10000. Sigma (see EARSM variables) at the middle plane using Alya with 250 Million DoF and P2 polynomials.

Revision as of 13:41, 1 December 2021

Flow in a 3D diffuser

Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format

Abstract

The 3D (Stanford) Diffuser is identified as a case with complex internal corner flow and 3D separation, is well documented (see, e.g., UFR_4-16_Test_Case) and has a relatively simple geometry. The diffuser has an inlet section, an expansion section and an outlet section (see Figure 1). The flow at the inlet is assumed to be a fully developed rectangular channel flow. At the outlet, standard Dirichlet condition for the pressure is prescribed. An inflow Reynolds number of 10000 is considered based on the duct height and the flow is considered to be incompressible.

This test case concerns a highly resolved Direct Numerical Simulation (DNS) using the in-house Finite Element Method (FEM) code Alya developed at BSC. The parallel multi-physics code Alya has been designed for large-scale parallel applications and uses a non-incremental fractional-step method to stabilize pressure, with an explicit fourth-order Runge-Kutta time integration scheme. The final scheme preserves momentum, angular momentum and kinetic energy at the discrete level.

The primary objective of this contribution is to provide a rich database of flow and turbulent statistics for verification and validation on subsequent computational campaigns. The provided statistical quantities in the database are:

  • mean pressure and velocities;
  • Reynolds stress components;
  • pressure autocorrelation;
  • Taylor microscale;
  • Kolmogorov length and time scales;
  • Reynolds stress equations budget terms;
  • triple velocity correlation;
  • pressure-velocity correlation.
DNS 1 3 fig1.png
Figure 1: Stanford double diffuser, Re=10000. Non-dimensional instantaneous velocity at the two middle planes using Alya with 250 Million DoF and P2 polynomials.


DNS 1 3 sigma.png
Figure 2: Stanford double diffuser, Re=10000. Sigma (see EARSM variables) at the middle plane using Alya with 250 Million DoF and P2 polynomials.




Contributed by: Oriol Lehmkuhl, Arnau Miro — Barcelona Supercomputing Center (BSC)

Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format


© copyright ERCOFTAC 2024