Skip to content

Algorithm

基本理论区

408 主题 2.1k 帖子
  • 这个方程怎么求?

    1
    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
    9 帖子
    8k 浏览
    L

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

  • 絮体的聚并和破碎

    1
    1 帖子
    2k 浏览
    李东岳

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

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

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

    3
    3 帖子
    4k 浏览

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

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

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

    1
    1 帖子
    2k 浏览
    李东岳

    IP访问

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

  • 群平衡模型,研究过的大神求指教呀!!!!

    5
    5 帖子
    6k 浏览
    李东岳

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

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

  • 最近固态硬盘空间不够

    1
    1 帖子
    2k 浏览
    李东岳

    把ubuntu装在了机械硬盘。

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

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

  • 水流的空化问题

    1
    1 帖子
    2k 浏览
    W

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

  • 1 帖子
    2k 浏览
    李东岳

    http://www.sciencedirect.com/science/article/pii/S0021999117306289

    0_1505210410009_捕获.JPG

  • 求助!

    1
    1 帖子
    2k 浏览
    K

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

  • 5 帖子
    5k 浏览
    J

    @李东岳 好的,谢谢!

  • 4 帖子
    4k 浏览
    李东岳

    你的理解是正确的

  • 多重参考系的疑问

    3
    3 帖子
    4k 浏览
    Z

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

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

    6
    6 帖子
    7k 浏览
    H

    0_1501229627502_figure1.PNG

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

    1
    1 帖子
    2k 浏览
    Q

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

  • 并行计算入门求解

    6
    6 帖子
    6k 浏览
    Q

    @Kingkong 搜如下教材

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

  • 2 帖子
    3k 浏览
    Q

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

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

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

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

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

    2
    2 帖子
    3k 浏览
    W

    @李东岳 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
    15 帖子
    18k 浏览

    @浪迹天大 是这样的,我最近也在读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 帖子
    3k 浏览
    李东岳

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

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