Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Hasegawa, Yuta; Onodera, Naoyuki; Asahi, Yuichi; Idomura, Yasuhiro
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 5 Pages, 2022/12
This study implemented and tested the ensemble data assimilation (DA) of turbulent flows using the lattice Boltzmann method and the local ensemble transform Kalman filter (LBM-LETKF). The computational code was implemented fully on GPUs. The test was carried out for the 3D turbulent flow around a square cylinder with 
meshes and 32 ensemble members using 32 GPUs. The time interval of the DA in the test was a half of the period of the Kalman vortex shedding. The normalized mean absolute errors (NMAE) of the lift coefficient were 132%, 148%, and 13.2% for the non-DA case, the nudging case (a simpler DA algorithm), and the LETKF case, respectively. It was found that the LETKF achieved good DA accuracy even though the observation was not frequent enough for the small scale turbulence, while the nudging showed systematic delays in its solution, and could not keep the DA accurately.
Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Nakayama, Hiromasa
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 3 Pages, 2022/12
We have developed a wind simulation code named CityLBM to realize wind digital twins. Mesoscale wind conditions are given as boundary conditions in CityLBM by using a nudging data assimilation method. It is found that conventional approaches with constant nudging coefficients fail to reproduce turbulent intensity in long time simulations, where atmospheric stability conditions change significantly. We propose a dynamic parameter optimization method for the nudging coefficient based on an ensemble Kalman filter. CityLBM was validated against plume dispersion experiments in the complex urban environment of Oklahoma City. The nudging coefficient was updated to reduce the error of the turbulent intensity between the simulation and the observation. The mean error of velocity variance is reduced by 
10% compared to the conventional nudging method with a constant nudging coefficient.
Asahi, Yuichi; Padioleau, T.*; Latu, G.*; Bigot, J.*; Grandgirard, V.*; Obrejan, K.*
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 8 Pages, 2022/12
We implement a kinetic plasma simulation code with multiple performance portable frameworks and evaluated its performance on Intel Icelake, NVIDIA V100 and A100 GPUs, and AMD MI100 GPU. Relying on the language standard parallelism stdpar and proposed language standard multi-dimensional array support mdspan, we demonstrate a performance portable implementation without harming the readability and productivity. With stdpar, we obtain a good overall performance for a kinetic plasma mini-application in the range of 
20% to the Kokkos version on Icelake, V100, A100 and MI100. We conclude that stdpar can be a good candidate to develop a performance portable and productive code targeting Exascale era platforms, assuming this programming model will be available on AMD and/or Intel GPUs in the future.
Sugihara, Kenta; Onodera, Naoyuki; Idomura, Yasuhiro; Yamashita, Susumu
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 5 Pages, 2022/12
The conventional Allen-Cahn type multi-phase field method was modified to conserve not only the sum of the masses of all phases but also the mass of each phase. The interface advection calculations within a two-dimensional rotational velocity field were performed as a verification problem, and the conservation was successfully achieved. The proposed method was used to calculate the horizontally aligned pair of bubbles rising, and it was found that the bouncing phenomenon between bubbles can be calculated at 1/50 resolution of the high-resolution calculation by Zhang et al. using the volume of fluid method.
Kawamura, Takuma; Sakamoto, Naohisa*
Dai-36-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 3 Pages, 2022/12
Volume rendering is useful for visualizing computer fluid dynamics (CFD) data, and its VR visualization helps to understand complex 3D data. Volume rendering of large scale data in remote locations in VR space is an important issue of the visualization field. A remote visualization application CS-PBVR can interactively visualize the large-scale datasets in remote locations with volume rendering. In order to extend CS-PBVR into VR-PBVR which is applicable to a head mount display (HMD) Oculus rift S, we added a stereo image generation function, a gesture control function, and a renewed processing flow. VR-PBVR achieved interactive visualization of remotely located test dataset (2M cells) with 90 fps.
Kamiya, Tomohiro; Yoshida, Hiroyuki
no journal, ,
We proposed equations of conservative form using a VOF function from a conservation law and the discretization method of these equations for liquid-gas two-phase compressible fluid simulations based on a sharp-interface model. The outline of these proposals and results of one-dimensional numerical tests will be reported.