Influence of seeded impurities on the fusion reactor operation

In order to study the influence of seeded impurities on ITER like reactor operation the COREDIV code has been extended to include the transport of several sputtered and/or injected impurities. In the COREDIV code the core plasma is treated in the frame of 1D radial transport model whereas in the edge the 2D multifluid code EPIT is used. The EPIT code solves in the slab geometry the MHD equations for densities and velocities for all ions species as well as for electron and ion temperatures. The iteration scheme in the code leads to steady state solution of coupled core and SOL system. The numerical results for Carbon and Nickel plate and seeded impurities Silicon and Neon are presented.     

Pressure Stabilization Method for Steady Viscous Flows in a System of Pipes

This paper deals with a singular perturbation of the stationary Stokes and Navier-Stokes systems. The term ε²Δp is added to the continuity equation, where ε is a small parameter. For a domain with cylindrical outlets to infinity and exponentially decaying data, existence and uniqueness of solutions under flux conditions at infinity are established for the linear problem and also for the nonlinear problem in the case of small data. Asymptotically exact estimates are proved for ε tending to zero. For sufficiently regular data, these estimates imply the convergence in H _loc ^5/2−δ for the velocity parts and in H _loc ^3/2−δ for the pressure parts, respectively. Bibliography: 17 titles.Dedicated to V. A. Solonnikov on the occasion of his 70th birthday__________Published in Zapiski Nauchnykh Seminarov POMI, Vol. 306, 2003, pp. 107–133.     

Variational Derivative for Buckling Loads of Beams and Frames

The finite-difference method is a numerical technique for obtaining approximate solutions to differential equations. The main objective of the present study is to give a new aspect to the finite-difference method by using a variational derivative. By applying this formulation, accurate values of the buckling loads of beams and frames with various end supports are obtained. The performance of this formulation is verified by comparison with numerical examples in the literature __________ Published in Prikladnaya Mekhanika, Vol. 41, No. 7, pp. 139–144, July 2005.     

Monotone Nonincreasing Random Fields on Partially Ordered Sets. I

For an arbitrary poset H and measure ρ on H × _R (where _R is the real axis), we construct a monotone decreasing stochastic field η_ρ and compute its finite-dimensional distributions. In the case where H is a Λ-semilattice and the measure ρ satisfies additional conditions, we compute various characteristics of the field η_ρ such as the expectation of the field value at a point, variance of the field value at a point, and correlation function of the field. The described construction of random fields gives a new method for constructing positive definite functions on posets. Bibliography: 6 titles.__________Translated from Zapiski Nauchnykh Seminarov POMI, Vol. 301, 2003, pp. 92–143.     

Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.

     

Thermal diffusion of envelope solitons on anharmonic atomic chains

We study the motion of envelope solitons on anharmonic atomic chains in the presence of dissipation and thermal fluctuations. We consider the continuum limit of the discrete system and apply an adiabatic perturbation theory which yields a system of stochastic integro-differential equations for the collective variables of the ansatz for the perturbed envelope soliton. We derive the Fokker-Planck equation of this system and search for a statistically equivalent system of Langevin equations, which shares the same Fokker-Planck equation. We undertake an analytical analysis of the Langevin system and derive an expression for the variance of the soliton position Var[x _s ] which predicts a stronger than linear time dependence of Var[x _s ] (superdiffusion). We compare these results with simulations for the discrete system and find they agree well. We refer to recent studies where the diffusion of pulse solitons were found to exhibit a superdiffusive behaviour on longer time scales.Received: 28 June 2004, Published online: 26 November 2004PACS: 05.10.Gg Stochastic analysis methods - 05.45.Yv Solitons - 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.50. + q Lattice theory and statistics