feiqunali 发表于 2010-12-5 20:40:27

请问openfoam下如何算升阻力系数

这个问题困扰我很久了,一直都得不到满意的答复。我在算一个圆柱绕流问题,想要得到升阻力及系数的情况,之前有人看到网上有人说用function object功能,然后我照着别人的办法设置controlDick文件,但是输出结果只有U和P,并没有我要的升力和阻力,请问有人能告诉我这是什么情况么?我该如何着手?或者说计算升阻力有没有别的办法?希望有达人能详细的介绍一下,不胜感激!

OpenFOAM 发表于 2010-12-7 17:14:32

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

FoamFile
{
    version         2.0;
    format          ascii;
    class         dictionary;
    location      system;
    object          sampleDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Set output format : choice of
//      xmgr
//      jplot
//      gnuplot
//      raw
setFormat raw;

// Surface output format. Choice of
//      null      : suppress output
//      foamFile    : separate points, faces and values file
//      dx          : DX scalar or vector format
//      vtk         : VTK ascii format
//      raw         : x y z value format for use with e.g. gnuplot 'splot'.
//
// Note:
// other formats such as obj, stl, etc can also be written (by proxy)
// but without any values!
surfaceFormat raw;

// interpolationScheme. choice of
//      cell          : use cell-centre value only; constant over cells (default)
//      cellPoint   : use cell-centre and vertex values
//      cellPointFace : use cell-centre, vertex and face values.
// 1] vertex values determined from neighbouring cell-centre values
// 2] face values determined using the current face interpolation scheme
//    for the field (linear, gamma, etc.)
interpolationScheme cellPoint;

// Fields to sample.
fields
(
    p
);


// Set sampling definition: choice of
//      uniform             evenly distributed points on line
//      face                one point per face intersection
//      midPoint            one point per cell, inbetween two face intersections
//      midPointAndFace   combination of face and midPoint
//
//      curve               specified points, not nessecary on line, uses
//                        tracking
//      cloud               specified points, uses findCell
//
// axis: how to write point coordinate. Choice of
// - x/y/z: x/y/z coordinate only
// - xyz: three columns
//(probably does not make sense for anything but raw)
// - distance: distance from start of sampling line (if uses line) or
//             distance from first specified sampling point
//
// type specific:
//      uniform, face, midPoint, midPointAndFace : start and end coordinate
//      uniform: extra number of sampling points
//      curve, cloud: list of coordinates
sets
(

);


// Surface sampling definition: choice of
//      plane : values on plane defined by point, normal.
//      patch : values on patch.
//
// 1] patches are not triangulated by default
// 2] planes are always triangulated
// 3] iso-surfaces are always triangulated
surfaces
(
    airfoil
    {
      type            patch;    // always triangulated
    patchName    Airfoil;
    }
   slat
    {
    type      patch;
    patchName    Slat;
    }
);


// *********************************************************************** //
页: [1]
查看完整版本: 请问openfoam下如何算升阻力系数