# scalarTransportFoam for arbitrary number of transported scalars

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

 June 4, 2011, 07:49 scalarTransportFoam for arbitrary number of transported scalars #1 New Member   Andreas P Join Date: Sep 2010 Posts: 26 Rep Power: 13 Hi, I am just taking my first steps programming OF, so please pardon my awkwardly foolish question. The only thing I want to do is extend scalarTransportFoam so that it can solve for an arbitrary number of transported (and independent) scalars. That means, T should not be a volScalarField, but a field of vectors of length N, and accordingly the user should provide the diffusivities also as a vector of length N. Of course it would be possible to solve the different transport equations one after another. But I don't want to do that, since I might need coupling between the different scalars some time later. So how could I achieve that? I tried to define T as a volVectorField, but that seems to always expect a vector of length 3. Is that true? I'd appreciate any help to get me on the right track. Thanks in advance! Andreas

 June 4, 2011, 11:52 #2 Senior Member   David Gaden Join Date: Apr 2009 Location: Winnipeg, Canada Posts: 437 Rep Power: 19 I don't think it would work as you hope. First of all, the vector idea is not suitable: vectors have direction, and the X vector of T will coincide with the X direction of velocity. What you are suggesting would require defining a new Type for fvScalarMatrix, which would be an involved process (a nightmare, actually... you'd have to add to all the fvc, fvm namespaces, and create all the field operators for starters...). Furthermore, I suspect that the algorithm would ultimately need to solve each component individually, so I don't think you'd be saving any computation time by defining this new Type. I wrote a scalarTransportFoam version that solves for an arbitirary number of scalars (individually), listed in a dictionary file. Let me know if you want more details about that.

June 4, 2011, 12:35
#3
New Member

Andreas P
Join Date: Sep 2010
Posts: 26
Rep Power: 13

Quote:
 Originally Posted by marupio I don't think it would work as you hope. First of all, the vector idea is not suitable: vectors have direction, and the X vector of T will coincide with the X direction of velocity.
OK, I kind of suspected that.

Quote:
 Originally Posted by marupio I wrote a scalarTransportFoam version that solves for an arbitirary number of scalars (individually), listed in a dictionary file. Let me know if you want more details about that.
Yes, I would definitely like to have more details about your implementation!

Thanks!
Andreas

June 4, 2011, 12:51
#4
Senior Member

Join Date: Apr 2009
Posts: 437
Rep Power: 19
Here's some of the code. I wrote it before I was very comfortable with OF, so it's a little awkward in places. It uses a custom class called "admScalar", which I can't seem to find anymore, but it was just a container - I'm sure you can guess what it was supposed to do.

Also, I grouped all the scalars together to use the same fvSchemes - this simplifies the fvSchemes file, but you don't have to do that. (This is what all that admConvection, admLaplacian, etc.. stuff is all about.)

This should at least give you an idea of what to do - especially the part about reading the variables in from a dictionary, and creating a pointer list of the volScalarFields.
Attached Files
 scalarTransportBits.tar.gz (794 Bytes, 41 views)

June 4, 2011, 16:36
#5
Senior Member

Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 33
Quote:
 Originally Posted by andreasp Yes, I would definitely like to have more details about your implementation!
You could also take a look at the "diesel" or "reacting" solvers. Species are solved in a sequence, and their number is arbitrary (related to components).

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541)
OpenQBMM - An open-source implementation of quadrature-based moment methods.

To obtain more accurate answers, please specify the version of OpenFOAM you are using.

 June 4, 2011, 18:34 #6 Senior Member   David Gaden Join Date: Apr 2009 Location: Winnipeg, Canada Posts: 437 Rep Power: 19 Now that you mention that, Alberto, I think that's what I used for inspiration when I wrote that code.

 June 5, 2011, 15:58 #7 New Member   Andreas P Join Date: Sep 2010 Posts: 26 Rep Power: 13 Thanks for your code snippets! I'll take a look at them and at the other foam solvers mentioned, and I'll come back if I run into problems...

 Tags extend, multiple scalars, scalartransportfoam