Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. thw1021
    T
    • Profile
    • Following
    • Followers
    • Topics
    • Posts
    • Groups

    thw1021

    @thw1021

    7
    Posts
    126
    Profile views
    0
    Followers
    0
    Following
    Joined Last Online

    thw1021 Follow

    Best posts made by thw1021

    This user hasn't posted anything yet.

    Latest posts made by thw1021

    • RE: 如何隐式离散应变率张量

      @李东岳 至少论文里面是这样说的。The strain rate tensor S is treated implicitly in the modified flow solver.

      posted in OpenFOAM
      T
      thw1021
    • 如何隐式离散应变率张量

      请教大家一个问题。
      如果将雷诺应力按如下方式分解
      $$
      \boldsymbol{\tau}^{m}=\nu_{t}^{L} \mathbf{S}+\left(\boldsymbol{\tau}-\boldsymbol{\tau}^{L}\right)+\operatorname{tr}(\boldsymbol{\tau})
      $$
      其中$\nu_{t}^{L}$可以理解为湍流涡黏系数(但有所不同)。$\mathbf{S}=\frac{1}{2}\left[\nabla \mathbf{U}+(\mathbf{\nabla} \mathbf{U})^{T}\right]$是应变率张量(strain-rate tensor)。
      如果要求$\mathbf{S}$隐式离散,在OpenFOAM中应该如何实现?(The strain rate tensor $\mathbf{S}$ is treated implicitly in the modified flow solver.)

      我个人的理解:
      以不可压缩流动为例,OpenFOAM-7中对应的雷诺应力计算方式为:

      template<class BasicTurbulenceModel>
      Foam::tmp<Foam::fvVectorMatrix>
      Foam::linearViscousStress<BasicTurbulenceModel>::divDevRhoReff
      (
          volVectorField& U
      ) const
      {
          return
          (
            - fvc::div((this->alpha_*this->rho_*this->nuEff())*dev2(T(fvc::grad(U))))
            - fvm::laplacian(this->alpha_*this->rho_*this->nuEff(), U)
          );
      }
      

      因此,只有$\nabla \mathbf{U}$是隐式离散,而$(\mathbf{\nabla} \mathbf{U})^{T}$其实是显式离散的。并且也没有直接使用$\mathbf{S}$,相当于是把 $2\mathbf{S}$ 做了分解。

      如果要对$\mathbf{S}$进行隐式离散,是不是应使用fvm::div() ?
      $\nu_{t}^{L} \mathbf{S}$ 应该如何表示 ?

      posted in OpenFOAM
      T
      thw1021
    • RE: C币运行机制测试

      好久没登录,想请教个问题,突然提示没有C币,只好来水一下了:mihu:
      水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水水

      posted in C斯达克
      T
      thw1021
    • RE: foam-extend 如何对特定区域采用特定的分区方法?

      @东岳 李老师,不是同一个人。

      posted in OpenFOAM
      T
      thw1021
    • RE: foam-extend 如何对特定区域采用特定的分区方法?

      @溯光 icoDyMIbFoam

      posted in OpenFOAM
      T
      thw1021
    • RE: foam-extend 如何对特定区域采用特定的分区方法?

      参考资料:
      【1】 http://www.cfd-china.com/topic/231/openfoam并行计算如何手动划分区域

      【2】 https://openfoam.com/releases/openfoam-v1712/parallel.php

      posted in OpenFOAM
      T
      thw1021
    • foam-extend 如何对特定区域采用特定的分区方法?

      问题来源:
      foam-extend中使用IBM求解器时,并行总是出现问题,似乎需要对部分区域进行人工分区。

      已解决:
      参考【1】,使 setFields,对边界存在的区域分配指定的核。

      存在的问题:

      首先使用 decomposePar 对整个计算域进行分区,然后再使用 setFields 对特定区域进行分区,导致不同的核处理的网格量有很大差别,降低了计算效率。

      待解决的问题:

      能否对特定的区域指定分区方法,然后运行 decomposePar 同时分区,而不是必须先对整个计算域进行分区。如附图所示,对区域1使用人工分区,对区域2使用 scotch 分区,以保证每个核处理的网格量相近。我找到了一些参考资料【2】,但没有在 foam-extend4.0 中试验成功。

      图片1.png

      posted in OpenFOAM
      T
      thw1021