怎样统计多相流计算中物体的受力
-
这个是我自己的代码
/*--------------------------------*- 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 ); }
-
我是这么搞的:
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的
-
@tidedrinker 我算的是双流体模型,这个方法好像不行(▼皿▼#)