24小时热门版块排行榜    

查看: 342  |  回复: 1

狼牙_nuaa

银虫 (小有名气)

[求助] les 的 channel flow validation 中添加体积力的相关问题

各位大神好,我在进行LES程序模块的验证,遇到点问题,非常头疼。
个人代码能力水水的,实在不知道这问题咋处理了,特在此把问题贴出来,望得到各位虫友的帮助,谢谢了
文献的公式和代码都在附件
特注:我的代码是有量纲的
代码的主题贴在这里 ,全部代码在附件:

alpha_f =   2._precision/dtstep
beta_f = -0.2_precision/dtstep
bfr = bfr + (alpha_f*(qfp - qf0) + beta_f*(qf  - qf0))*dtstep/Ln/Lp    !bfr表示体积力,这一步是体积力的更新,
                                                                                                         !qf表示质量通量,qfp是一阶预估步,qf0是预设的预估步,我直接取为第一步的值。
                                                                                                          ! Ln 和Lp表示物体在展向和法相的体长
一阶预估步为:
qfp = qf - dtstep*(Ln*Lp*bfr + 2.0_precision * Lp*vis_wall(1)*dudx(1))  !vis_wall 和  dudx(1)表示壁面处的剪切速度

les 的 channel flow validation 中添加体积力的相关问题
可压缩槽道流.png
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : channel.F90
  • 2015-11-19 16:46:35, 9.46 K

» 猜你喜欢

哈哈哈
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

狼牙_nuaa

银虫 (小有名气)

这么多人,就给我个助贴也行哇
哈哈哈
2楼2015-11-19 17:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 狼牙_nuaa 的主题更新
信息提示
请填处理意见