A domain decomposition scheme linking linearized kinetic and aerodynamic equations is considered. Convergence of the alternating scheme is shown. Moreover the physical correctness of the obtained coupled solutions is discussed.
Linear half-space problems can be used to solve domain decomposition problems between Boltzmann and aerodynamic equations. A new fast numerical method computing the asymptotic states and outgoing distributions for a linearized BGK half-space problem is presented. Relations with the so-called variational methods are discussed. In particular, we stress the connection between these methods and Chapman-Enskog type expansions.
A nonequilibrium situation governed by kinetic equations with strongly contrasted Knudsen numbers in different subdomains is discussed. We consider a domain decomposition problem for Boltzmann- and Euler equations, establish the correct coupling conditions and prove the validity of the obtained coupled solution. Moreover numerical examples comparing different types of coupling conditions are presented.
As it is well known in statistical physics the stationary distribution can be obtained by maximizing entropy. We show how one can reconstruct the formula for entropy knowing the formula for the stationary distribution. A general case is discussed and some concrete physical examples are considered.
We consider the numerical computation of nonlinear functionals of distribution functions approximated by point measures. Two methods are described and estimates for the speed of convergence as the number of points tends to infinity are given. Moreover numerical results for the entropy functional are presented.
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.
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 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.
Based on state-space formulas for coprime factorizations over ... and an algebraic characterization of J-inner functions, normalized doubly-coprime factorizations for different classes of continuous- and discrete-time transfer functions are derived by using a single general construction method. The parametrization of the factors is in terms of the stabilizing solutions of general degenerate continuous- respectively discrete-time Riccati equations, which are obtained by examining state-space representations of J-normalized factor matrices.
This paper deals with domain decomposition methods for kinetic and drift diffusion semiconductor equations. In particular accurate coupling conditions at the interface between the kinetic and drift diffusion domain are given. The cases of slight and strong nonequilibrium situations at the interface are considered and some numerical examples are shown.