High resolution numerical study of Rayleigh-Taylor turbulence using a thermal lattice Boltzmann scheme

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

Numerical solution of a singular integral equation with Cauchy kernel in the plane contact problem

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.

Simplified particulate model for coarse-grained hemodynamics simulations

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.

Numerical simulations of compressible Rayleigh-Taylor turbulence in stratified fluids

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.

Exponential Lawson integration for nearly Hamiltonian systems arising in optimal control

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.

Thermal and hydrodynamic effects in the ordering of lamellar fluids

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.

Polynomial functors and opetopes

We give an elementary and direct combinatorial definition of opetopes in terms of trees, well-suited for graphical manipulation and explicit computation. To relate our definition to the classical definition, we recast the Baez-Dolan slice construction for operads in terms of polynomial monads: our opetopes appear naturally as types for polynomial monads obtained by iterating the Baez-Dolan construction, starting with the trivial monad. We show that our notion of opetope agrees with Leinster's. Next we observe a suspension operation for opetopes, and define a notion of stable opetopes.