BAIL 2006 Book of Abstracts - Institut für Numerische und ...

BAIL 2006 Book of Abstracts - Institut für Numerische und ... BAIL 2006 Book of Abstracts - Institut für Numerische und ...

num.math.uni.goettingen.de
from num.math.uni.goettingen.de More from this publisher
04.02.2013 Views

F. ALIZARD, J.-CH. ROBINET: Two-dimensional temporal modes in nonparallel flows ✬ ✫ Figure 1: Discrete global spectrum from a convectively instable boundary layer. The real part of the eigenvalue is represented on the Xaxis and the imaginary part on the Y-axis. Figure 2: Real part of the horizontal velocity eigenfunction of a discrete eigenvalue from the spectrum ( Fig. 1). with U(x, y) =[U, V ] t and û(x, y) =[û, ˆv] t . The boundary conditions on velocity disturbances and constraints on pressure (T. N. Phillips & G. W. Roberts (1993), C. Canuto et al. (1988)) complete the eigenvalue value problem. The partial differential stability equations (2) are discretized using an algorithm based on the collocation method based on Chebyshev Gauss-Lobatto grid. The algebraic eigenvalue problem (A − ΩB) X = 0 is solved by the QZ algorithm. 4. Two-Dimensional Temporal Modes, results and perspectives The figure 1 represents a global spectrum of the boundary layer. The stable discrete eigenvalues appear in concordance with the fact a boundary layer is globally stable. These discrete values represent spatio-temporal convective modes as it can been shown fig. 2 (structure similar to a spatial exponential growth). At the conference, this analysis will be also applied on a separated incompressible boundary layer which is not absolutely unstable but only convectively unstable. References [1] M. Gaster , “The development of a two-dimensional wavepacket in a growing boundary layer ”, Proc. R. Soc. London, 3, 317–332 (1982). [2] Lin R.S. and Malik M.R , “On the stability of attachment-line boundary layers. ”, J. Fluid Mech., 3,239–255 (1995). [3] Theofilis V. Hein S and Dallmann U , “On the origins of unsteadiness and three dimensionality in a laminar separation bubble ”, Proc. R. Soc. London, 3,3229–3246 (2000). [4] U. Ehrenstein and F. Gallaire , “On two dimensional temporal modes in spatially evolving open flows: the flat-plate boundary layer ”, J. Fluid Mech., (2005). [5] Briley WR., “A numerical study of laminar separation bubbles using the Navier-Stokes equations.”, J. Fluid Mech., 3,713–736 (1971). [6] Timothy N. Phillips and Gareth W.Roberts, “The Treatment of Spurious Pressure Modes in Spectral Incompressible Flow Calculations.”, Journal of Computational Physics, 3,150–164 (1993). [7] C. Canuto and M. Y. Hussaini and A. Quarteroni and T. A. Zang, “Spectral Methods in Fluids Dynamics”, Springer, (1988). Speaker: ALIZARD, F. 68 BAIL 2006 2 ✩ ✪

TH. ALRUTZ, T. KNOPP: Near-wall grid adaptation for wall functions ✬ ✫ Near-wall grid adaptation for wall-functions Th. Alrutz, T. Knopp Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Bunsenstr. 10, 37073 Göttingen, Germany Thomas.Alrutz@dlr.de A near-wall grid adaptation method for RANS turbulence modelling using wall-functions is proposed. Universal wall functions allow a considerable acceleration of the flow solver, but their predictions may deviate from the results with near-wall integration in flow situations which are beyond the underlying modelling assumptions of wall-functions. In aerodynamic flows these include (i) stagnation points (flow impingement) and subsequent transition from laminar to turbulent flow, (ii) large surface curvature in conjunction with strong pressure gradients (as can be observed around leading edge of an airfoil), (iii) strong adverse pressure gradients leading to separation, and (iv) regions of separated flow. Boundary layer theory shows that the deviation from the universal wall-law becomes more significant as the walldistance of the first node above the wall is increased. But close agreement with the universal wall-law can be achived if the first node is close enough to the wall. To reduce the modelling error of universal wall-functions, a near-wall grid-adaptation technique is pro- Ýposed. �Ý Regions of strong surface curvature, large pressure gradients and recirculating flow are detected by a flow and geometry based sensor. Then in critical regions, the nodes of the prismatic nearwall grid are shifted towards the wall so that a user-specified target value forÝ is ensured, where Ù���with wall-distance of the first nodeÝ , friction velocityÙ�and viscosity�. This approach is applied to a transonic airfoil flow with shock induced separation [5] and to a subsonic highlift airfoil close to stall [6]. Thereby, the predictions around the leading edge (suction peak) and of the separation point and the recirculation region can be improved significantly. y + 100 80 60 40 20 y + (1) = 20 y + (1) = 40 y + (1) = 60 0 0 0.25 0.5 0.75 1 x/c y + 60 40 20 y + (1) = 20 y + (1) = 40 y + (1) = 60 0 0 0.25 0.5 0.75 1 x/c Figure 1: RAE2822 case 10 [5]: Distribution ofÝ for Menter SST model [4] without adaptation (left) and withÝ-adaptation of the structured (prismatic) near-wall grid (right). References [1] T. Knopp, T. Alrutz and D. Schwamborn, “A grid and flow adaptive wall-function method for RANS turbulence modelling”, Journal of Computational Physics (submitted). Speaker: ALRUTZ, TH 69 BAIL 2006 1 ✩ ✪

F. ALIZARD, J.-CH. ROBINET: Two-dimensional temporal modes in nonparallel<br />

flows<br />

✬<br />

✫<br />

Figure 1: Discrete global spectrum from a convectively<br />

instable bo<strong>und</strong>ary layer. The real<br />

part <strong>of</strong> the eigenvalue is represented on the Xaxis<br />

and the imaginary part on the Y-axis.<br />

Figure 2: Real part <strong>of</strong> the horizontal velocity<br />

eigenfunction <strong>of</strong> a discrete eigenvalue from the<br />

spectrum ( Fig. 1).<br />

with U(x, y) =[U, V ] t and û(x, y) =[û, ˆv] t . The bo<strong>und</strong>ary conditions on velocity disturbances<br />

and constraints on pressure (T. N. Phillips & G. W. Roberts (1993), C. Canuto et al. (1988))<br />

complete the eigenvalue value problem. The partial differential stability equations (2) are discretized<br />

using an algorithm based on the collocation method based on Chebyshev Gauss-Lobatto<br />

grid. The algebraic eigenvalue problem (A − ΩB) X = 0 is solved by the QZ algorithm.<br />

4. Two-Dimensional Temporal Modes, results and perspectives<br />

The figure 1 represents a global spectrum <strong>of</strong> the bo<strong>und</strong>ary layer. The stable discrete eigenvalues<br />

appear in concordance with the fact a bo<strong>und</strong>ary layer is globally stable. These discrete values<br />

represent spatio-temporal convective modes as it can been shown fig. 2 (structure similar to a<br />

spatial exponential growth). At the conference, this analysis will be also applied on a separated<br />

incompressible bo<strong>und</strong>ary layer which is not absolutely unstable but only convectively unstable.<br />

References<br />

[1] M. Gaster , “The development <strong>of</strong> a two-dimensional wavepacket in a growing bo<strong>und</strong>ary<br />

layer ”, Proc. R. Soc. London, 3, 317–332 (1982).<br />

[2] Lin R.S. and Malik M.R , “On the stability <strong>of</strong> attachment-line bo<strong>und</strong>ary layers. ”, J. Fluid<br />

Mech., 3,239–255 (1995).<br />

[3] The<strong>of</strong>ilis V. Hein S and Dallmann U , “On the origins <strong>of</strong> unsteadiness and three dimensionality<br />

in a laminar separation bubble ”, Proc. R. Soc. London, 3,3229–3246 (2000).<br />

[4] U. Ehrenstein and F. Gallaire , “On two dimensional temporal modes in spatially evolving<br />

open flows: the flat-plate bo<strong>und</strong>ary layer ”, J. Fluid Mech., (2005).<br />

[5] Briley WR., “A numerical study <strong>of</strong> laminar separation bubbles using the Navier-Stokes<br />

equations.”, J. Fluid Mech., 3,713–736 (1971).<br />

[6] Timothy N. Phillips and Gareth W.Roberts, “The Treatment <strong>of</strong> Spurious Pressure Modes in<br />

Spectral Incompressible Flow Calculations.”, Journal <strong>of</strong> Computational Physics, 3,150–164<br />

(1993).<br />

[7] C. Canuto and M. Y. Hussaini and A. Quarteroni and T. A. Zang, “Spectral Methods in<br />

Fluids Dynamics”, Springer, (1988).<br />

Speaker: ALIZARD, F. 68 <strong>BAIL</strong> <strong>2006</strong><br />

2<br />

✩<br />

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!