UFR 2-14 Test Case: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page SilverP:UFR 2-14 Test Case to UFR 2-14 Test Case)
 
(63 intermediate revisions by 2 users not shown)
Line 1: Line 1:


= Fluid-structure interaction II =
= Fluid-structure interaction in turbulent flow past cylinder/plate configuration II (Second swiveling mode) =
{{UFRHeader
{{UFRHeader
|area=2
|area=2
Line 8: Line 8:
== Flows Around Bodies ==
== Flows Around Bodies ==
=== Underlying Flow Regime 2-14 ===
=== Underlying Flow Regime 2-14 ===


= Test Case Study =
= Test Case Study =
Line 14: Line 13:
== Description of the geometrical model and the test section ==
== Description of the geometrical model and the test section ==


FSI-PfS-2a consists of a flexible rubber structure with an attached steel weight clamped behind a fixed rigid non-rotating cylinder installed in a water channel (see Fig. 1). The experiments use the same set-up as used in [http://uriah.dedi.melbourne.co.uk/w/index.php/UFR_2-13 FSI-PfS-1a]. As a consequence all channel related parameters suchs as test section geometries, blocking ratio and the working conditions including the inflow profile remain the same as described in FSI-PfS-1a. The difference in the investigations are the changed structure definitions of FSI-PfS-2a. The deformable structure used in the experiment behind the cylinder has a slightly shorter length with <math>l_1 \operatorname{=} 0.050\,m ~ (l_1/D \approx 2.27)</math>. The attached steel weight has a length of <math>l_2 \operatorname{=} 0.010\,m ~ (l_2/D \approx 0.46)</math> and the width w that the addition of <math>l_1</math> and <math>l_2</math> yields the length <math>l\operatorname{=} 0.060\,m</math> identical to the plate of FSI-PfS-1a. The whole structure including the rigid cylinder, the flexible plate and the steel weight has a width <math>w \operatorname{=} 0.177\,m ~ (w/D \approx 8.05)</math>. Again the small gap of about <math>1.5 \times 10^{-3}\,m</math> between the side of the structure and both lateral channel walls is present. In contrast to the rubber material applied in FSI-PfS-1a the rubber used in FSI-PfS-2a has an almost constant thickness <math>h \operatorname{=} 0.002\,m ~ (h/D \approx 0.09)</math>. All parameters of the geometrical configuration of the FSI-PfS-2a benchmark are summarized as follows:
FSI-PfS-2a consists of a flexible rubber plate with an attached steel weight clamped behind a fixed rigid non-rotating cylinder installed in a water channel (see Fig. 1). The experiments use the same set-up as applied in [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13 FSI-PfS-1a]. As a consequence all channel related parameters such as the test section geometry, the blocking ratio and the working conditions including the inflow profile remain the same as described in FSI-PfS-1a. The differences of the present investigation are related to the modified structure definition of FSI-PfS-2a. The deformable structure used in the experiment behind the cylinder has a slightly shorter length with <math>l_1 \operatorname{=} 0.050\,m ~ (l_1/D \approx 2.27)</math>. The attached steel weight has a length of <math>l_2 \operatorname{=} 0.010\,m ~ (l_2/D \approx 0.46)</math>. Consequently, the summation of <math>l_1</math> and <math>l_2</math> yields the length <math>l\operatorname{=} 0.060\,m</math> identical to the plate of FSI-PfS-1a. The whole structure including the rigid cylinder, the flexible plate and the steel weight has a width <math>w \operatorname{=} 0.177\,m ~ (w/D \approx 8.05)</math>. Again small gaps of about <math>1.5 \times 10^{-3}\,m</math> are present between the lateral sides of the structure and the lateral channel walls. In contrast to the rubber material applied in FSI-PfS-1a the rubber used in FSI-PfS-2a has an almost constant thickness <math>h \operatorname{=} 0.002\,m ~ (h/D \approx 0.09)</math>. All parameters of the geometrical configuration of the FSI-PfS-2a benchmark are summarized as follows:


[[Image:qnet_FSI-PfS-2a_measures.png]]
[[Image:qnet_FSI-PfS-2a_measures.png]]
Line 20: Line 19:
[[Image:qnet_FSI-PfS-2a_geo3.png]]
[[Image:qnet_FSI-PfS-2a_geo3.png]]


Fig. 1: Geometrical configuration of the FSI-PfS-2a Benchmark.
Fig. 1: Geometrical configuration of the FSI-PfS-2a benchmark.
 
== Description of the water channel ==
 
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Description_of_the_water_channel description of UFR 2-13].


== Flow parameters ==
== Flow parameters ==
Line 26: Line 29:
Several preliminary tests were performed to find the best working conditions in terms of maximum structure displacement, good reproducibility and measurable structure frequencies within the turbulent flow regime.  
Several preliminary tests were performed to find the best working conditions in terms of maximum structure displacement, good reproducibility and measurable structure frequencies within the turbulent flow regime.  


Figure 3 depicts the measured extrema of the structure displacement versus the inlet velocity and Fig. 4 gives the frequency and Strouhal number as a function of the inlet velocity. These data were achieved by measurements with the laser distance sensor explained in Section Laser distance sensor. The entire diagrams are the result of a measurement campaign in which the inflow velocity was consecutively increased from 0 to 2.4 m/s. At an inflow velocity of <math>u_\text{inflow} = 1.385 m/s</math> (same inflow velocity as in [http://uriah.dedi.melbourne.co.uk/w/index.php/UFR_2-13 FSI-PfS-1a]) the displacements are symmetrical, reasonably large and well reproducible. Based on the inflow velocity chosen and the cylinder diameter, the Reynolds number is equal to Re = 30,470.  
Figure 2 depicts the measured extrema of the structure displacement versus the inlet velocity and Fig. 3 gives the frequency and Strouhal number as a function of the inlet velocity. These data were achieved by measurements with the laser distance sensor explained in Section [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Laser_distance_sensor Laser distance sensor]. The diagrams are the result of a measurement campaign in which the inflow velocity was consecutively increased from 0 to 2.4 m/s. At an inflow velocity of <math>u_\text{inflow} = 1.385 \, m/s </math> (same inflow velocity as in [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13 FSI-PfS-1a]) the displacements are symmetrical, reasonably large and well reproducible. Based on the inflow velocity chosen and the cylinder diameter, the Reynolds number is equal to Re = 30,470.  


[[Image:qnet_FSI-PfS-2a_dis_vel.png]]
[[Image:qnet_FSI-PfS-2a_dis_vel.png]]


Fig. 3: Experimental displacements of the structure extremity versus the inlet velocity.  
Fig. 2: Experimental displacements of the structure extremity versus the inlet velocity.  




[[Image:qnet_FSI-PfS-2a_f_vel.png]]
[[Image:qnet_FSI-PfS-2a_f_vel.png]]


Fig. 4: Experimental measurements of the frequency and the corresponding Strouhal number of the FSI phenomenon versus the inlet velocity.  
Fig. 3: Experimental measurements of the frequency and the corresponding Strouhal number of the FSI phenomenon versus the inlet velocity.  
 


Regarding the flow around the front cylinder, at this inflow velocity the flow is again in the sub-critical regime. That means the boundary layers are still laminar, but transition to turbulence takes place in the free shear layers evolving from the separated boundary layers behind the apex of the cylinder.
Regarding the flow around the front cylinder, at this inflow velocity the flow is again in the sub-critical regime. That means the boundary layers are still laminar, but transition to turbulence takes place in the free shear layers evolving from the separated boundary layers behind the apex of the cylinder.
Line 41: Line 45:
== Material Parameters ==
== Material Parameters ==


The density of the rubber material is found to be <math>\rho_{rubber plate}=1090 ~kg/m^3</math> for a thickness of the plate h = 0.002 m. This permits the accurate modeling of inertia effects of the structure and thus static and dynamic test cases can be used to calibrate the material constants (see Kalmbach and Breuer, 2013). Again the St. Venant-Kirchhoff constitutive law is chosen as the simplest hyper-elastic material model. Similar to FSI-PfS-1a, there are only two parameters to be defined: the Young's modulus E and the Poisson's ratio ν. Complementary experimental/numerical structure test studies (static, dynamic and decay test scenarios) indicate that the Young's modulus is E=3.15 MPa and the Poisson's ratio is ν=0.48 (a detailed description of the structure tests is available in Kalmbach and Breuer, 2013). The density of the steel weight is given by <math>\rho_{steel}=7850~kg/m^3</math>.
The density of the rubber material is found to be <math>\rho_{rubber}=1090 ~kg/m^3</math> for a thickness of the plate h = 0.002 m. This permits the accurate modeling of inertia effects of the structure and thus static and dynamic test cases can be used to calibrate the material constants (see Kalmbach and Breuer, 2013; De Nayer and Breuer, 2014). Again the St. Venant-Kirchhoff constitutive law is chosen as the simplest hyper-elastic material model. Similar to FSI-PfS-1a, there are only two parameters to be defined: the Young's modulus <math> E_{rubber} </math> and the Poisson's ratio <math> \nu_{rubber} </math>.  
In Kalmbach and Breuer (2013) the Young’s modulus of the rubber was determined based on a purely static test case to
<math> E_{rubber} = 3.15 ~ MPa</math>. In De Nayer and Breuer (2014) a purely dynamic test case was carried out and led to a slightly different value of <math> E_{rubber} = 4.10 ~ MPa</math>. The Poisson's ratio of the rubber is <math> \nu_{rubber} = 0.48</math>. A detailed description of all structural tests is available in Kalmbach and Breuer (2013) and De Nayer and Breuer (2014). The properties of the steel weight are given in the table below.
 
 
Structural parameters (flexible structure: para-rubber):
 
{| class="wikitable"
| Density || <math>\rho_{rubber} ~ \operatorname{=} ~ 1090 ~kg/m^3</math>
|-
| Young's modulus  || <math>E_{rubber} ~\operatorname{=} ~4.1 ~MPa</math>
|-
| Poisson's ratio || <math>\nu_{rubber} ~\operatorname{=} ~0.48</math>
|}
 


Structure Parameters:  
Structural parameters (rear mass: steel):  


{| class="wikitable"
{| class="wikitable"
| Density || <math>\rho_{rubber plate} ~ \operatorname{=} ~ 1090 ~kg/m^3</math> ||
| Density || <math>\rho_{steel} ~ \operatorname{=} ~ 7850 ~kg/m^3</math>
|-  
|-  
| Young's modulus  || <math>E ~\operatorname{=} ~3.15 ~MPa</math> ||
| Young's modulus  || <math>E_{steel} ~\operatorname{=} ~ 210 \times 10^3 ~MPa</math>  
|-
|-
| Poisson's ratio ||<math>\nu ~\operatorname{=} ~0.48</math> ||
| Poisson's ratio || <math>\nu_{steel} ~\operatorname{=} ~0.3</math>
|}
|}


Line 57: Line 75:
== Particle-image velocimetry ==
== Particle-image velocimetry ==


[http://uriah.dedi.melbourne.co.uk/w/index.php/UFR_2-13_Test_Case#Particle-image_velocimetry]
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Particle-image_velocimetry description for UFR 2-13].
 
== Laser distance sensor ==
 
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Laser_distance_sensor description for UFR 2-13].
 
= Numerical Simulation Methodology =


== Computational fluid dynamics (CFD) ==


== Laser distance sensor ==
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Computational_fluid_dynamics_.28CFD.29 description for UFR 2-13].
 
== Computational structural dynamics (CSD) ==
 
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Computational_structural_dynamics_.28CSD.29 description for UFR 2-13].
 
== Coupling algorithm ==
 
See [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Coupling_algorithm description for UFR 2-13].
 
== Numerical CFD Setup ==
 
Preliminary tests were conducted for the geometrically similar FSI-PfS-1a test case to identify the appropriate CFD domain size (see [http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13 FSI-PfS-1]). The outcome was that a consideration of a part of the whole test section combined with periodic boundary conditions on the lateral boundaries is sufficient to get reliable
flow results. This configuration is denoted ''subset case''. The geometry of the current FSI-PfS-2a test case is similar to
FSI-PfS-1a. Owing to the rear mass in the present case avoiding three-dimensional deformations of the flexible structure, the
restriction to the subset case is even more justified than for FSI-PfS-1a. Therefore, the subset case is directly applied here: The
depth of the computational domain is not the entire test section width <math>W/D \approx 8.18</math>, but a subset of the total width equal to the
entire length of the splitter plate <math>\left( l_1 + l_2 \right)/D \approx 2.73</math> yielding a quadratic plate. A block-structured grid is
generated consisting of 91 blocks involving about 14 million control volumes (CVs). Note that this resolution is chosen based on the
experiences made for the FSI-PfS-1a test case (de Nayer et al., 2014). Figure 4 shows the x-y cross-section of the grid. Since only every fourth grid line of
the mesh is shown in Fig. 4, the angles between the grid lines and the transitions between the blocks appear to be worse than they are in the full grid. In order to fully resolve
the viscous sublayer on the elastic structure, the first cell center is located at a distance of <math>\Delta y / D \approx 6.8 \times 10^{-4}</math> from the flexible structure, which leads to average <math>y^+</math> values below 0.8, mostly even below 0.5 near the body. The geometrical stretching factors are kept below 1.1. In spanwise
direction the grid consists of 72 equidistant cells. The inflow side is rounded in order to use a C-grid. Consequently, the computational
domain in front of the cylinder is slightly larger than in the test section depicted in Fig. 1. On the right side
a convective outlet boundary condition is applied. On the cylinder and on the flexible structure no-slip walls are defined. The top and the
bottom of the domain are relatively far away from the flexible structure and are thus taken into account via slip walls. The lateral
sides are assumed to be periodic as mentioned above.
 
The numerical method developed is based on an explicit time-marching scheme involving the low-storage Runge Kutta algorithm described
above. Hence a small time step of <math>\Delta t = 10^{-5}\,s</math> is used. This time step size corresponds to a CFL number in the order of
unity.
 
The subgrid-scale (SGS) model chosen for the main simulation is the Smagorinsky model (Smagorinsky, 1963). It is applied with the well established
standard constant <math>C_s = 0.1</math> and a Van Driest damping function near solid walls. Owing to the moderate Reynolds number considered and the
fine grid applied, the SGS model is expected to have a limited influence on the results. Nevertheless, in order to investigate and
verify this issue, FSI simulations using the dynamic procedure as suggested by Germano et al. (1991) and Lilly (1992) and the WALE model
by Nicoud and Ducros (1999) are carried out and analyzed in Section [http://qnet-ercoftac.cfms.org.uk/w/index.php/Lib:UFR_2-14_Evaluation#Sensitivity_study Sensitivity study]. That includes a study on the
model parameters and their influence on the results.
 
[[File:Benchmark_FSI-PfS-2a_mesh_coarser_4.png]]
 
Fig. 4: x-y cross-section of the grid used for the simulation with 197,136 CVs (Note that only every fourth grid line in each direction is displayed here).
 
== Numerical CSD Setup ==
 
The thin flexible structure and the rear mass are modeled by
quadrilateral four-nodes shell elements. These elements are special
seven-parameters shell elements (Bischoff, 1999). Preliminary
structural tests show that <math>30 \times 10</math> shell elements for
the flexible para-rubber and <math>5 \times 10</math> shell elements for
the stainless rear mass are sufficient to deliver accurate
results. The well-known Assumed Natural Strain (ANS)
(Hughes and Tezduyar (1981), Park and Stanley (1986)) is deactivated in the current case and the
Enhanced Assumed Strain (EAS) method is set to (M,B,T,Q,S)=(4,0,4,4,0)
avoiding locking phenomena as proposed by Bischoff et al. (2004).
 
The flexible structure is clamped to the cylinder: All the nodes at
the shell edge in contact with the cylinder have zero degree of
freedom (DoF). On the opposite downstream trailing-edge, the rear mass
is free to move and all nodes have the full set of six degrees of
freedom. The nodes on the lateral sides have to fit to the periodic
CFD boundary conditions: The x-displacements, the y-displacements and
the rotations have to be the same for a lateral node and its opposite
counterpart. Moreover, the periodic boundary conditions imply that the
z-displacement of the nodes on the sides are forced to be zero. This
special treatment of the lateral nodes is explained in detail and
validated in De Nayer et al. (2014) for the FSI-PfS-1a test case.
 
A FSI sub-cycling algorithm is not implemented. Therefore, the same time
step as for the CFD solver is set.
 
== Coupling conditions ==
 
In accordance with previous laminar and turbulent test cases in Breuer et al. (2012) and De Nayer et al. (2014) a linear extrapolation of the structural
deformation is used at the beginning of any new time step to
accelerate the convergence. A second-order extrapolation is not taken
into account, since it was observed that it does not improve the
convergence for the current small time step size. Furthermore, a
constant underrelaxation factor of <math>\omega = 0.5</math> is considered for
the displacements during the FSI-subiterations to stabilize the
coupling. The forces are exchanged without underrelaxation. The FSI
convergence criterion is set to <math>\varepsilon_\text{FSI} = 10^{-4}</math> for
the <math>L_2</math> norm of the displacement differences between two
FSI-subiterations normalized by the displacement differences between
the current and the old time step. With these settings about 5
FSI-subiterations are required to reach the convergence
criterion.
 
For the sake of clarity, the whole computational setup is summarized below.
 
[[File:FSI_setup.png]]
 
 
For the current FSI simulation the flow field is initialized by assuming that the entire structure
is non-deformable. For this rigid configuration the turbulent flow is predicted until it reaches
a quasi-periodic state. In this case the shell attached to the backside of the cylinder acts
like a splitter plate attenuating the generation of a von Karman vortex street behind the
cylinder. Nevertheless, quasi-periodic vortex shedding is still observed with a Strouhal number
of <math>\text{St}_\text{fixed} \approx 0.175 </math>. Then, the structure is released and the FSI simulation begins.
A time interval of about 2 seconds physical time is computed. The structure requires a few cycles until a
statistically steady-state is reached. After that a phase-averaging is carried out.
 
= Generation of Phase-resolved Data =
 
The results obtained by the experiment and the simulation cannot be directly compared.
Indeed, due to turbulence, the data contain a random chaotic part. This contribution has
to be filtered out in order to obtain phase-averaged data which can be compared directly. A
special procedure, called phase-averaging method, has been explained in detail in Kalmbach
and Breuer (2013) to treat both, the experimental and the numerical data sets. For this procedure
a representative signal of the FSI problem is needed. In the present case the deformations of
the plate are quasi two-dimensional due to the rear mass. Therefore, the y-displacement of
the flexible structure at the monitoring point close to the trailing edge is chosen for this purpose.
The current FSI phenomenon is quasi-periodic. This quasi-periodicity can also be confirmed by looking at the predicted
phase plane which shows a typical form of a ”8” for the second swiveling
mode.


[[http://uriah.dedi.melbourne.co.uk/w/index.php/UFR_2-13_Test_Case#Laser_distance_sensor]]
Each period found is divided into n equidistant sub-parts (here n = 23). The sub-part j of the
period <math>T_i </math> corresponds to the sub-part j of the period <math>T_{i+1}</math> and so on. Each data set found in
a sub-part j is averaged with the other data sets found in the sub-parts j of all other periods.
Finally, data sets of n = 23 phase-averaged positions for the representative reference period
are achieved. For the experiment, about 47 particle-image velocimetry (PIV) data sets are
averaged for each of the 23 structure positions. For the simulation, about 40 numerical data
sets are averaged for each of the 23 structure positions. Note that the numerical data sets are
also averaged in spanwise direction, in order to get better phase-averaged mean values. That
is the main reason why the phase-averaging procedure for the measurements comprises about
150 periods, whereas for the numerical results only 15 periods are taken into account.


For a full description of the phase averaging procedure we refer to the description provided for
[http://qnet-ercoftac.cfms.org.uk/w/index.php/UFR_2-13_Test_Case#Generation_of_Phase-resolved_Data UFR 2-13].
----
----
{{ACContribs
{{ACContribs

Latest revision as of 12:13, 12 February 2017

Fluid-structure interaction in turbulent flow past cylinder/plate configuration II (Second swiveling mode)

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References

Flows Around Bodies

Underlying Flow Regime 2-14

Test Case Study

Description of the geometrical model and the test section

FSI-PfS-2a consists of a flexible rubber plate with an attached steel weight clamped behind a fixed rigid non-rotating cylinder installed in a water channel (see Fig. 1). The experiments use the same set-up as applied in FSI-PfS-1a. As a consequence all channel related parameters such as the test section geometry, the blocking ratio and the working conditions including the inflow profile remain the same as described in FSI-PfS-1a. The differences of the present investigation are related to the modified structure definition of FSI-PfS-2a. The deformable structure used in the experiment behind the cylinder has a slightly shorter length with . The attached steel weight has a length of . Consequently, the summation of and yields the length identical to the plate of FSI-PfS-1a. The whole structure including the rigid cylinder, the flexible plate and the steel weight has a width . Again small gaps of about are present between the lateral sides of the structure and the lateral channel walls. In contrast to the rubber material applied in FSI-PfS-1a the rubber used in FSI-PfS-2a has an almost constant thickness . All parameters of the geometrical configuration of the FSI-PfS-2a benchmark are summarized as follows:

Qnet FSI-PfS-2a measures.png

Qnet FSI-PfS-2a geo3.png

Fig. 1: Geometrical configuration of the FSI-PfS-2a benchmark.

Description of the water channel

See description of UFR 2-13.

Flow parameters

Several preliminary tests were performed to find the best working conditions in terms of maximum structure displacement, good reproducibility and measurable structure frequencies within the turbulent flow regime.

Figure 2 depicts the measured extrema of the structure displacement versus the inlet velocity and Fig. 3 gives the frequency and Strouhal number as a function of the inlet velocity. These data were achieved by measurements with the laser distance sensor explained in Section Laser distance sensor. The diagrams are the result of a measurement campaign in which the inflow velocity was consecutively increased from 0 to 2.4 m/s. At an inflow velocity of (same inflow velocity as in FSI-PfS-1a) the displacements are symmetrical, reasonably large and well reproducible. Based on the inflow velocity chosen and the cylinder diameter, the Reynolds number is equal to Re = 30,470.

Qnet FSI-PfS-2a dis vel.png

Fig. 2: Experimental displacements of the structure extremity versus the inlet velocity.


Qnet FSI-PfS-2a f vel.png

Fig. 3: Experimental measurements of the frequency and the corresponding Strouhal number of the FSI phenomenon versus the inlet velocity.


Regarding the flow around the front cylinder, at this inflow velocity the flow is again in the sub-critical regime. That means the boundary layers are still laminar, but transition to turbulence takes place in the free shear layers evolving from the separated boundary layers behind the apex of the cylinder.

Material Parameters

The density of the rubber material is found to be for a thickness of the plate h = 0.002 m. This permits the accurate modeling of inertia effects of the structure and thus static and dynamic test cases can be used to calibrate the material constants (see Kalmbach and Breuer, 2013; De Nayer and Breuer, 2014). Again the St. Venant-Kirchhoff constitutive law is chosen as the simplest hyper-elastic material model. Similar to FSI-PfS-1a, there are only two parameters to be defined: the Young's modulus and the Poisson's ratio . In Kalmbach and Breuer (2013) the Young’s modulus of the rubber was determined based on a purely static test case to . In De Nayer and Breuer (2014) a purely dynamic test case was carried out and led to a slightly different value of . The Poisson's ratio of the rubber is . A detailed description of all structural tests is available in Kalmbach and Breuer (2013) and De Nayer and Breuer (2014). The properties of the steel weight are given in the table below.


Structural parameters (flexible structure: para-rubber):

Density
Young's modulus
Poisson's ratio


Structural parameters (rear mass: steel):

Density
Young's modulus
Poisson's ratio

Measuring Technique

Particle-image velocimetry

See description for UFR 2-13.

Laser distance sensor

See description for UFR 2-13.

Numerical Simulation Methodology

Computational fluid dynamics (CFD)

See description for UFR 2-13.

Computational structural dynamics (CSD)

See description for UFR 2-13.

Coupling algorithm

See description for UFR 2-13.

Numerical CFD Setup

Preliminary tests were conducted for the geometrically similar FSI-PfS-1a test case to identify the appropriate CFD domain size (see FSI-PfS-1). The outcome was that a consideration of a part of the whole test section combined with periodic boundary conditions on the lateral boundaries is sufficient to get reliable flow results. This configuration is denoted subset case. The geometry of the current FSI-PfS-2a test case is similar to FSI-PfS-1a. Owing to the rear mass in the present case avoiding three-dimensional deformations of the flexible structure, the restriction to the subset case is even more justified than for FSI-PfS-1a. Therefore, the subset case is directly applied here: The depth of the computational domain is not the entire test section width , but a subset of the total width equal to the entire length of the splitter plate yielding a quadratic plate. A block-structured grid is generated consisting of 91 blocks involving about 14 million control volumes (CVs). Note that this resolution is chosen based on the experiences made for the FSI-PfS-1a test case (de Nayer et al., 2014). Figure 4 shows the x-y cross-section of the grid. Since only every fourth grid line of the mesh is shown in Fig. 4, the angles between the grid lines and the transitions between the blocks appear to be worse than they are in the full grid. In order to fully resolve the viscous sublayer on the elastic structure, the first cell center is located at a distance of from the flexible structure, which leads to average values below 0.8, mostly even below 0.5 near the body. The geometrical stretching factors are kept below 1.1. In spanwise direction the grid consists of 72 equidistant cells. The inflow side is rounded in order to use a C-grid. Consequently, the computational domain in front of the cylinder is slightly larger than in the test section depicted in Fig. 1. On the right side a convective outlet boundary condition is applied. On the cylinder and on the flexible structure no-slip walls are defined. The top and the bottom of the domain are relatively far away from the flexible structure and are thus taken into account via slip walls. The lateral sides are assumed to be periodic as mentioned above.

The numerical method developed is based on an explicit time-marching scheme involving the low-storage Runge Kutta algorithm described above. Hence a small time step of is used. This time step size corresponds to a CFL number in the order of unity.

The subgrid-scale (SGS) model chosen for the main simulation is the Smagorinsky model (Smagorinsky, 1963). It is applied with the well established standard constant and a Van Driest damping function near solid walls. Owing to the moderate Reynolds number considered and the fine grid applied, the SGS model is expected to have a limited influence on the results. Nevertheless, in order to investigate and verify this issue, FSI simulations using the dynamic procedure as suggested by Germano et al. (1991) and Lilly (1992) and the WALE model by Nicoud and Ducros (1999) are carried out and analyzed in Section Sensitivity study. That includes a study on the model parameters and their influence on the results.

Benchmark FSI-PfS-2a mesh coarser 4.png

Fig. 4: x-y cross-section of the grid used for the simulation with 197,136 CVs (Note that only every fourth grid line in each direction is displayed here).

Numerical CSD Setup

The thin flexible structure and the rear mass are modeled by quadrilateral four-nodes shell elements. These elements are special seven-parameters shell elements (Bischoff, 1999). Preliminary structural tests show that shell elements for the flexible para-rubber and shell elements for the stainless rear mass are sufficient to deliver accurate results. The well-known Assumed Natural Strain (ANS) (Hughes and Tezduyar (1981), Park and Stanley (1986)) is deactivated in the current case and the Enhanced Assumed Strain (EAS) method is set to (M,B,T,Q,S)=(4,0,4,4,0) avoiding locking phenomena as proposed by Bischoff et al. (2004).

The flexible structure is clamped to the cylinder: All the nodes at the shell edge in contact with the cylinder have zero degree of freedom (DoF). On the opposite downstream trailing-edge, the rear mass is free to move and all nodes have the full set of six degrees of freedom. The nodes on the lateral sides have to fit to the periodic CFD boundary conditions: The x-displacements, the y-displacements and the rotations have to be the same for a lateral node and its opposite counterpart. Moreover, the periodic boundary conditions imply that the z-displacement of the nodes on the sides are forced to be zero. This special treatment of the lateral nodes is explained in detail and validated in De Nayer et al. (2014) for the FSI-PfS-1a test case.

A FSI sub-cycling algorithm is not implemented. Therefore, the same time step as for the CFD solver is set.

Coupling conditions

In accordance with previous laminar and turbulent test cases in Breuer et al. (2012) and De Nayer et al. (2014) a linear extrapolation of the structural deformation is used at the beginning of any new time step to accelerate the convergence. A second-order extrapolation is not taken into account, since it was observed that it does not improve the convergence for the current small time step size. Furthermore, a constant underrelaxation factor of is considered for the displacements during the FSI-subiterations to stabilize the coupling. The forces are exchanged without underrelaxation. The FSI convergence criterion is set to for the norm of the displacement differences between two FSI-subiterations normalized by the displacement differences between the current and the old time step. With these settings about 5 FSI-subiterations are required to reach the convergence criterion.

For the sake of clarity, the whole computational setup is summarized below.

FSI setup.png


For the current FSI simulation the flow field is initialized by assuming that the entire structure is non-deformable. For this rigid configuration the turbulent flow is predicted until it reaches a quasi-periodic state. In this case the shell attached to the backside of the cylinder acts like a splitter plate attenuating the generation of a von Karman vortex street behind the cylinder. Nevertheless, quasi-periodic vortex shedding is still observed with a Strouhal number of . Then, the structure is released and the FSI simulation begins. A time interval of about 2 seconds physical time is computed. The structure requires a few cycles until a statistically steady-state is reached. After that a phase-averaging is carried out.

Generation of Phase-resolved Data

The results obtained by the experiment and the simulation cannot be directly compared. Indeed, due to turbulence, the data contain a random chaotic part. This contribution has to be filtered out in order to obtain phase-averaged data which can be compared directly. A special procedure, called phase-averaging method, has been explained in detail in Kalmbach and Breuer (2013) to treat both, the experimental and the numerical data sets. For this procedure a representative signal of the FSI problem is needed. In the present case the deformations of the plate are quasi two-dimensional due to the rear mass. Therefore, the y-displacement of the flexible structure at the monitoring point close to the trailing edge is chosen for this purpose. The current FSI phenomenon is quasi-periodic. This quasi-periodicity can also be confirmed by looking at the predicted phase plane which shows a typical form of a ”8” for the second swiveling mode.

Each period found is divided into n equidistant sub-parts (here n = 23). The sub-part j of the period corresponds to the sub-part j of the period and so on. Each data set found in a sub-part j is averaged with the other data sets found in the sub-parts j of all other periods. Finally, data sets of n = 23 phase-averaged positions for the representative reference period are achieved. For the experiment, about 47 particle-image velocimetry (PIV) data sets are averaged for each of the 23 structure positions. For the simulation, about 40 numerical data sets are averaged for each of the 23 structure positions. Note that the numerical data sets are also averaged in spanwise direction, in order to get better phase-averaged mean values. That is the main reason why the phase-averaging procedure for the measurements comprises about 150 periods, whereas for the numerical results only 15 periods are taken into account.

For a full description of the phase averaging procedure we refer to the description provided for UFR 2-13.



Contributed by: Andreas Kalmbach, Guillaume De Nayer, Michael Breuer — Helmut-Schmidt Universität Hamburg

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References


© copyright ERCOFTAC 2024