CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

blockCoupledScalarTransport solver

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

Like Tree1Likes
  • 1 Post By pweber

Reply
 
LinkBack Thread Tools Display Modes
Old   March 22, 2014, 15:54
Exclamation blockCoupledScalarTransport solver
  #1
New Member
 
Join Date: Mar 2014
Posts: 20
Rep Power: 3
alisina-s is on a distinguished road
Hi everybody,
I am trying to solve 2 transport equations that are coupled each other and I modified blockCoupledScalarTransport.C in the way of my equations and the same as mentioned here
http://sourceforge.net/p/openfoam-ex...TransportFoam/
I added all this files with makefile in my directory but when I want to compile my solver(the modified version of blockCoupledScalarTransport) I have this:
In file included from tensor2Field.H:40:0,
from tensor2Field.C:27:
tensor2.H:39:21: fatal error: TensorN.H: No such file or directory
compilation terminated.
Does anybody knows what is this and what should I do?
thanks in advance
alisina-s is offline   Reply With Quote

Old   March 22, 2014, 18:23
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,308
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings alisina-s and welcome to the forum!

I'm not sure I understand your question. Are you trying to build this solver with or without OpenFOAM?

Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   March 22, 2014, 18:39
Default
  #3
New Member
 
Join Date: Mar 2014
Posts: 20
Rep Power: 3
alisina-s is on a distinguished road
thanks,
I am trying to use this solver by a little bit modification int the way of my equations in openfoam.
alisina-s is offline   Reply With Quote

Old   March 22, 2014, 18:53
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,308
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by alisina-s View Post
thanks,
I am trying to use this solver by a little bit modification int the way of my equations in openfoam.
OK, then there are two important details that you should keep in mind:
  1. When creating your own customized solver, follow the same strategy outlined on this tutorial: http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam
  2. This particular solver you're using, namely "blockCoupledScalarTransport", should be used only with FOAM-Extend 3.0.
    In other words, do not try to use it in OpenFOAM 2.3 or 2.2 or older, as it is not compatible. The reason is simple: equation coupling in FOAM-Extend is done differently from OpenFOAM. More specifically, I'm not sure if OpenFOAM even has this feature.
wyldckat is online now   Reply With Quote

Old   March 23, 2014, 06:11
Default
  #5
New Member
 
Join Date: Mar 2014
Posts: 20
Rep Power: 3
alisina-s is on a distinguished road
thanks bruno for your answer
I knew about how to make our own solver but I did not know that I have to use this specific version of openfoam you mentioned.my equations are like below
A) dT/dt +d/dx(UT)=0
B) dT/dt +d/dx(Uk)=0
dk/dt +d/dx(UT)=0
and k is another quantity that want to be transfered like T.
I know these are not notation of openfoam in writing the equations I just want to mention my equations.
do you have any other suggestion for this except using "blockCoupledScalarTransport"?

regards,
sina
alisina-s is offline   Reply With Quote

Old   March 23, 2014, 13:28
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,308
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Unfortunately, I don't know the answer.
But I do have the vague idea that this has been asked here in the forum a few times in the past...

Good luck! Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   March 24, 2014, 06:45
Default
  #7
New Member
 
Paul Weber
Join Date: Mar 2014
Posts: 3
Rep Power: 3
pweber is on a distinguished road
Hello sina,

I think you do not need a block coupled solver since your equations decouple.

Quote:
Originally Posted by alisina-s View Post
my equations are like below
A) dT/dt +d/dx(UT)=0
B) dT/dt +d/dx(Uk)=0
C) dk/dt +d/dx(UT)=0
Eq. A) just uses T, eq. C) becomes

Code:
0 = dk/dt + d/dx(UT)
  = dk/dt - dT/dt
  = dk/dt + d/dx(Uk)
by first using eq. A) and then eq. B). Thus you have

Code:
dT/dt + d/dx(UT) = 0
dk/dt + d/dx(Uk) = 0
Best,
Paul
wyldckat likes this.
pweber is offline   Reply With Quote

Old   March 27, 2014, 12:47
Default
  #8
New Member
 
Join Date: Mar 2014
Posts: 20
Rep Power: 3
alisina-s is on a distinguished road
Quote:
Originally Posted by pweber View Post
Hello sina,

I think you do not need a block coupled solver since your equations decouple.



Eq. A) just uses T, eq. C) becomes

Code:
0 = dk/dt + d/dx(UT)
  = dk/dt - dT/dt
  = dk/dt + d/dx(Uk)
by first using eq. A) and then eq. B). Thus you have

Code:
dT/dt + d/dx(UT) = 0
dk/dt + d/dx(Uk) = 0
Best,
Paul
thanks paul
but what if I wanted to solve these two quantity together,I mean if I want k be under effect of Quantity T,Do you have any idea?
alisina-s is offline   Reply With Quote

Old   April 3, 2014, 07:44
Default
  #9
New Member
 
Paul Weber
Join Date: Mar 2014
Posts: 3
Rep Power: 3
pweber is on a distinguished road
Dear Sina,

the nature of your equations is such that T and k decouple, i.e. they have no influence on each other. Reformulations of the system and it's simulation will never show dependence of T and k on each other, simply because it does not exist. However, if the underlying physics modeled by your system indeed imply a coupling, the only way to get it included is to redo the modeling in a proper way and thus derive a new, then coupled, equation system.

In that case, the coupled system can be written as
Code:
dT/dt + d/dx(UT) = f(k)
dk/dt + d/dx(Uk) = g(T)
where f and g are functions (in general, operators) depending on k and T, respectively. Having that, proceed to define non-zero parts of f and g which you would like to be treated implicitly. Then you are prepared to make use of block coupled simulation techniques, but as long as you use the model given to us here, you will never see any coupling of the two variables.

Best,
Paul
pweber is offline   Reply With Quote

Old   April 3, 2014, 11:12
Default
  #10
New Member
 
Join Date: Mar 2014
Posts: 20
Rep Power: 3
alisina-s is on a distinguished road
Quote:
Originally Posted by pweber View Post
Dear Sina,

the nature of your equations is such that T and k decouple, i.e. they have no influence on each other. Reformulations of the system and it's simulation will never show dependence of T and k on each other, simply because it does not exist. However, if the underlying physics modeled by your system indeed imply a coupling, the only way to get it included is to redo the modeling in a proper way and thus derive a new, then coupled, equation system.

In that case, the coupled system can be written as
Code:
dT/dt + d/dx(UT) = f(k)
dk/dt + d/dx(Uk) = g(T)
where f and g are functions (in general, operators) depending on k and T, respectively. Having that, proceed to define non-zero parts of f and g which you would like to be treated implicitly. Then you are prepared to make use of block coupled simulation techniques, but as long as you use the model given to us here, you will never see any coupling of the two variables.

Best,
Paul
let consider I have just
dT/dt + d/dx(Uk) = 0
dk/dt + d/dx(UT) = 0
how can I solve this?
I was thinking to define a vector a(T k2) and a'(k -T) and then solve da/dt + da'/dx = 0
do you know how can I implement this in openfoam?
Thanks in advance,
Sina
alisina-s is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Fluent Adjoint Solver? ex10148 FLUENT 15 July 21, 2015 01:13
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
Divergence problem Smaras FLUENT 13 February 21, 2013 06:03
3d vof Smaras FLUENT 2 February 19, 2013 07:58
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 11:48.