2,034 research outputs found
A non-hybrid method for the PDF equations of turbulent flows on unstructured grids
In probability density function (PDF) methods of turbulent flows, the joint
PDF of several flow variables is computed by numerically integrating a system
of stochastic differential equations for Lagrangian particles. A set of
parallel algorithms is proposed to provide an efficient solution of the PDF
transport equation, modeling the joint PDF of turbulent velocity, frequency and
concentration of a passive scalar in geometrically complex configurations. An
unstructured Eulerian grid is employed to extract Eulerian statistics, to solve
for quantities represented at fixed locations of the domain (e.g. the mean
pressure) and to track particles. All three aspects regarding the grid make use
of the finite element method (FEM) employing the simplest linear FEM shape
functions. To model the small-scale mixing of the transported scalar, the
interaction by exchange with the conditional mean model is adopted. An adaptive
algorithm that computes the velocity-conditioned scalar mean is proposed that
homogenizes the statistical error over the sample space with no assumption on
the shape of the underlying velocity PDF. Compared to other hybrid
particle-in-cell approaches for the PDF equations, the current methodology is
consistent without the need for consistency conditions. The algorithm is tested
by computing the dispersion of passive scalars released from concentrated
sources in two different turbulent flows: the fully developed turbulent channel
flow and a street canyon (or cavity) flow. Algorithmic details on estimating
conditional and unconditional statistics, particle tracking and particle-number
control are presented in detail. Relevant aspects of performance and
parallelism on cache-based shared memory machines are discussed.Comment: Accepted in Journal of Computational Physics, Feb. 20, 200
Wavelet transforms and their applications to MHD and plasma turbulence: a review
Wavelet analysis and compression tools are reviewed and different
applications to study MHD and plasma turbulence are presented. We introduce the
continuous and the orthogonal wavelet transform and detail several statistical
diagnostics based on the wavelet coefficients. We then show how to extract
coherent structures out of fully developed turbulent flows using wavelet-based
denoising. Finally some multiscale numerical simulation schemes using wavelets
are described. Several examples for analyzing, compressing and computing one,
two and three dimensional turbulent MHD or plasma flows are presented.Comment: Journal of Plasma Physics, 201
Synchronization and optimization of Large Eddy Simulation using an online Ensemble Kalman Filter
An online Data Assimilation strategy based on the Ensemble Kalman Filter
(EnKF) is used to improve the predictive capabilities of Large Eddy Simulation
(LES) for the analysis of the turbulent flow in a plane channel, . The algorithm sequentially combines the LES prediction with
high-fidelity, sparse instantaneous data obtained from a Direct Numerical
Simulation (DNS). It is shown that the procedure provides an augmented state
which exhibits higher accuracy than the LES model and it synchronizes with the
time evolution of the high-fidelity DNS data if the hyperparameters governing
the EnKF are properly chosen. In addition, the data-driven algorithm is able to
improve the accuracy of the subgrid-scale model included in the LES, the
Smagorinsky model, via the optimization of a free coefficient. However, while
the online EnKF strategy is able to reduce the global error of the LES
prediction, a discrepancy with the reference DNS data is still observed because
of structural flaws of the subgrid-scale model used
Effect of turbulent fluctuations on the drag and lift forces on a towed sphere and its boundary layer
The impact of turbulent fluctuations on the forces exerted by a fluid on a
towed spherical particle is investigated by means of high-resolution direct
numerical simulations. The measurements are carried out using a novel scheme to
integrate the two-way coupling between the particle and the incompressible
surrounding fluid flow maintained in a high-Reynolds-number turbulent regime.
The main idea consists in combining a Fourier pseudo-spectral method for the
fluid with an immersed-boundary technique to impose the no-slip boundary
condition on the surface of the particle. Benchmarking of the code shows a good
agreement with experimental and numerical measurements from other groups. A
study of the turbulent wake downstream the sphere is also reported. The mean
velocity deficit is shown to behave as the inverse of the distance from the
particle, as predicted from classical similarity analysis. This law is
reinterpreted in terms of the principle of "permanence of large eddies" that
relates infrared asymptotic self-similarity to the law of decay of energy in
homogeneous turbulence.
The developed method is then used to attack the problem of an upstream flow
that is in a developed turbulent regime. It is shown that the average drag
force increases as a function of the turbulent intensity and the particle
Reynolds number. This increase is significantly larger than predicted by
standard drag correlations based on laminar upstream flows. It is found that
the relevant parameter is the ratio of the viscous boundary layer thickness to
the dissipation scale of the ambient turbulent flow. The drag enhancement can
be motivated by the modification of the mean velocity and pressure profile
around the sphere by small scale turbulent fluctuations.Comment: 24 pages, 22 figure
- …