Skip to content

Algorithm

基本理论区

416 主题 2.1k 帖子
  • 将多相流模型完全双曲化都有什么好处呢?

    1
    1 帖子
    1k 浏览
    yhdthuY

    如题。

  • 一个朋友的答疑,关于动量交换

    1
    1 帖子
    1k 浏览
    李东岳

    考虑最简单的动量交换,液相:
    \begin{equation}
    \frac{\rd \bfU_\rc}{\rd t}=-\bfA
    \end{equation}
    气相:
    \begin{equation}
    \frac{\rd \bfU_\rd}{\rd t}=\bfA
    \end{equation}
    乍一看,如果$\bfA$是负的,那么会导致$\bfU_\rd$向下走(为负),$\bfU_\rc$(为正)。但是这并不符合物理,考虑一个管子的颗粒,如果颗粒向下走,必然会带动空气同时向下走。但是方程缺不是这种体现,为什么呢?

    实际上,上述方程并没有写完整,完整形式应该是这样:
    液相:
    \begin{equation}
    \frac{\rd \bfU_\rc}{\rd t}=-\bfA=-\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
    \end{equation}
    气相:
    \begin{equation}
    \frac{\rd \bfU_\rd}{\rd t}=\bfA=\frac{1}{\tau}(\bfU_\rc-\bfU_\rd)
    \end{equation}
    如果开始的时候$\bfU_\rc=0$,则为

  • 1 帖子
    1k 浏览
    FluidGaoF

    请教各位老师,外力驱动的流动,常用哪些算例做benchmark,比如Poiseuille流动,非定常的Womersley流动,还有吗?

  • 2 帖子
    2k 浏览
    bestucanB

    速度的物质倒数两端乘以密度,再化简,就是时间项和对流项加在一起。

    密度就是单位体积的质量,速度的物质导数就是加速度。
    所以时间项和对流项加在一起就是单位体积内的$F=ma$,即$F=\rho a$。这就是惯性力吧?

    力分为体积力和表面力,本来表面力应该是产生相对于体心的力矩。但是无限小体积的情况下把力臂忽略了(我猜的),力矩造成的转动靠物理规则演化出来而不故意设计出来。所以左边是流体微元变化需要的力,右边是流体微元真实受到的力,两者相等,就推出来流体微元怎么变化的了。

    物理意义有助于了解方程怎么来的,但是后来方程形式因为代换化简早没有了当初的模样,只能从大概上理解。像广义的源项,这个和运输方程放在一起好理解。NS方程化简到最后是三个运输方程,运输的物理量分别是密度、速度、熵还是啥。源项,就像一个无源封闭体内进来多少就出去多少,但是有了源,就兴许进来的少出去的多。就是对运输产生了扰动。

    要是搁到动量方程上,比如加个源项,磁流体在磁场中受电磁力,或者离心机里受额外的等效重力。

    另外粘性造成的应力是和速度有关系的吧?叫本构关系,用来区别牛顿和非牛顿流体的。系数小不打紧,只要湍流涡够多,速度够快。耗散能量的能力就很强。

  • 3 帖子
    3k 浏览
    FluidGaoF

    @东岳 那请问东岳老师,在t=0,和t=1时刻,这个函数返回的值是什么意义呢?

  • 5 帖子
    4k 浏览
    Z

    只算过解析解,不知道你用的是什么方法算flowmap?

  • 请问如何使用N-S方程反求粘度?

    1
    1 帖子
    1k 浏览
    S

    如题,这段时间在做一个分子模拟的案例,需要应用到N-S方程反求粘度,想知道这里边的一些基本思路。

  • 波方程

    1
    1 帖子
    1k 浏览
    李东岳

    \begin{equation}
    \zeta =x-ct,\eta=x+ct
    \end{equation}
    \begin{equation}
    \frac{\p\eta}{\p x}=1,\frac{\p\zeta}{\p x}=1
    \end{equation}
    \begin{equation}
    \frac{\p\eta}{\p t}=c,\frac{\p\zeta}{\p t}=-c
    \end{equation}
    对任意变量关于$\zeta,\eta$的函数
    \begin{equation}
    \frac{\p}{\p x}=\frac{\p}{\p\eta}\frac{\p\eta}{\p x}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p x}=\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}
    \end{equation}
    \begin{equation}
    \frac{\p}{\p t}=\frac{\p}{\p\eta}\frac{\p\eta}{\p t}+\frac{\p}{\p\zeta}\frac{\p\zeta}{\p t}=c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}
    \end{equation}
    \begin{equation}
    \frac{\p^2}{\p x^2}=\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)\left(\frac{\p}{\p\eta}+\frac{\p}{\p\zeta}\right)=\frac{\p^2}{\p\eta^2}+\frac{\p^2}{\p\zeta^2}+2\frac{\p^2}{\p\eta\p\zeta}
    \end{equation}
    \begin{equation}
    \frac{\p^2}{\p t^2}=\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)\left(c\frac{\p}{\p\eta}-c\frac{\p}{\p\zeta}\right)=c^2\frac{\p^2}{\p\eta^2}+c^2\frac{\p^2}{\p\zeta^2}-2c^2\frac{\p^2}{\p\eta\p\zeta}
    \end{equation}

  • 多变量分布矩

    1
    1 帖子
    1k 浏览
    李东岳

    二变量高斯分布:
    \begin{equation}
    f(u,v)=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left(-\frac{1}{2(1-\rho^2)}\left[\frac{(u-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(u-\mu_1)(v-\mu_2)}{\sigma_1\sigma_2}+\frac{(v-\mu_2)^2}{\sigma_2^2}\right]\right)
    \end{equation}
    MGF为:
    \begin{equation}
    m_{i,j}=\exp\left(i\mu_1+j\mu_2+0.5(\sigma_1^2i^2+2\rho\sigma_1\sigma_2ij+\sigma_2^2j^2)\right)
    \end{equation}
    纯矩计算方法为
    \begin{equation}
    \begin{split}
    m_{0,0}&=1\\
    m_{1,0}&=\mu\\
    m_{2,0}&=\mu^2+\sigma^2\\
    m_{3,0}&=\mu^3+3\mu\sigma^2\\
    \end{split}
    \end{equation}
    假设$\mu_1=10,\mu_2=20,\sigma_1=\sigma_2=2,\rho=0.5$,有纯矩:
    \begin{split}
    m_{0,0}&=1\\
    m_{1,0}&=10\\
    m_{2,0}&=104\\
    m_{3,0}&=1120\\
    m_{0,1}&=20\\
    m_{0,2}&=404\\
    m_{0,3}&=8240\\
    \end{split}
    同时有混合矩
    \begin{equation}
    m_{i,j}=\exp\left(36\right)
    \end{equation}
    不行,混合矩计算方法不对

  • 17 帖子
    13k 浏览
    H

    @J 把你最近学习的成果总结一下,发上来嘛

  • 在两相流代码中添加空化模型的一些问题

    2
    2 帖子
    2k 浏览
    李东岳

    看起来很复杂的样子.... :143:

  • 求关于openfoam矩阵求解器的相关论文

    5
    5 帖子
    4k 浏览

    @队长别开枪 感谢感谢 😊

  • 一个非守恒的连续性方程

    3
    3 帖子
    3k 浏览
    B

    @东岳 谢谢李老师的回复!
    我看了下跟公式(7)到(8)的过程很相似,但是为什么公式7到公式8的过程中将公式7左边第四项给略去呢?

  • HR算法应用于求解PBE

    2
    2 帖子
    2k 浏览
    李东岳

    这是用离散法求解PBE的增长过程。如果没有增长过程,如果不是离散法,就不需要这么做。a是波速,Minmod是limiter。如果Minmod目前你还不熟悉,我大体知道你的CFD基础在哪,你需要补CFD基本理论,再来看PBE:138:

  • 流动分离后减少的总压去哪里了?

    5
    5 帖子
    3k 浏览

    @东岳 总压减小是因为速度变小,还是静压变小,或二者都有?无粘不可压的时候就是动能与压力势能的转化,是守恒的。

  • 8 帖子
    5k 浏览
    C

    @Gordonaero 感谢,这两篇论文之前没有看到过。
    这个动态填充的研究是想预测填充以及降温时间吗?关于这方面管道预冷的研究,之前看过国内的西交大有做这方面的研究,包括实验与数值计算。
    按照我的理解,这个均相模型将混合物作为一种赝流体,封闭的关键其实是混合物的状态方程以及如何将混合物拆解成气液两相的比例。

  • 一维不可压缩流体如何计算压力变化?

    5
    5 帖子
    4k 浏览
    李东岳

    SIMPLEC在这里:http://dyfluid.com/simplefoam.html
    另外,@Gordonaero 可以参考这个 http://www.cfd-china.com/topic/3395

  • entropy function的一个问题

    4
    4 帖子
    4k 浏览
    李东岳

    $\epsilon \rightarrow 0$ 并不是$\epsilon = 0$,否则不会有...if $u$ is smooth at $x_1$ and $x_2$...。在粘度趋向于很小的时候,不连续变成具备一定厚度的光滑解,同样承认有厚度的激波。所以
    $$
    \epsilon\int_\Omega\frac{\p}{\p u}\left(\frac{\p \eta}{\p u}\right)\left(\frac{\p u}{\p x}\right)^2\rd x\rd t \geq 0
    $$
    另外,
    \begin{equation}
    \int_{x_2}^{x_1}\left(\epsilon(\eta_q q_x)_x\rd x -\epsilon\eta(\eta_q)_q q_x^2\right)\rd x=\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)-\epsilon\eta(\eta_q)_q q_x^2\Delta x
    \end{equation}
    考虑一个非常小的$\epsilon=1e-10$,在控制体内$\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)\rightarrow 0$,$\epsilon\eta(\eta_q)_q q_x^2\Delta x$还是大于0.

  • 关于《icoFOAM解析》的一些疑问

    4
    4 帖子
    3k 浏览

    @东岳 感谢回复。建议原文适当增加配图,这将提高可阅读性。

  • 什么是薄层粘性假设?

    3
    3 帖子
    2k 浏览

    @东岳 Thin-Layer Navier—Stokes Approximation