Rey DeLeon

Large-Eddy Simulations of Turbulent Incompressible Flows on GPU Clusters

A dual-level parallel incompressible flow solver accelerates turbulent flow computations on GPU clusters. This approach solves the pressure Poisson equation with a full-depth, amalgamated parallel geometric multigrid method, and implements a Lagrangian dynamic model for subgrid-scale turbulence modeling.

An Immersed Boundary Geometric Preprocessor for Arbitrarily Complex Terrain and Geometry

AbstractThere is a growing interest to apply the immersed boundary method to compute wind fields over arbitrarily complex terrain. The computer implementation of an immersed boundary module into an existing flow solver can be accomplished with minor modifications to the rest of the computer program. However, a versatile preprocessor is needed at the first place to extract the essential geometric information pertinent to the immersion of an arbitrarily complex terrain inside a 3D Cartesian mesh. Errors in the geometric information can negatively impact the correct implementation of the immersed boundary method as part of the solution algorithm. Additionally, the distance field from the terrain is needed to implement various subgrid-scale turbulence models and to initialize wind fields over complex terrain. Despite the popularity of the immersed boundary method, procedures used in the geometric preprocessing stage have received less attention. The present study found that concave and convex regions of compl...

Toward a GPU-Accelerated Immersed Boundary Method for Wind Forecasting Over Complex Terrain

A short-term wind power forecasting capability can be a valuable tool in the renewable energy industry to address load-balancing issues that arise from intermittent wind fields. Although numerical weather prediction models have been used to forecast winds, their applicability to micro-scale atmospheric boundary layer flows and ability to predict wind speeds at turbine hub height with a desired accuracy is not clear. To address this issue, we develop a multi-GPU parallel flow solver to forecast winds over complex terrain at the micro-scale, where computational domain size can range from meters to several kilometers. In the solver, we adopt the immersed boundary method and the Lagrangian dynamic large-eddy simulation model and extend them to atmospheric flows. The computations are accelerated on GPU clusters with a dual-level parallel implementation that interleaves MPI with CUDA. We evaluate the flow solver components against test problems and obtain preliminary results of flow over Bolund Hill, a coastal hill in Denmark.Copyright

GPU-Accelerated Large-Eddy Simulation of Turbulent Channel Flows

High performance computing clusters that are augmented with cost and power efficient graphics processing unit (GPU) provide new opportunities to broaden the use of large-eddy simulation technique to study high Reynolds number turbulent flows in fluids engineering applications. In this paper, we extend our earlier work on multi-GPU acceleration of an incompressible Navier-Stokes solver to include a large-eddy simulation (LES) capability. In particular, we implement the Lagrangian dynamic subgrid scale model and compare our results against existing direct numerical simulation (DNS) data of a turbulent channel flow at Reτ = 180. Overall, our LES results match fairly well with the DNS data. Our results show that the Reτ = 180 case can be entirely simulated on a single GPU, whereas higher Reynolds cases can benefit from a GPU cluster.

Direct Numerical Simulation of Turbulent Katabatic Slope Flows with an Immersed-Boundary Method

We investigate a Cartesian-mesh immersed-boundary formulation within an incompressible flow solver to simulate laminar and turbulent katabatic slope flows. As a proof-of-concept study, we consider four different immersed-boundary reconstruction schemes for imposing a Neumann-type boundary condition on the buoyancy field. Prandtl’s laminar solution is used to demonstrate the second-order accuracy of the numerical solutions globally. Direct numerical simulation of a turbulent katabatic flow is then performed to investigate the applicability of the proposed schemes in the turbulent regime by analyzing both first- and second-order statistics of turbulence. First-order statistics show that turbulent katabatic flow simulations are noticeably sensitive to the specifics of the immersed-boundary formulation. We find that reconstruction schemes that work well in the laminar regime may not perform as well when applied to a turbulent regime. Our proposed immersed-boundary reconstruction scheme agrees closely with the terrain-fitted reference solutions in both flow regimes.

Simulations of Turbulent Flow Over Complex Terrain Using an Immersed-Boundary Method

We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin–Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.