|
[Sponsors] |
July 23, 2012, 12:04 |
non-stationary scalarTransportFoam
|
#1 |
Member
Join Date: Jul 2012
Posts: 31
Rep Power: 13 |
hi everyone,
i am trying to modify the scalarTransportFoam solver so that it can be applied to time-varying velocity fields (that i have computed with another solver). My idea was to insert the following lines from createFields.H inside the main loop of scalarTransportFoam.C: ------------------- volVectorField U ( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); # include "createPhi.H" ------------------------------- i thought this would update the U-field for every time step. The U-fields in the respective folders are, however, all changed to the initial U-field ( /0/U). Does anyone have an idea what's the problem here? Thx in advance. |
|
August 2, 2012, 08:52 |
|
#2 |
Member
Join Date: Jul 2012
Posts: 31
Rep Power: 13 |
anyone has an idea?!
|
|
August 2, 2012, 15:48 |
|
#3 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings kpax and welcome to the forum!
Perhaps you should think the other way around, namely use an existing solver as a basis for a non-stationary scalarTransportFoam. Better yet, you'll find at openfoamwiki.net the following two pages:
Best regards, Bruno
__________________
|
|
August 3, 2012, 10:12 |
|
#4 |
Member
Join Date: Jul 2012
Posts: 31
Rep Power: 13 |
hey wyldckat,
the other way around works perfectly! thx alot! |
|
October 30, 2013, 16:59 |
|
#5 |
Member
Vishal Achasrya
Join Date: Nov 2011
Posts: 38
Rep Power: 14 |
Hi,
But if icoFoam is used and a scalar transport is added to it, then both velocity and T would be solved for simultaneously coupled. What if I want to use scalarTransportFoam, but instead have a time-varying velocity field that is not solved but prescribed by the user? For example, what if my velocity is such that u=A+B*cos(t) and I want to use this time-varying field as the velocity with scalarTransportFoam, then how is this done? |
|
November 1, 2013, 19:12 |
|
#6 | ||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings Vishal,
Quote:
One example of coupled vs segregated flow: http://www.cfd-online.com/Wiki/Fluen...ated_solver.3F Quote:
Best regards, Bruno
__________________
|
|||
November 5, 2013, 12:01 |
|
#7 |
Member
Vishal Achasrya
Join Date: Nov 2011
Posts: 38
Rep Power: 14 |
Bruno,
Thanks for the reply. Indeed, for a passive scalar, the U and T are independent. However what I have is a case where my U is a function of space and time and I have a mathematical form for it such as A+B*cos(t). I wanted to hard-code this into the solver. I looked at the code and U (and phi) are defined in createFields.H before the simple.loop(). Which means that the default scalarTransport can only use a U that is defined at the beginning and kept the same throughout. I changed this by removed the block of code from "volVectorField U ( " until #include "createPhi.H" and included it in my simple.loop(). But since I am not reading or writing U to a file, I need not define U as an IOobject. Instead, I need to have a loop over all cells that calculates U as A(x,y)+B(x,y)*cos(t). If I use the following inside the simple.loop(): volVectorField U = (U0+uamp*sin(2.0*pi*f0*runTime))*vector(1, 0, 0) + (V0+vamp*sin(2.0*pi*f0*runTime))*vector(0, 1, 0); It does not work ! Here, U0, uamp, V0 and vamp are dimensionedScalars that have units m/s and f0 is a dimensionedScalar with unit 1/s, runTime is the time value I have used, pi is a const scalar. All of these values are read from the transportProperties dictionary file in the createFields.H file. |
|
November 23, 2013, 05:19 |
|
#8 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Vishal,
I've finally managed to get some time and clear mind to look into your post. The first thing that comes to mind is potentialFoam: http://openfoamwiki.net/index.php/PotentialFoam - this solver is as simple or simpler than scalarTransportFoam and allows for you to do tests of modifying "U" more easily with this solver. The second thing that comes to mind: can you please share the source code of the modified solver you're working on? This way I could test for myself and figure out what's wrong/missing. Best regards, Bruno
__________________
|
|
January 21, 2017, 19:44 |
Same kind of Problem!!!
|
#9 | |
Senior Member
Himanshu Sharma
Join Date: Jul 2012
Posts: 101
Rep Power: 13 |
Quote:
Where you able to solve your problem, and implement it in OpenFOAM. I am currently doing same kind of problem where my U =B*cos(t). Any help will be really appreciated. Thank you |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
ScalarTransportFoam for RTD calculations | santoo_cfd | OpenFOAM Running, Solving & CFD | 39 | July 12, 2021 01:15 |
Can CFX do CHT simulations with a solid domain rotating in a stationary fluid domain? | acro | CFX | 15 | September 23, 2016 11:16 |
Problem with scalarTransportFoam illistrated using pitzDaily tutorial | mlawson | OpenFOAM | 2 | January 18, 2011 13:39 |
Stationary Vs. Moving Mesh | startingcfd | Main CFD Forum | 1 | July 18, 2010 11:05 |
How to output data from stationary part only? | Aerolex | FLUENT | 0 | November 16, 2009 22:46 |