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 survey contains a description of different types of mathematical models used for the simulation of vehicular traffic. It includes models based on ordinary differential equations, fluid dynamic equations and on equations of kinetic type. Connections between the different types of models are mentioned. Particular emphasis is put on kinetic models and on simulation methods for these models.
In the present paper a general criticism of kinetic equations for vehicular traffic is given. The necessity of introducing an Enskog-type correction into these equations is shown. An Enskog-line kinetic traffic flow equation is presented and fluid dynamic equations are derived. This derivation yields new coefficients for the standard fluid dynamic equations of vehicular traffic. Numerical simulations for inhomogeneous traffic flow situations are shown together with a comparison between kinetic and fluid dynamic models.
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 some adaptive load balance techniques for the simulation of rarefied gas flows on parallel computers. It is shown that a static load balance is insufficient to obtain a scalable parallel efficiency. Hence, two adaptive techniques are investigated which are based on simple algorithms. Numerical results show that using heuristic techniques one can achieve a sufficiently high efficiency over a wide range of different hardware platforms.
We study a model for learning periodic signals in recurrent neural networks proposed by Doya and Yoshizawa  that can be considered as a model for temporal pattern memory in animal motoric systems. A network receives an external oscillatory input and adjusts its weights so that this signal can be reproduced approximately as the network output after some time. We use tools from adaptive control theory to derive an algorithm for weight matrices with special structure. If the input is generated by a network of the same structure the algorithm converges globally and does not exhibit the deficiencies of the back-propagation based approach of Doya and Yoshizawa under a persistency of excitation condition. This simple algorithm can also be used for open loop identification under quite restructive assumptions. The persistency of excitation condition cannot be proven even for the matrices with special structure but for a 3d system. For higher dimensional systems we give connections to the theory of linear time-varying systems where this condition is generically true (under assumption which are also needed in the time-invariant case). However, we cannot show that the linearized system related to the nonlinear neural network fulfills these generic assumptions.
In the modeling of biological phenomena, in living organisms whether the measurements are of blood pressure, enzyme levels, biomechanical movements or heartbeats, etc., one of the important aspects is time variation in the data. Thus, the recovery of a "smooth" regression or trend function from noisy time-varying sampled data becomes a problem of particular interest. Here we use non-linear wavelet thresholding to estimate a regression or a trend function in the presence of additive noise which, in contrast to most existing models, does not need to be stationary. (Here, nonstationarity means that the spectral behaviour of the noise is allowed to change slowly over time.). We develop a procedure to adapt existing threshold rules to such situations, e.g., that of a time-varying variance in the errors. Moreover, in the model of curve estimation for functions belonging to a Besov class with locally stationary errors, we derive a near-optimal rate for the L2-risk between the unknown function and our soft or hard threshold estimator, which holds in the general case of an error distribution with bounded cumulants. In the case of Gaussian errors, a lower bound on the asymptotic minimax rate in the wavelet coefficient domain is also obtained. Also it is argued that a stronger adaptivity result is possible by the use of a particular location and level dependent threshold obtained by minimizing Stein's unbiased estimate of the risk. In this respect, our work generalizes previous results, which cover the situation of correlated, but stationary errors. A natural application of our approach is the estimation of the trend function of nonstationary time series under the model of local stationarity. The method is illustrated on both an interesting simulated example and a biostatistical data-set, measurements of sheep luteinizing hormone, which exhibits a clear nonstationarity in its variance.