Skip to main content Accessibility help

A Scalable Numerical Method for Simulating Flows Around High-Speed Train Under Crosswind Conditions

  • Zhengzheng Yan (a1), Rongliang Chen (a1), Yubo Zhao (a1) and Xiao-Chuan Cai (a2)


This paper presents a parallel Newton-Krylov-Schwarz method for the numerical simulation of unsteady flows at high Reynolds number around a high-speed train under crosswind. With a realistic train geometry, a realistic Reynolds number, and a realistic wind speed, this is a very challenging computational problem. Because of the limited parallel scalability, commercial CFD software is not suitable for supercomputers with a large number of processors. We develop a Newton-Krylov-Schwarz based fully implicit method, and the corresponding parallel software, for the 3D unsteady incompressible Navier-Stokes equations discretized with a stabilized finite element method on very fine unstructured meshes. We test the algorithm and software for flows passing a train modeled after China’s high-speed train CRH380B, and we also compare our results with results obtained from commercial CFD software. Our algorithm shows very good parallel scalability on a supercomputer with over one thousand processors.


Corresponding author



Hide All
[1]Ahmed, S. R., Ramm, G., and Faltin, G., Some salient features of the time-averaged ground vehicle wake, SAE paper, 840300, 1984.
[2]Balay, S., Buschelman, K., Gropp, W. D., Kaushik, D., Knepley, M., Mc, L. C.Innes, Smith, B. F., and, Zhang, H., PETSc Users Manual, Technical Report, Argonne National Laboratory, 2012.
[3]Cai, X.-C. and Sarkis, M., A restricted additive Schwarz preconditioner for general sparse linear systems, SIAM J. Sci. Comput., 21 (1999), 792797.
[4]Chen, R. and Cai, X.-C., Parallel one-shot Lagrange-Newton-Krylov-Schwarz algorithms for shape optimization of steady incompressible flows, SIAM J. Sci. Comput., 34 (2012), 584605.
[5]Chen, R., Wu, Y., Yan, Z., Zhao, Y., and Cai, X.-C., A parallel domain decomposition method for 3D unsteady incompressible flows at high Reynolds number, J. Sci. Comput., (to appear).
[6]Diedrichs, B., Studies of Two Aerodynamic Effects on High-Speed Trains: Crosswind Stability and Discomforting Car Body Vibrations Inside Tunnels, Ph.D. Thesis, Royal Institute of Technology (KTH), 2006.
[7]Favre, T., Diedrichs, B., and Efraimsson, G., Detached-Eddy simulations applied to unsteady crosswind aerodynamics of ground vehicles, Progress in Hybrid RANS-LES Modelling, Springer, 2010, 167177.
[8]Franca, L. P. and Frey, S. L., Stabilized finite element method: II. The incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Engrg., 99 (1992), 209233.
[9]Guilmineau, E., Computational study of flow around a simplified car body, J. Wind. Eng. Ind. Aerod., 96 (2008), 12071217.
[10]Hemida, H. and Baker, C., Large-eddy simulation of the flow around a freight wagon subjected to a crosswind, Comput. Fluids, 39 (2010), 19441956.
[11]Hemida, H. and Krajnovic, S., LES study of the impact of the wake structures on the aerodynamics of a simplified ICE2 train subjected to a side wind, presented at the 4th International Conference on Computational Fluid Dynamics, Ghent, Belgium, July 1014, 2006.
[12]Hu, Q. and Zou, J., Nonlinear inexact Uzawa algorithms for linear and nonlinear saddle-point problems, SIAM J. Optim., 16 (2006), 798825.
[13]Hwang, F.-N. and Cai, X.-C., A parallel nonlinear additive Schwarz preconditioned inexact Newton algorithm for incompressible Navier-Stokes equations, J. Comput. Phys., 204 (2005), 666691.
[14]Hwang, F.-N., Wu, C.-Y., and Cai, X.-C., Numerical simulation of three-dimensional blood flows using domain decomposition method on parallel computer, J. CSME., 31 (2010), 199208.
[15]Karypis, G., Aggarwal, R., Schloegel, K., Kumar, V., and Shekhar, S., ParMETIS home page,
[16]Khier, W., Breuer, M., and Durst, F., Flow structure around trains under side wind conditions: a numerical study, Comput. Fluids, 29 (2000), 179195.
[17]Krajnovic, S., Numerical simulation of the flow around an ICE2 train under the influence of a wind gust, presented at the International Conference on Railway Engineering, Hong Kong, China, March 2528, 2008.
[18]Saad, Y., Iterative Method for Sparse Linear Systems, SIAM, 2003.
[19]Smith, B. F., Bjørstad, P. E., and Gropp, W. D., Domain Decomposition Parallel Multilevel Methods for Elliptic Partial Differential Equations, Cambridge university press, Cambridge, 1996.
[20]Wu, Y. and Cai, X.-C., A parallel two-level method for simulating blood flows in branching arteries with the resistive boundary condition, Comput. Fluids, 45 (2011), 92102.
[21]Yang, G. W., Guo, D. L., Yao, S. B., and Liu, C. H., Aerodynamic design for China new highspeed trains, Sci. China Tech. Sci., 55 (2012), 19231928.
[22]Zhou, D., Tian, H. Q., and Lu, Z. J., Influence of the strong crosswind on aerodynamic performance of passenger train running on embankment, J. Traffic Transp. Eng., 7 (2007), 69.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed