|
[Sponsors] |
March 21, 2016, 09:27 |
Coupling interFoam and DPMFoam
|
#1 |
New Member
Li Linmin
Join Date: Nov 2015
Location: China
Posts: 27
Rep Power: 10 |
Hi,
I am trying to incorporate the interFoam and DPMFoam to simulate the air-water-particles systems. I want to create the viscosity field of continuous phases using the mixture law. volScalarField muc ( IOobject ( "muc", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), alpha1*rho1*"nu1" + alpha2*rho2*"nu2", alpha1.boundaryField().types() ); but I am in trouble with reading "nu" from "transportProperties" file. Anyone can help me? if anyone has similar codes for the Multiphase-DPM coupled problem and share here, it will be very helpful. Thank you very much! Linmin |
|
March 22, 2016, 07:19 |
|
#2 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Quote:
From your code, you can just use: Code:
twoPhaseProperties.mu();
__________________
My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html |
||
March 22, 2016, 07:30 |
|
#3 |
New Member
Li Linmin
Join Date: Nov 2015
Location: China
Posts: 27
Rep Power: 10 |
twoPhaseProperties.mu(); It returns the "mu" of mixture? And if it is calculated according to "alpha1*rho1*nu1+alpha2*rho2*nu2"?
|
|
March 22, 2016, 07:36 |
|
#4 | |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
Quote:
Code:
Foam::tmp<Foam::volScalarField> Foam::incompressibleTwoPhaseMixture::mu() const { const volScalarField limitedAlpha1 ( min(max(alpha1_, scalar(0)), scalar(1)) ); return tmp<volScalarField> ( new volScalarField ( "mu", limitedAlpha1*rho1_*nuModel1_->nu() + (scalar(1) - limitedAlpha1)*rho2_*nuModel2_->nu() ) ); }
__________________
My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html |
||
March 22, 2016, 07:44 |
|
#5 | |
New Member
Li Linmin
Join Date: Nov 2015
Location: China
Posts: 27
Rep Power: 10 |
Quote:
volScalarField muc ( IOobject ( "muc", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), twoPhaseProperties.mu(), alpha1.boundaryField().types() ); But it shows the error: " ‘twoPhaseProperties’ was not declared" Very sorry to interrupt again! And you documents are very helpful! Thank you very much ! |
||
March 22, 2016, 08:04 |
|
#6 |
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Beijing, China
Posts: 838
Rep Power: 17 |
hi
i will give u a more explaination after i finish my dinner.
__________________
My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html |
|
March 22, 2016, 09:12 |
|
#7 |
New Member
Li Linmin
Join Date: Nov 2015
Location: China
Posts: 27
Rep Power: 10 |
||
March 22, 2016, 22:30 |
|
#8 | |
New Member
Li Linmin
Join Date: Nov 2015
Location: China
Posts: 27
Rep Power: 10 |
Quote:
incompressibleTwoPhaseMixture twoPhaseProperties(U, phi); and used the code: twoPhaseProperties.mu(), to read the viscosity of mixture and successfully complied the code. Thank you very much again for your help and your web is very helpful. Linmin |
||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Multiphase Solid Particle Tracking | alexlupo | OpenFOAM Running, Solving & CFD | 114 | March 17, 2022 21:52 |
Issue with pressure boundary and gravity in DPMFoam | eric | OpenFOAM Running, Solving & CFD | 6 | May 8, 2019 09:28 |