Computational experiments are one of the most used and flexible investigation tools in fluid dynamics. The Lattice Boltzmann Equation is a well established computational method particularly promising for multi-phase flows at micro and macro scales. Here we present preliminary results on performances of the LBE method on the Cell Broadband Engine platform.

We present a study of Eulerian and Lagrangian statistics from a high-resolution numerical simulation of isotropic and homogeneous turbulence using the FLASH code, with an estimated Taylor microscale Reynolds number of around 600. Statistics are evaluated over a data set with 18563 spatial grid points and with 2563 = 16.8 million particles, followed for about one large-scale eddy turnover time. We present data for the Eulerian and Lagrangian structure functions up to the tenth order. We analyze the local scaling properties in the inertial range.

We present the results of a high resolution numerical study of two-dimensional (2D) Rayleigh-Taylor turbulence using a recently proposed thermal lattice Boltzmann method The goal of our study is both methodological and physical We assess merits and limitations concerning small- and large-scale resolution/accuracy of the adopted integration scheme We discuss quantitatively the requirements needed to keep the method stable and precise enough to simulate stratified and unstratified flows driven by thermal active fluctuations at high Rayleigh and high Reynolds numbers We present data with spatial

This paper describes a collocation method for solving numerically a singular integral equation with Cauchy and Volterra operators, associated with a proper constraint condition. The numerical method is based on the transformation of the given integral problem into a hypersingular integral equation and then applying a collocation method to solve the latter equation. Convergence of the resulting method is then discussed, and optimal convergence rates for the collocation and discrete collocation methods are given in suitable weighted Sobolev spaces.

Human blood flow is a multiscale problem: in first approximation, blood is a dense suspension of plasma and deformable red cells. Physiological vessel diameters range from about one to thousands of cell radii. Current computational models either involve a homogeneous fluid and cannot track particulate effects or describe a relatively small number of cells with high resolution but are incapable to reach relevant time and length scales. Our approach is to simplify much further than existing particulate models.

We present the results of our numerical simulations of the Rayleigh-Taylor turbulence, performed using a recently proposed (Sbragaglia et al 2009 J. Fluid Mech. 628 299, Scagliarini et al 2010 Phys. Fluids 22 055101) lattice Boltzmann method that can describe consistently a thermal compressible flow subjected to an external forcing. The method allowed us to study the system in both the nearly Boussinesq regime and the strongly compressible regime. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.

We are concerned with the discretization of optimal control problems when a Runge-Kutta scheme is selected for the related Hamiltonian system. It is known that Lagrangian's first order conditions on the discrete model, require a symplectic partitioned Runge-Kutta scheme for state-costate equations. In the present paper this result is extended to growth models, widely used in Economics studies, where the system is described by a current Hamiltonian. We prove that a correct numerical treatment of the state-current costate system needs Lawson exponential schemes for the costate approximation.

Phase separation in a complex fluid with lamellar order has been studied in the case of cold thermal fronts propagating diffusively from external walls. The velocity hydrodynamic modes are taken into account by coupling the convection-diffusion equation for the order parameter to a generalized Navier-Stokes equation. The dynamical equations are simulated by implementing a hybrid method based on a lattice Boltzmann algorithm coupled to finite difference schemes.