scalarTransportFoam does not reach 100% of the input scalar

RobertHB · December 9, 2016, 07:03

Hello Foamers,

i'm running a scalar Transport case in OpenFOAM 3.0.1 (Ubuntu 14.04) aiming to do an analysis of the residence time distribution. Since my input velocity field comes from a turbulent pimpleFoam simulation i followed the advice given here by D. Combest to add the necessary support for turbulent diffusivity and the additional term to ensure boundedness. I'm using the fvSchemes and fvSolution entries from the tutorial case. I'm running the mapFields utility to map the simulated velocity field onto my transport case. The simulation runs fine, without any errors.

I'm aware of the problem with the field phi calculation in scalarTransportFoam but this should be avoided by the addition suggested by D. Combest. mapFields sadly does not map my original phi field to the transport case but writes a phi.unmapped file.

But the concentration at my outlet does not reach 100% of the input and caps at 0.51.

Can anyone help me to fix this problem?

Regards,

Robert

/edit:
fvSchemes
fvSolution

RobertHB · December 12, 2016, 06:52

I tested the standard scalarTransportFoam as my custom TransportFoam with the pitzDaily tutorial case. For the test case both solvers do get anywhere close to 1.0 reaching the highest concentration at the outlet with 2 - 3e-5 of my input scalar. For my own case both solvers reach up to 0.5 of my input scalar.

Do you have any idea how this problem can be fixed?

RobertHB · December 12, 2016, 10:51

Ran the transport simulation with the above mentioned fvSchemes and fvSolution files on a case without a stl file as solid boundary. The solver does reach 100% of the input scalar under these circumstances. But as soon as the same model is applied to a case with a stl boundary the solution does not reach the 100%.

RobertHB · December 13, 2016, 04:28

Hello together,
i ran a second model with a different stl topography yesterday. The code runs fine, but as soon as i add a stl file instead of a flat wall BC the solution does not reach 1.0. It instead converges at a value of 5.0 this time.

Igor141 · February 1, 2021, 09:53

Hi RobertHB,

Did you find a solution to your problem ? I'm experiencing the same problem with the interFoam solver on OF6 by using the scalarTransport function in my controlDict file (my scalar reach 60% of the total concentration in a mixing model, although there is no diffusion/source or well term in my transport equation). Also if possible, could you please share again your fvSolutions/fvSchemes/controlDict files ? Maybe my problem is as simple as an error in these files.

Thanks

December 9, 2016, 07:03	scalarTransportFoam does not reach 100% of the input scalar	#1
RobertHB Senior Member Robert Join Date: May 2015 Location: Bremen, GER Posts: 292 Rep Power: 12	Hello Foamers, i'm running a scalar Transport case in OpenFOAM 3.0.1 (Ubuntu 14.04) aiming to do an analysis of the residence time distribution. Since my input velocity field comes from a turbulent pimpleFoam simulation i followed the advice given here by D. Combest to add the necessary support for turbulent diffusivity and the additional term to ensure boundedness. I'm using the fvSchemes and fvSolution entries from the tutorial case. I'm running the mapFields utility to map the simulated velocity field onto my transport case. The simulation runs fine, without any errors. I'm aware of the problem with the field phi calculation in scalarTransportFoam but this should be avoided by the addition suggested by D. Combest. mapFields sadly does not map my original phi field to the transport case but writes a phi.unmapped file. But the concentration at my outlet does not reach 100% of the input and caps at 0.51. Can anyone help me to fix this problem? Regards, Robert /edit: fvSchemes fvSolution Last edited by RobertHB; December 9, 2016 at 08:53. Reason: added fvSchemes and fvSolution files

December 12, 2016, 10:51	Successfully ran the simulation.	#3
RobertHB Senior Member Robert Join Date: May 2015 Location: Bremen, GER Posts: 292 Rep Power: 12	Ran the transport simulation with the above mentioned fvSchemes and fvSolution files on a case without a stl file as solid boundary. The solver does reach 100% of the input scalar under these circumstances. But as soon as the same model is applied to a case with a stl boundary the solution does not reach the 100%.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
problem during mpi in server: expected Scalar, found on line 0 the word 'nan'	muth	OpenFOAM Running, Solving & CFD	3	August 27, 2018 05:18
Issue symmetryPlane 2.5d extruded airfoil simulation	281419	OpenFOAM Running, Solving & CFD	5	November 28, 2015 14:09
Diverging solution in transonicMRFDyMFoam	tsalter	OpenFOAM Running, Solving & CFD	30	July 7, 2014 07:20
compressible flow in turbocharger	riesotto	OpenFOAM	50	May 26, 2014 02:47
[OpenFOAM] Paraview 3.98 - errors when saving geometry file	pajot	ParaView	1	September 28, 2013 11:45

December 12, 2016, 06:52		#2
RobertHB Senior Member Robert Join Date: May 2015 Location: Bremen, GER Posts: 292 Rep Power: 12	I tested the standard scalarTransportFoam as my custom TransportFoam with the pitzDaily tutorial case. For the test case both solvers do get anywhere close to 1.0 reaching the highest concentration at the outlet with 2 - 3e-5 of my input scalar. For my own case both solvers reach up to 0.5 of my input scalar. Do you have any idea how this problem can be fixed?

December 13, 2016, 04:28		#4
RobertHB Senior Member Robert Join Date: May 2015 Location: Bremen, GER Posts: 292 Rep Power: 12	Hello together, i ran a second model with a different stl topography yesterday. The code runs fine, but as soon as i add a stl file instead of a flat wall BC the solution does not reach 1.0. It instead converges at a value of 5.0 this time.

February 1, 2021, 09:53		#5
Igor141 New Member ap141 Join Date: Sep 2015 Posts: 2 Rep Power: 0	Hi RobertHB, Did you find a solution to your problem ? I'm experiencing the same problem with the interFoam solver on OF6 by using the scalarTransport function in my controlDict file (my scalar reach 60% of the total concentration in a mixing model, although there is no diffusion/source or well term in my transport equation). Also if possible, could you please share again your fvSolutions/fvSchemes/controlDict files ? Maybe my problem is as simple as an error in these files. Thanks