CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Convert denseParticleFoam to a conserved solver

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 8, 2021, 11:15
Default Convert denseParticleFoam to a conserved solver
New Member
Mingze Gao
Join Date: Sep 2020
Posts: 7
Rep Power: 3
mgao1 is on a distinguished road
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!

mgao1 is offline   Reply With Quote


dense dpm, solver, solver modification

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
fluent divergence for no reason sufjanst FLUENT 2 March 23, 2016 16:08
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
3d vof Smaras FLUENT 2 February 19, 2013 06:58
Working directory via command line Luiz CFX 4 March 6, 2011 20:02
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08

All times are GMT -4. The time now is 00:13.