Algorithm

392 主题 2.0k 帖子
  • High resolution scheme离散对流项

    14 帖子
    8k 浏览

    @李东岳 您好,抱歉,很久没来了。当时在学习FVM,最近在尝试利用openfoam解决NS方程。

  • 这个方程怎么求?

    1 帖子
    2k 浏览

    \begin{equation}
    m_0=w_1+w_2
    \end{equation}
    \begin{equation}
    m_1=z_1w_1\xi_1+z_2w_2\xi_2
    \end{equation}
    \begin{equation}
    m_2=z_1^4w_1\xi_1^2+z_2^4w_2\xi_2^2
    \end{equation}
    \begin{equation}
    m_3=z_1^9w_1\xi_1^3+z_2^9w_2\xi_2^3
    \end{equation}
    \begin{equation}
    m_4=z_1^{16}w_1\xi_1^4+z_2^{16}w_2\xi_2^4
    \end{equation}
    \begin{equation}
    m_5=z_1^{25}w_1\xi_1^5+z_2^{25}w_2\xi_2^5
    \end{equation}
    左边的$m$都已知,右边的6个未知数$z_1$, $z_2$, $w_1$, $w_2$, $\xi_1$, $\xi_2$需要求。

    除了Newton-Raphson迭代方法?

  • 论坛有做有限元的么?

    9 帖子
    5k 浏览

    毕超的《计算流体力学有限元方法及其编程详解》这本书里面详细介绍了方程各项的离散,不过使用MATLAB编写的,我现在主要就是参考这本书的

  • 絮体的聚并和破碎

    1 帖子
    1k 浏览

    目前絮体的聚并和破碎用的什么模型预测的?
    是否有用到PBM?

    http://lib.dr.iastate.edu/cbe_pubs/118/
    这篇文章用PBM模拟聚并行为,不知道和絮体这个是不是很相关。

  • 狄拉克δ函数的积分特性函数

    3 帖子
    3k 浏览

    @李东岳
    两种方式,简单的是走极限,Gauss分布函数当方差趋于0时可以逼近Dirac函数,所以Gauss分布和函数的积分在方差趋于0时的极限就是Dirac的积分特性函数。

    复杂的就是泛函分析,丫们直接像有限元的加权余量法一样定义函数弱导数有满足某种加权积分关系的性质,反正就是定义都改了,想怎么玩怎么玩。

  • 新版本CFD百科,欢迎测试

    1 帖子
    2k 浏览

    IP访问

    注册用户可以自行创建、修订词条,试试看?谢谢

  • 5 帖子
    4k 浏览

    欧拉拉格朗日模型中离散型用ODE方程描述,追踪每个粒子,这和PBM模型是冲突的,PBM是数量密度函数的模型。

    目前还没见到过文献讨论这个 :cheeky:

  • 最近固态硬盘空间不够

    1 帖子
    1k 浏览

    把ubuntu装在了机械硬盘。

    用了大约4年的固态硬盘,忽然发现难以忍受机械硬盘。用了1星期之后,又装在了固态硬盘。只有40G大小。

    现在固态硬盘便宜了,250G的640元。

  • 水流的空化问题

    1 帖子
    2k 浏览

    实验室有个试验台,里面充满水,水流速度最大60m/s,无叶轮,最近要做单相流实验,想确定增压泵打多少压力水才不会空化。按照我的想法,100℃水温时,水的汽化压力为0.1MPa,要保证不空化,需要线速度最大的地方局部压力达到0.1MPa,按照伯努利方程,需要加上动压,这个时候发现系统压力居然要达到2.8MPa。仔细想想应该是不适用伯努利方程的,所以想问一下,有没有什么好的方法可以估算系统压力呢?

  • 1 帖子
    2k 浏览
  • 求助!

    1 帖子
    1k 浏览

    最近在看Turbulence Modeling for CFD,里面有个CD,谁有CD的内容呀!

  • 5 帖子
    3k 浏览

    @李东岳 好的,谢谢!

  • 4 帖子
    3k 浏览

    你的理解是正确的

  • 多重参考系的疑问

    3 帖子
    3k 浏览

    首先感谢李老师的热情帮助~~
    这个PPT中的第9页给出公式(1.5.9)。如果这个公式是正确的,则我的疑问就没有了,因为我原图中(4)式的右端项,第二项与第三项抵消了。但可疑的,这个公式(1.5.9)似乎并不正确。贴上来再麻烦李老师帮我看看:
    0_1502119393399_7fd74579-5b34-4c15-ac3c-1cc3aae56236-image.png

  • 一阶迎风格式数值扩散的原因

    6 帖子
    5k 浏览

    0_1501229627502_figure1.PNG

  • 滑移边界条件和对称边界条件有什么区别?

    9 帖子
    8k 浏览

    @Qing 谢谢啦:laughing:

  • 求 Von Neumann Stability Analysis 相关的教材?

    1 帖子
    1k 浏览

    Von Neumann Stability Analysis 在普通的CFD教材里面讲的太简单了
    但是到复杂非线性问题(比如Euler方程)中分析时间推进格式的阻尼特性 (damping characteritics) 的时候就比较复杂
    而且一般期刊文章中只给应用结果而且缺乏细节
    它背后的逻辑是什么,我一直没有搞清楚。。。
    希望有大神能推荐一些详细一点的教材,能把思路讲的比较透彻的

  • 并行计算入门求解

    6 帖子
    4k 浏览

    @Kingkong 搜如下教材

    高性能计算并行编程技术_MPI并行程序设计_[都志辉]:非常适合入门
    Using MPI :结合上一个教材一起看
    Iterative Methods for Sparse Linear Systems [ Yousef Saad ] :介绍 deomain decomposition 的数学理论

  • 2 帖子
    2k 浏览

    @yuan_neu 建议看看《张量分析》的相关课程
    二阶 【各向同性】的 【张量函数】 ( Dij(r) )
    只可能是上述 (6.25) 的表达形式

    所以,你首先要明白“张量”,其次再搞清楚“张量函数”的定义,然后需要理解“各向同性”的概念,最后才能明白 6.25 式是怎么推导出来的

    实际上,NS 方程的 Newton 本构关系,也是 “可以” 基于 【线性】和【各向同性】的假设直接写出来的 !

    推荐看看清华黄克智的《张量分析》教材

  • 有关压力速度耦合算法的文章

    2 帖子
    2k 浏览

    @李东岳 coupled solver比segregated solver能更快收敛,但在处理复杂几何形状物体时不太容易收敛[1]。

    [1] http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2012/KlasJareteg/KlasJareteg-OF2012-Project-3-Version-2.0.pdf

  • LES和RANS在代码中实现过程是怎样的?

    15 帖子
    13k 浏览

    @浪迹天大 是这样的,我最近也在读OpenFOAM中LES模型,发现Lund提到
    The nearly universal approach is to simply write down the filtered Navier-Stokes equations together with an assumed model for the subgrid-scale stresses and then apply the desired spatial discretization to this “filtered” system. Although it is rarely mentioned, what one is doing by adopting this procedure is to imagine that the finite support of the computational mesh together with the low-pass characteristics of the discrete differentiating operators act as an effective filter. One then directly associates the computed velocity field with the filtered velocity. This procedure will be referred to as implicit filtering since an explicit filtering operation never appears in the solution procedure.
    也就是说OpenFOAM程序除了动态模型外都没有对N-S方程进行显式的滤波操作,那么我们在方程植入中只要修改$\nu_{sgs}$就OK了,我在后面会做几个滤波函数的选择的对比,看看最后结果一样不一样

  • 2 帖子
    2k 浏览

    不是很清楚你说的那种涡脱落现象。不过用户指南中我记得有一章曾经表示,大体意思为:“如果用稳态求解器模拟比较强的非稳态问题会出现残差的震荡”。

    能否附稳态求解涡街的文献?谢谢。

  • 1 帖子
    1k 浏览

    无意间发现这个帖子,里面提到了normal turbulent model,
    https://www.cfd-online.com/Forums/openfoam/125473-low-reynolds-turbulence-models.html

  • 边界条件引起的震荡

    2 帖子
    2k 浏览

    0_1497945059409_2017-06-20 15-50-42屏幕截图.png

    震荡已经非常非常小了,

  • 关于气液两相流边界层y+的确定

    2 帖子
    2k 浏览

    能否提供一个关于气液两相流关于y+的文献?可能研究侧重点不同,目前我关注的领域对y+要求不是很严格。

    一般情况下,多相流动湍流要依据连续相而定,因此我的回答是应该根据液体来算。在某些极端情况下,气体的湍流场可以认为和液相湍流场相同。

  • 请教雷诺方程的问题?

    3 帖子
    2k 浏览

    @赵一铭 润滑理论的雷诺方程就是NS方程简化了,不考虑厚度方向的压力变化。因为在润滑中动静环之间的液膜很薄,可以认为在厚度方向没有压力变化,是雷诺提出的,就叫雷诺方程了。不知道和那个雷诺是不是一个雷诺。

  • 涡致振动的问题

    5 帖子
    3k 浏览

    @赵一铭 目前想通过现有商业软件做一做···openfoam有点难,项目要求又很急····

  • Python

    11 帖子
    10k 浏览

    @arctan 非常感谢!

  • 请教几个关于壁面函数的问题,谢谢

    8 帖子
    8k 浏览

    @cfd-china link text!![alt text]0_1492654992588_QQ截图20170420102158.png 这是这个文献上说的在固壁上的条件;还有一篇文献直接在壁面设置为0了link text0_1492655090308_QQ截图20170420101901.png

  • SIMPLE 算法中密度如何求解?

    7 帖子
    5k 浏览

    @李东岳已经非常感谢了,我再继续努努力吧!:expressionless:

  • 气体密度如何更新?

    3 帖子
    2k 浏览

    @赵一铭 整个程序我用的SIMPLE算法,所以连续性方程改成了压力修正方程来求解压力,不能直接求解出密度,怎么办?

  • 6 帖子
    4k 浏览

    谢谢,我刚把这个导入到程序中,但结果仍然不理想,解出得P'偏大,在修正速度时会导致速度变成负的,而且是进口速度得负几十倍,所以总感觉压力方程和速度方程不匹配,肿么破??:crying:

  • 关于FVM的几何修正

    3 帖子
    3k 浏览

    不是很了解这个技术。

  • 3 帖子
    2k 浏览

    @赵一铭 不是用软件,是编程计算,现在的问题是一个不收敛,和出现负值

  • 11 帖子
    7k 浏览

    3个月前沟通过,基金会无心出理论指南,除非有重大的赞助。基金会的风格和Henry个人风格很有关系,毕竟Henry是老大,并且Henry和学术界无任何关系。考虑一些商业事情,要从人的本性出发。

    国内CFDer如果想要官方版TG,需要募集大量资金赞助基金会。

  • 气液固多相流

    2 帖子
    2k 浏览

    没看过三相的,两相的很多。又分为欧拉欧拉,欧拉拉格朗日。Ishii2010年出了一个新书,名称不记得了。这个人和Drew是欧拉欧拉模型的先驱。

  • 一阶迎风的耗散性

    2 帖子
    2k 浏览

    0_1489821091713_vanLeer.png

    上图是细画网格后采用TVD格式计算的,型线保持的可以了。

  • 焓和温度的关系

    2 帖子
    2k 浏览

    Cp不是常数的情况下,焓和温度不是线性的关系吧。还是说你要找某个特定流体的焓随温度变化的关系?

  • 如何理解PDE的弱解?

    2 帖子
    4k 浏览

    @random_ran

    不满足微分方程,但是满足积分方程的解就是弱解吧。本来俩方程应该等价的。

    其实我觉得关键单不在于什么是弱解,而是这帮数学家为啥要引入弱解这个概念,我的理解是丫们发现物理上一些有意义的解不严格满足他们推导出来的微分方程,但是似乎又差得不多,然后整吧整吧搞出一个弱解的概念来。

    有点感觉想和稀泥。。。。

  • 8 帖子
    5k 浏览

    @random_ran :upset: 老师让做这方面的,没办法啊

  • 最近看流体力学数学导论中到的一段推导

    6 帖子
    4k 浏览

    @李东岳 哈哈,其实是从nodebb的帮助中抄来的:cheeky:

  • 求矩阵的逆一定要经过阶数的3次幂操作?

    3 帖子
    2k 浏览

    @李东岳

    常见的问题是求$\mathbf A \cdot \mathbf x=\mathbf b$中的$\mathbf x$,可以有直接解法(如LU),也可以用迭代法(如Gauss-Seidal),求逆矩阵肯定是下下之选,原因有几点: 操作太多,$~O(N^3)$; 稀疏矩阵的逆矩阵不是稀疏的; 你要求得$\mathbf x$还要再乘一次,$~O(N^2)$ 逆矩阵求了之后没有用,因为CFD方程组是非线性的,每次求的矩阵$\mathbf A$都不一样,所以你求得的逆矩阵只用一次,下次迭代又得重新求解。如果是针对线性问题,每次针对不同的$\mathbf b$求对应的$\mathbf x$,用逆矩阵还勉强说得过去。 现代的CFD程序似乎迭代法用得更多,JFNK, IDR(s),最不济还有CG和BiCGSTAB.有的连矩阵都不需要组装,何苦去求逆...
  • 请教各种系数的定义及其作用

    5 帖子
    5k 浏览

    @random_ran 谢谢,很有帮助,祝好

  • 如果存在一个波,如何获取波的周期?

    3 帖子
    2k 浏览

    恩,对时域的数据用傅里叶变换到频率域上是工程用得最多的方法。

    知乎上3年前Heinrich写了一篇科普文,觉得特别浅显易懂。

    我在对升力系数做处理求karman vortex频率的时候,就用的Matlab。如果需要的话,我会贴在这。

    另外,PSD图,也就是频率域上的能量分布曲线,也是经常在和周期性相关的问题中出现。

  • 有关高斯积分积分点的确定

    3 帖子
    4k 浏览

    @kasitefy
    非常对。

    另外一个问题是这种加和的数值积分忒慢了

  • 时间离散格式和时间推进格式有何不同?

    2 帖子
    2k 浏览

    时间离散格式应该是计算时间导数的格式,模板没有限制;而推进格式应该要求模板是单侧偏心的。http://www.scholarpedia.org/article/General_linear_methods是一般的时间推进格式,偏心模板最多也就是隐式格式。

  • 2016目录

    2 帖子
    3k 浏览

    李博继续多出好文!

  • 2 帖子
    3k 浏览

    Fluent我们不太清楚开发人员是如何做的数值稳定性工作。不过如果你看OpenFOAM源代码会发现OpenFOAM里面对湍流方程的守恒、有界、减少数值扩散、以及针对网格质量的自适应数值调整是非常非常多的。没有一个方程在openfoam中是直接求解的。So~

  • 关于计算的收敛性

    2 帖子
    2k 浏览

    很普适性的问题。如果你说的是如何判定非稳态的收敛性的话,可以监控一些参数,如果稳定波动一定周期,就可以认为收敛了。

  • 关于求导精度是否会降阶

    3 帖子
    2k 浏览

    用有限差分法得到的二阶精度流场变量和用有限体积法得到的二阶精度流场变量,再分别求涡量场,用两种方法得到的涡量场是否都具有二阶精度或者一阶精度?我不太清楚这其中的原理。你说不会降阶,能否具体说一下原因?