Refine
Year of publication
Document Type
- Preprint (29)
Language
- English (29)
Has Fulltext
- yes (29)
Keywords
Faculty / Organisational entity
In the paper we discuss the transition from kinetic theory to macroscopic fluid equations, where the macroscopic equations are defined as aymptotic limits of a kinetic equation. This relation can be used to derive computationally efficient domain decomposition schemes for the simulaion of rarefied gas flows close to the continuum limit. Moreover, we present some basic ideas for the derivation of kinetic induced numerical schemes for macroscopic equations, namely kinetic schemes for general conservation laws as well as Lattice-Boltzmann methods for the incompressible Navier-Stokes equations.
In the present paper we investigate the Rayleigh-Benard convection in rarefied gases and demonstrate by numerical experiments the transition from purely thermal conduction to a natural convective flow for a large range of Knudsen numbers from 0.02 downto 0.001. We address to the problem how the critical value for the Rayleigh number defined for incompressible vsicous flows may be translated to rarefied gas flows. Moreover, the simulations obtained for a Knudsen number Kn=0.001 and Froude number Fr=1 show a further transition from regular Rayleigh-Benard cells to a pure unsteady behavious with moving vortices.
The paper presents some new estimates on the gain term of the Boltzmann collision operator. For Maxwellian molecules, it is shown that the L -norm of the gain term can be bounded in terms of the L1 and L -norm of the density function f. In the case of more general collision kernels, like the hard-sphere interaction potential, the gain term is estimated pointwise by the L -norm of the density function and the loss term of the Boltzmann collision operator.
The paper presents numerical results on the simulation of boundary value problems for the Boltzmann equation in one and two dimensions. In the one-dimensional case, we use prescribed fluxes at the left and diffusive conditions on the right end of a slab to study the resulting steady state solution. Moreover, we compute the numerical density function in velocity space and compare the result with the Chapman-Enskog distribution obtained in the limit for continuous media. The aim of the two-dimensional simulations is to investigate the possibility of a symmetry break in the numerical solution.
This report contains the following three papers about computations of rarefied gas flows:; ; a) Rarefied gas flow around a disc with different angles of attack, published in the proceedings of the 17th RGD Symposium, Aachen, 1990.; ; b) Hypersonic flow calculations around a 3D-deltawing at low Knudsen numbers, published in the proceedings of the 17th RGD Symposium,; Aachen, 1990.; ; c) Rarefied gas flow around a 3D-deltawing, published in the proceedings of the Workshop on Hypersonic Flows for Reentry Problems,; Part 1, Antibes, France, January 22-25, 1990.; ; All computations are part of the HERMES Research and Development Program.
Second Order Scheme for the Spatially Homogeneous Boltzmann Equation with Maxwellian Molecules
(1995)
In the standard approach, particle methods for the Boltzmann equation are obtained using an explicit time discretization of the spatially homogeneous Boltzmann equation. This kind of discretization leads to a restriction of the discretization parameter as well as on the differential cross section in the case of the general Boltzmann equation. Recently, it was shown, how to construct an implicit particle scheme for the Boltzmann equation with Maxwellian molecules. The present paper combines both approaches using a linear combination of explicit and implicit discretizations. It is shown that the new method leads to a second order particle method, when using an equiweighting of explicit and implicit discretization.
The paper presents a parallelization technique for the finite pointset method, a numerical method for rarefied gas flows.; First we give a short introduction to the Boltzmann equation, which describes the behaviour of rarefied gas flows. The basic ideas of the finite pointset method are presented and a strategy to parallelize the algorithm will be explained. It is shown that a static processor partition leads to an insufficient load-balance of the processors. Therefore an optimized parallelization technique based on an adaptive processor partition will be introduced, which improves the efficiency of the simulation code over the whole region of interesting flow situation. Finally we present a comparison of the CPU-times between a parallel computer and a vector computer.
The asymptotic analysis of IBVPs for the singularly perturbed parabolic PDE ... in the limit epsilon to zero motivate investigations of certain recursively defined approximative series ("ping-pong expansions"). The recursion formulae rely on operators assigning to a boundary condition at the left or the right boundary a solution of the parabolic PDE. Sufficient conditions for uniform convergence of ping-pong expansions are derived and a detailed analysis for the model problem ... is given.
The system of shallow water waves is one of the classical examples for nonlinear, twodimensional conservation laws. The paper investigates a simple kinetic equation depending on a parameter e which leads for e to 0 to the system of shallow water waves. The corresponding equilibrium distribution function has a compact support which depends on the eigenvalues of the hyperbolic system. It is shown that this kind of kinetic approach is restricted to a special class of nonlinear conservation laws. The kinetic model is used to develop a simple particle method for the numerical solution of shallow water waves. The particle method can be implemented in a straightforward way and produces in test examples sufficiently accurate results.
Numerical Simulation of the Stationary One-Dimensional Boltzmann Equation by Particle Methods
(1995)
The paper presents a numerical simulation technique - based on the well-known particle methods - for the stationary, one-dimensional Boltzmann equation for Maxwellian molecules. In contrast to the standard splitting methods, where one works with the instationary equation, the current approach simulates the direct solution of the stationary problem. The model problem investigated is the heat transfer between two parallel plates in the rarefied gas regime. An iteration process is introduced which leads to the stationary solution of the exact - space discretized - Boltzmann equation, in the sense of weak convergence.
As an alternative to the commonly used Monte Carlo Simulation methods for solving the Boltzmann equation we have developed a new code with certain important improvements. We present results of calculations on the reentry phase of a space shuttle. One aim was to test physical models of internal energies and of gas-surface interactions.
The paper presents some approximation methods for the Boltzmann equation. In the first part fully implicit discretization techniques for the spatially homogeneous Boltzmann equation are investigated. The implicit equation is solved using an iteration process. It is shown that the iteration converges to the correct solution for the moments of the distribution function as long as the mass conservation is strictly fulfilled. For a simple model Boltzmann equation some unexpected features of the implicit scheme and the corresponding iteration process are clarified. In the second part a new iteration algorithm is proposed which should be used for the stationary Boltzmann equation. The realization of the method is very similar to the standard splitting algorithms except some new stochastic elements.
Monte-Carlo methods are widely used numerical tools in various fields of application, like rarefied gas dynamics, vacuum technology, stellar dynamics or nuclear physics. A central part in all applications is the generation of random variates according to a given probability law. Fundamental techniques to generate non-uniform random variates are the inversion principle or the acceptance-rejection method. Both procedures can be quite time-consuming if the given probability law has a complicated structure.; In this paper we consider probability laws depending on a small parameter and investigate the use of asmptotic expansions to generate random variates. The results given in the paper are restrictedto first order expansions. We show error estimates for the discrepancy as well as for the bounded Lipschitz distance of the asymptotic expansion. Furthermore the integration error for some special classes of functions is given. The efficiency of the method is proved by a numerical example from rarefied gas flows.
The paper presents a fast implementation of a constructive method to generate a special class of low-discrepancy sequences which are based on Van Neumann-Kakutani tranformations. Such sequences can be used in various simulation codes where it is necessary to generate a certain number of uniformly distributed random numbers on the unit interval.; From a theoretical point of view the uniformity of a sequence is measured in terms of the discrepancy which is a special distance between a finite set of points and the uniform distribution on the unit interval.; Numerical results are given on the cost efficiency of different generators on different hardware architectures as well as on the corresponding uniformity of the sequences. As an example for the efficient use of low-discrepancy sequences in a complex simulation code results are presented for the simulation of a hypersonic rarefied gas flow.
Based on general partitions of unity and standard numerical flux functions, a class of mesh-free methods for conservation laws is derived. A Lax-Wendroff type consistency analysis is carried out for the general case of moving partition functions. The analysis leads to a set of conditions which are checked for the finite volume particle method FVPM. As a by-product, classical finite volume schemes are recovered in the approach for special choices of the partition of unity.