《Water》：Numerical Investigation of Tank Geometry and Hydraulic Forcing to Increase Ammonia Water Renewal Efficiency in Stratified Aquaculture Earth Tanks
编辑推荐：
在水产养殖中，高效的水更新系统对于保障水质和防止氨积累至关重要。本研究采用二维垂向（2DV）有限元模型中的流函数-涡量（streamfunction-vorticity）公式，模拟基于SISdATA项目中Figueira半集约化鱼池的实验室尺度养殖水槽内浮
在水产养殖中，高效的水更新系统对于保障水质和防止氨积累至关重要。本研究采用二维垂向（2DV）有限元模型中的流函数-涡量（streamfunction-vorticity）公式，模拟基于SISdATA项目中Figueira半集约化鱼池的实验室尺度养殖水槽内浮力驱动的水更新与氨输运过程。该模型通过Boussinesq近似将不可压缩流动与氨的对流-扩散耦合，并采用流线迎风/Petrov-Galerkin（SUPG）稳定化与算子分裂（operator-splitting）格式，以确保在对流主导输运条件下的数值稳定性。为分离几何形状对更新性能的影响，研究了五种水槽构型，其长度相同（m）而水深由0.0227 m逐步增至0.350 m；相应进水口高度为0.00454～0.070 m，同时保持雷诺数与Péclet数不变。在给定流量下，较深水槽导致较低的进水口流速、环流对污染层穿透的延迟，以及环流控制型更新相对于射流直接冲刷的相对重要性上升。在基准条件下，达到83%氨去除所需时间由Tank1的48 min至Tank4的208 min不等，而最深的Tank5因后期再循环增强，91 min即实现83%去除。水力缩放进一步表明，维持恒定进水口流速可使Tank5的更新时间缩短至75 min，而固定进、出水口开度条件总体上有利于最浅构型。定容分析表明，将水槽长深比由110.13降至0.46，会使83%更新时间由48 min增至353 min。结果表明，水更新效率受水槽几何、水力边界条件与长深比共同控制。所提出的模型为优化养殖水槽设计与水力运行以改善氨去除与水质提供了数值工具。
## 论文解读
水产养殖是增长最快的食品生产部门之一，对满足全球日益增长的水生蛋白需求具有重要意义。然而集约化进程引发了关于用水、废物管理与环境可持续性的担忧。循环水养殖系统（RAS）被视为应对上述挑战的有效生产技术，可最大限度减少用水并实现受控养殖条件。Martins等总结了欧洲RAS的最新进展，指出反硝化反应器、污泥处理技术与营养物循环利用是提升系统可持续性的重要因素；Davidson等证明了陆基淡水RAS的可行性并成功养殖上市规格大西洋鲑；其后续研究将膜生物反应器引入淡水RAS，使淡水需求降低6倍以上而鱼类性能保持稳定；Zhang等构建了包含生物塘、人工湿地与生态沟渠的综合RAS系统，获得更佳水质与更高产量。尽管技术不断进步，废物积累仍是集约化与半集约化养殖中的重大运行难题：Dauda等指出未食饲料、鱼类排泄物与代谢废物是导致水质恶化的主要养殖废物来源；Liu等讨论了养殖尾水中的污染物，强调实施更优尾水处理与更新策略的必要性。水质有效监测与管理同样必要，Lindholm-Lehto综述了现代传感器技术，Shubairi等指出预测性水动力与水质模型是未来数字孪生与决策支持应用的关键组成部分，Zuhaer等提出了基于物联网的实时水质监测系统。数值模拟因而成为理解养殖设施内流动与输运过程的重要工具：Behroozi与Couturier建立圆形水槽CFD模型，表明环流对水力性能有显著影响；Gorle等研究浮式封闭养殖场，证明内部流动结构显著影响养殖单元环境条件；Covas等通过实验与数值模拟量化储水槽中的混合与更新，强调水槽几何与进/出水口配置的重要性。进一步CFD研究聚焦于水力优化：Lopez-Rebollar等报道环形沉淀器对圆形水槽混合特性与固体去除效率有显著影响；Zhao等数值考察直径-水深比对循环水养殖水槽流动环流与颗粒排放的影响，发现圆形水槽存在使水力性能与颗粒排放效率较优的长深比；Zhang等报道水槽几何与排水配置是影响双排水养殖系统环流模式的强因子。近期CFD日益用于水质分析：Yu等将CFD与多相流及组分输运模型耦合以最大化曝气效率并降低能耗；Haris等表明外驱动流动与浮力的相互作用可显著改变二维空腔系统中的环流与分层。尽管水动力与水质建模已有进展，明确考虑氨输运与更新的研究仍很少：Rattanakul等证明了平流-扩散公式在水生环境中模拟输运过程的适用性；Jiang等建立了氨动力学数学模型并讨论了理解氮积累与输运的预测工具需求。流函数-涡量公式相较原始变量公式是模拟不可压缩流动的高效途径，Erturk等确立了其充分捕捉高雷诺数下复杂涡结构的潜力，Bonaventura等进一步采用半拉格朗日与有限元方法改进该公式，使其对输运主导流动更稳定、更精确。水槽与工业反应器中溶解气体输运的数值模拟为水动力、气液传质与组分输运之间的耦合提供了重要认识：Matko等将两相CFD与组分输运耦合模拟氧化沟中气液流动、氧传质与溶解氧分布；Shah等开发并实验验证了全尺度工业废水曝气池的多相CFD模型；Yu等将Euler-Euler多相公式与组分输运模型耦合模拟养殖水槽中溶解氧的时空演化；Moilanen等研究曝气发酵罐中的气液水动力与传质；Mallick等开发了基于界面面积浓度模型的高效CFD框架预测曝气搅拌槽反应器中的氧传质。这些研究表明水中溶解组分的输运强烈受制于底层流场与传质机制，为本研究提供了方法学背景。然而，水更新过程中水动力环流、分层、氨输运与水槽几何的耦合效应尚未被充分理解，水深与长深比等几何特征对环流模式与氨冲洗效率的修正程度尚未被系统量化。为此，研究人员建立了基于流函数-涡量公式的二维垂向（2DV）数值模型以研究分层养殖水槽中耦合的水动力与氨输运，其方法学贡献在于在开源FreeFEM-cs软件中，将浮力效应的Boussinesq近似、对流主导氨输运的SUPG稳定化有限元格式以及半隐式算子分裂策略相结合，并特别强调水槽几何作为控制参数。研究提出四个研究问题：（1）在雷诺数与Péclet数保持恒定时，水深如何影响环流结构与氨更新？（2）不同水力边界条件缩放方式（恒定进水口流速、固定进水口或出水口开度）如何改变氨去除性能？（3）在定容条件下，长深比如何影响达到特定氨更新水平所需时间？（4）有利环流模式能在多大程度上补偿水深增加对氨冲洗的负面影响？
在技术方法上，研究采用二维垂向矩形计算域，尺寸依据SISdATA项目所研究的Figueira鱼池确定；下部区域初始为氨污染水，上部为淡水并设置进、出水边界；进口层厚度取为0.2倍总水深，恒定单宽流量驱动污染层冲刷。流动由连续性方程与含浮力的Navier-Stokes方程描述，密度变化采用Boussinesq线性化；氨输运由对流-扩散方程控制，其中预设的表观标量扩散系数取为水温约室温条件下淡水运动黏度的5倍，属未经率定的建模假设。为精确满足连续性并提高计算效率，方程重写为流函数-涡量公式，得到联系流函数与涡量的Poisson方程、含黏性扩散与浮力强迫的涡量输运方程以及氨输运方程。初始条件为静止流体与以平滑双曲正切剖面给出的氨浓度分层，界面厚度取一个有限元高度；进水口给定均匀水平流速与零氨浓度，出水口给定氨零梯度条件，底部与侧壁下部为不可渗透无滑移壁面，自由面采用刚盖近似。数值上采用非结构三角形网格有限元离散，Poisson方程用标准Galerkin法，涡量输运方程用半隐式后向Euler格式（扩散隐式、非线性对流与浮力显式），氨输运方程采用SUPG稳定化以抑制陡峭浓度梯度附近的非物理振荡，稳定化参数按单元特征尺寸定义；空间离散对流函数采用二次元（）、对涡量与氨浓度采用线性元（），时间推进采用一阶半隐式后向Euler格式结合算子分裂，每步依次求解涡量输运方程、流函数Poisson方程、由流函数计算速度场、求解SUPG稳定化氨输运方程，并施加正性修正以消除数值离散产生的小负浓度。案例研究考察Tank1–Tank5五种实验室尺度构型，由SISdATA原型（约m量级）几何缩放至适配实验水槽的尺寸，涵盖浅而高长深比（Tank1、Tank2）、中等水深（Tank3、Tank4）与深而低长深比（Tank5）三类几何；所有构型长度、单宽流量、运动黏度、氨扩散系数、浮力系数与初始氨浓度相同，进水口开度恒为总水深的20%，故水深增加使进水口高度增加而进水口流速降低；基于进水口的雷诺数与Péclet数在各构型间保持不变。研究还对Tank3开展了网格敏感性分析，对Tank1开展了表观标量扩散系数与浮力系数敏感性分析。
研究结果方面：网格敏感性分析表明，随网格加密氨浓度预测变化逐渐减小，最细两套网格间差异仅约0.14%，最粗网格在后段模拟中出现数值不稳定，说明较细离散下解趋于网格无关。表观标量扩散系数敏感性分析表明，扩散系数取,,时下层浓度诊断分别约在55、48、46 min达到设定阈值，相对参考情形变化约与；两个较低扩散系数值产生非物理的氨量增加，可靠性不足，说明预测依赖于所设标量输运系数，但该系数并未经物理率定。浮力系数敏感性分析表明，对Tank1，浮力系数取与0时阈值时间均为71 min，而取原始系数时为48 min，相对零浮力情形缩短约32.4%，表明弱浓度依赖浮力影响甚微，而预测对假定浮力幅值高度敏感，该结果依赖浮力参数化而非对氨致密度分层的物理确认。基准工况分析以1、10、30 min三个时刻的浓度场快照比较五种构型：Case 1中Tank1在1 min时氨主要滞留于下部、存在清晰分层界面、氨衰减约，10 min时环流发展将氨带向出口、衰减约，30 min时污染显著减少、衰减约，更新接近完成；Case 2中Tank2在1 min时界面基本保持、衰减约，10 min时环流深入、衰减约，30 min时大部分区域浓度降低但底部残留污染区大于Tank1、衰减约；Case 3中Tank3在1 min时存在陡峭过渡带、衰减约，10 min时污染区收缩、衰减约，30 min时显著区域降低但底部仍存污染、衰减约；Case 4中Tank4在1 min时界面清晰、衰减约，10 min时混合增强、衰减约，30 min时大部分区域耗竭、下部中等浓度、衰减约；Case 5中Tank5在1 min时氨仅进入下部、衰减约，10 min时污染区收缩、衰减约，30 min时大部分区域大幅下降、仅底部局地残留、衰减约。平均氨浓度时序分析表明，Tank1因浅水深与最高进水口流速而衰减最快，Tank3、Tank4衰减曲线更平缓，Tank5初始衰减虽慢但后期超越中等水深构型；定量比较显示Tank1约16 min达50%去除、48 min达83%去除，Tank2约40 min与133 min，Tank3约55 min与204 min，Tank4约43 min与208 min（最慢），而更深的Tank5仅需13 min与91 min，显示其后期更新显著更快，说明给定流量下氨更新强烈受几何影响但并非随水深单调变化。水力缩放参数分析设置恒定进水口流速（S1）、固定进水口高度（S2）与固定出水口高度（S3）三种情景：S1下各构型流量随进水口高度改变，Tank5更新时间最短为75 min、Tank1为116 min，而Tank2、Tank3、Tank4分别约需207、206、196 min；S2下Tank1最快（48 min）、Tank2约145 min，而Tank3、Tank4、Tank5在360 min模拟期内均未达到更新标准；S3下Tank1最短（48 min）、Tank2最长（139 min），Tank3与Tank4均约227 min，而更深的Tank5以76 min成为第二快构型；表明更新性能对进、出水条件极为敏感，恒定进水口流速使最深构型最有利，而固定进或出水口开度使最浅构型始终更新最短，中等水深构型对水力缩放最为敏感。定容研究通过随水深增加而缩短槽长保持总水量恒定，长深比由Tank1的110.13降至Tank5的0.46，而进水口高度、出水口高度、流量、进水口流速、输运系数与数值参数均不变：Tank1约48 min达目标浓度，Tank2、Tank3、Tank4分别约需113、128、186 min，Tank5最慢约需353 min，超过Tank1的七倍；衰减曲线间距随长深比减小而增大，表明长深比对氨更新有显著影响，且独立于总水量。
讨论部分指出，尽管基准构型的基于进水口雷诺数与Péclet数相同，更新时间仍差异显著，说明相同的入口无量纲数并不意味着相同的内部流动条件，观测差异应归因于几何、入口动量、流动发展与所设输运参数的共同作用，这与Behroozi与Couturier、Zhang等关于几何与流场结构影响水交换效率的认识一致，但属定性一致而非可直接比较。浓度场随水深的变化表明浅水槽入口流动较快到达污染下层，而深水槽因进水口流速降低而直接穿透减弱；Tank5虽进水口流速最低但在初始期后形成较大尺度环流使氨再分布并输向出口，这种非单调响应表明水深本身不足以预测更新性能，与Zhao等有关水槽比例影响流动分配的认识相呼应。水力缩放分析进一步表明进、出水条件须与几何共同考虑：恒定进水口流速下较强水力驱动可部分补偿水深增加的影响，而固定进水口开度限制了随水深变化的水力强迫改变；因此研究结果不支持单一普适的最优水深，几何与进/出水设计须在明确的水力缩放框架内评估。定容分析将几何比例变化与水量变化分离，表明长深比是影响输运路径与更新过程的重要几何因素，但因同时改变槽长与水深，观测趋势反映二者共同作用而非长深比的独立影响。研究还澄清氨去除的物理含义：模型不含硝化、挥发、吸附或生物消耗，模拟的浓度降低源于水力交换、对流与所设表观标量扩散引起的再分布，而非化学或生物降解；与Jiang等的氨氮动力学模型互补，后者需用于预测实际养殖运行中的氨浓度。与Yu等关于溶解氧分布、Matko等关于氧化沟氧输运的研究一致于流动影响组分分布的一般原理，但溶解氧涉及气液传质与生化消耗，本模型未予表征，故不构成对氨衰减预测的直接验证。敏感性分析揭示了更多的定量解释局限：不同表观标量扩散系数改变预测阈值时间且低值下存在数值可靠性问题，说明所报道更新时间应视为条件性预测而非独立验证的物理去除时间；浮力敏感性表明原始系数所预测的加速更新强烈依赖假定的浮力幅值，初始条件应解释为浓度分层而非氨致密度分层的实验证据；浮力测试仅针对Tank1，其他构型在弱浮力或无浮力下排序是否改变尚待解决。从工程角度看，基准、水力缩放与定容研究提供互补认识：水深增加不必然单调改变更新效率，入口动量与开口配置重要，水槽比例的大幅改变可强烈影响更新过程；本研究的贡献在于采用耦合溶解氨输运模型考察这些几何与水力因素，而非仅依赖流速或颗粒排放指标。但若干局限限制了向全尺度设施的直接应用：2DV公式无法解析三维环流、鱼类诱导混合及真实氨产生与转化的全部过程，模型仍需标量质量守恒、浓度阈值定义、网格与时间步无关性及输运系数的进一步验证，并需速度与浓度实测、溶质质量守恒检验及独立支持的密度-浓度关系。
结论部分译文为：本研究建立了二维垂向（2DV）有限元模型，以考察水槽几何与水力驱动对氨水更新的影响。基准结果表明更新效率取决于环流发展而非仅取决于水深：Tank1在48 min达到设定更新阈值，而Tank4需208 min；尽管Tank5为最深构型，却以91 min达到同一阈值，表明水深与更新性能之间呈非单调关系。水力缩放与几何比例也强烈影响预测更新时间：恒定进水口流速下Tank5更新时间最短，为75 min；定容分析中长深比由110.13降至0.46使更新时间由48 min增至353 min。这些结果表明，评估水更新时须将进水口条件、水深与几何比例共同考虑。敏感性分析进一步发现预测更新行为依赖于所设表观标量扩散系数与浮力系数，特别是去除浓度诱导浮力后Tank1的浓度阈值时间由48 min增至71 min；因此定量预测仍以所选模型参数为条件，需进一步验证与确认。所提出的数值框架为在实体建造或改造前评估替代水槽几何与水力运行条件提供了有用基础。结果强调水深不应脱离进、出水口设计单独选取，因为二者均影响环流与富氨水的输运。未来研究应将模型扩展至三维水槽构型，并利用实验室与现场测量验证预测流场与氨浓度历程；进一步工作还应建立物理上合理的输运与密度参数，并在将模型用于定量养殖设计建议之前纳入相关生物与水质过程。
