Parameter optimization for turbulent boundary layer generation using ensemble Kalman filter

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.

Development of a surface heat flux model for urban wind simulation using locally mesh-refined lattice Boltzmann method

Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Nakayama, Hiromasa

Dai-35-Kai Suchi Ryutai Rikigaku Shimpojiumu Koen Rombunshu (Internet), 3 Pages, 2021/12

A detailed wind simulation is very important for designing smart cities. Since a lot of tall buildings and complex structures make the air flow turbulent in urban cities, large-scale CFD simulations are needed. We develop a GPU-based CFD code based on a Lattice Boltzmann Method (LBM) with a block-based Adaptive Mesh Refinement (AMR) method. In order to reproduce real wind conditions, the wind condition and ground temperature of the mesoscale weather forecasting model are given as boundary conditions. In this research, a surface heat flux model based on the Monin-Obukhov similarity theory was introduced to improve the calculation accuracy. We conducted a detailed wind simulation in Oklahoma City. By executing this computation, wind conditions in the urban area were reproduced with good accuracy.

Large-eddy simulation of plume dispersion in the central district of Oklahoma City by coupling with a mesoscale meteorological simulation model and observation

Nakayama, Hiromasa; Takemi, Tetsuya*; Yoshida, Toshiya

Atmosphere (Internet), 12(7), p.889_1 - 889_15, 2021/07

https://doi.org/10.3390/atmos12070889

Contaminant gas dispersion within urban area resulting from accidental or intentional release is of great concern to public health and social security. When estimating plume dispersion in a built-up urban area under real meteorological conditions by computational fluid dynamics (CFD), a crucial issue is how to prescribe the model input conditions. There are typically two approaches: using the outputs of a meso-scale meteorological simulation (MMS) model and meteorological observations (OBS). However, the influence of the different approaches on the simulation results have not been fully demonstrated. In this study, we conducted large-eddy simulations (LESs) of plume dispersion in the urban central district under real meteorological conditions by coupling with a MMS model and OBS obtained at a single stationary point, and evaluated the two different coupling simulations in comparison with the field experimental data. The LES-MMS coupling showed better performance than the LES-OBS one. However, the latter one also showed a reasonable performance comparable to the acceptance criteria on the model prediction within a factor of two of the experimental data. These facts indicate that the approach of using observations at a single stationary point still has enough potential to drive CFD models for plume dispersion under real meteorological conditions.

Real-time tracer dispersion simulations in Oklahoma City using the locally mesh-refined lattice Boltzmann method

Onodera, Naoyuki; Idomura, Yasuhiro; Hasegawa, Yuta; Nakayama, Hiromasa; Shimokawabe, Takashi*; Aoki, Takayuki*

Boundary-Layer Meteorology, 179(2), p.187 - 208, 2021/05

https://doi.org/10.1007/s10546-020-00594-x

A plume dispersion simulation code named CityLBM enables a real time simulation for several km by applying adaptive mesh refinement (AMR) method on GPU supercomputers. We assess plume dispersion problems in the complex urban environment of Oklahoma City (JU2003). Realistic mesoscale wind boundary conditions of JU2003 produced by a Weather Research and Forecasting Model (WRF), building structures, and a plant canopy model are introduced to CityLBM. Ensemble calculations are performed to reduce turbulence uncertainties. The statistics of the plume dispersion field, mean and max concentrations show that ensemble calculations improve the accuracy of the estimation, and the ensemble-averaged concentration values in the simulations over 4 km areas with 2-m resolution satisfied factor 2 agreements for 70% of 24 target measurement points and periods in JU2003.

Development of a data assimilation method using vibration equation for Large-Eddy Simulations of turbulent boundary layer flows

Nakayama, Hiromasa; Takemi, Tetsuya*

Journal of Advances in Modeling Earth Systems (Internet), 12(8), p.e2019MS001872_1 - e2019MS001872_18, 2020/08

https://doi.org/10.1029/2019MS001872

We proposed a data assimilation method using a vibration equation which can incorporate turbulence winds toward target mean winds while maintaining small-scale turbulent fluctuations. First, we conducted test simulations in which nudging is applied in a basic turbulent boundary layer (TBL) flow toward a target one. It was shown that the basic TBL flow can be reasonably nudged toward the target one while maintaining the turbulent fluctuations well, when prescribing the natural frequency in the vibration equation smaller than the spectral peak frequency in the TBL flow. Then, we applied the proposed nudging method by incorporating data obtained from meteorological observations located in the southern part of the actual city of Kyoto. The mean wind velocity profiles are reasonably nudged toward the target observed profile and the turbulence statistics are also favourably maintained. It is concluded that the LES data assimilation method successfully nudges toward the turbulence winds well.

Development of the data assimilation method suitable for large-eddy simulation model using the vibration equation

Nakayama, Hiromasa; Takemi, Tetsuya*

Proceedings of 19th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes (HARMO-19) (USB Flash Drive), 5 Pages, 2019/06

We proposed the data assimilation method using the vibration equation suitable for LES models. First, we generated the basic turbulent boundary layer (TBL) flow with a power law exponent of 0.14 at the upstream area and then attempted to nudge toward the TBL flow with a power law exponent of 0.25 by our proposed data assimilation. For a case that the natural frequency set in the vibration equation was larger than the peak frequency of the basic TBL flow, the turbulent fluctuations were rapidly damped although the mean velocity profile matched with the target one. On the other hand, for a case that the natural frequency was smaller than the peak frequency, the mean wind velocity profile was reasonably nudged toward the target profile while maintaining the turbulent fluctuations. It is concluded that our proposed data assimilation method using the vibration equation successfully nudges toward the target mean winds with small-scale turbulent fluctuations.

LES studies for predicting plume concentrations around nuclear facilities using an overlapping technique

Nakayama, Hiromasa; Takemi, Tetsuya*

Proceedings of 18th International Conference on Harmonisation within Atmospheric Dispersion Modelling for Regulatory Purposes (HARMO-18) (USB Flash Drive), p.843 - 847, 2017/10

We have attempted to develop a practical and quick local-scale atmospheric dispersion calculation method using an overlapping technique for plume concentration distributions in an emergency response to nuclear accidents. In order to evaluate the overlapping approach, we performed LESs of turbulent flows and plume dispersion under varying mean wind directions using the meteorological data as the model input. When compared with the simulated results under the real meteorological condition, it is shown that the concentration distribution patterns are reasonably simulated by the overlapping method. It can be concluded that the atmospheric dispersion calculation method using the overlapping technique has potential performance for emergency responses to nuclear accidents.

Local-scale high-resolution atmospheric dispersion model using large-eddy simulation; LOHDIM-LES

Nakayama, Hiromasa; Nagai, Haruyasu

JAEA-Data/Code 2015-026, 37 Pages, 2016/03

JAEA-Data-Code-2015-026.pdf:2.48MB

We developed LOcal-scale High-resolution atmospheric DIspersion Model using Large-Eddy Simulation (LOHDIM-LES). This dispersion model is designed based on LES which is effective to reproduce unsteady behaviors of turbulent flows and plume dispersion. The basic equations are the continuity equation, the Navier-Stokes equation, and the scalar conservation equation. Buildings and local terrain variability are resolved by high-resolution grids with of a few meters and these turbulent effects are represented by immersed boundary method. In simulating atmospheric turbulence, boundary layer flows are generated by a recycling turbulent inflow technique in a driver region set up at the upstream of the main analysis region. This turbulent inflow data are imposed at the inlet of the main analysis region. By this approach, the LOHDIM-LES can provide detailed information on wind velocities and plume concentration in the investigated area.

Large-Eddy Simulation of turbulent winds during the Fukushima Daiichi Nuclear Power Plant accident by coupling with a meso-scale meteorological simulation model

Nakayama, Hiromasa; Takemi, Tetsuya*; Nagai, Haruyasu

Advances in Science & Research (Internet), 12(1), p.127 - 133, 2015/00

https://doi.org/10.5194/asr-12-127-2015

One of the important issues to be solved to couple LES and MM models is to impose time-dependent turbulent inflow data for LESs from the MM model outputs, because the MM models are not able to reproduce high-frequency turbulent fluctuations appropriate to drive LES models. In this study, we apply the recycling technique to couple the CFD and MM models. We conduct an LES of turbulent winds during the Fukushima Daiichi Nuclear Power Plant accident and report the usefulness of our approach by comparing the simulated results with the observations.

Development of local-scale high-resolution atmospheric dispersion model using Large-Eddy Simulation, 4; Turbulent flows and plume dispersion in an actual urban area

Nakayama, Hiromasa; Leitl, B.*; Harms, F.*; Nagai, Haruyasu

Journal of Nuclear Science and Technology, 51(5), p.626 - 638, 2014/02

https://doi.org/10.1080/00223131.2014.885400

We apply the LES model to turbulent flows and plume dispersion in an actual urban area. Although some of the turbulence and dispersion characteristics are quantitatively different from the wind tunnel experimental data, the distribution patterns are generally similar to those of the experiments. It is concluded that our LES model simulates reasonably the unsteady behavior of turbulent flows and plume dispersion even for complex heterogeneous urban areas.