43 research outputs found

    Three-dimensional simulations of Bingham plastic flows with the multiple-relaxation-time lattice Boltzmann model

    Get PDF
    This paper presents a three-dimensional (3D) parallel multiple-relaxation-time lattice Boltzmann model (MRT-LBM) for Bingham plastics which overcomes numerical instabilities in the simulation of non-Newtonian fluids for the Bhatnagar–Gross–Krook (BGK) model. The MRT-LBM and several related mathematical models are briefly described. Papanastasiou’s modified model is incorporated for better numerical stability. The impact of the relaxation parameters of the model is studied in detail. The MRT-LBM is then validated through a benchmark problem: a 3D steady Poiseuille flow. The results from the numerical simulations are consistent with those derived analytically which indicates that the MRT-LBM effectively simulates Bingham fluids but with better stability. A parallel MRT-LBM framework is introduced, and the parallel efficiency is tested through a simple case. The MRT-LBM is shown to be appropriate for parallel implementation and to have high efficiency. Finally, a Bingham fluid flowing past a square-based prism with a fixed sphere is simulated. It is found the drag coefficient is a function of both Reynolds number (Re) and Bingham number (Bn). These results reveal the flow behavior of Bingham plastics

    Anomalous kinetics and transport from 1D self--consistent mode--coupling theory

    Full text link
    We study the dynamics of long-wavelength fluctuations in one-dimensional (1D) many-particle systems as described by self-consistent mode-coupling theory. The corresponding nonlinear integro-differential equations for the relevant correlators are solved analytically and checked numerically. In particular, we find that the memory functions exhibit a power-law decay accompanied by relatively fast oscillations. Furthermore, the scaling behaviour and, correspondingly, the universality class depends on the order of the leading nonlinear term. In the cubic case, both viscosity and thermal conductivity diverge in the thermodynamic limit. In the quartic case, a faster decay of the memory functions leads to a finite viscosity, while thermal conductivity exhibits an even faster divergence. Finally, our analysis puts on a more firm basis the previously conjectured connection between anomalous heat conductivity and anomalous diffusion

    Lattice Boltzmann simulations of two-phase flow with high density ratio in axially symmetric geometry

    Get PDF
    In this paper, a two-phase lattice Boltzmann (LB) model, developed for simulating fluid flows on a Cartesian grid at high liquid-to-gas density ratios, is adapted to an axisymmetric coordinate system. This is achieved by incorporating additional source terms in the planar evolution equations for the density and pressure distribution functions such that the axisymmetric mass and momentum conservation equations are recovered in the macroscopic limit. Appropriate numerical treatment of the terms is performed to obtain stable computations at high density ratio for this axisymmetric model. The particle collision is modeled by employing multiple relaxation times to attain stability at low viscosity. The model is evaluated by verifying the Laplace-Young relation for a liquid drop, comparing computed frequency of oscillations of an initially ellipsoidal drop with analytical values and comparing the behavior of a spherical drop impinging on a wet wall with prior results. The time evolution of the radial distance of the tip of the corona, formed when the drop impinges, agrees well with prior data.Shiladitya Mukherjee and John Abraha

    Simulation of static critical phenomena in non-ideal fluids with the Lattice Boltzmann method

    Full text link
    A fluctuating non-ideal fluid at its critical point is simulated with the Lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order-parameter.Comment: 23 pages, 16 figure

    Energetic α-particle sources produced through proton-boron reactions by high-energy high-intensity laser beam

    No full text
    In an experiment performed with a high-intensity and high-energy laser system, α-particle production in proton-boron reaction by using a laser-driven proton beam was measured. α particles were observed from the front and also from the rear side, even after a 2-mm-thick boron target. The data obtained in this experiment have been analyzed using a sequence of numerical simulations. The simulations clarify the mechanisms of α-particle production and transport through the boron targets. α-particle energies observed in the experiment and in the simulation reach 10–20 MeV through energy transfer from 20–30 MeV energy incident protons. Despite the lower cross sections for protons with energy above the sub-MeV resonances in the proton-boron reactions, 10^8–10^9 α particles per steradian have been detected
    corecore