code-saturne / code-saturne/code_saturne

Bug: Symmetric tensor diffusivity

未关闭
#64 3 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看

还没有人认领这个 Issue。

主要语言
C++
星标
307
派生
102
PR 合并指标
30 天内没有已合并 PR

描述

The combination of turb_model->iturb = CS_TURB_NONE; and iturt(iscalt) = 10 (i.e. GGDH) leads to a segfault. I would expect the code to ignore iturt for laminar flow so the user can simply change one parameter when switching on the turbulence model. IMO, we should at least throw a warning that the settings are incompatible, rather than segfault.

It seems the issue stems from cpltur.f90 where a test for GGDH leads to an assumption of a symmetric tensor diffusivity for the thermal variable, when in reality the diffusivity is a scalar (my settings have diffusivity as variable in space and time, if it makes a difference).

贡献指南

打开贡献指南

从这里开始

  1. 先读完整个 Issue,再读项目的贡献指南。
  2. 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
  3. Fork 仓库,在一个分支上完成修改。
  4. 提交 Pull Request,并在描述里引用这个 Issue 编号。

调研方向

从 cpltur.f90 开始,使用 CS_TURB_NONE 和 iturt(iscalt)=10 重现 segfault,同时包括空间和时间变化扩散率的设置。跟踪热变量的 GGDH 分支,并将其对对称张量的假设与层流标量情况进行比较。当不兼容的设置不再导致 segfault,并且能够被安全忽略或通过警告报告时,即表示完成。

由索引模型根据 Issue 内容生成。

评估

技术栈
fortran
领域
hpc
Issue 类型
缺陷
难度
4/5
预计耗时
3-5 天
活跃度
停滞
描述清晰度
基本清楚
新手友好度
42/100

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。