@东岳 也是一个比较强大的软件,不过主要在化工方面用的多,像一些化工单元流程,如精馏,吸收,蒸发,萃取什么的
lsprxd
帖子
-
-
这个不是用Aspen比较好吗?
-
有的,可以参考毕超的《计算流体力学有限元方法及其编程详解 》,上面有详细的利用FEM法对NS方程离散的过程,还有MATLAB程序
-
@东岳 是的,我看的文献大多数是关注速度,温度和涡量场变化,几乎没有涉及压力场变化,于是我编程计算的时候也没有考虑,结果老师问我为什么,我都不知道如何回答了
-
是有什么原因吗?
-
看了许多文献,在方腔自然对流中大家都分析速度,温度,怒塞尔数的变化,为什么没有考虑压力的变化?
-
是的,我编程做一个三维后台阶的例子,如果按照h , L, k求得话,那这个h对流换热系数如何确定?
-
@东岳 在资料上看的有平均努塞尔数和局部努塞尔数,在求解的时候有用微分方程表达(但都是无量纲形式),如果是原始变量的方式应该如何计算?
-
也就是说如何根据原始变量方程得出来的解求得努塞尔数
-
我是用MATLAB编程的
-
努赛尔数为无量纲的一个准则数,为壁面法向无量纲过余温度梯度,最后的定义我们可以看到是这个形式,是用无量纲形式表达;
但我在编程计算的时候采用的是原始变量,没有采用无量纲方程,这样的话如何计算努塞尔数? -
@李东岳 谢谢李老师
-
@李东岳 温度和热流密度应该是标量吧?那热流密度为0就是对称边界了
-
对称边界条件是法向梯度为0,热流密度是法向为q,如果q=0,不是和对称一样吗?
-
我想问一下,如果热流密度为0,是不是就变成对称边界条件了?
-
@李东岳 哦,明白,谢谢李老师
-
我要模拟一个简单的二维强迫对流换热,边界条件为速度入口,压力出口0,入口温度为0,上下壁面温度为1,出口温度不知道,在压力出口时该如何设置,出口温度应该是热流密度为0啊!
他这个回流温度是什么? -
毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散,不过使用MATLAB编写的,我现在主要就是参考这本书的
-
时间项的处理基本上都是差分发,然后引入sita因子
-
@李东岳 对于时间项的处理也是利用差分格式比如速度时间项的偏导:(u^n+1-u^n)/delta t
-
如果说三节点单元和四节点单元的话,没有四节点单元的精度高
-
@Qing 谢谢啦:laughing:
-
@Qing 如果壁面是水平的就可以用对称边界条件代替滑移
-
@Qing 像对称边界条件我只要设置法向为0,切向就不用管了;在这个滑移不仅要给出法向为0,设置出切向速度?如何插值得到,在没有求解之前怎么知道壁面附近的流场信息?PS:我看的一篇文献,上面在求湍流的时候设置了三种边界条件,其中一种就是壁面uw设置为free slip,然后我就不知道这个free slip是如何设置的了!!
-
@Qing 你好,不太明白第二句话什么意思:confused:
-
具体比如一根管子,轴处是对称边界条件,ux不知道,uy对y的偏导为0;管壁如果是无滑移就都是uw=0,如果是滑移的话该怎么表达?
-
@李东岳 没事,麻烦李老师了!一个?比如这个图
这些都是p点吧?那不同的p点的kp应该不同吧? -
@李东岳 李老师,您这个简单计算y+的时候带入的kp是p点的值,如何确定的?还有这个p点是距离壁面y处的位置,但距离壁面y处的点应该有好多啊,而且每个应该都不一样吧?那求出来的y+也不止一个吧?
-
@李东岳 谢谢李老师,OpenFOAM这个不熟悉,您贴的这个我在底下先仔细看看
-
在网上看的有人这样说的“最简单的办法是用对第一个节点的K,E直接赋值,由U+,Y+计算ut(摩擦速度),K=ut*ut/sqrt(0.09),E=ut**3/y/0.42”,这个赋值是赋在壁面上?这个p点是距壁面y处的值,那这些赋值的位置应该是距壁面y处的所有值吗?
-
有的文献是这样说的
有点不明白的是“处在对数率层,可用牛顿迭代法求出u*”这个是如何迭代?这个公式里面的u是哪个地方的u?是全场的还是y出的u?得到的k和e也是y处的u吗?
这个壁面函数是影响主流区吗?为什么最后他说壁面切应力对NS方程的影响可用用等效粘性系数表示? -
@random_ran 好多文献就是这样的,给定初值,求解N-S,得到速度带入k-e,更新湍流粘度。就是第二张图片那个连接里面的,而且他直接在wall上设置k和e为0
-
@李东岳 上面的第5步更新wall粘度,这个wall粘度是不是整个区域的wall粘度?
像一般的边界都是自然边界或者本质边界,在这个k-e方程中,k的边界可以是第二类边界,在y方向上的偏导为0,Kp按照k方程计算;
按照书上说的p点的epsilon是按照代数方程计算就是式(28);
那这个在p点的k和e是不是应该是一组数,而不是一个确定的数?
还有一个就是在壁面处的计算和整个区域的计算是不是一个控制方程? -
@李东岳 wall粘度和流场中的湍流粘度是不是一个粘度?
-
@李东岳 李老师我是采用无网格伽辽金法来做的,整个方程的求解策略是这样的
这个壁面函数是如何引入到这个策略中的? -
@李东岳 李老师,在获得全场的k值,是不是要提取壁面节点处的kP,然后带入式(28)确定壁面节点处的εP?
-
实在是不懂了,求解释:crying:
-
-
各位老师,这个壁变函数一直搞得浑浑噩噩的,看文献说“近 壁 区 高 Re数 采 用 的紊 动 各 向同 性 假 设 不 再 适 用,因而k-e紊 流 模式在 此 不 能 应 用,故 应 选 取 距 固 壁某 一 距 离 的 边 界作 为 计 算 边界。”这句有话有点不懂,k和e是取“距 固 壁某 一 距 离的值”作为边界条件吗?那固壁处怎么办,就不处理了吗
-
@李东岳 前辈你好,这个壁变函数一直搞得浑浑噩噩的,看文献说“近 壁 区 高 Re数 采 用 的紊 动 各 向同 性 假 设 不 再 适 用,因而k-e紊 流 模式在 此 不 能 应 用,故 应 选 取 距 固 壁某 一 距 离 的 边 界作 为 计 算 边界。”这句有话有点不懂,k和e是取“距 固 壁某 一 距 离的值”作为边界条件吗?那固壁处怎么办,就不处理了吗
-
@赵一铭 不是用软件,是编程计算,现在的问题是一个不收敛,和出现负值
-
k和e出现了负值
-
源文件链接丢了
-
@random_ran :upset: 老师让做这方面的,没办法啊
-
@random_ran 如果只是求解NS方程没问题,如果在加上k-e模型就困难了:upset:
-
@cfd-china 是有限单元法的,知道CFD好多是有限差分和有限容积法,我是要用无网格伽辽金法数值分析的,因为这个方法和有限单元法类似,就找这方面的内容。湍流粘度是根据k-e算出的,但这个动量方程也和湍流粘度有关的,在这两部分耦合始终有问题
-
@random_ran 是有限单元法的,知道CFD好多是有限差分和有限容积法,我是要用无网格伽辽金法数值分析的,因为这个方法和有限单元法类似,就找这方面的内容。
CFD在反应精馏塔中的应用
CFD在反应精馏塔中的应用
利用有限元方法做计算流体
在方腔自然对流中为什么大家都不分析压力的变化
在方腔自然对流中为什么大家都不分析压力的变化
在方腔自然对流中为什么大家都不分析压力的变化
关于努塞尔数问题请教
关于努塞尔数问题请教
关于努塞尔数问题请教
关于努塞尔数问题请教
关于努塞尔数问题请教
关于温度对称边界条件
关于温度对称边界条件
关于温度对称边界条件
关于温度对称边界条件
关于压力出口温度设置
关于压力出口温度设置
论坛有做有限元的么?
论坛有做有限元的么?
论坛有做有限元的么?
论坛有做有限元的么?
滑移边界条件和对称边界条件有什么区别?
滑移边界条件和对称边界条件有什么区别?
滑移边界条件和对称边界条件有什么区别?
滑移边界条件和对称边界条件有什么区别?
滑移边界条件和对称边界条件有什么区别?
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
关于壁面函数有些不懂
请教几个关于壁面函数的问题,谢谢
请教几个关于壁面函数的问题,谢谢
湍流KE方程数值模拟k和e出现负值应该如何处理?
湍流KE方程数值模拟k和e出现负值应该如何处理?
OpenFOAM与fluent的比较
在用有限单元法求解k-e模型时,湍流粘度怎么处理
在用有限单元法求解k-e模型时,湍流粘度怎么处理
在用有限单元法求解k-e模型时,湍流粘度怎么处理
在用有限单元法求解k-e模型时,湍流粘度怎么处理