5,579 research outputs found
Fluid Simulations with Localized Boltzmann Upscaling by Direct Simulation Monte-Carlo
In the present work, we present a novel numerical algorithm to couple the
Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann
equation with a finite volume like method for the solution of the Euler
equations. Recently we presented in [14],[16],[17] different methodologies
which permit to solve fluid dynamics problems with localized regions of
departure from thermodynamical equilibrium. The methods rely on the
introduction of buffer zones which realize a smooth transition between the
kinetic and the fluid regions. In this paper we extend the idea of buffer zones
and dynamic coupling to the case of the Monte Carlo methods. To facilitate the
coupling and avoid the onset of spurious oscillations in the fluid regions
which are consequences of the coupling with a stochastic numerical scheme, we
use a new technique which permits to reduce the variance of the particle
methods [11]. In addition, the use of this method permits to obtain estimations
of the breakdowns of the fluid models less affected by fluctuations and
consequently to reduce the kinetic regions and optimize the coupling. In the
last part of the paper several numerical examples are presented to validate the
method and measure its computational performances
A thermal lattice Boltzmann model for micro/nano-flows
The dynamic behavior of charged micro and nanofluids plays a crucial role in a large variety of industrial and biological processes. Such dynamic behavior is characterized by the simultaneous occurrence of several competing mechanisms, such as electrostatic interactions, viscous dissipation and hydrodynamic effects, often taking place in complex geometries. This paper focuses on a thermal lattice Boltzmann model for micro/nano-flows
Mesoscopic simulation of diffusive contaminant spreading in gas flows at low pressure
Many modern production and measurement facilities incorporate multiphase
systems at low pressures. In this region of flows at small, non-zero Knudsen-
and low Mach numbers the classical mesoscopic Monte Carlo methods become
increasingly numerically costly. To increase the numerical efficiency of
simulations hybrid models are promising. In this contribution, we propose a
novel efficient simulation approach for the simulation of two phase flows with
a large concentration imbalance in a low pressure environment in the low
intermediate Knudsen regime. Our hybrid model comprises a lattice-Boltzmann
method corrected for the lower intermediate Kn regime proposed by Zhang et al.
for the simulation of an ambient flow field. A coupled event-driven
Monte-Carlo-style Boltzmann solver is employed to describe particles of a
second species of low concentration. In order to evaluate the model, standard
diffusivity and diffusion advection systems are considered.Comment: 9 pages, 8 figure
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
Relaxation time simulation method with internal-energy exchange for perfect gas at near continuum
This paper presents an internal energy exchange scheme for the relaxation time simulation method (RTSM) which solves the BGK equation for the perfect gas flow at near-continuum region discrete rotational energies are introduced to model the relaxation of internal energy modes. This development improved the agreements between RTSM and DSMC with little additional computational cost. The result shows a possibility of an improved hybrid RTSM/DSMC code for the continuum/rarefied gas flow
Employing per-component time step in DSMC simulations of disparate mass and cross-section gas mixtures
A new approach to simulation of stationary flows by Direct Simulation Monte
Carlo method is proposed. The idea is to specify an individual time step for
each component of a gas mixture. The approach consists of modifications mainly
to collision phase and recommendation on choosing time step ratios. It allows
softening the demands on the computational resources for cases of disparate
collision diameters of molecules and/or disparate molecular masses. These are
the cases important in vacuum deposition technologies. Few tests of the new
approach are made. Finally, the usage of new approach is demonstrated on a
problem of silver nanocluster diffusion in carrier gas argon in conditions of
silver deposition experiments.Comment: The goal of submission is to find native English speaker willing to
help me polish the paper. This is paper draft sent to Communications in
Computational Physics. It is recommended to publication. The need of
polishing was one of editors decision. See Additional data for MS Word source
fil
Unified Gas-kinetic Wave-Particle Methods III: Multiscale Photon Transport
In this paper, we extend the unified gas-kinetic wave-particle (UGKWP) method
to the multiscale photon transport. In this method, the photon free streaming
and scattering processes are treated in an un-splitting way. The duality
descriptions, namely the simulation particle and distribution function, are
utilized to describe the photon. By accurately recovering the governing
equations of the unified gas-kinetic scheme (UGKS), the UGKWP preserves the
multiscale dynamics of photon transport from optically thin to optically thick
regime. In the optically thin regime, the UGKWP becomes a Monte Carlo type
particle tracking method, while in the optically thick regime, the UGKWP
becomes a diffusion equation solver. The local photon dynamics of the UGKWP, as
well as the proportion of wave-described and particle-described photons are
automatically adapted according to the numerical resolution and transport
regime. Compared to the -type UGKS, the UGKWP requires less memory cost
and does not suffer ray effect. Compared to the implicit Monte Carlo (IMC)
method, the statistical noise of UGKWP is greatly reduced and computational
efficiency is significantly improved in the optically thick regime. Several
numerical examples covering all transport regimes from the optically thin to
optically thick are computed to validate the accuracy and efficiency of the
UGKWP method. In comparison to the -type UGKS and IMC method, the UGKWP
method may have several-order-of-magnitude reduction in computational cost and
memory requirement in solving some multsicale transport problems.Comment: 27 pages, 15 figures. arXiv admin note: text overlap with
arXiv:1810.0598
- …