Reduced-order Kalman-filtered hybrid simulation combining particle tracking velocimetry and direct numerical simulation

2012 ◽  
Vol 709 ◽  
pp. 249-288 ◽  
Author(s):  
Takao Suzuki

AbstractThe capability of state-of-the-art techniques integrating experimental and computational fluid dynamics has been expanding recently. In our previous study, we have developed a hybrid unsteady-flow simulation technique combining particle tracking velocimetry (PTV) and direct numerical simulation (DNS) and demonstrated its capability at low Reynolds numbers. Similar approaches have also been proposed by a few groups; however, applying algorithms of this type generally becomes more challenging with increasing Reynolds number because the time interval of the frame rate for particle image velocimetry (PIV) becomes much greater than the required computational time step, and the PIV/PTV resolution tends to be lower than that necessary for computational fluid dynamics. To extend the applicability to noisy time-resolved PIV/PTV data, the proposed algorithm optimizes the data input temporally and spatially by introducing a reduced-order Kalman filter. This study establishes a framework of the Kalman-filtered hybrid simulation and proves the concept by tackling a planar-jet flow at $\mathit{Re}\approx 2000$ as an example. We evaluate the filtering functions as well as convergence of the proposed algorithm by comparing with the existing PTV–DNS hybrid simulation, and show some techniques available to hybrid velocity fields by analysing vortical motion in the shear layers of the jet.

2010 ◽  
Vol 655 ◽  
pp. 344-379 ◽  
Author(s):  
TAKAO SUZUKI ◽  
HUI JI ◽  
FUJIO YAMAMOTO

Instability waves in a laminar planar jet are extracted using hybrid unsteady-flow simulation combining particle tracking velocimetry (PTV) and direct numerical simulation (DNS). Unsteady velocity fields on a laser sheet in a water tunnel are measured with time-resolved PTV; subsequently, PTV velocity fields are rectified in a least squares sense so that the equation of continuity is satisfied, and they are transplanted to a two-dimensional incompressible Navier–Stokes solver by setting a multiple of the computational time step equal to the frame rate of the PTV system. As a result, the unsteady hybrid velocity field approaches that of the measured one over time, and we can simultaneously acquire the unsteady pressure field. The resultant set of flow quantities satisfies the governing equations, and their resolution is comparable to that of numerical simulation with the noise level much lower than the original PTV data. From hybrid unsteady velocity fields, we extract eigenfunctions using bi-orthogonal decomposition as a spatial problem for viscous instability. We also investigate stability/convergence characteristics of the hybrid simulation referring to linear stability analysis.


ASAIO Journal ◽  
2020 ◽  
Vol 66 (2) ◽  
pp. 173-182 ◽  
Author(s):  
Bente Thamsen ◽  
Utku Gülan ◽  
Lena Wiegmann ◽  
Christian Loosli ◽  
Marianne Schmid Daners ◽  
...  

Author(s):  
Utku Gülan ◽  
Diego Gallo ◽  
Raffaele Ponzini ◽  
Beat Lüthi ◽  
Markus Holzner ◽  
...  

The complex hemodynamics observed in the human aorta make this district a site of election for an in depth investigation of the relationship between fluid structures, transport and pathophysiology. In recent years, the coupling of imaging techniques and computational fluid dynamics (CFD) has been applied to study aortic hemodynamics, because of the possibility to obtain highly resolved blood flow patterns in more and more realistic and fully personalized flow simulations [1]. However, the combination of imaging techniques and computational methods requires some assumptions that might influence the predicted hemodynamic scenario. Thus, computational modeling requires experimental cross-validation. Recently, 4D phase contrast MRI (PCMRI) has been applied in vivo and in vitro to access the velocity field in aorta [2] and to validate numerical results [3]. However, PCMRI usually requires long acquisition times and suffers from low spatial and temporal resolution and a low signal-to-noise ratio. Anemometric techniques have been also applied for in vitro characterization of the fluid dynamics in aortic phantoms. Among them, 3D Particle Tracking Velocimetry (PTV), an optical technique based on imaging of flow tracers successfully used to obtain Lagrangian velocity fields in a wide range of complex and turbulent flows [4], has been very recently applied to characterize fluid structures in the ascending aorta [5].


2017 ◽  
Vol 10 (8) ◽  
pp. 3145-3165 ◽  
Author(s):  
Chiel C. van Heerwaarden ◽  
Bart J. H. van Stratum ◽  
Thijs Heus ◽  
Jeremy A. Gibbs ◽  
Evgeni Fedorovich ◽  
...  

Abstract. This paper describes MicroHH 1.0, a new and open-source (www.microhh.org) computational fluid dynamics code for the simulation of turbulent flows in the atmosphere. It is primarily made for direct numerical simulation but also supports large-eddy simulation (LES). The paper covers the description of the governing equations, their numerical implementation, and the parameterizations included in the code. Furthermore, the paper presents the validation of the dynamical core in the form of convergence and conservation tests, and comparison of simulations of channel flows and slope flows against well-established test cases. The full numerical model, including the associated parameterizations for LES, has been tested for a set of cases under stable and unstable conditions, under the Boussinesq and anelastic approximations, and with dry and moist convection under stationary and time-varying boundary conditions. The paper presents performance tests showing good scaling from 256 to 32 768 processes. The graphical processing unit (GPU)-enabled version of the code can reach a speedup of more than an order of magnitude for simulations that fit in the memory of a single GPU.


2014 ◽  
Vol 47 (12) ◽  
pp. 3149-3155 ◽  
Author(s):  
Diego Gallo ◽  
Utku Gülan ◽  
Antonietta Di Stefano ◽  
Raffaele Ponzini ◽  
Beat Lüthi ◽  
...  

2017 ◽  
Author(s):  
Chiel C. van Heerwaarden ◽  
Bart J. H. van Stratum ◽  
Thijs Heus ◽  
Jeremy A. Gibbs ◽  
Evgeni Fedorovich ◽  
...  

Abstract. This paper describes MicroHH 1.0, a new and open source (www.microhh.org) computational fluid dynamics code for the simulation of turbulent flows in the atmosphere. It is primarily made for direct numerical simulation, but also supports large-eddy simulation (LES). The paper covers the description of the governing equations, their numerical implementation, and the parametrizations included in the code. Furthermore, the paper presents the validation of the dynamical core in the form of convergence and conservation tests, and comparison of simulations of channel flows and slope flows against well-established test cases. The full numerical model, including the associated parametrizations for LES, has been tested for a set of cases under stable and unstable conditions, under the Boussinesq and anelastic approximation, and with dry and moist convection under stationary and time-varying boundary conditions. The paper presents performance tests showing good scaling from 256 to 32,768 processes. The Graphical Processing Unit-enabled version of the code reaches speedups of more than an order of magnitude with respect to the conventional code for a variety of cases.


1996 ◽  
Vol 33 (9) ◽  
pp. 163-170 ◽  
Author(s):  
Virginia R. Stovin ◽  
Adrian J. Saul

Research was undertaken in order to identify possible methodologies for the prediction of sedimentation in storage chambers based on computational fluid dynamics (CFD). The Fluent CFD software was used to establish a numerical model of the flow field, on which further analysis was undertaken. Sedimentation was estimated from the simulated flow fields by two different methods. The first approach used the simulation to predict the bed shear stress distribution, with deposition being assumed for areas where the bed shear stress fell below a critical value (τcd). The value of τcd had previously been determined in the laboratory. Efficiency was then calculated as a function of the proportion of the chamber bed for which deposition had been predicted. The second method used the particle tracking facility in Fluent and efficiency was calculated from the proportion of particles that remained within the chamber. The results from the two techniques for efficiency are compared to data collected in a laboratory chamber. Three further simulations were then undertaken in order to investigate the influence of length to breadth ratio on chamber performance. The methodology presented here could be applied to complex geometries and full scale installations.


Sign in / Sign up

Export Citation Format

Share Document