[gschaider]

Originally Posted by tomf

Hi Bernhard,

Thank you for the clarification. That will save me some trial and error. Maybe I can find another solution.

Regards,
Tom

I haven't tried that (but would be interested if you succeed): use the sample-utility to write out a sampledSurface type sampledPatch and read that with timeVaryingMappedFixedValue (you'll have to find out which writer-format is working)

[tomf]

Thanks for the suggestion, did not yet consider that option, but it sounds promising. I'll keep you updated once I find a solution.

[atoof]

Originally Posted by gschaider

- Plugin functions to access fields from the radiation mode or the combustion model (for instance the reaction rates)

Dear Bernahrd,

Is there any function in the new version of swak4foam to evoke reaction rate (combustion->R(Yi) in YEqn.H) of species? I need to write reaction rate of some soecies in a file.

Thank you in advance for any hint

Hossein

[gschaider]

Originally Posted by atoof

Dear Bernahrd,

Is there any function in the new version of swak4foam to evoke reaction rate (combustion->R(Yi) in YEqn.H) of species? I need to write reaction rate of some soecies in a file.

Thank you in advance for any hint

Hossein

Yes. That is what I was saying with that line.

Just add "libswakChemistryModelFunctionPlugin.so" to the libs-entry in controlDict. Now if you use any swak-thing (function-object or also funkySetFields) that uses the Field-parser for the first time you see an output like this:

Creating expression field muEff ...
"Loaded plugin functions for 'FieldValueExpressionDriver':"
psiChem_RR:
"volScalarField psiChem_RR(primitive/word speciesName)"
psiChem_RRError:
"volScalarField psiChem_RRError()"
psiChem_RRSumPositive:
"volScalarField psiChem_RRSumPositive()"
psiChem_Sh:
"volScalarField psiChem_Sh()"
psiChem_dQ:
"volScalarField psiChem_dQ()"
psiChem_deltaTChem:
"volScalarField psiChem_deltaTChem()"
psiChem_tc:
"volScalarField psiChem_tc()"
psiChem_updateChemistry:
"volScalarField psiChem_updateChemistry(primitive/scalar timestep)"

which is a list of the additional functions including their signatures. The one you're interested in is psiChem_RR which takes a species-Name as an argument.

All this uses the results from the last run-through of the chemistry-solver. In funkySetFields you've got to trigger this calculation with the psiChem_updateChemistry-function (only use it once):
"psiChem_updateChemistry(0.001)+psiChem_RR(CO)+psi Chem_RR(CO2)"
for Reaction rates of CO and CO2.

Of course you can use it to check whether the time-step influences the reaction rates
"(psiChem_updateChemistry(0.001)+psiChem_RR(CO ))-(psiChem_updateChemistry(0.002)+psiChem_RR(CO))"
but that of course calculates everything twice

[atoof]

Originally Posted by gschaider

Yes. That is what I was saying with that line.

Thanks Bernhard. It's very useful but when I try implementing the following sentences in my solver, I encounter the following error that states ‘psiChem_RR’ was not declared in this scope.

if (Y[i].name() == "O2")
{
volScalarField wO2 = psiChem_RR(Y[i].name());
wO2.write();
}


I inserted the #include "swakPsiChemistryModelPluginFunction.H" header in my .C file and the address of it in "option" which is in following:

EXE_INC = \
.
.
.
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/swak4FoamParsers/lnInclude/ \
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/swakSourceFields/lnInclude/ \
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/functionPlugins/swakChemistryModelFunctionPlugin/lnInclude

EXE_LIBS = \
.
.
.
-lsampling \
-lmeshTools \
-L$(FOAM_USER_LIBBIN) \
-lswak4FoamParsers \
-lswakChemistryModelFunctionPlugin \
-lswakSourceFields

How can I solve above problem?

Best Regards,

Hossein

[gschaider]

Originally Posted by atoof

Thanks Bernhard. It's very useful but when I try implementing the following sentences in my solver, I encounter the following error that states ‘psiChem_RR’ was not declared in this scope.

if (Y[i].name() == "O2")
{
volScalarField wO2 = psiChem_RR(Y[i].name());
wO2.write();
}


I inserted the #include "swakPsiChemistryModelPluginFunction.H" header in my .C file and the address of it in "option" which is in following:

EXE_INC = \
.
.
.
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/swak4FoamParsers/lnInclude/ \
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/swakSourceFields/lnInclude/ \
-I$(HOME)/Swak4Foam_2.x_release_0.2.3/Libraries/functionPlugins/swakChemistryModelFunctionPlugin/lnInclude

EXE_LIBS = \
.
.
.
-lsampling \
-lmeshTools \
-L$(FOAM_USER_LIBBIN) \
-lswak4FoamParsers \
-lswakChemistryModelFunctionPlugin \
-lswakSourceFields

How can I solve above problem?

Best Regards,

Hossein

Plugin-functions are designed to be used in parsers. It could be possible to use them on their own, but the calling would be plain weird. The two best options you have are
- create an instance of FieldValueExpressionParser. Feed it the expression you want to get and get the result from that. I don't explain this in detail because for the purpose of what you want to do this would be like using dynamite to go after the cockroaches in your drain: lot of fun, but the neighbours might think it is a little bit excessive.
- have a look at the source of the plugin and try to lift the parts you need (If I remember it correctly: get from the registry the object-registry the object with the name "thermophysicalProperties". Cast it down to the combustion-thermo and then call the RR there). That might be the easiest (for your purposes)

[atoof]

Originally Posted by gschaider

- have a look at the source of the plugin and try to lift the parts you need

Thanks for quick reply. I looked up the "class swakPsiChemistryModelPluginFunction_RR" in swakPsiChemistryModelPluginFunction.C and found my question.

Thanks again

[gschaider]

Originally Posted by atoof

Thanks for quick reply. I looked up the "class swakPsiChemistryModelPluginFunction_RR" in swakPsiChemistryModelPluginFunction.C and found my question.

I was hoping you might find the answer there

[gschaider]

A new version of swak4Foam has been released. It is the version found on the USB-stick of the 8th OpenFOAM workshop. New features include (but for all see the README which has a complete list):

- functionObjects now store their data in the sub-directory postProcessing according to the new convention OpenFOAM has (sorry: this may break existing scripts relying on the old location)
- slight change in the format written may break the restart of cases where swak-specific restart data was written with older swak-version (should not hit many people)
- Example cases are adapted to run with OF 2.2. This may prohibit them from running with 2.1 or 2.0

But now the good parts:

- build scripts have been adapted to run on more machines, allow tuning of the compilation process, compile on non-gcc compilers and give more information for finding problems
- A new library with additional searchableSurfaces (to be used in snappyHexMesh) that allow boolean operations and coordinate transformations of searchable surfaces
- Function objects to calculate the distributions of variables (and expressions)
- specification of aliases for field names that contain special characters (like thermo:rho)
- functionObject to list all registered objects (basically for finding out what is there)
- functionObjects to handle global variables
- functionObject that executes a list of other functionObjects only if the run is parallel/serial
- A new type of variable (stacked) that allows recording multiple values
- More flexible output control for functionObjects based on simpleFunctionObjects
- swakDataEntry is improved
- other stuff and bug fixes (check the README)

Also to make life easier for the non-Emacs people who have problems converting the (still incomplete ) reference guide to a more readable form it can be downloaded here: http://openfoamwiki.net/images/d/db/...mReference.pdf

[JR22]

Hi Bernhard,

Will swak4Foam 0.2.4 install correctly on the new OpenFOAM 2.2.1 ?

Also, when will swak4Foam 0.2.4 be posted to the openfoamwiki.net site?

Thank you very much.

[gschaider]

Originally Posted by JR22

Hi Bernhard,

Will swak4Foam 0.2.4 install correctly on the new OpenFOAM 2.2.1 ?

I think so. But the only test is trying to compile it

Originally Posted by JR22

Also, when will swak4Foam 0.2.4 be posted to the openfoamwiki.net site?

I like to think that the SVN-download at http://openfoamwiki.net/index.php/Co...am#Downloading is sufficient. Providing the tarballs was only a workaround while the SVN was down but I don't see the point creating/testing/uploading tarballs if there is already a way to get the sources

[gschaider]

8 months ago the last release of swak4Foam happened. A week later Edward Snowden blew the whistle on the NSA. Coincidence? Unlikely! (who is back there saying "Correlation does not imply causality"? Obviously someone who doesn't appreciate a good conspiracy theory). Today the next version of swak4Foam is announced and the world is expecting the next Spy-Agency to be exposed. The only question is: will it be China, Russia or Austria.

Most of the features have already leaked to the public through the public development repository and the twitter-acount @swakPyFoam but for all those whose intelligence-branch failed to interpret these informations correctly here is an incomplete list of the new features

• the version of the release has jumped to 0.3.0 as this release has now a parser for lagrangian particles (these were the last major OF "things" that swak4Foam couldn't access. It is now possible to let swak4Foam do calculations on clouds of particles
• the new release now supports Foam-3 (this is in the 2.x branch). The 1.x branch supports OF 1.7 and 1.6-ext, all newer Versions are supported in the 2.x-branch
• a number of new accumulations have been added (median, quantile and weighted counterparts for most of them)
• Python-integration has been improved (IPython inclusion, stabler etc)
• new utilities: fieldReport to quickly get information about fields, funkyPythonPostproc allows loading OF-data into a Python-workspace and analyze it with the NumPy/SciPy/Matplotlib-ecosystem of scientific software, funkySetLagrangianFields allows setting values on the particles of a lagrangian cloud
• a number of new function objects (calculation of potential flow, recalculation of thermophyiscs etc). The one with the longest name is swakExpressionAverageDistribution (applications for this are for instance the calculation of the average pressure as a function of the channel length)
• several other enhancements to utilities and function objects (all found in the README)
• a number of bugs was mercilessly squashed. Among them a problem with the correct summation of expressions on faceSets/faceZones in parallel runs
• attempts to improve compilation

A few highlights of the release are demonstrated here http://openfoamwiki.net/staticPages/...K_again.html#/

For further information see
- the Wiki page http://openfoamwiki.net/index.php/Contrib/swak4Foam
- the README that comes with the sources
- the Documentation-folder of the sources has an incomplete reference guide. It documents all the expression-grammars

Have fun and look out for the men in the trenchcoats

[gschaider]

This release has been sitting on my hard-disk for several weeks but now I've gotten around to releasing it officially (it is technically identical to the version on the USB-sticks distributed at the Workshop and also the version that is pulled when you compile it as one of the ThirdParty-packages with foam-extend-3.1)

This Release supports the latest versions of OpenFOAM (2.3) and Foam-Extend (3.1). It has its fair share of bug-fixes and also some new features:

• Compiling a "private" bison-version. The reason for this is that a growing number of distributions install bison 3.0 which is incompatible with swak4Foam. On such machines a script downloads a compatible version of bison (2.7), compiles it and installs it locally so that it will be only used for the compilation of swak4Foam
• New utility funkyWarpMesh allows deforming the mesh according to an expression
• A number of function-objects that trigger the writing of time-steps according to conditions
• Function object writeOldTimesOnSignal that writes the current state of the simulation when a floating-point-exception etc happens
• Additional library with fvOptions for OF-versions that support it
• Macro expansion in expressions

For a full list see the section about 0.3.1 in the README

[atoof]

Dear Bernhard,

How can I store the reaction rate of species O2 in each time step using psiChem_RR? Which expressions should be added to controlDict file after adding "libswakChemistryModelFunctionPlugin.so" in the libs or do something else?
I'm using swak4Foam_2.x_release_0.2.3 now.

Thank you in advance,

Hossein

[gschaider]

Originally Posted by atoof

Dear Bernhard,

How can I store the reaction rate of species O2 in each time step using psiChem_RR? Which expressions should be added to controlDict file after adding "libswakChemistryModelFunctionPlugin.so" in the libs or do something else?
I'm using swak4Foam_2.x_release_0.2.3 now.

Thank you in advance,

Hossein

Not sure what you mean with "store at each timestep" but maybe you may find the slides 125ff of http://openfoamwiki.net/images/7/7a/...icTraining.pdf helpful

[atoof]

Originally Posted by gschaider

Not sure what you mean with "store at each timestep" but maybe you may find the slides 125ff of http://openfoamwiki.net/images/7/7a/...icTraining.pdf helpful

Thanks.
That is what I need. Sorry for another question. For example If I set the deltaT value in controlDict to 1e-07, Does psichem_RR function consider that value to calculate the reaction rate?

Thanks

[gschaider]

Originally Posted by atoof

Thanks.
That is what I need. Sorry for another question. For example If I set the deltaT value in controlDict to 1e-07, Does psichem_RR function consider that value to calculate the reaction rate?

Thanks

It reuses the RR that was calculated by the chemistry and used to solve the species transport equations. If you want to check the influence of the timestep there is a separate function for it (it is explained in slides 138ff. Please get into the habit of checking more than just a few pages before asking again. You're one "under-researched" question away from being permanently ignored)

[Masgar]

Dear foamers,

I am trying to install swak4Foam on OF2.3, and I already tried the OF2.2x version suggested by Bruno but still there are some errors. I am wondering if there is a package developed for OF2.3. I am looking forward to your answers,

Thanks in advance and Best,
Mahdi

[gschaider]

Originally Posted by Masgar

Dear foamers,

I am trying to install swak4Foam on OF2.3, and I already tried the OF2.2x version suggested by Bruno but still there are some errors. I am wondering if there is a package developed for OF2.3. I am looking forward to your answers,

Thanks in advance and Best,
Mahdi

The current 2.x-release should compile with 2.3 (at least it did in June when the release was done). Please provide a link to the release you mean (Bruno is very active so looking for something he might have said takes a long time)

[Masgar]

Originally Posted by gschaider

The current 2.x-release should compile with 2.3 (at least it did in June when the release was done). Please provide a link to the release you mean (Bruno is very active so looking for something he might have said takes a long time)

Hi Bernhard,

Thanks a lot for your quick response as always. I tried the 2.x-release downloaded from:

svn checkout svn://svn.code.sf.net/p/openfoam-extend/svn/trunk/Breeder_2.0/libraries/swak4Foam/ swak4Foam_2.x

I attached the make log file output from this compile. Then I also tried this one which was suggested by Bruno (https://github.com/wyldckat/swak4foam/tree/OF22X), but I am not sure if this principally should work for me or not; anyway it did not

Thanks and Best,
Mahdi