|
[Sponsors] |
September 27, 2021, 12:18 |
Proper way to add a settling scalar field
|
#1 |
New Member
Roie Ezraty
Join Date: Oct 2020
Posts: 4
Rep Power: 6 |
Hello all,
I want to add to my solver a diffusion equation for the field "C" which settles in a constant velocity V. In order to assign a constant velocity to the field I go by the thread: How a scalar FIeld can be moved at constant velocity So in my case: Code:
fvScalarMatrix CEqn ( fvm::ddt(C) // changed - C + fvm::div(phiC, C) // changed - phiC & C - fvm::laplacian(kappaC, C) // changed - kappaC & C == fvOptions(C) // changed - C ); Code:
surfaceScalarField phiC ( IOobject ( "phiC", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), fvc::flux(U+V) ); My problem is this: If I define phiC, as a check, to be exactly like phi (meaning the last line only has fvc::flux(U)) I should get the same results as if I take phi instead of phiC in the diffusion equation itself - but I actually get different results! Did I miss something on the way? Thanks in advance, Roie |
|
September 28, 2021, 07:43 |
|
#2 | |
Member
Hosein
Join Date: Nov 2011
Location: Germany
Posts: 94
Rep Power: 14 |
Quote:
|
||
September 28, 2021, 08:14 |
|
#3 |
New Member
Roie Ezraty
Join Date: Oct 2020
Posts: 4
Rep Power: 6 |
Thanks a lot Hosein!
So yes, I am using a modified solver built upon buoyantBoussinesqPimpleFoam - Does buoyantBoussinesqPimpleFoam use CorrectPhi? and if so, how can I add phiC to the correction procedure as well? Thanks again, Roie |
|
Tags |
constant velocity, diffusion equation, phi, scalar field |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Foam::error::PrintStack | almir | OpenFOAM Running, Solving & CFD | 92 | May 21, 2024 07:56 |
Velocity inlet boundary condition | kit607 | OpenFOAM Pre-Processing | 1 | November 20, 2019 05:04 |
Steady Advection-diffusion equation, with scalar velocity field, openfoam | andrea5 | OpenFOAM Programming & Development | 0 | September 17, 2019 04:10 |
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' | muth | OpenFOAM Running, Solving & CFD | 3 | August 27, 2018 04:18 |
Issue symmetryPlane 2.5d extruded airfoil simulation | 281419 | OpenFOAM Running, Solving & CFD | 5 | November 28, 2015 13:09 |