|
[Sponsors] |
New phaseScalarTransport function object for a single phase in interFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 14, 2021, 11:28 |
New phaseScalarTransport function object for a single phase in interFoam
|
#1 |
New Member
rasoul
Join Date: Mar 2015
Posts: 24
Rep Power: 11 |
Hi everyone!
As you know OpenFOAM announced a new function object for solving scalar transport in only one phase of multi-phase flow solvers especially interFoam. (https://twitter.com/cfddirect/status...167872?lang=en) I believe here is the best place to discuss about phaseScalarTransport features, advantages and limitations. when OpenFOAM applies this part of code to the algorithm, basically solute (e.g. s.water) diffuses into the second phase but this algorithm will multiply it to the alpha.water (volume fraction in interFoam solver) and then solve the next time step. Code:
fvScalarMatrix fieldEqn ( fvm::ddt(alpha, s_) + fvm::div(alphaPhi, s_, divScheme) - fvm::laplacian ( fvc::interpolate(alpha)*fvc::interpolate(D), s_, laplacianScheme ) == fvOptions_(alpha, s_) - fvm::ddt(residualAlpha_, s_) + fvc::ddt(residualAlpha_, s_) ); and when I tried a simple test case: (interFoam + phaseTransportFoam) initial conditions of alpha.water and s.water: (link below) https://drive.google.com/file/d/1Jth...ew?usp=sharing indeed when solute spreads to the blue phase it propagates and gives artificial diffusion https://drive.google.com/file/d/1ja9...ew?usp=sharing so the second question is: how can we prevent s.water to diffuse to the blue phase causing this artificial diffusion. (maybe related to the first question) actually with multiplying s.water to alpha it reduces the whole initial amount of solute in the main phase. any explanation or help would be highly appreciated! Rasoul |
|
March 5, 2021, 14:53 |
|
#2 |
Senior Member
Reviewer #2
Join Date: Jul 2015
Location: Knoxville, TN
Posts: 141
Rep Power: 11 |
Rasoul,
I think I have the answer to your questions in this post. interFoam with scalar transport in single phase Thanks, Rdf |
|
May 26, 2021, 07:31 |
|
#3 | |
New Member
Sudarshan
Join Date: Sep 2020
Posts: 3
Rep Power: 6 |
Hi,
I'm also trying to simulate a case (similar to sand in water + air) however I do not understand how to implement the phaseTransportFoam (in particular could you please explain what you meant by interfoam + phaseTransportFoam). I tried to implement this in the dambreak case. In the controlDict, file I included the following lines before the end. phaseScalarTransport1 { type phaseScalarTransport; libs ("libsolverFunctionObjects.so"); field s.water; p p_rgh; } phaseScalarTransport2 { type phaseScalarTransport; libs ("libsolverFunctionObjects.so"); field alpha.water * s.water; p p_rgh; } Then I run the command "interfoam" which ends well but when I enter the command "phaseTransportFoam" it is not found. Could you please tell the steps clearly on how to implement phaseTransportFoam. Thanks. Sudarshan. Quote:
|
||
Tags |
interfoam, only one phase, phasetransportfoam, solute transport |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Applying strain rate function to fluid phase only | pquinn | OpenFOAM Post-Processing | 0 | March 6, 2020 20:03 |
[snappyHexMesh] Problem with parallel run of snappyHexMesh | Lorenzo92 | OpenFOAM Meshing & Mesh Conversion | 5 | April 15, 2016 04:12 |
[snappyHexMesh] How to define to right point for locationInMesh | Mirage12 | OpenFOAM Meshing & Mesh Conversion | 7 | March 13, 2016 14:07 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 07:36 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 14:00 |