8,736 research outputs found
Lattice Boltzmann study on Kelvin-Helmholtz instability: the roles of velocity and density gradients
A two-dimensional lattice Boltzmann model with 19 discrete velocities for
compressible Euler equations is proposed (D2V19-LBM). The fifth-order Weighted
Essentially Non-Oscillatory (5th-WENO) finite difference scheme is employed to
calculate the convection term of the lattice Boltzmann equation. The validity
of the model is verified by comparing simulation results of the Sod shock tube
with its corresponding analytical solutions. The velocity and density gradient
effects on the Kelvin-Helmholtz instability (KHI) are investigated using the
proposed model. Sharp density contours are obtained in our simulations. It is
found that, the linear growth rate for the KHI decreases with
increasing the width of velocity transition layer but increases with
increasing the width of density transition layer . After the
initial transient period and before the vortex has been well formed, the linear
growth rates, and , vary with and
approximately in the following way, and
, where , ,
and are fitting parameters and is the effective
interaction width of density transition layer. When
the linear growth rate does not vary significantly any more.
One can use the hybrid effects of velocity and density transition layers to
stabilize the KHI. Our numerical simulation results are in general agreement
with the analytical results [L. F. Wang, \emph{et al.}, Phys. Plasma
\textbf{17}, 042103 (2010)].Comment: Accepted for publication in PR
A hybrid method for hydrodynamic-kinetic flow - Part I -A particle-gridmethod for reducing stochastic noise in kinetic regimes
In this work we present a hybrid particle-grid Monte Carlo method for the Boltzmann equation, which is characterized by a significant reduction of the stochastic noise in the kinetic regime. The hybrid method is based on a first order splitting in time to separate the transport from the relaxation step. The transport step is solved by a deterministic scheme, while a hybrid DSMC-based method is used to solve the collision step. Such a hybrid scheme is based on splitting the solution in a collisional and a non-collisional part at the beginning of the collision step, and the DSMC method is used to solve the relaxation step for the collisional part of the solution only. This is accomplished by sampling only the fraction of particles candidate for collisions from the collisional part of the solution, performing collisions as in a standard DSMC method, and then projecting the particles back onto a velocity grid to compute a piecewise constant reconstruction for the collisional part of the solution. The latter is added to a piecewise constant reconstruction of the non-collisional part of the solution, which in fact remains unchanged during the relaxation step. Numerical results show that the stochastic noise is significantly reduced at large Knudsen numbers with respect to the standard DSMC method. Indeed in this algorithm, the particle scheme is applied only on the collisional part of the solution, so only this fraction of the solution is affected by stochastic fluctuations. But since the collisional part of the solution reduces as the Knudsen number increases, stochastic noise reduces as well at large Knudsen number
Modeling asteroid collisions and impact processes
As a complement to experimental and theoretical approaches, numerical
modeling has become an important component to study asteroid collisions and
impact processes. In the last decade, there have been significant advances in
both computational resources and numerical methods. We discuss the present
state-of-the-art numerical methods and material models used in "shock physics
codes" to simulate impacts and collisions and give some examples of those
codes. Finally, recent modeling studies are presented, focussing on the effects
of various material properties and target structures on the outcome of a
collision.Comment: Chapter to appear in the Space Science Series Book: Asteroids IV.
Includes minor correction
A Particle-based Multiscale Solver for Compressible Liquid-Vapor Flow
To describe complex flow systems accurately, it is in many cases important to
account for the properties of fluid flows on a microscopic scale. In this work,
we focus on the description of liquid-vapor flow with a sharp interface between
the phases. The local phase dynamics at the interface can be interpreted as a
Riemann problem for which we develop a multiscale solver in the spirit of the
heterogeneous multiscale method, using a particle-based microscale model to
augment the macroscopic two-phase flow system. The application of a microscale
model makes it possible to use the intrinsic properties of the fluid at the
microscale, instead of formulating (ad-hoc) constitutive relations
Comparison of macro- and microscopic solutions of the Riemann problem I. Supercritical shock tube and expansion into vacuum
The Riemann problem is a fundamental concept in the development of numerical methods for the macroscopic flow equations. It allows the resolution of discontinuities in the solution, such as shock waves, and provides a powerful tool for the construction of numerical flux functions. A natural extension of the Riemann problem involves two phases, a liquid and a vapour phase which undergo phase change at the material boundary. For this problem, we aim at a comparison with the macroscopic solution from molecular dynamics simulations. In this work, as a first step, the macroscopic solution of two important Riemann problem scenarios, the supercritical shock tube and the expansion into vacuum, were compared to microscopic solutions produced by molecular dynamics simulations. High fidelity equations of state were used to accurately approximate the material behaviour of the model fluid. The results of both scenarios compare almost perfect with each other. During the vacuum expansion, the fluid obtained a state of non-equilibrium, where the microscopic and macroscopic solutions start to diverge. A limiting case was shown, where liquid droplets appeared in the expansion fan, which was approximated by the macroscopic solution, assuming an undercooled vapour.DFG, 84292822, TRR 75: Tropfendynamische Prozesse unter extremen Umgebungsbedingunge
Recommended from our members
Advances and Challenges in Computational Research of Micro and Nano Flows
This paper was presented at the 4th Micro and Nano Flows Conference (MNF2014), which was held at University College, London, UK. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute, ASME Press, LCN London Centre for Nanotechnology, UCL University College London, UCL Engineering, the International NanoScience Community, www.nanopaprika.eu.This paper presents a collective overview of recent studies regarding the computational modelling
of micro- and nano-fluidic systems. The review provides an introduction to atomistic, mesoscale and hybrid
methods for simulating micro and nano-flows, as well as discusses recent applications and results from the
application of such methods
- …