Applying Ensemble Kalman Filter to Transonic Flows Through a Two-Dimensional Turbine Cascade

2021 ◽  
Vol 143 (12) ◽  
Author(s):  
Sasuga Ito ◽  
Masato Furukawa ◽  
Kazutoyo Yamada ◽  
Kaito Manabe

Abstract Turbulence is one of the most important phenomena in fluid dynamics. Large eddy simulation (LES) generally allows us to analyze smaller eddies than when using simulations based on unsteady Reynolds-averaged Navier–Stokes equations (URANS). In addition, the numerical solutions of LES show good agreements with experiments and numerical solutions based on direct numerical simulation. URANS simulations are, however, frequently used in academia and industry because LES computations are much more expensive compared with URANS simulations. In this investigation, an optimization of unsolved coefficients of the k–ω two equations model is performed on the transonic flow around T106A low-pressure turbine cascade to improve the accuracy of turbulence prediction with URANS. For the optimization approach, two-dimensional URANS is combined with ensemble Kalman filter which is one of the data assimilation techniques. In the assimilation process, a time- and spanwise-averaged LES result is used as pseudo-experimental data. Three-dimensional URANS simulations are performed for the evaluation of the optimization effect. URANS simulations are also applied to a different turbine cascade flow for the evaluation of the robustness of the optimized coefficients. These URANS results confirmed that the optimized coefficients improve the accuracy of turbulence prediction.

1995 ◽  
Vol 291 ◽  
pp. 369-392 ◽  
Author(s):  
Ronald D. Joslin

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier–Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic-source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in flat-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.


2018 ◽  
Author(s):  
Jiajun Chen ◽  
Yue Sun ◽  
Hang Zhang ◽  
Dakui Feng ◽  
Zhiguo Zhang

Mixing in pipe junctions can play an important role in exciting force and distribution of flow in pipe network. This paper investigated the cross pipe junction and proposed an improved plan, Y-shaped pipe junction. The numerical study of a three-dimensional pipe junction was performed for calculation and improved understanding of flow feature in pipe. The filtered Navier–Stokes equations were used to perform the large-eddy simulation of the unsteady incompressible flow in pipe. From the analysis of these results, it clearly appears that the vortex strength and velocity non-uniformity of centerline, can be reduced by Y-shaped junction. The Y-shaped junction not only has better flow characteristic, but also reduces head loss and exciting force. The results of the three-dimensional improvement analysis of junction can be used in the design of pipe network for industry.


Inventions ◽  
2020 ◽  
Vol 5 (3) ◽  
pp. 33
Author(s):  
Horng-Wen Wu ◽  
Tang-Hong Chen ◽  
Nugroho-Putra Kelana ◽  
De-An Huang

This study analyzes transient turbulent modeling of three-dimensional multiple dimpled fin array using large eddy simulation (LES). The Navier–Stokes equations as well as the energy equation were constructed by the finite volume method and then discretized to form algebraic equations, which were solved by semi-implicit method for pressure-linked equation (SIMPLE). The solutions of temperature and velocity were obtained by iterating computation until it converged within each step. This simulation places nine fins on the bottom surface of a channel and changes the height of the dimple (0.4, 0.8, and 1.2 mm) with three different levels of Reynolds number (Re) (3500, 5000, and 6500) to investigate the temperature and flow field without gravity in forced convection. The results indicate that the dimpled fin array can generate vortices between the convex/concave dimples and the fin base and increase the influences of the height of the dimple on the flow field around the fin array. The averaged time-mean of the Nusselt number (Nu) for the dimple height of 0.8 mm is higher than that of the no-dimple case up to 14.4%, while the averaged time-mean Nu for the dimple height of 1.2 mm is lower than that of the no-dimple case up to 11.6%.


2006 ◽  
Vol 128 (6) ◽  
pp. 1394-1399 ◽  
Author(s):  
Donghyun You ◽  
Meng Wang ◽  
Rajat Mittal ◽  
Parviz Moin

A novel structured grid approach which provides an efficient way of treating a class of complex geometries is proposed. The incompressible Navier-Stokes equations are formulated in a two-dimensional, generalized curvilinear coordinate system complemented by a third quasi-curvilinear coordinate. By keeping all two-dimensional planes defined by constant third coordinate values parallel to one another, the proposed approach significantly reduces the memory requirement in fully three-dimensional geometries, and makes the computation more cost effective. The formulation can be easily adapted to an existing flow solver based on a two-dimensional generalized coordinate system coupled with a Cartesian third direction, with only a small increase in computational cost. The feasibility and efficiency of the present method have been assessed in a simulation of flow over a tapered cylinder.


1990 ◽  
Vol 220 ◽  
pp. 397-411 ◽  
Author(s):  
Laura L. Pauley ◽  
Parviz Moin ◽  
William C. Reynolds

The separation of a two-dimensional laminar boundary layer under the influence of a suddenly imposed external adverse pressure gradient was studied by time-accurate numerical solutions of the Navier–Stokes equations. It was found that a strong adverse pressure gradient created periodic vortex shedding from the separation. The general features of the time-averaged results were similar to experimental results for laminar separation bubbles. Comparisons were made with the ‘steady’ separation experiments of Gaster (1966). It was found that his ‘bursting’ occurs under the same conditions as our periodic shedding, suggesting that bursting is actually periodic shedding which has been time-averaged. The Strouhal number based on the shedding frequency, local free-stream velocity, and boundary-layer momentum thickness at separation was independent of the Reynolds number and the pressure gradient. A criterion for onset of shedding was established. The shedding frequency was the same as that predicted for the most amplified linear inviscid instability of the separated shear layer.


2011 ◽  
Vol 666 ◽  
pp. 506-520 ◽  
Author(s):  
F. DOMENICHINI

The vortex formation behind an orifice is a widely investigated phenomenon, which has been recently studied in several problems of biological relevance. In the case of a circular opening, several works in the literature have shown the existence of a limiting process for vortex ring formation that leads to the concept of critical formation time. In the different geometric arrangement of a planar flow, which corresponds to an opening with straight edges, it has been recently outlined that such a concept does not apply. This discrepancy opens the question about the presence of limiting conditions when apertures with irregular shape are considered. In this paper, the three-dimensional vortex formation due to the impulsively started flow through slender openings is studied with the numerical solution of the Navier–Stokes equations, at values of the Reynolds number that allow the comparison with previous two-dimensional findings. The analysis of the three-dimensional results reveals the two-dimensional nature of the early vortex formation phase. During an intermediate phase, the flow evolution appears to be driven by the local curvature of the orifice edge, and the time scale of the phenomena exhibits a surprisingly good agreement with those found in axisymmetric problems with the same curvature. The long-time evolution shows the complete development of the three-dimensional vorticity dynamics, which does not allow the definition of further unifying concepts.


1987 ◽  
Vol 184 ◽  
pp. 207-243 ◽  
Author(s):  
Ralph W. Metcalfe ◽  
Steven A. Orszag ◽  
Marc E. Brachet ◽  
Suresh Menon ◽  
James J. Riley

The three-dimensional stability of two-dimensional vortical states of planar mixing layers is studied by direct numerical integration of the Navier-Stokes equations. Small-scale instabilities are shown to exist for spanwise scales at which classical linear modes are stable. These modes grow on convective timescales, extract their energy from the mean flow and exist at moderately low Reynolds numbers. Their growth rates are comparable with the most rapidly growing inviscid instability and with the growth rates of two-dimensional subharmonic (pairing) modes. At high amplitudes, they can evolve into pairs of counter-rotating, streamwise vortices, connecting the primary spanwise vortices, which are very similar to the structures observed in laboratory experiments. The three-dimensional modes do not appear to saturate in quasi-steady states as do the purely two-dimensional fundamental and subharmonic modes in the absence of pairing. The subsequent evolution of the flow depends on the relative amplitudes of the pairing modes. Persistent pairings can inhibit three-dimensional instability and, hence, keep the flow predominantly two-dimensional. Conversely, suppression of the pairing process can drive the three-dimensional modes to more chaotic, turbulent-like states. An analysis of high-resolution simulations of fully turbulent mixing layers confirms the existence of rib-like structures and that their coherence depends strongly on the presence of the two-dimensional pairing modes.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Donald P. Rizzetta ◽  
Miguel R. Visbal

Numerical calculations were carried out to investigate control of transition on a flat plate by means of local dynamic surface deformation. The configuration and flow conditions are similar to a previous computation which simulated transition mitigation. Physically, the surface modification may be produced by piezoelectrically driven actuators located below a compliant aerodynamic surface, which have been employed experimentally. One actuator is located in the upstream plate region and oscillated at the most unstable frequency of 250 Hz to develop disturbances representing Tollmien–Schlichting instabilities. A controlling actuator is placed downstream and oscillated at the same frequency, but with an appropriate phase shift and modified amplitude to decrease disturbance growth and delay transition. While the downstream controlling actuator is two-dimensional (spanwise invariant), several forms of upstream disturbances were considered. These included disturbances which were strictly two-dimensional, those which were modulated in amplitude and those which had a spanwise variation of the temporal phase shift. Direct numerical simulations were obtained by solution of the three-dimensional compressible Navier–Stokes equations, utilizing a high-fidelity computational scheme and an implicit time-marching approach. A previously devised empirical process was applied for determining the optimal parameters of the controlling actuator. Results of the simulations are described, features of the flowfields elucidated, and comparisons made between solutions of the uncontrolled and controlled cases for the respective incoming disturbances. It is found that the disturbance growth is mitigated and the transition is delayed for all forms of the upstream perturbations, substantially reducing the skin friction.


Author(s):  
Sasuga Ito ◽  
Masato Furukawa ◽  
Yamada Kazutoyo ◽  
Kaito Manabe

Abstract Turbulence is one of the most important phenomena in fluid dynamics. In general, turbulent phenomena can be resolved more clearly with Large Eddy Simulation (LES) compared with Unsteady Reynolds Averaged Navier-Stokes (URANS), and the numerical solution shows good agreements with that based on Direct Numerical Simulation (DNS). However, more time and computational power are needed on LES than those on URANS. If possible, the ideal simulation method is that the method is able to resolve the turbulent phenomena same quality as the results based on DNS and LES with less time and less computational power same as that on URANS. This paper shows an adaptive simulation method based on URANS and Ensemble Kalman Filter (Enkf) to reproduce the flow fields based on LES. In this study, a two-dimensional turbine cascade flow has been solved with URANS and LES. The adaptive simulation method has been also applied to the cascade flow. As the results, in the flow field of URANS with the assimilated turbulence model’s parameters, the separation phenomenon and the boundary layer thickness was close to that of the time averaged LES.


Author(s):  
Mustafa Koz ◽  
Serhat Yesilyurt

Microorganisms such as bacteria use their rotating helical flagella for propulsion speeds up to tens of tail lengths per second. The mechanism can be utilized for controlled pumping of liquids in microchannels. In this study, we aim to analyze the effects of control parameters such as axial span between helical rounds (wavelength), angular velocity of rotations (frequency), and the radius of the helix (amplitude) on the maximum time-averaged flow rate, maximum head, rate of energy transfer, and efficiency of the micropump. The analysis is based on simulations obtained from the three-dimensional time-dependent numerical model of the flow induced by the rotating spiral inside a rectangular-prism channel. The flow is governed by Navier-Stokes equations subject to continuity in time-varying domain due to moving boundaries of the spiral. Numerical solutions are obtained using a commercial finite-element package which uses arbitrary Lagrangian-Eulerian method for mesh deformations. Results are compared with asymptotic results obtained from the resistive-force-theory available in the literature.


Sign in / Sign up

Export Citation Format

Share Document