|
[Sponsors] |
[swak4Foam] Efficient run-time postprocessing: combination of swak4foam and solver modification? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 21, 2013, 12:16 |
Efficient run-time postprocessing: combination of swak4foam and solver modification?
|
#1 |
Member
Marcus Letzel
Join Date: Sep 2012
Location: Bremen
Posts: 35
Rep Power: 13 |
Dear Foamers,
similar to this thread on how to calculate a turbulent heat flux <u'T'>, where <...> denotes a local time average, I plan to calculate the vertical turbulent momentum fluxes <w'u'> and <w'v'>. Plus, the vertical total momentum fluxes <wu> and <wv>. This is in order to monitor OpenFOAM LES of the atmospheric boundary layer. I want to output these data not for the total domain but only as one-dimensional vertical profiles at say ten or twenty locations within the simulation domain, e.g. for postprocessing with XmGrace. My question is which is the best way to efficiently calculate and output these data. Since I have a modified solver anyway, it would be possible to code the time averaging according to this post for U as well as the products wu and wv. The time-averages <w'u'> and <w'v'> could be calculated in a similar manner in the solver. But I would expect that this wastes CPU time as well as RAM since the output is actually only needed for ten to twenty one-dimensional vertical profiles. So, is there a way to calculate these data more efficiently only where output is needed? I have a vague feeling that swak4foam may accomplish this. If so, how? And is there any elegant way to combine the definition of ten, twenty vertical profiles? If I use ten, twenty different sampledSet, my controlDict would become quite long. Looking forward to hear your advices, Marcus |
|
October 21, 2013, 19:21 |
|
#2 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
October 22, 2013, 06:57 |
|
#3 | |
Member
Marcus Letzel
Join Date: Sep 2012
Location: Bremen
Posts: 35
Rep Power: 13 |
Quote:
Is there any chance that the compulsory condensation of swakExpression results to single values could become optional at some point in the future? I must admit that I am a bit hesitant regarding the dumpSwakExpression approach. IMHO, if the columns in the output file are not labelled, it seems to be prone to human errors during further postprocessing. For example, columns may get mixed up later. So, I am now considering to resort to the computationally less efficient combination of solver modifications (averaging and flux calculation) and the native OpenFOAM sampling functionality by declaring ten, twenty sampledSets. |
||
October 22, 2013, 10:06 |
|
#4 | ||
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
Also having the "dump" as the default functionality is asking for trouble. Approx 2.7 days after release somebody will use it on an internalField, write it every timestep and complain that swak filled up his disk Quote:
One of the reasons why there are no labels with dump is a) that this must be implemented for every valueType separately b) What should the label be? Position? Cell number? c) the hardest part: what should be done if the entity changes its size, ordering? (for instance due to mesh motion or a sampledSurface of type isoSurface) Start a new file? Just write a new heading?
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] How to run the case in swak4Foam | ce73stargazer | OpenFOAM Community Contributions | 7 | December 19, 2017 11:38 |
InterDyMFoam+simpleFunctionObject | Elham | OpenFOAM Running, Solving & CFD | 5 | July 10, 2017 11:59 |
[swak4Foam] and snappyHexMesh gives erratic mesh result | Coke Rivas Ordenes | OpenFOAM Community Contributions | 3 | April 17, 2017 13:06 |
[swak4Foam] How to run pulsedPitzdaily example in swak4Foam Example ? (oodles solver) | nanavati | OpenFOAM Community Contributions | 1 | January 27, 2015 04:06 |