在N-S微分方程中加入狄拉克函数（脉冲函数）源项，在OF中如何实现？

[hongfu2233]

在N-S微分方程中加入狄拉克函数（脉冲函数）源项，在OF中如何实现？请教一下各位有什么思路吗？这个函数怎么用编程语言实现？

[CFD中文网]

@hongfu2233

狄拉克函数，就是在某点存在值，其他为0. 为何不尝试用if语句？

[hongfu2233]

多谢啊，我是小白，不知道在哪个位置加入这个代码，以icoFoam为例，能给我指点一二吗？

[CFD中文网]

@hongfu2233

源项中添加这个函数，还得具体问题具体分析。如果单说这个函数，可以这样写（C++）：

if(x == 1.0)
{
    y = 1;
}
else
{
    y = 0;
}

就是这样..

[hongfu2233]

@cfd-china 多谢啦，我试试！