자료실

안녕하십니까,

현재 OpenFOAM-3.0.1을 이용하여 최적화를 수행중에 있습니다.

기존 adjointShapeOptimizationFoam을 수정한 UNAdjointShapeOptimizationFoam 솔버를 Compile 하는 과정에서 아래와 같은 에러가 계속 발생하고 있습니다.

UNAdjointShapeOptimizationFoam.C:50:22: fatal error: RASModel.H: No such file or directory

 #include "RASModel.H"

                      ^

compilation terminated.

make: *** [Make/linux64GccDPInt32Opt/UNAdjointShapeOptimizationFoam.o] Error 1

왜 RASModel.H 파일을 찾아가지 못하는지 궁금합니다.

솔버파일을 첨부하였습니다.

안녕하세요?

scalarTransportFoam 관련해서 문의를 드립니다.

현재 하려고 하는 일이 scalarTransportFoam 의 수송 방정식을 simpleFoam 솔버에 추가해서 농도장을 계산 하는 코드를 구성하는게 목표인데요

(icoFoam 솔버에 온도장을 푸는 방정식을 추가하는 부분은 수행해봤습니다)

simpleFoam 솔버 폴더를 복사해서 새로운 폴더(mysolver)를 만들고

my 솔버에(아래)

-----------------------------------------------------

    while (simple.loop())

    {

        Info<< "Time = " << runTime.timeName() << nl << endl;

        // --- Pressure-velocity SIMPLE corrector

        {

            #include "UEqn.H"

            #include "pEqn.H"

        }

        laminarTransport.correct();

        turbulence->correct();

/*------------추가-------------*/    

       fvScalarMatrix cEqn

        (

               fvm::div(phi,C)

            - fvm::laplacian(turbulence->nuEff()/0.9 , C)

            );     

        cEqn.solve();

/*------------추가-------------*/  

        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"

            << "  ClockTime = " << runTime.elapsedClockTime() << " s"

            << nl << endl;

    }

---------------------------------------------------------- 

위와같이 식을 추가했습니다

코드를 컴파일한 후(에러없이 컴파일 되었습니다)  

0폴더에 C 파일을 추가했습니다 그리고 constant 폴더의 transportProperties 에 concentration diffusivity(kc) 값도 추가하였습니다 

createFields.H 파일 내 transportProperties 내용에 농도장 C를 불러오는 volScalarField C 와 diffusivity 값이 불러 와지도록 구문을 추가했습니다.

위와같이 경계조건 파일들을 만든 후에 솔버를 실행하면 속도장은 농도장 방정식을 추가하기 전처럼 계산이 되지만 농도장은 계산이 되지 않습니다

(농도장 초기조건은 velocity inlet 부분에 무차원 1값을주고 나머지 경계조건에는 zerogradient 조건으로 줬습니다)

혹시 제가 지금까지 했던 과정중에 빠뜨리거나 잘못한 부분이 있으면 조언 부탁드립니다.

 ----------------------------------------------------------------

아 그리고 추가질문으로  

simplefoam.c 의 UEqn 부분에

      + turbulence->divDevReff(U)

부분이 있는데 이 부분을 보려고 하면 어디로 들어가면 되는지 혹시 아시나요?

divDevReff(U) 수식이 책에는 

divDevReff(U) = - fvm :: laplacian ( nuEff () , U )

                       - fvc :: div ( nuEff () * dev ( T ( fvc :: grad ( U ) ) ) )

라고 나오는데, 이부분을 확인하고싶어서  

src/TurbulenceModels/...... 하위 폴더에 있는 소스파일들을 뒤져도 찾지를 못하겠네요;;

 혹시 이부분도 알고 계신다면 조언 부탁드립니다.

읽어주셔서 감사합니다.  

안녕하세요.

얼마 전에 글을 올렸었는데 아직 문제가 해결되지를 않아 다시 올립니다.

지금 사용하고 있는 2.3.0 버전은 이후 버전과 다르게 vorticity.C 파일이 openfoam230/applications/utilities/postProcessing/velocityField/vorticity에 위치하고 있고 vorticity.H 파일이 없는 것 같습니다.

지금 source term에 식을 이용하여 새로운 값을 계산하는데 vorticity, magvorticity의 값이 계산의 일부로 들어갑니다.

전에 말씀하신 방법은 postprocessing을 위해 서만 function object를 사용해 시간에 따른 vorticity 값을 계산하는 방법이라 생각하는데요, 지금 경우같이 postprocessing에서 필요한 값을 얻기 위한 것이 아닌 시간마다 source term에 식으로 새로운 값을 추가 시 식의 부분 값으로 vorticity의 값이 필요한 경우 vorticity.C가 소속된 디렉토리를 작업디렉토리로 copy해서 관련 application source file과 같이 Make 디렉토리에 있는 file, option에서 수정을 한 후 같이 새롭게 compile해야 되는 것이 아닌가 생각합니다. 그래서, 이러한 방법이나, 아니면 vorticity.C가 소속된 디렉토리를 원래 위치에 두고 관련 application 파일이 있는 디렉토리만을 가지고 compile을 해도 계속해서   

pisoVCFoam.C:95:37: error: ‘vorticity’ was not declared in this scope 

pisoVCFoam.C:118:36: error: ‘vorticity’ was not declared in this scope

pisoVCFoam.C:118:47: error: ‘magVorticity’ was not declared in this scope

위와 같은 에러가 발생하고 있습니다. 참고로 pisoVCFoam은 제가 새로운 방법을 implement하기위해 만든 새로운 pisoFoam 입니다.

만약 파일을 직접 살펴보시기를 원하시면 이메일주소를 알려주시면 첨부하겠습니다. 

감사합니다. 

안녕하세요

지금 OpenFOAM 2.3.0 을 사용하고 있는데 OpenFOAM extended 3.2를 우분투 12.04에 설치해서 사용하여도 기존 OpenFOAM 2.3.0 과 충돌같은 것이 없는지 알고 싶고 OpenFOAM extended 3.2를 설치하여 사용 시 기존의 설치한 pyfoam, swak4foam을 사용하는데도 다른 조치없이 사용 할 수있는지 알고 싶습니다. 그리고, OpenFOAM extended 3.2가 기존의 OpenFOAM과 비교해서 무엇이 추가되고 장점은 무엇인지 알려주시면 감사하겠습니다.  

안녕하세요.

지금 Ubuntu 12.04에서 OpenFOAM 2.3.0을 계속 주로 사용해 왔는데, 최신 버전을 사용 할 필요가 있어 OpenFOAM 3.0.1 을 설치하여 사용해 보려고 합니다. 이런 경우 Ubuntu 12.04에서 14.04로 update 해서(지금 하드디스크를 두개의 PARTITION으로 나누어 Ubuntu와 Window을 사용하고 있습니다) OpenFOAM 3.0.1을 설치해 사용해도 문제가 없는지 (기존의 2.3.0에서 작업한 피일이나 iimplement 한 scheme등에) 알고 싶습니다. 아니면 12.04에서 Source Pack installation을 사용하는 방법이 더 나은지 알려주시면 감사하겠습니다. 

안녕하세요.

지금 제가 개발하고 있는 방법에서 하나의 과정으로 source term에 vorticity의 time step마다 계산된 값을 이용하려고 합니다.

postprocess로 vorticity를 계산하는 방법은 계산 후에  "vorticity" 라는 명령어를 이용해서 간단히 얻을 수 있는데,  

저와 같은 경우처럼 time step마다 vorticity 를 계산하기 위해서는 어떻게 해야 하는지 문의 드립니다. 지금 v2.3을 주로 이용하고 있는데  

이 버전까지는  "vorticity" directory 가 /opt/openfoam230/applications/utilities/postProcessing/velocityField/ 에 위치하고 있고 "vorticity.H" 이라는 헤더 파일도 없는데 버전 2.3 이후에는 해당 directory도 다른 곳에 위치하고 "vorticity.H" 도 존재하는 것 같습니다.

저와 같은 경우 $(FOAM_USER_APPBIN)/newApp 와 같이 새로운 application의 디렉토리에 기존의 vorticity에 관련된 전체 디렉토리를 복사해서 "file", "option"에서 설정 후 같이 compile해서 값을 얻으려고 하는데 문제가 발생하고 있습니다.

비슷한 경험을 가진 분께서는 조언 주시면 감사하겠습니다. 

안녕하십니까. 작년에 기본교육을 받은 교육생입니다.

OpenFOAM을 이용하여서 thermostat을 비정상상태로 모델링 하려고 합니다.

chtMultiRegionFoam을 이용하여서 비정상 해석을 하며,

system/solid1/fvOption을 이용하여서 heat source를 제어합니다.

(서모스텟의 region이름은 solid1입니다.)

문제는 fvOption에서 heat source를 정의할 때 제어를 어떻게 해야 될 지 모르겠습니다.

가령 (x,y,z) = (1, 1, 0)에서 온도가 600K이 넘으면 source를 끄고,

600K이하면 source를 켜는 방식으로 조절하고 싶습니다.

아래 소스를 첨부드리오니 확인 부탁드립니다.

/*--------------------------------*- C++ -*----------------------------------*

| =========                 |                                                 |

| \      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \    /   O peration     | Version:  2.3.x                                 |

|   \  /    A nd           | Web:      www.OpenFOAM.org                      |

|    \/     M anipulation  |                                                 |

*---------------------------------------------------------------------------*/

FoamFile

{

    version    2.0;

    class      dictionary;

    format     ascii;

    object     fvOptions;

}

heatSource

{

    type            scalarCodedSource;

    active          true;

    selectionMode   all;

    scalarCodedSourceCoeffs

    {

        fieldNames      (h);

        redirectType    sourceTime;

        codeInclude

        #{

        #};

        codeCorrect

        #{

            Pout<< "**codeCorrect**" << endl;

        #};

        codeAddSup

        #{

            const Time& time = mesh().time();

            const scalarField& V = mesh_.V();

            const vectorField& C = mesh_.C();

            scalarField& heSource = eqn.source();

            forAll(C, i)

            {

                heSource[i] -= 3000*V[i]; // source is 3000W/m3

            }

            Pout << "***codeAddSup***" << endl;

            Pout << P << endl;

        #};

        codeSetValue

        #{

            Pout<< "**codeSetValue**" << endl;

        #};

        // Dummy entry. Make dependent on above to trigger recompilation

        code

        #{

            $codeInclude

            $codeCorrect

            $codeAddSup

            $codeSetValue

        #};

    }

    sourceTimeCoeffs

    {

        // Dummy entry

    }

안녕하세요.

기존 오픈폼에서 제공하지 않는... 사용자 정의 wall function을 새로 만들고 싶은데, 

솔버나, 난류 라이브러리와 다르게.. 너무 막연하네요..

혹시 wall function 관련 참고자료나 소스 코드 분석(?)자료 같은 것이 존재하는지 궁금합니다ㅠ

안녕하세요. 난류모델을 implement하는 가운데 질문이 있습니다.

우선 제가 심으려는 모델은 

Lam-Bremhorst나 Abe-Nagano-Kondoh 등의 Low-Reynolds number RANS모델입니다.

난류 상수들이 다른 모델처럼 주어져있습니다. 예를들어 Cmu 0.09, C1 1.44 등...

이때 부력을 고려하는 경우,

C3 상수를 고려해줘야하는데, 

이때 C3 = tanh |v/u| 입니다. 

저는 다른 값들 처럼 lookupOr.. 로 상수값을 선언하는 방식이 아니라 

equation안에 집어넣으려고 합니다..

이때 vertical velocity와 horizontal velocity 를 어떻게 정의해야하는지 궁금합니다.

U.component(0) 이런 식으로 정의를 하니

Foam::incompressible::turbulenceModel::U’ does not have class type

안녕하세요.

소스코드를 살펴보다가 궁금한점이 있어서 질문 올립니다.

난류를 계산할때, 

예를 들어 kEpsilon 모델이라 한다면,

k equation과 e equation을 풀어서 nut 를 return하는 것으로 알고있습니다.

제가 관심이 있는 영역은 열을 포함하는 문제인데, 

그럴 경우 k와 e 식에 Gk(buoyancy production term)이 포함이 되어야 할것같은데, 

제가 찾아본 바로는 없는 것 같습니다.

일반 오픈폼 난류 모델 라이브러리나 솔버에서는 이를 어떻게 접근하고 있는지 궁금합니다.