Computational Fluid Dynamics Laboratory
Civil Engineering, National Ilan University, Taiwan.
Welcome to the research group of the CFD Lab under the direction of Prof. CHING-SEN WU.
The research of the CFD Lab focuses on the study of fundamental phenomena in complex flows. Our studies are motivated by environmental hydraulics and water resource engineering carried out with experimental observations and with developing and running Direct Numerical Simulation (DNS) or Large Eddy Simulation (LES) solvers. Our most recent research has focused on:
-
Density currents
At the present time, the members in our group include 4 master students and a few undergraduate students.
For more detailed information, please visit our research pages using the menu above.
27th National Computational Fluid Dynamics Conference
Aug. 11-12 2022
46th Conference on Theoretical and Applied Mechanics, CTAM 2022
Nov. 18-19 2022
The 44th Ocean Engineering Conference in Taiwan
Nov. 17-18 2022
75th Annual Meeting of the APS Division of Fluid Dynamics
Nov. 20-22 2022
Indianapolis, USA
UPCOMING ACADEMIC EVENTS
LATEST RESEARCH
-
Gravity Current over an Array of Densified Obstacles
The evolution of bottom-propagating gravity currents with the presence of an array of densified obstacles submerged in a channel is investigated using large-eddy simulations. Our attention is particularly focused on the flow transition of gravity currents over rough surfaces with extra resistance that provokes significant dissipative processes. Two geometric parameters of the roughness elements, namely, the submergence ratio of the obstacle and the array density ratio, govern their kinematic and dynamic effects on the propagation of gravity currents. Physically, the submergence ratio plays a significant role in the control of the current diversion, and the array density regulates the flow pathway of gravity current propagation. The integrated measures show that two distinct flow morphologies are identified.
-
Interactions Between a Gravity Current and an Immersed Obstacle
The interaction of bottom-propagating gravity currents with immersed obstacles in the path of the lock-exchange configuration was investigated based on the large-eddy simulations. The three-dimensional Navier-Stokes solver was quantitatively employed to resolve the flow structure of gravity currents and their dynamics during impact. The integral measure of analysis comprises the front condition, the energy budget, the turbulent mixing, and the force response. Notably, the integrated analysis of gravity currents provides insights into physical mechanisms by identifying distinct propagation stages during transitions, including impact stage, transient stage, and quasi-steady stage.
-
Stratification of Gravity Currents in the Slumping Phase
Experiments on gravity currents produced from full-depth two-layer stratified buoyancy sources propagating in the slumping phase are presented. The flow morphology of two-layer stratfied gravity currents in the slumping phase can be categorized into two different regimes.
-
Unstable Cylindrical Gravity Currents
High-resolution simulations of unstable cylindrical gravity currents when wandering and splitting motions occur in a rotating system are reported. In this study, our attention is focused on the situation of unstable rotating cylindrical gravity currents when the ratio of Coriolis to inertia forces is larger, namely, 0.5 ≤ C ≤ 2.0. A striking feature of the unstable rotating cylindrical gravity currents is the persistent upwelling and downwelling motions, which are observed for both the vortex-wandering and vortex splitting motions.
-
Cylindrical Gravity Currents in a Rotating System
Cylindrical gravity currents, produced by a full-depth lock release, in a rotating system are investigated by means of three-dimensional high-resolution simulations of the incompressible variable-density Navier–Stokes equations with the Coriolis term and using the Boussinesq approximation for a small density difference. Our attention is focused on the situation when the ratio of Coriolis to inertia forces is not large, namely 0.1<C<0.3, and the non-rotating case, namely C=0, is also briefly considered. Depending on the ratio of Coriolis to inertia forces, such a contraction–relaxation motion may be initiated after or before the attainment of a maximum radius of propagation.