Skip to content
  • 0 赞同
    1 帖子
    2k 浏览
    李东岳

    现在很多AI号,翻译sci基本就是10秒钟一篇,就比如那种jcp昨天见刊的文章,今天国内的马上出来一篇中文爆火全网。这些都是ai写的。感兴趣的可以把类似这些文章的文字复制到这里然后识别一下。

    我测试了一半ai写的,一半人工写的,准确率非常高,太厉害了。

    https://smodin.io/zh-cn/ai内容检测器

    捕1获.JPG

  • 0 赞同
    2 帖子
    5k 浏览
    李东岳

    VOF+动网格是interDyMFoam,用这个求解器就行

  • 1 赞同
    3 帖子
    6k 浏览

    @李东岳

    常见的问题是求$\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.有的连矩阵都不需要组装,何苦去求逆...
  • 0 赞同
    6 帖子
    7k 浏览
    X

    @benf 不客气不客气,欢迎引用我们组的论文~这个模型有外国组已经验证过了,所以如果你有兴趣编应该可以原样复现

  • 有关泡沫金属与石蜡传热传质问题

    Fluent
    1
    0 赞同
    1 帖子
    2k 浏览
    Z

    替代文字 我模拟的是一个矩形泡沫铜浸润石蜡的传热传质过程,现在模型验证一直没弄好。首先是编写了一个简单的UDF,采用热平衡模型,导入fluent后提醒我定义比热的那一块有错误,求各位前辈帮我看一下;其次,我添加了Boussinesq项后,优化了一点点液相云图中熔融线一直平行的问题,但是与论文中的还是有差距,且我的熔融线上下厚度保持一致(实际上应该下面厚很多),UDF中还需要添加什么项来解决这个问题呢?以下三张图分别是我的UDF、我跑出来的液相云图以及论文的液相云图,求各位帮我解惑,不甚感激。2.png 3.pn替代文字g 4.png

  • 关于划分网格的无量纲数的定义

    Meshy
    2
    0 赞同
    2 帖子
    5k 浏览
    D

    @hhj0219 和y+定义一样。其实就是:流向网格间距除以壁面第一层距离,然后再乘以y+。

  • 请问FLUENT能不能模拟真空

    Fluent
    8
    0 赞同
    8 帖子
    10k 浏览
    F

    我觉得应该搞清楚液滴在真空状态下的变化,然后针对性的设置环境条件

  • 圆柱相切的划分策略

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

    谢谢老师~~

  • 计算流体力学开源软件推荐

    CFD彩虹条
    3
    0 赞同
    3 帖子
    5k 浏览

    PeleC(AMReX)

  • 气液固多相流

    Algorithm
    2
    0 赞同
    2 帖子
    5k 浏览
    C

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

  • CFD结果渲染

    Fluent
    6
    0 赞同
    6 帖子
    7k 浏览
    huntersongH

    x3d
    用这个导出,导入到blender

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

    Algorithm
    3
    0 赞同
    3 帖子
    8k 浏览
    C

    @kasitefy
    非常对。

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

  • CFD青年成长支持计划(2021)

    公告
    137
    0 赞同
    137 帖子
    165k 浏览
    F

    想问一下李老师这个计划还有吗!!很想要这本湍流模型:mianmo:

  • 0 赞同
    1 帖子
    3k 浏览

    工科末流985硕士毕业,想在武汉找CFD相关的工作,家人们有没有推荐的公司。

  • 0 赞同
    1 帖子
    2k 浏览

    a80fed28-45ef-4507-a08d-13abf19c7ccc-image.png 这个udf具体该怎么编写和使用?

  • fluent数据处理成二维数组

    Fluent
    3
    0 赞同
    3 帖子
    5k 浏览
    N

    不知道当时的问题是否解决了,我去年做的一部分工作里面也需要相关的处理。当时采用的Matlab对导出的数据进行处理,不过在数据处理之前需要将 x,y当中一些极小的值变成0,比如10的负十几次方的值,因为这些fluent输出的值可能会对处理产生影响,然后导入程序就可以变成矩阵形式了,程序大概如下。供以后有需要的朋友们参考。

    data1 = xlsread('D:\User.....');

    % 获取X、Y坐标和密度值的向量
    X = data1(:, 1);
    Y = data1(:, 2);
    density = data1(:, 3);

    % 确定X和Y的唯一值,并按升序排列
    unique_x = unique(X);
    unique_y = unique(Y);

    % 创建行列编号映射表
    x_mapping = containers.Map(unique_x, 1:length(unique_x));
    y_mapping = containers.Map(unique_y, 1:length(unique_y));

    % 初始化矩阵
    num_rows = length(unique_y);
    num_cols = length(unique_x);
    densityMatrix1 = zeros(num_rows, num_cols);

    % 将密度值填充到矩阵中
    for i = 1:length(X)
    row = y_mapping(Y(i));
    col = x_mapping(X(i));
    densityMatrix1(row, col) = density(i);
    end

  • OpenFoam DNS计算

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

    谢谢 学到啦 还有一个问题 采用高分辨率网格可以弥补精度不够高的不足嘛

  • paraview打开sprayFoam的案例文件报错

    OpenFOAM
    1
    0 赞同
    1 帖子
    2k 浏览
    mingyangM

    用paraview打开sparyFoam求解的aachenBomb案例文件,结果有这样的报错:
    ERROR: In C:\bbd\8fbebceb\build\superbuild\paraview\src\VTK\IO\Geometry\vtkOpenFOAMReader.cxx, line 7891
    vtkOpenFOAMReaderPrivate (0000022DE32F20E0): Error reading line 20 of E:\OpenFOAM\OpenFOAM-8\sprayFoam\aachenBomb\5e-05/lagrangian/sprayCloud/positions: Expected punctuation token ')', found f
    我用的是Windows环境下的paraview5.11.0。路过的大佬,帮忙看看。

  • 0 赞同
    9 帖子
    10k 浏览
    G

    在foam-extend4.1的ibm里折腾了两三个月,发现用foam-extend4.1的浸没边界法计算升阻力系数产生震荡是其本身的问题。具体原因可以参考这篇论文Sources of spurious force oscillations from an immersed boundary
    method for moving-body problems
    这篇论文提到了三个解决方法:1、在力源项后面再增加一个质量源项以抑制动边界的力震荡 2、增大时间步长 3、将动边界所处位置的网格加密
    目前我只是用了后面两种方法可以有效的减少其震荡,但震荡还是会存在,如果想获得一条平滑曲线,可能需要与第一种方法结合。第一种方法涉及到一些算法问题,鉴于目前的水平有限还不会在openfoam中实现,希望有大佬能够实现并给大伙解惑。

  • mvConvection->fvmDiv

    OpenFOAM
    5
    0 赞同
    5 帖子
    5k 浏览
    李东岳
    tmp<fv::convectionScheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)") ) ); { reaction->correct(); dQ = reaction->dQ(); label inertIndex = -1; volScalarField Yt(0.0*Y[0]); forAll(Y, i) { if (Y[i].name() != inertSpecie) { volScalarField& Yi = Y[i]; fvScalarMatrix YiEqn ( fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi)

    这个看起来并不会导致区别。有可能cfd-online的用户没做特别仔细的判定吧。