CFD中文网

最近在编写程序的时候遇到并行边界上网格面心处数值通量(φ_f)的同步问题，单核运行不存在并行边界，所以没有这个问题。以下图为例：

并行计算时，求解区域分为了4个区域，两个区域之间产生了并行边界面，对于内部网格面元，网格面心处的数值通量(φ_f)的计算跟单核运行时是一样的，但是对于并行边界处的网格面元，数值通量(φ_f)在两个相邻并行边界上应该是一样的，因为并行边界的存在不应该改变原来的解。因此，我用了如下的syncFaceList()函数来进行同步，alpha_f是定义在网格面心上的surfaceScalarField，这个函数在头文件"syncTools.H"中有定义，我用了以下三种同步方式，发现流体穿过并行边界时还是产生了误差，单核运行时误差(数值解与解析解之差)可以达到10E-12次方量级，但是并行运行时流体一穿过并行边界，误差就上升到了10E-4量级，我觉得这是因为并行边界上数据同步没有做好的原因，请问有小伙伴遇到过这个问题嘛？使用syncFaceList()函数的正确方式是什么？如果使用syncFaceList()函数进行同步网格面心处的数值通量(φ_f)还是有问题的话，打算自己写一下同步过程。

// syncTools::syncFaceList(mesh(), alpha_f, maxEqOp<scalar>(), false);
syncTools::syncFaceList(mesh(), alpha_f, plusEqOp<scalar>(), false);
// syncTools::syncFaceList(mesh(), alpha_f, eqOp<scalar>(), false);

最近在收集一些计算流体力学招聘的信息，看到华为有算法这面的招聘，当然还有热设计工程师，有老铁了解里面的工作内容和加班情况吗？我想了解一下，虽然不一定能进去

@浪迹天大 在 运行时动态加载库 中说：

@Micro 我有个朋友写过相关的博文，感兴趣可以看看
https://zmeng.tech/2020/10/DynamicLibraries/

https://zmeng.tech/2020/10/DynamicLibraries/

支持！李老师是国内CFD界为数不多有想法而且行动力强的人，希望你创业成功，开创一个新的环境。

@huicfd 并行时每个进程完成该进程的任务，这是只输出了0号进程（openfoam好像默认该进程为master prosessor），而其他进程的结果没有输出。要获取其他进程计算得到的数据，需要同时使用open foam中的Pstrean::gatherList函数和Pstream::scatterList函数进行数据同步。
贴一段我的测试代码，希望对你有帮助。

if ( Pstream::parRun() ) 
{   
    List<scalarField> test(Pstream::nProcs());   
    test[Pstream::myProcNo()].setSize(3, Pstream::myProcNo());   
    Pstream::gatherList<scalarField>(test); 
    Pstream::scatterList(test); 
    // Pout << test << endl; 

    List<scalar> numcell(Pstream::nProcs());  
    numcell[Pstream::myProcNo()] = mesh.nCells();  
    Pstream::gatherList<scalar>(numcell); 
    Pstream::scatterList<scalar>(numcell); 
    // reduce( numcell, sumOp<List<scalar> >() );  

    reduce( nCells, sumOp<scalar>() );  
    // Pout << "CPU = " << Pstream::myProcNo() << ", masterProcNo = " 
    // << Pstream::masterNo() << endl;   

    if ( Pstream::master() ) 
    { 
        scalar size = 0.0; 
        for ( label procI=0; procI<Pstream::nProcs(); procI++ ) 
        {
            size += numcell[procI]; 
        }
        // Pout << "grid size = " << size << endl; 
        // Pout << "Cell size in each Process is " << numcell 
        // << ", Total cell size = " << nCells << endl; 
    }   

}   

@zp气固两相流 博主这个问题解决了吗？想知道后续

老铁找到了吗？可以分享一下面试了哪些CFD公司吗？

我把fixed_value文件夹放在$(WM_PROJECT_DIR)这个目录下就可以编译成功了，看来openfom中不允许把库文件放在任意目录下。链接文本

@新喻庸
file中的配置：

oscillatingfixedvalueFvPatchFields.C
LIB = $(FOAM_USER_LIBBIN)/liboscillating    

options中的配置：

EXE_INC = \
    -I~/Desktop/MyFile/NewBC3/fixed_value/lnInclude \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude

LIB_LIBS = \
    -lfiniteVolume \
    -lmeshTools \
    -L$(FOAM_USER_LIBBIN) -lfixed_value

最近在编写一个边界条件的类，我打算把这个边界条件的类做成库函数。但是在编译的过程中遇到一些问题，感觉是Make中的file和options没有配置好，有没有CFDer遇到过类似的问题，请指教一下！下面贴上我的file和options配置：
file中的配置：

oscillatingfixedvalueFvPatchFields.C
LIB = $(FOAM_USER_LIBBIN)/liboscillating

options中的配置：

EXE_INC = \
    -I~/Desktop/MyFile/NewBC3/fixed_value/lnInclude \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude

LIB_LIBS = \
    -lfiniteVolume \
    -lmeshTools \
    -L$(FOAM_USER_LIBBIN) -lfixed_value

以上为Make中的配置内容，用命令(wmake libso)编译时提示的错误为：

错误说找不到文件“fixedvalueFvPatchFields.H”，但是在options中不是已经配置了吗？这里有一些疑惑，请各位CFDer指教！