It has been shown that the equation of diffusion, linear and nonlinear, can be obtained in a suitable scaling limit by a two-velocity model of the Boltzmann equation [7] . Several numerical approximations were introduced in order to discretize the corresponding multiscale hyperbolic systems [8, 1, 4]. In the present work we consider relaxed approximations for multiscale kinetic systems with asymptotic state represented by nonlinear diffusion equations. The schemes are based on a relaxation approximation that permits to reduce the second order diffusion equations to first order semi-linear hyperbolic systems with stiff terms. The numerical passage from the relaxation system to the nonlinear diffusion equation is realized by using semi-implicit time discretization combined with ENO schemes and central differences in space. Finally, parallel algorithms are developed and their performance evaluated. Application to porous media equations in one and two space dimensions are presented.
Parallel algorithms for nonlinear diffusion by using relaxation approximation / F. Cavalli, G. Naldi, M. Semplice - In: Numerical mathematics and advanced applications : Proceedings of ENUMATH 2005, the 6th European Conference on Numerical Mathematics and Advanced Applications, Santiago de Compostela, Spain, July 2005 / [a cura di] A. Bermúdez de Castro, D. Gómez, P. Quintela, P. Salgado. - Berlin : Springer, 2006. - ISBN 978-3-540-34287-8. - pp. 404-411 (( Intervento presentato al 6. convegno ENUMATH tenutosi a Santiago de Compostela nel 2005 [10.1007/978-3-540-34288-5_35].
Parallel algorithms for nonlinear diffusion by using relaxation approximation
F. CavalliPrimo
;G. NaldiSecondo
;M. SempliceUltimo
2006
Abstract
It has been shown that the equation of diffusion, linear and nonlinear, can be obtained in a suitable scaling limit by a two-velocity model of the Boltzmann equation [7] . Several numerical approximations were introduced in order to discretize the corresponding multiscale hyperbolic systems [8, 1, 4]. In the present work we consider relaxed approximations for multiscale kinetic systems with asymptotic state represented by nonlinear diffusion equations. The schemes are based on a relaxation approximation that permits to reduce the second order diffusion equations to first order semi-linear hyperbolic systems with stiff terms. The numerical passage from the relaxation system to the nonlinear diffusion equation is realized by using semi-implicit time discretization combined with ENO schemes and central differences in space. Finally, parallel algorithms are developed and their performance evaluated. Application to porous media equations in one and two space dimensions are presented.
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
