Grant: Discovery Grants Program - Individual

Canada ranks as the second country in the world with the most tornadoes per year. This disastrous weather may inflict colossal property damage and even take numerous human lives. While computer technology advancement promises to aid computational analysis of tornado dynamics, several challenges are currently found in this research area. First, when using conventional numerical methods to solve the Navier-Stokes equations based tornado model, the common drawback lies in the macroscopic nature of this flow model, leading to difficulties in capturing finer flow behaviors necessary for a reliable tornado dynamics study. Second, commercial computational fluid dynamics (CFD) tools were found unable to attain satisfactory simulations of complex tornado phenomena, because their limited interface often prevents the user from investigating tornado details in a desirable manner. Third, parallel processing is considered compulsory in dealing with large-scale computations while this technology has been little employed in available in-house code specifically used for the tornado dynamics study. To overcome these major challenges, this research work is built upon the gas kinetics based lattice Boltzmann model (LBM) owing to its inherent parallelizability for algorithmic development and noticeable extensibility towards its coupling with other advanced numerical tools. OpenLB, a free open-source LBM software package with strong parallel execution capabilities, will serve as the basis code, and the following major new research tasks aimed at economical and more reliable tornado simulations will be conducted: (1) Develop a finite volume (FV) discretization of the lattice Boltzmann model, so that many constraints in the LBM context can be removed, leading to the use of unstructured mesh and implicit solution methods, which are essentially demanded when a highly turbulent tornadic wind is numerically investigated. (2) Develop an innovative immersed boundary (IB) approach based tornado-building interaction model in order to avoid updating the velocity at the outer boundary of computational domain when the tornado evolves with time. (3) Develop new IB-FV-LBM adaptive meshing tools powered by a novel a posteriori error estimator corresponding to the LES based turbulence model in the 3-D context, and embed it in the code to improve the reliability of simulation results. (4) Develop a pertinent discrete element model to enable the air-debris interaction for a more realistic study of tornado dynamics. With the above features embedded in OpenLB, the resulting adaptive hybrid CFD tools will be able to perform a large series of tornado-building interaction tests. Valuable insight and guidance will be gained for better understanding tornado dynamics and improving the wind-resistant capabilities in building design. The training of HQP through this program will also benefit Canada with future wind engineering innovations.