|
[Sponsors] |
Convert denseParticleFoam to a conserved solver |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 8, 2021, 12:15 |
Convert denseParticleFoam to a conserved solver
|
#1 |
New Member
Mingze Gao
Join Date: Sep 2020
Posts: 7
Rep Power: 6 |
Hi there,
I am recently modifying the denseParticleFoam (DPFoam) and I want to convert it to a conserved solver. The original code is: 1 fvVectorMatrix UcEqn 2 ( 3 fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc) 4 - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) 5 + continuousPhaseTurbulence->divDevTau(Uc) 6 == 7 (1.0/rhoc)*cloudSU 8 ); In my opinion, the 4th line will make the momentum equation a non-consered form. Therefore, when I modified the solver, I think I should remove this line, and the code will look like: alphaRhoc = rhoc*alphac; alphaRhocf = fvc::interpolate(alphaRhoc); alphaRhoPhic = alphaRhocf*phic; 9 fvm::ddt(alphaRhoc, Uc) + fvm::div(alphaRhoPhic, Uc) 10 + turbulence->divDevTau(alphaRhoc,Uc) 11 == 12 cloudSU However, this will return me different results in all the fields (such as velocity, pressure, etc...) due to the 4th will not be a zero filed. My question is, should I remove the 4th line of code to make it conserved form? If yes, is there any other code I should add to get a similar result? Any help will be much appreciated! Regards, Mingze |
|
Tags |
dense dpm, solver, solver modification |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
fluent divergence for no reason | sufjanst | FLUENT | 2 | March 23, 2016 17:08 |
thobois class engineTopoChangerMesh error | Peter_600 | OpenFOAM | 4 | August 2, 2014 10:52 |
3d vof | Smaras | FLUENT | 2 | February 19, 2013 07:58 |
Working directory via command line | Luiz | CFX | 4 | March 6, 2011 21:02 |
why the solver reject it? Anyone with experience? | bearcat | CFX | 6 | April 28, 2008 15:08 |