位置: IT常识 - 正文

【OpenFOAM】-olaFlow-算例1- baseWaveFlume(openfoam中文教程)

编辑:rootadmin
【OpenFOAM】-olaFlow-算例1- baseWaveFlume

推荐整理分享【OpenFOAM】-olaFlow-算例1- baseWaveFlume(openfoam中文教程),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:openfoam中文教程,openfoam lbm,open foam,openfoam sample,openfoam sph,openfoam sample,openfoam-a little user-manual,openfoam reconstruct,内容如对您有帮助,希望把文章链接给更多的朋友!

算例路径: olaFlow\tutorials\baseWaveFlume 算例描述: 一个基础的二维波浪水槽 学习目标: olaFlow 求解器的造波、消波方法和算例设置 算例快照:

图1 波浪模拟结果

图2 算例网格

文件结构:

├── 0.org│ ├── U │ ├── alpha.water│ ├── alpha.water.org│ └── p_rgh├── cleanCase├── constant│ ├── dynamicMeshDict│ ├── g│ ├── transportProperties│ ├── turbulenceProperties│ └── waveDict --> 设置波浪要素├── runCase└── system ├── blockMeshDict ├── controlDict ├── decomposeParDict ├── fvSchemes ├── fvSolution └── setFieldsDict --> 设置水深

算例文件解析:

【0.org\U】

dimensions [0 1 -1 0 0 0 0]; // 量纲 m/sinternalField uniform (0 0 0); // 内部速度场 均一场 0 0 0boundaryField //边界场{ inlet // 造波边界 { type waveVelocity; // 波浪速度 waveDictName waveDict; // 读取constant\waveDict中的波浪要素 value uniform (0 0 0); //初值为 0 0 0 } outlet // 消波边界 { type waveAbsorption2DVelocity; // 使用了二维消波理论,olaFlow采用主动消波法 value uniform (0 0 0); } bottom // 底部边界为固壁边界,边界上速度为零 { type fixedValue; value uniform (0 0 0); } atmosphere // 大气边界,允许空气流出和流入 { type pressureInletOutletVelocity; value uniform (0 0 0); } frontAndBack // 前后面,empty指示模型为二维模型 { type empty; }}

【0.org\p_rgh】

// p_rgh = p - rgh,实际压力减去静水压力dimensions [1 -1 -2 0 0 0 0]; // M(1) L(-1) T(-2)internalField uniform 0;boundaryField{ frontAndBack { type empty; } outlet { type fixedFluxPressure; //将压力梯度设置为0,边界上的通量由速度边界条件指定 value uniform 0; } inlet { type fixedFluxPressure; value uniform 0; } bottom { type fixedFluxPressure; value uniform 0; } atmosphere { type totalPressure; //总压条件:流出 p = p0; 流入 p = p0 - 0.5|U|^2 U U; phi phi; rho rho; psi none; gamma 1; p0 uniform 0; value uniform 0; }}【OpenFOAM】-olaFlow-算例1- baseWaveFlume(openfoam中文教程)

【0.org\alpha.water.org】

// 设置流体体积分数dimensions [0 0 0 0 0 0 0]; internalField uniform 0; boundaryField{ inlet { type waveAlpha; // 根据波浪条件设置 waveDictName waveDict; value uniform 0; } frontAndBack { type empty; } outlet { type zeroGradient; // 零梯度 } bottom { type zeroGradient; } atmosphere { type inletOutlet; // 当流体流出时,α的梯度为0;如果流入,α为0,100%的空气流入 inletValue uniform 0; value uniform 0; }}

【constant\dynamicMeshDict】

dynamicFvMesh staticFvMesh;

【constant\g】

dimensions [0 1 -2 0 0 0 0];value ( 0 0 -9.81 );

【constant\transportProperties】

phases (water air);water{ transportModel Newtonian; nu [0 2 -1 0 0 0 0] 1e-06; // 流体运动粘度 rho [1 -3 0 0 0 0 0] 1000; // 流体密度}air{ transportModel Newtonian; nu [0 2 -1 0 0 0 0] 1.48e-05; rho [1 -3 0 0 0 0 0] 1;}sigma [1 0 -2 0 0 0 0] 0.07; // 水和空气之间的表面张力参数

【constant\turbulenceProperties】

simulationType laminar; // 设置为层流模型

【constant\waveDict】

waveType regular; // 规则波waveTheory cnoidal; // 椭圆余弦波genAbs 1; // 考虑造波边界的消波性能 1/0absDir 0.0; // 造波边界的消波方向nPaddles 1; // 主动消波的Paddles数量设置waveHeight 0.10; // 波高wavePeriod 3; // 波周期waveDir 0.0; // 波向wavePhase 1.57079633; // 初始相位// Change both entries to true to re-read this dictionary upon restart.rereadAlpha false;rereadU false;

【system\blockMeshDict】

scale 1;vertices ( (0.0 -0.02 0.0) (10.0 -0.02 0.0) (10.0 -0.02 0.7) (0.0 -0.02 0.7) (0.0 0.0 0.0) (10.0 0.0 0.0) (10.0 0.0 0.7) (0.0 0.0 0.7));blocks ( hex (0 1 5 4 3 2 6 7) (500 1 70) simpleGrading (1 1 1));edges ();patches ( patch inlet // 造波边界 ( (0 4 7 3) ) patch outlet // 消波边界 ( (1 5 6 2) ) wall bottom // 水槽底部边界 ( (0 1 5 4) ) patch atmosphere // 大气边界 ( (3 2 6 7) ) empty frontAndBack // 水槽侧面边界 ( (0 1 2 3) (4 5 6 7) ));mergePatchPairs();

【system\controlDict】

application olaFlow; // olaFlow求解器startFrom latestTime;startTime 0;stopAt endTime;endTime 60;deltaT 0.001; // 计算时间步writeControl adjustableRunTime;writeInterval 0.05; // 写出时间步purgeWrite 0;writeFormat ascii;writePrecision 6;compression off; // 是否压缩格式写出,可节约硬盘空间, on/offtimeFormat general;timePrecision 6;runTimeModifiable yes;adjustTimeStep yes; // 采用自适应时间步,可能会加速计算,也可能造成时间步极小maxCo 0.5; // CFL条件的Courant数, 一般<1, 设置一个小值会使计算结果更精确,但也减小了时间步长,增加了计算成本maxAlphaCo 0.5; // 两相交界面上的最大Courant数maxDeltaT 0.025;

【system\decomposeParDict】

numberOfSubdomains 2; // 并行区域分解数目method scotch; // 区域分解方法...

【system\fvSchemes】

// 指定控制方程中各项的有限体积法的离散格式ddtSchemes // 指定时间离散格式{ default Euler; // Euler法,一阶精度,条件稳定} gradSchemes // 梯度项离散格式{ default Gauss linear; // 高斯定理,将网格中心的量插值到网格面上}// olaFlow 的算例中给出了几乎所有可能出现的项,具体算例可能不会包含全部项divSchemes // 对流项与散度项的离散格式, 将网格中心的量插值到网格面上,因此实际上选用的是interpolationSchemes{ div(rhoPhi,U) Gauss limitedLinearV 1; // Guass limitedLinear(一种TVD格式,使同时满足精度和有界) V类(采用限制器时考虑了流动方向) div(U) Gauss linear; // 二阶精度,无界 div((rhoPhi|interpolate(porosity)),U) Gauss limitedLinearV 1; div(rhoPhiPor,UPor) Gauss limitedLinearV 1; div(rhoPhi,UPor) Gauss limitedLinearV 1; div(rhoPhiPor,U) Gauss limitedLinearV 1; div(phi,alpha) Gauss vanLeer; // Gauss vanLeer(一种TVD格式,使同时满足精度和有界) div(phirb,alpha) Gauss interfaceCompression; // 界面压缩格式,基于一般限制格式 div((muEff*dev(T(grad(U))))) Gauss linear; div(phi,k) Gauss upwind; // 一阶迎风格式,有界 div(phi,epsilon) Gauss upwind; div((phi|interpolate(porosity)),k) Gauss upwind; div((phi|interpolate(porosity)),epsilon) Gauss upwind; div(phi,omega) Gauss upwind; div((phi|interpolate(porosity)),omega) Gauss upwind;}laplacianSchemes // 拉普拉斯项离散格式{ default Gauss linear corrected; // Guass线性插值,corrected(显式的非正交网格修正)}interpolationSchemes{ default linear; // 线性插值格式}snGradSchemes // 面法向梯度格式{ default corrected;}fluxRequired // { default no; p_rgh; pcorr; alpha.water;}

【system\fvSolution】

// 指定方程组矩阵求解器、残差以及其他算法控制solvers{ "alpha.water.*" { nAlphaCorr 1; // nAlphaSubCycles 2; alphaOuterCorrectors yes; cAlpha 1; MULESCorr no; nLimiterIter 3; solver smoothSolver; // 求解器:光顺求解器。对称和非对称矩阵均适用 smoother symGaussSeidel; // 光顺器:对称Gauss-Seidel方法 tolerance 1e-8; // 标准化残差 relTol 0; // 相对残差,表征残差的递减量 } "pcorr.*" // 压强校正量 { solver PCG; // 求解器:预处理共轭梯度法。仅适用于对称矩阵 preconditioner DIC; // 预处理器:基于对角的不完全Cholesky预处理器 tolerance 1e-5; relTol 0; } p_rgh { solver PCG; preconditioner DIC; tolerance 1e-07; relTol 0.05; } p_rghFinal // 在求解p_rgh时可设置较大的残差,在最后一步设置严格的误差 { $p_rgh; relTol 0; } U { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-06; relTol 0; }}// 离散方程组的分离式解法PIMPLE // SIMPLE 与 PISO 算法的耦合{ momentumPredictor no; nOuterCorrectors 1; // 值不大于1时,算法变为纯PISO算法;若大于1,为非定常SIMPLE算法 nCorrectors 3; nNonOrthogonalCorrectors 0;}relaxationFactors // 松弛方法,为降低数值波动{ fields // 表示需要使用松弛方法(显式)的变量 { } equations // 表示需要使用松弛方法(隐式)的方程 { ".*" 1; }}

【system\setFieldsDict】

defaultFieldValues( volScalarFieldValue alpha.water 0 // 初始化流体体积分数);regions( boxToCell { box (-10 -1 -1) (30 1 0.4); fieldValues ( volScalarFieldValue alpha.water 1 // 将box范围内的流体设置为a = 1的相,即水深 ); });
本文链接地址:https://www.jiuchutong.com/zhishi/300854.html 转载请保留说明!

上一篇:基于OC端的Bridge-API组件化应用(oc底层原理)

下一篇:vue导入处理Excel表格详解(导入vue.js)

  • 苹果se3是双卡双待吗(苹果se3是双卡双待吗..2022)

    苹果se3是双卡双待吗(苹果se3是双卡双待吗..2022)

  • 淘宝直播不能上下滑动(淘宝直播不能上下切换)

    淘宝直播不能上下滑动(淘宝直播不能上下切换)

  • 三星2016+和三星2017有什么区别(三星w2016对比2017)

    三星2016+和三星2017有什么区别(三星w2016对比2017)

  • 蓝牙power off什么意思(蓝牙说power off是什么意思)

    蓝牙power off什么意思(蓝牙说power off是什么意思)

  • 爱奇艺vip可以同时几个人用(爱奇艺vip可以两人共用吗)

    爱奇艺vip可以同时几个人用(爱奇艺vip可以两人共用吗)

  • 勾选框怎么打出来(勾选方框怎么打)

    勾选框怎么打出来(勾选方框怎么打)

  • 华为点开图片退不出去(华为手机图片点开后不能返回)

    华为点开图片退不出去(华为手机图片点开后不能返回)

  • 苹果11没有挂电话键(苹果11没有挂电话功能吗)

    苹果11没有挂电话键(苹果11没有挂电话功能吗)

  • 开飞行模式能收到短信吗(开飞行模式能收到验证码吗)

    开飞行模式能收到短信吗(开飞行模式能收到验证码吗)

  • 哪款ipad能插电话卡(ipad支持充电)

    哪款ipad能插电话卡(ipad支持充电)

  • 密保手机和绑定手机有什么区别(密保手机和绑定手机区别)

    密保手机和绑定手机有什么区别(密保手机和绑定手机区别)

  • 华为m6怎么设置分辨率(华为M6怎么设置手写笔防误触)

    华为m6怎么设置分辨率(华为M6怎么设置手写笔防误触)

  • 苹果11下面的横线怎么去掉(苹果11下面的横线怎么没震动了?)

    苹果11下面的横线怎么去掉(苹果11下面的横线怎么没震动了?)

  • 唯品会怎么申请退货(唯品会怎么申请差价)

    唯品会怎么申请退货(唯品会怎么申请差价)

  • 设备信息权限是什么意思(设备信息权限怎么开启)

    设备信息权限是什么意思(设备信息权限怎么开启)

  • nova4发热怎么解决(nova4手机发热是怎么回事)

    nova4发热怎么解决(nova4手机发热是怎么回事)

  • 小米9多少倍变焦(小米九几倍光学变焦)

    小米9多少倍变焦(小米九几倍光学变焦)

  • 按键精灵怎么设置20分钟按一个键(按键精灵怎么设置循环)

    按键精灵怎么设置20分钟按一个键(按键精灵怎么设置循环)

  • oppor15x怎么开启快充(oppor15x怎么开启游戏模式)

    oppor15x怎么开启快充(oppor15x怎么开启游戏模式)

  • 怎么让港版x支持电信(港版怎么用支付宝支付)

    怎么让港版x支持电信(港版怎么用支付宝支付)

  • Windows7系统更改家庭组的密码的方法(win7系统修改)

    Windows7系统更改家庭组的密码的方法(win7系统修改)

  • element table数据量太大导致网页卡死崩溃(elements table)

    element table数据量太大导致网页卡死崩溃(elements table)

  • python怎么获取键盘监听(python获取值对应的键)

    python怎么获取键盘监听(python获取值对应的键)

  • 爱上源码,重学Spring IoC深入

    爱上源码,重学Spring IoC深入

  • 未按规定缴纳保险
  • 个人所得税年度汇算时间
  • 税款的汇兑损益如何做账
  • 建筑简易征收需要成本发票吗
  • 股东参与经营可以分红吗
  • 国税地税合并后叫什么名称
  • 收到赔偿款怎么做分录
  • 非现金资产债务怎么算
  • 城市维护建设税怎么做分录
  • 资本公积金转增股本是利好吗
  • 买什么样的房子可以贷款
  • 多交附加税怎么做会计分录
  • 当期允许抵扣的农产品进项税额
  • 超出发票使用范围
  • 仓储行业税率是多少
  • 工会没有税号怎么开普票
  • 公司当年福利费没有用到限额,可以按照14%计提税前扣除吗?
  • 非直接结算是什么意思
  • 电子税务局里的利润表,本月金额是填累计数吗
  • 国际货运代理免税怎么做账
  • 建筑安装预缴增值税
  • 外企投资应该怎么投资
  • 关联公司代付款协议
  • 短期借款利息的处理
  • 会员充值增送如何取消
  • 调整以前年度应付工资怎么调
  • win10 bios设置
  • 会计记账凭证摘要模板
  • 高新企业开办费是多少
  • 补交以前年度的城建税会计分录
  • php模块
  • 劳务费要申报个税吗怎么申报
  • 出租房屋收取的水电费的税率
  • IE浏览器如何取消全屏
  • 增值税按简易计价方式
  • 解聘员工补偿金能拖欠吗
  • PHP:JulianToJD()的用法_日历函数
  • 单臂路由实验目的
  • 职工意外保险比例
  • Spring Security(七)-- AuthenticationProvider
  • vue2路由跳转页面不刷新问题
  • web自动化selenium实战项目
  • 劳务派遣开具发票
  • 进项税额销项税额
  • phpcms文档
  • 库存亏损要缴纳增值税吗
  • 什么是进项税转出成本
  • vue使用百度地图md5报错
  • 织梦使用手册
  • 固定资产可以一次性折旧吗
  • 建筑劳务公司的税率是多少
  • 股东权益是资产吗
  • 地下车库的成本包括哪些
  • 固定资产未提完折旧
  • 实收资本挂在其他应收款现在股权转让
  • 固定资产清理账户借方的核算内容包括
  • win10文件夹设置
  • win10预览版选哪个
  • xp怎么安装framework
  • mac系统操作的小技巧
  • Mac怎么强制关闭窗口
  • centos下载命令
  • windows任务管理器怎么用
  • ubuntu sudo not found
  • 进程 cmd
  • 使用windows防火墙禁止软件联网
  • ubuntu18.04没有声音
  • win10各个版本的桌面
  • linux常用桌面环境
  • python字典常用操作以及字典的嵌套
  • perl-v
  • shell脚本的fi
  • dos命令怎么写
  • cmd开机启动命令
  • 有关的拼音
  • 技术咨询服务费税率
  • 个体税务怎么查询
  • 合并申报后老申报表还可以继续使用吗
  • 西安税务总局
  • 地税局和税务局有什么区别
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设