Cases for "A flow decomposition method based on computational fluid dynamics for rock weir head-discharge relationship"

This folder contains all the cases in the paper titled \"A Flow Decomposition Method based on Computational Fluid Dynamics
for Rock Weir Head-Discharge Relationship\".

To demonstrate how to run the cases, the RL10 case is used as an example:
1. The mesh should be generated by executing the script \"runMesh\".
2. The simulation should be run by executing the script \"runCFD\".

All simulations used OpenFOAM v5.x. For this research, all simulations were run in parallel with 160 cores. If run
with less cores, the simulation time might be long.

For other simulation cases, the bed geometries (bed.stl) should be changed. All bed geomerites can be found in the folder
named \"bedSTL\". The names of the files are defined in the paper. The simulated discharge Q can be changed in velocity file \"U\",
specifically the \"volumetricFlowRate\" for the \"inlet\" boundary condition.

Postprocessing:
1. The simulated water surface elevation (WSE.csv) was extracted in ParaView with contour alpha.water=0.5.
2. The simulated discharge was calculated by the integration of U normal to the streamwise direction with the weight of alpha.water.

This folder also contains the python script for the use of the flow decomposition method. In addition, the Matlab code to generate
the figure comparing new results and existing recommended Cd values is also included.

本文件夹包含论文《基于计算流体动力学（Computational Fluid Dynamics）的岩堰水头-流量关系流分解方法》中的所有算例。

为演示算例的运行流程，本文以RL10算例为例：
1. 需通过执行脚本`runMesh`生成计算网格；
2. 需通过执行脚本`runCFD`启动数值模拟。

所有数值模拟均基于OpenFOAM v5.x版本开展，本次研究中所有模拟均采用160核心并行运行；若使用更少核心，模拟耗时可能较长。

对于其他模拟算例，需修改床面几何模型（bed.stl）。所有床面几何模型均存于名为`bedSTL`的文件夹中，文件名已在论文中定义。可通过修改速度文件`U`中的边界条件调整模拟流量Q，具体为`inlet`边界的`volumetricFlowRate`参数。

后处理流程：
1. 利用ParaView软件提取模拟水面高程数据（WSE.csv），提取时选取等值面alpha.water=0.5；
2. 模拟流量通过对速度场U沿流向的法向分量，并以alpha.water为权重积分计算得到。

本文件夹同时包含用于流分解方法的Python脚本。此外，还附带了用于生成对比新结果与现有推荐Cd值图表的Matlab代码。