怎样统计多相流计算中物体的受力

[dzw05]

各位大佬，在统计物体受力时，如果是单相流可以统计，那么如果是多相流（multiPhaseEulerFoam）呢，怎样给定rho、U和P的变量？

[李东岳]

统计物体受力？什么意思。

rho等存在一个rhom，是相平均密度，其他类似。

[dzw05]

@李东岳 就是和单相流一样，统计一个旋转部件的受力。openfoam有force的functionObject可以统计单相流动力部件受力，但需要指定rho，U和P，分别计算粘性力和压力对受力的贡献。目前好像没找到多相流怎么统计。

[李东岳]

目前我知道的也只是按照单相流动计算。

[队长别开枪]

http://openfoamwiki.net/index.php/Contrib/swak4Foam/Example_sectionalForces

[队长别开枪]

这个是我自己的代码

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1706                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

forcesSectional
{	
	type                    patchExpression;
	//outputControlMode       outputTime;
	enabled                 true;
        verbose                 true;
	patches
        (
		foil                // change to your name or names of patch
	);

        aliases{alpha alpha.water;}
 
	variables
        (
	    "rho_water=997.4;"
	    "nu_water=9.274e-7;"
            "rho_vapor=0.02099;"
	    "nu_vapor=4.68e-04;"
            "u_inf=5.0;"
            "c=0.15;"
	    "pressure_force=p_rgh*normal()*area();"
	    "viscous_force=-((alpha*rho_water + (1.0 - alpha)*rho_vapor))*((alpha*nu_water + (1.0 - alpha)*nu_vapor))*snGrad(U)*area();"
            "dynamic_force=0.5*rho_water*u_inf*u_inf;"
            "area_ref=0.1*c;"
	    "forces_coeff=(pressure_force+viscous_force)/(dynamic_force*area_ref);"
	);
 
        expression "forces_coeff";
        accumulations
        (
            sum
        );
}

[tidedrinker]

我是这么搞的：

    force1
    {
        type forces;
        libs ( "libforces.so" );
        patches (BLADE HUB);
        log     on;
   	rho         rhoInf;      // Indicates incompressible
  	rhoInf      1000;
        writeControl   timeStep;
        writeInterval  4;
        pitchAxis   (1 0 0);
        CofR (0 0 0);
    }

不用管那个rho，二相流里有它没它一个样，程序会自动计算平均rho的

[dzw05]

@队长别开枪 多谢啦😊，我试试

[dzw05]

@tidedrinker 我算的是双流体模型，这个方法好像不行(▼皿▼#)

[小龙]

@队长别开枪 你好，这个是需要安装swak4Foam才可以用的吗？

[队长别开枪]

@小龙 是的。

[小龙]

@队长别开枪 OF5.0版本支持swak4Foam吗？

[队长别开枪]

@小龙 我用的openfoam plus，其他的没试过。你可以查看一下这个
https://openfoamwiki.net/index.php/Installation/swak4Foam/Downloading

[小龙]

@队长别开枪 好的，多谢啦

[王金成]

@队长别开枪 请问您用的是openfoam plus的哪一个版本呢? 好像v1606+也不能用swak4foam

[队长别开枪]

@王金成 我用的1812，可以编译SWAK，不过需要自己做些配置，不然会报错，过段时间我写个教程。