Efficient GPU implementation of a two waves TVD-WAF method for the two-dimensional one layer shallow water system on structured meshes. de la Asunción, M., Castro Díaz, M. J., Fernández Nieto, E. D., M. Mantas, J., Ortega-Acosta, S., & González-Vida, J. 80:441–452.
Efficient GPU implementation of a two waves TVD-WAF method for the two-dimensional one layer shallow water system on structured meshes [link]Paper  abstract   bibtex   
The numerical solutions of shallow water equations are useful for applications related to geophysical flows that usually take place in large computational domains and could require real time calculation. Therefore, parallel versions of accurate and efficient numerical solvers for high performance platforms are needed to be able to deal with these simulation scenarios in reasonable times. In this paper we present an efficient CUDA implementation of a first and second order HLL methods and a two-waves TVD-WAF one. We propose to write all these methods under a common framework, such as, their CUDA implementations share the same structure. In particular, the reformulation of TVD-WAF numerical flux and the improved definition of the flux limiter allows us to obtain a more robust solver in situations like wet/dry fronts. Finally, some numerical tests are presented showing that the TVD-WAF method is slightly slower that the first order HLL method and two times faster than the second order HLL method, but it provides numerical results almost as accurate as the second order HLL scheme.
@article{de_la_asuncion_efficient_2013,
	title = {Efficient {GPU} implementation of a two waves {TVD}-{WAF} method for the two-dimensional one layer shallow water system on structured meshes},
	volume = {80},
	url = {http://dx.doi.org/10.1016/j.compfluid.2012.01.012},
	abstract = {The numerical solutions of shallow water equations are useful for applications related to geophysical flows that usually take place in large computational domains and could require real time calculation. Therefore, parallel versions of accurate and efficient numerical solvers for high performance platforms are needed to be able to deal with these simulation scenarios in reasonable times. In this paper we present an efficient {CUDA} implementation of a first and second order {HLL} methods and a two-waves {TVD}-{WAF} one. We propose to write all these methods under a common framework, such as, their {CUDA} implementations share the same structure. In particular, the reformulation of {TVD}-{WAF} numerical flux and the improved definition of the flux limiter allows us to obtain a more robust solver in situations like wet/dry fronts. Finally, some numerical tests are presented showing that the {TVD}-{WAF} method is slightly slower that the first order {HLL} method and two times faster than the second order {HLL} method, but it provides numerical results almost as accurate as the second order {HLL} scheme.},
	pages = {441--452},
	journaltitle = {Computers \& Fluids},
	author = {de la Asunción, Marc and Castro Díaz, Manuel J. and Fernández Nieto, E. D. and M. Mantas, José and Ortega-Acosta, Sergio and González-Vida, J.-M.},
	date = {2013},
}

Downloads: 0