煙の拡散について
560 views
Skip to first unread message
loid
Jun 21, 2010, 12:43:48 AM6/21/10
to OpenFOAM
loidと申します。
OpenFOAMは、チュートリアルを試した程度の初心者です。
煙や匂いの拡散をシミュレーションしてみたいと思っています。
手始めに、線香やタバコのような煙の拡散をシミュレーションしたいのですがどなたか、
参考になるサンプルプログラムやチュートリアルをご紹介いただけると幸です。
OpenFOAMは、チュートリアルを試した程度の初心者です。
煙や匂いの拡散をシミュレーションしてみたいと思っています。
手始めに、線香やタバコのような煙の拡散をシミュレーションしたいのですがどなたか、
参考になるサンプルプログラムやチュートリアルをご紹介いただけると幸です。
koba
Jul 11, 2010, 11:23:30 PM7/11/10
to OpenFOAM
kobaと申します。
非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
ピタッと当てはまるものがないと思います。
最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
濃度拡散は可能と思います。
具体的には、chemkinフォルダー内にあるchem.inpファイルの
AL => BL 6.80E+06 0.0 33500.0! 1
の頻度因子の6.8e6を０にすれば、反応しません。
上記はver1.6の場合です。 ver1.7はまだよく見ていません。
詳細は下記などを参考にしてください。 サンプルファイルもあります。
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/reactingFoam.pdf
非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
ピタッと当てはまるものがないと思います。
最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
濃度拡散は可能と思います。
具体的には、chemkinフォルダー内にあるchem.inpファイルの
AL => BL 6.80E+06 0.0 33500.0! 1
の頻度因子の6.8e6を０にすれば、反応しません。
上記はver1.6の場合です。 ver1.7はまだよく見ていません。
詳細は下記などを参考にしてください。 サンプルファイルもあります。
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/reactingFoam.pdf
loid
Jul 13, 2010, 9:04:28 AM7/13/10
to OpenFOAM
koba様
loidと申します。
対応できそうなソルバーや具体的な対応法方、参考文献までご教授いただきありがとうございます。
OpenFOAM初心者でしかも反応系の知識も乏しいのでうまくいくかわかりませんが、ご教授いただ
いたreactingFOAMで試してみようと思います。
また、ご相談させていただきと思いますのでよろしくお願いいたします。
On 7月12日, 午後12:23, koba <koba100...@gmail.com> wrote:
> kobaと申します。
> 非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
> ピタッと当てはまるものがないと思います。
> 最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
> 濃度拡散は可能と思います。
> 具体的には、chemkinフォルダー内にあるchem.inpファイルの
> AL => BL 6.80E+06 0.0 33500.0! 1
> の頻度因子の6.8e6を０にすれば、反応しません。
> 上記はver1.6の場合です。 ver1.7はまだよく見ていません。
> 詳細は下記などを参考にしてください。 サンプルファイルもあります。http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/...
loidと申します。
対応できそうなソルバーや具体的な対応法方、参考文献までご教授いただきありがとうございます。
OpenFOAM初心者でしかも反応系の知識も乏しいのでうまくいくかわかりませんが、ご教授いただ
いたreactingFOAMで試してみようと思います。
また、ご相談させていただきと思いますのでよろしくお願いいたします。
On 7月12日, 午後12:23, koba <koba100...@gmail.com> wrote:
> kobaと申します。
> 非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
> ピタッと当てはまるものがないと思います。
> 最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
> 濃度拡散は可能と思います。
> 具体的には、chemkinフォルダー内にあるchem.inpファイルの
> AL => BL 6.80E+06 0.0 33500.0! 1
> の頻度因子の6.8e6を０にすれば、反応しません。
> 上記はver1.6の場合です。 ver1.7はまだよく見ていません。
loid lisp
Jul 16, 2010, 9:43:27 AM7/16/10
to OpenFOAM
koba様
loidです。
koba様にご紹介いただいた参考文献とサンプルファイルをダウンロードし試してみました。
サンプルファイルの一部（chemistryProperties とcombustionProperties)にミスがありましたの修正して実行してみました。
blockMeshはうまくいきましたが、reactingFoamを実行すると下記のようなエラーで止まってしまいます。
アドバイスを頂けますと幸です。
loid@dolphin:~/OpenFOAM/loid-1.6/run/reactingFoamTutorial$ reactingFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : reactingFoam
Date : Jul 16 2010
Time : 22:20:18
Host : dolphin
PID : 7803
Case : /home/loid/OpenFOAM/loid-1.6/run/reactingFoamTutorial
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
// using new solver syntax:
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 0;
}
// using new solver syntax:
Yi
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
Reading chemistry properties
Reading g
cannot open file
file: /home/loid/OpenFOAM/loid-1.6/run/reactingFoamTutorial/constant/g at line 0.
From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 62.
FOAM exiting
loidです。
koba様にご紹介いただいた参考文献とサンプルファイルをダウンロードし試してみました。
サンプルファイルの一部（chemistryProperties とcombustionProperties)にミスがありましたの修正して実行してみました。
blockMeshはうまくいきましたが、reactingFoamを実行すると下記のようなエラーで止まってしまいます。
アドバイスを頂けますと幸です。
loid@dolphin:~/OpenFOAM/loid-1.6/run/reactingFoamTutorial$ reactingFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : reactingFoam
Date : Jul 16 2010
Time : 22:20:18
Host : dolphin
PID : 7803
Case : /home/loid/OpenFOAM/loid-1.6/run/reactingFoamTutorial
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
// using new solver syntax:
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-09;
relTol 0;
}
// using new solver syntax:
Yi
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
// using new solver syntax:
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}
Reading chemistry properties
Reading g
cannot open file
file: /home/loid/OpenFOAM/loid-1.6/run/reactingFoamTutorial/constant/g at line 0.
From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 62.
FOAM exiting
2010年7月13日22:04 loid <lisp...@gmail.com>:
--
このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
このグループから退会するには、openfoam+u...@googlegroups.com にメールを送信してください。
詳細については、http://groups.google.com/group/openfoam?hl=ja からこのグループにアクセスしてください。
koba
Jul 17, 2010, 4:33:56 AM7/17/10
to OpenFOAM
loid様
kobaです。
サンプルはver1.4だと思います。 ver1.6とフォーマットが変わっています。
このことを先にお伝えしなければいけなかったですが、ご容赦ください。
いろいろと変更されたと思いますが、
system/fvSolutionsのsolversは画面で表示された
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
：
：
のとおりに書き直してください。
そのほかに、
重力加速度は、environmentalProperties -> ｇ
乱流モデルは、turbulenceProperties -> turbulenceProperties RASproperties
濃度拡散のみの計算の場合は
chemistryPropertiesは
psiChemistryModel ODEChemistryModel<gasThermoPhysics>;
chemistry off;
thermophysicalPropertiesは
thermoType hPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
CHEMKINFile "$FOAM_CASE/chemkin/chem.inp";
CHEMKINThermoFile "$FOAM_CASE/chemkin/therm.dat";
のように変更しました。
On 7月16日, 午後10:43, loid lisp <lispl...@gmail.com> wrote:
> koba様
>
> loidです。
>
> koba様にご紹介いただいた参考文献とサンプルファイルをダウンロードし試してみました。
> サンプルファイルの一部（chemistryProperties
> とcombustionProperties)にミスがありましたの修正して実行してみました。
> blockMeshはうまくいきましたが、reactingFoamを実行すると下記のようなエラーで止まってしまいます。
> アドバイスを頂けますと幸です。
>
> loid@dolphin:~/OpenFOAM/loid-1.6/run/reactingFoamTutorial$ reactingFoam
> /*----------------------------------------------------------------------------*\
> 2010年7月13日22:04 loid <lispl...@gmail.com>:
> > ..
> > 詳細については、http://groups.google.com/group/openfoam?hl=jaからこのグループにアクセスしてください。- 引用テキストを表示しない -
>
> - 引用テキストを表示 -
kobaです。
サンプルはver1.4だと思います。 ver1.6とフォーマットが変わっています。
このことを先にお伝えしなければいけなかったですが、ご容赦ください。
いろいろと変更されたと思いますが、
system/fvSolutionsのsolversは画面で表示された
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}
：
のとおりに書き直してください。
そのほかに、
重力加速度は、environmentalProperties -> ｇ
乱流モデルは、turbulenceProperties -> turbulenceProperties RASproperties
濃度拡散のみの計算の場合は
chemistryPropertiesは
psiChemistryModel ODEChemistryModel<gasThermoPhysics>;
chemistry off;
thermophysicalPropertiesは
thermoType hPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
CHEMKINFile "$FOAM_CASE/chemkin/chem.inp";
CHEMKINThermoFile "$FOAM_CASE/chemkin/therm.dat";
のように変更しました。
On 7月16日, 午後10:43, loid lisp <lispl...@gmail.com> wrote:
> koba様
>
> loidです。
>
> koba様にご紹介いただいた参考文献とサンプルファイルをダウンロードし試してみました。
> サンプルファイルの一部（chemistryProperties
> とcombustionProperties)にミスがありましたの修正して実行してみました。
> blockMeshはうまくいきましたが、reactingFoamを実行すると下記のようなエラーで止まってしまいます。
> アドバイスを頂けますと幸です。
>
> loid@dolphin:~/OpenFOAM/loid-1.6/run/reactingFoamTutorial$ reactingFoam
> | =========
> | |
> | \\ / F ield | OpenFOAM: The Open Source CFD
> Toolbox |
> | \\ / O peration | Version:
> 1.6 |
> | \\ / A nd | Web: www.OpenFOAM.org
> |
> | \\/ M anipulation
> | |
> \*----------------------------------------------------------------------------*/
> | |
> | \\ / F ield | OpenFOAM: The Open Source CFD
> Toolbox |
> | \\ / O peration | Version:
> 1.6 |
> | \\ / A nd | Web: www.OpenFOAM.org
> |
> | \\/ M anipulation
> | |
>
>
>
> > koba様
>
> > loidと申します。
>
> > 対応できそうなソルバーや具体的な対応法方、参考文献までご教授いただきありがとうございます。
> > OpenFOAM初心者でしかも反応系の知識も乏しいのでうまくいくかわかりませんが、ご教授いただ
> > いたreactingFOAMで試してみようと思います。
> > また、ご相談させていただきと思いますのでよろしくお願いいたします。
>
> > On 7月12日, 午後12:23, koba <koba100...@gmail.com> wrote:
> > > kobaと申します。
> > > 非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
> > > ピタッと当てはまるものがないと思います。
> > > 最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
> > > 濃度拡散は可能と思います。
> > > 具体的には、chemkinフォルダー内にあるchem.inpファイルの
> > > AL => BL 6.80E+06 0.0 33500.0! 1
> > > の頻度因子の6.8e6を０にすれば、反応しません。
> > > 上記はver1.6の場合です。 ver1.7はまだよく見ていません。
> > > 詳細は下記などを参考にしてください。 サンプルファイルもあります。
> >http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2007/AndreasLundstrom/.<http://www.tfd.chalmers.se/%7Ehani/kurser/OS_CFD_2007/AndreasLundstrom/.>
>
>
> > koba様
>
> > loidと申します。
>
> > 対応できそうなソルバーや具体的な対応法方、参考文献までご教授いただきありがとうございます。
> > OpenFOAM初心者でしかも反応系の知識も乏しいのでうまくいくかわかりませんが、ご教授いただ
> > いたreactingFOAMで試してみようと思います。
> > また、ご相談させていただきと思いますのでよろしくお願いいたします。
>
> > On 7月12日, 午後12:23, koba <koba100...@gmail.com> wrote:
> > > kobaと申します。
> > > 非圧縮性の濃度拡散問題になると思いますが、OpenFOAMに用意されているソルバーには、
> > > ピタッと当てはまるものがないと思います。
> > > 最も近いと思われるのはreactingFOAMで、反応させないか、もしくは着火しなければ
> > > 濃度拡散は可能と思います。
> > > 具体的には、chemkinフォルダー内にあるchem.inpファイルの
> > > AL => BL 6.80E+06 0.0 33500.0! 1
> > > の頻度因子の6.8e6を０にすれば、反応しません。
> > > 上記はver1.6の場合です。 ver1.7はまだよく見ていません。
> > > 詳細は下記などを参考にしてください。 サンプルファイルもあります。
> > ..
>
> > > On 6月21日, 午後1:43, loid <lispl...@gmail.com> wrote:
>
> > > > loidと申します。
> > > > OpenFOAMは、チュートリアルを試した程度の初心者です。
> > > > 煙や匂いの拡散をシミュレーションしてみたいと思っています。
> > > > 手始めに、線香やタバコのような煙の拡散をシミュレーションしたいのですがどなたか、
> > > > 参考になるサンプルプログラムやチュートリアルをご紹介いただけると幸です。
>
> > --
> > このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> > このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> > このグループから退会するには、openfoam+u...@googlegroups.com<openfoam%2Bunsub...@googlegroups.com>にメールを送信してください。
> > > On 6月21日, 午後1:43, loid <lispl...@gmail.com> wrote:
>
> > > > loidと申します。
> > > > OpenFOAMは、チュートリアルを試した程度の初心者です。
> > > > 煙や匂いの拡散をシミュレーションしてみたいと思っています。
> > > > 手始めに、線香やタバコのような煙の拡散をシミュレーションしたいのですがどなたか、
> > > > 参考になるサンプルプログラムやチュートリアルをご紹介いただけると幸です。
>
> > --
> > このメールは Google グループのグループ「OpenFOAM」の登録者に送られています。
> > このグループに投稿するには、open...@googlegroups.com にメールを送信してください。
> > 詳細については、http://groups.google.com/group/openfoam?hl=jaからこのグループにアクセスしてください。- 引用テキストを表示しない -
>
> - 引用テキストを表示 -
loid lisp
Jul 17, 2010, 9:03:41 AM7/17/10
to open...@googlegroups.com
koba様
loidです。
お忙しいにもかかわらず早速お返事をいただきありがとうございます。
詳細なご指導をいただき感謝申し上げます。早速試してみたいと思います。
loidです。
お忙しいにもかかわらず早速お返事をいただきありがとうございます。
詳細なご指導をいただき感謝申し上げます。早速試してみたいと思います。
2010年7月17日17:33 koba <koba1...@gmail.com>:
Reply all
Reply to author
Forward
0 new messages