Difference between revisions of "UFR 3-35 Best Practice Advice"

From KBwiki
Jump to navigation Jump to search
Line 10: Line 10:
 
The inflow condition is of major relevance for the vortex system. The horseshoe vortex dynamics are driven by the downflow in front of the cylinder. This downward directed flow is caused by a vertical pressure gradient, which in turn depends on the shape of the approaching inflow profile. Therefore, to study the horseshoe vortex in detail and in a generic way, the development of a fully developed turbulent open-channel or boundary layer flow approaching the cylinder has be ensured in the first place.
 
The inflow condition is of major relevance for the vortex system. The horseshoe vortex dynamics are driven by the downflow in front of the cylinder. This downward directed flow is caused by a vertical pressure gradient, which in turn depends on the shape of the approaching inflow profile. Therefore, to study the horseshoe vortex in detail and in a generic way, the development of a fully developed turbulent open-channel or boundary layer flow approaching the cylinder has be ensured in the first place.
  
A high spatial resolution is required for both the CFD and the experiment, especially to capture wall-shear. The near-wall velocity distribution in the region of interest (the wall jet) is hard to be modelled or approximated and thus has to be resolved. An evaluation of every single term of the stress balance reveals various terms to contribute significantly to the wall shear stress (Schanderl et al. 2017a).
+
A high spatial resolution is required for both the CFD and the experiment, especially to capture wall-shear. The near-wall velocity distribution in the region of interest (the wall jet) is hard to be modelled or approximated and thus has to be resolved. An evaluation of every single term of the stress balance reveals various terms to contribute significantly to the wall shear stress (Schanderl et al. 2017a). Further, at the investigated Reynolds number the dissipation of turbulent kinetic energy is not necessarily isotropic (Schanderl & Manhart 2018). Thus, the turbulence model should be chosen with care or its contribution to the momentum balance kept small.  
  
 
== Numerical Modelling Issues ==
 
== Numerical Modelling Issues ==
Line 29: Line 29:
  
 
== Physical Modelling ==
 
== Physical Modelling ==
The horseshoe vortex system is a complex three dimensional flow configuration. The two dimensional data acquisition method is, therefore, a limitation as the out-of-plane velocity component leads to a corresponding loss of particles. The number of valid samples suffered from this issue in combination of the low seeding density resulting from the large size of the flume. To overcome this issue, we additionally evaluated the PIV images with a <math>32\times32\mathrm{px}</math> grid. Whenever the instantaneous velocity fields based on a <math>16\times16\mathrm{px}</math> grid revealed a missing vector, the corresponding vector of the coarser evaluation was taken as a substitute, if possible. In this way, we could improve the number of valid samples and still keep the spatial resolution high.
+
The horseshoe vortex system is a complex three-dimensional flow configuration. Therefore, the two-dimensional data acquisition method is a limitation as the out-of-plane velocity component leads to a corresponding loss of particles. The number of valid samples suffered from this issue in combination of the low seeding density resulting from the large size of the flume. To overcome this issue, we additionally evaluated the PIV images with a <math>32\times32\mathrm{px}</math> grid. Whenever the instantaneous velocity fields based on a <math>16\times16\mathrm{px}</math> grid revealed a missing vector, the corresponding vector of the coarser evaluation was taken as a substitute, if possible. In this way, we could improve the number of valid samples and still keep the spatial resolution high.
 
However, the spatial resolution of the PIV data was too coarse to resolve the velocity gradient correctly. Therefore, a single pixel evaluation is recommended, in order to capture the wall-shear stress correctly.
 
However, the spatial resolution of the PIV data was too coarse to resolve the velocity gradient correctly. Therefore, a single pixel evaluation is recommended, in order to capture the wall-shear stress correctly.
  

Revision as of 18:52, 8 September 2019

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Cylinder-wall junction flow

Underlying Flow Regime 3-35

Best Practice Advice

Key Physics

The inflow condition is of major relevance for the vortex system. The horseshoe vortex dynamics are driven by the downflow in front of the cylinder. This downward directed flow is caused by a vertical pressure gradient, which in turn depends on the shape of the approaching inflow profile. Therefore, to study the horseshoe vortex in detail and in a generic way, the development of a fully developed turbulent open-channel or boundary layer flow approaching the cylinder has be ensured in the first place.

A high spatial resolution is required for both the CFD and the experiment, especially to capture wall-shear. The near-wall velocity distribution in the region of interest (the wall jet) is hard to be modelled or approximated and thus has to be resolved. An evaluation of every single term of the stress balance reveals various terms to contribute significantly to the wall shear stress (Schanderl et al. 2017a). Further, at the investigated Reynolds number the dissipation of turbulent kinetic energy is not necessarily isotropic (Schanderl & Manhart 2018). Thus, the turbulence model should be chosen with care or its contribution to the momentum balance kept small.

Numerical Modelling Issues

We conducted the Large Eddy Simulation using our in-house code MGLET as a finite volume code with the WALE model to represent the subgrid scale stresses without a damping function. In this model, the subgrid scale viscosity decreases towards the wall following the correct limiting behaviour of .

Discretisation method

The spatial and temporal discretization were done by a second order central differences and a thrid order Runge-Kutta method, respectively.

Grids and grid resolution

We used a staggered Cartesian grid fo the entire domain, and modelled the surface of the cylinder using a conservative second order immersed boundary method. The grid resolution in thevertical direction was using the wall-shear stress of the precursor simulation, and based on the local shear stress. The grid resolution was, therefore, sufficently fine for a LES. In case of further refinement, it would become a Direct Numerical Simulation (DNS).

Boundary conditions and computational domain

The computaitional domain containing the cylinder had a rectangular cross-section with the following size: . The bottom and side boundaries were treated with a no-slip condition, while the top boundary was modelled using a slip condition. As inflow condition a fully developed turbulent boundary layer was used stemming from a precursor simulation.

The precursor domain had a size of with periodic boundary conditions and no-slip condition at the bottom and the side walls. The surface was also modelled as slip condition.

Physical Modelling

The horseshoe vortex system is a complex three-dimensional flow configuration. Therefore, the two-dimensional data acquisition method is a limitation as the out-of-plane velocity component leads to a corresponding loss of particles. The number of valid samples suffered from this issue in combination of the low seeding density resulting from the large size of the flume. To overcome this issue, we additionally evaluated the PIV images with a grid. Whenever the instantaneous velocity fields based on a grid revealed a missing vector, the corresponding vector of the coarser evaluation was taken as a substitute, if possible. In this way, we could improve the number of valid samples and still keep the spatial resolution high. However, the spatial resolution of the PIV data was too coarse to resolve the velocity gradient correctly. Therefore, a single pixel evaluation is recommended, in order to capture the wall-shear stress correctly.

Application Uncertainties

The standard error of the mean value of the measured velocities was determined as follows:

,

the standard error of the higher central moments was obtained likewise:

.


The standard errors with respect to the standard deviation were quantified as follows:

Recommendations for Future Work

Performing a LES considering the surface waves by a level set method would improve the numerical results. Aiming to reduce the portion of modelled stresses, a further grid refinement could be implemented. The experiments could be improved by stereoscopic or tomographic PIV to acquire three dimensional data sets. Furthermore, the temporal resolution could be increased, in order to analyse the time scales of the horseshoe vortex system. The experimental set-up can be improved by providing the light sheet from below passing through the transparent bottom plate, while the PIV camera(s) are mounted at the side outside of the flume.





Contributed by: Ulrich Jenssen, Wolfgang Schanderl, Michael Manhart — Technical University Munich

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2019