Skip to content
  • 0 赞同
    10 帖子
    17k 浏览
    M

    @ibelief 请问是视频教程么,可否分享一下,最近在自学这本书,有视频学起来应该效果更好

  • 方柱绕流的周期性网格

    Meshy
    5
    0 赞同
    5 帖子
    9k 浏览
    Z

    是否为容差 matchTolerance 太大导致的衔接不上~

  • 拉格朗日粒子

    Algorithm
    1
    0 赞同
    1 帖子
    4k 浏览
    李东岳

    在1954年,IBM创立了Fortran语言。

    在1955年,多相流中的粒子法就开始进行了。在当时,其被取名为Particle In Cell (PIC)。

    PIC由Los Alamos国家实验室的F.H. Harlow提出。在PIC中,拉格朗日粒子通过其坐标(x, y, z)进行定位并且具有质量。有些情况下,拉格朗日粒子还具有多个物理特征。例如在PIC中,拉格朗日粒子可能还具有大小的区别等。

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    640k 浏览
    H

    我写 固我在,嘿嘿

  • 变比热容计算

    Algorithm
    4
    0 赞同
    4 帖子
    8k 浏览
    L

    之前有点问题,最终如下
    \par $C_p$采用4次多项式分段拟合
    \begin{equation}\label{equ:NS_Cp}
    C_p(T)=a_1+a_2T+a_3T^2+a_4T^3+a_5T^4
    \end{equation}
    \par静焓
    \begin{equation}\label{equ:NS_H}
    H(T)=\int_{T_{0}}^{T_{x}}C_p(T)dT=a_1T+\frac{a_2}{2}T^2+\frac{a_3}{3}T^3+\frac{a_4}{4}T^4+\frac{a_5}{5}T^5+a_6
    \end{equation}
    \par熵
    \begin{equation}\label{equ:NS_S}
    S(T)=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}=a_1\ln{T}+a_2T+\frac{a_3}{2}T^2+\frac{a_4}{3}T^3+\frac{a_5}{4}T^4+a_7
    \end{equation}
    \par$\bullet$求解静温(已知总温和马赫数)
    \begin{equation}\label{equNSUs}
    U_s^2=2\left(H(T_{tot})-H(T_{sta})\right)
    \end{equation}
    \begin{equation}\label{equNSmach}
    Mach^2=\frac{U^2}{\gamma(T)R_gT}=\frac{U^2}{\frac{C_p(T)}{C_p(T)-R_g}R_gT}
    \end{equation}
    \par由(\ref{equ:NSmach})和(\ref{equ:NSUs})得
    \begin{equation}\label{equNSTsta}
    T_{sta}=\frac{2\left(H(T_{tot})-H(T_{sta})\right)}{Mach^2\frac{C_p(T_{sta})R_g}{C_p(T_{sta})-R_g}}
    \end{equation}
    $\bullet$求解静压(已知总温、总压和静温)
    \par由p等熵过程
    \begin{equation}
    ds = C_p(T)\frac{dT}{T} -R_g\frac{d p}{p}=0
    \end{equation}
    \par两边同时积分有
    \begin{equation}
    \int_{T_{tot}}^{T_{sta}}C_p(T)\frac{dT}{T} =\int_{p_{tot}}^{p_{sta}} R_g\frac{d p}{p}
    \end{equation}
    \par记
    \begin{equation}
    S(T_{x})=\int_{T_{0}}^{T_{x}}C_p(T)\frac{dT}{T}
    \end{equation}
    \par则
    \begin{equation}
    S(T_{sta}) - S(T_{tot}) = R_g\ln\frac{p_{sta}}{p_{tot}}
    \end{equation}
    \par那么
    \begin{equation}\label{equNS_psta}
    p_{sta}=p_{tot}e^{\left(\frac{S(T_{sta})-S(T_{tot})}{R_g}\right)}
    \end{equation}

  • 0 赞同
    5 帖子
    8k 浏览
    肖艺

    @yfclark

  • 关于processor边界的问题

    OpenFOAM
    5
    0 赞同
    5 帖子
    6k 浏览
    A

    @东岳
    我现在是对所有网格做遍历,遍历每个网格上的面,判断是内部的面还是边界的面,如果是边界的面的话,通量值就会和另一块processor上对应的网格上的值相关,目前就是不知道怎么去找这个网格。

  • 请问哪种机械键盘适合敲字?

    CFD彩虹条
    4
    0 赞同
    4 帖子
    4k 浏览
    李东岳

    @桎梏 那青轴可以说是爆炸 :xinxin:

  • 0 赞同
    17 帖子
    26k 浏览
    liujunCFDL

    @TINGHAIK 谢谢前辈

  • 测试下新表情

    CFD彩虹条
    1
    0 赞同
    1 帖子
    2k 浏览
    李东岳

    :wolaile: :wolaile: :wolaile: :wolaile: :wolaile:

    :xiabanle: :xiabanle: :xiabanle: :xiabanle:

    :tishizi:

  • mesh转换前还要dos2unix吗?

    Meshy
    4
    0 赞同
    4 帖子
    7k 浏览
    cfd_xfxC

    dos2unix命令是用来将DOS格式的文本文件转换成UNIX格式下的文本文件,这两个系统下的文本文件的差异在于换行标志,DOS下是用\r\n断行,UNIX下采用\n作为断行。将windows下面写好的一个文本文件转移到linux下,无法直接运行,打开会发现出现很多^M的标志,可以采用dos2unix命令进行清除。dos2unix命令起且仅起这个作用。因此,对于windows下生成的mesh文件,如果转移到linux下能够正常运行跑通算例,则完全不需要dos2unix命令。
    同样,checkMesh也是,如果算例都能跑,则不需要这一步骤hhhh:papa:

  • [ICEM]如何画出如下的Blocks?

    Meshy
    5
    0 赞同
    5 帖子
    9k 浏览
    Z

    整体划分,然后删除多余的block,看看好用不

  • PBM之QMOM中最后一步。。。

    Algorithm
    11
    0 赞同
    11 帖子
    18k 浏览
    AstremonA

    @李东岳 http://dx.doi.org/10.1016/j.compchemeng.2015.09.008
    这篇文章对三种重构方法做了对比。

  • -div(phi,T) 和div(-phi,T)结果为何不同?

    OpenFOAM
    12
    0 赞同
    12 帖子
    19k 浏览
    李东岳

    这是div(phi,T)与-div(phi,T)等的解释

    微信图片_20200320103052.jpg

  • 用driftFluxFoam时出现报错

    OpenFOAM
    3
    0 赞同
    3 帖子
    4k 浏览
    T

    @东岳 好的,谢谢老师!

  • paraView-v1906安装问题

    OpenFOAM
    4
    0 赞同
    4 帖子
    6k 浏览
    李东岳

    http://www.cfd-china.com/topic/3182

    有一些介绍,你试试二楼的

  • 关于温度对称边界条件

    Algorithm
    6
    0 赞同
    6 帖子
    11k 浏览
    L

    @李东岳 谢谢李老师

  • 0 赞同
    2 帖子
    5k 浏览
    chanstormstoutC

    0_1519692297906_捕获.JPG
    是否是这样?应力张量我可以理解。我的疑问是速度梯度在笛卡尔坐标系下变换之后直接就变成在极坐标系中了速度梯度了?

  • 0 赞同
    2 帖子
    4k 浏览
    E

    大神如何让stl 导入 运行第一步blockMesh?
    我一运行就报错

  • C++11数组初始化

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    R

    @wwzhao

    简单的感受了一下两种赋值方式,效率确实很大提高。

    #include <stdio.h> #include <string> #include <chrono> #include <iostream> using namespace std ; class Timer { private: // Type aliases to make accessing nested type easier using clock_t = std::chrono::high_resolution_clock; using second_t = std::chrono::duration<double, std::ratio<1> >; std::chrono::time_point<clock_t> m_beg; public: Timer() : m_beg(clock_t::now()) { } void reset() { m_beg = clock_t::now(); } double elapsed() const { return std::chrono::duration_cast<second_t>(clock_t::now() - m_beg).count(); } }; class SomethingBefore11 { private: int m_array[5]; public: SomethingBefore11() // zero the member array { // If we want the array to have values, we'll have to use assignment here m_array[0] = 1; m_array[1] = 2; m_array[2] = 3; m_array[3] = 4; m_array[4] = 5; } }; class SomethingAfter11 { private: int m_array[5]; public: SomethingAfter11(): m_array { 1, 2, 3, 4, 5 } //zero the member array { } }; int main(){ Timer tBefore11; SomethingBefore11 m_array_before11; std::cout << "Time elapsed: " << tBefore11.elapsed() << ‘n’; Timer tAfter11; SomethingAfter11 m_array_after11; std::cout << "Time elapsed: " << tAfter11.elapsed() << ‘n’; return 0; } [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.506e-06 Time elapsed: 1.47e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 8.664e-06 Time elapsed: 1.9e-07 [xx OFtutorial0_helloWorld]$ whatAboutThisGuy Time elapsed: 7.646e-06 Time elapsed: 1.89e-07