Zalesak 问题

[李东岳]

用setFields不行，因为是不规则区域，尤其是有圆形存在。只能自己写代码，范例如：https://coding.net/u/dyfluid/p/Solvers_DYFLUID/git/tree/master/taylorGreenU

[浪迹天大]

@东岳
最近发现一个 setFields 的神奇用法，用 STL 来指定形状，就可以实现任意形状了呀。


完整文件见附件：
0_1535073080541_programming_OF_ICs.pdf

[李东岳]

@浪迹天大

看起来不错！

@myjimp

你可以试试做个Zalesak圆盘的STL？然后试试

[myjimp]

@浪迹天大 厉害，非常感谢！

[myjimp]

@东岳 请教一下，，谢谢，，速度场是函数u(x , y)=-π(y -0 .5), v(x , y)=π(x -0 .5)的话，OF怎么设置？？？

[李东岳]

forAll(U, cell)
{
    const scalar& x = mesh.C()[cell].x();
    const scalar& y = mesh.C()[cell].y();

    U[cell].x() = - M_PI*(y - 0.5);
    U[cell].y() = M_PI*(x - 0.5);
}

[myjimp]

谢谢，李老师指导。

[myjimp]

@东岳 李老师你说的这个是在U文件下用这种格式吗？

[myjimp]

@东岳 这段代码在哪里实现

[wwzhao]

Zalesak disk可以参考isoAdvector里面的算例配置：
https://github.com/isoAdvector/isoAdvector/tree/master/OpenFOAM/run/prescribedU/notchedDiscInSolidBodyRotation