|
[Sponsors] |
Printing out mixture-averaged diffusion coefficients in Openfoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 16, 2023, 11:43 |
Printing out mixture-averaged diffusion coefficients in Openfoam
|
#1 |
New Member
Harish Gopal
Join Date: Mar 2023
Posts: 4
Rep Power: 3 |
Hi everyone,
I hope you are having a nice day ! For my research, I am studying the dynamics of premixed flames forced by harmonic velocity disturbances. To achieve this goal, I am running some forced computations of a slot-stabilized flame in OpenFoam (reactingFoam). I am using mixture averaged Wilke Rule to compute the species diffusion coefficients (mixture averaged Dmix) from the binary diffusion coefficients (Dij). These binary diffusion coefficients, I am specifying as a table in my thermophysicalTransport file (Specify transport properties using function2 object in OpenFoam). To check if the mixture-averaged diffusion coefficient is calculated correctly (and to check if openfoam has correctly read the table values that I gave in), I want to print the mixture averaged diffusion coefficients (Dmix). The definition of Dmix is given in the file /opt/openfoam10/src/ThermophysicalTransportModels/laminar/Fickian/Fickian.H as the following //- Effective mass diffusion coefficient // for a given specie mass-fraction [kg/m/s] virtual tmp<volScalarField> DEff(const volScalarField& Yi) const; Please, can someone tell me how to print it each and every write time interval along with p,T, U and so on.. I know it is possible to print out the thermophysical properties by writing print statements in the solver file. For example, suppose one wants to print kappa (thermal conductivity), one has to add the following lines to the reactingFoam.c file k_output = thermo.kappa(); and then add the following in the createFields.H file volScalarField k_output ( IOobject ( "k_output", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), thermo.kappa() ); But how do I do the same for DEff ? The reason I ask this is, DEff is a quantity which varies for each species, no ? So should I write it species-by-species for example: thermophysicalTransport.DEff('O2'), thermophysicalTransport.DEff('CH4')... and so on.. I would be extremely grateful if one of you can give some insight on this please. Thanks, Harish |
|
February 16, 2024, 22:54 |
|
#2 |
New Member
Join Date: Nov 2022
Posts: 2
Rep Power: 0 |
Hi Harish,
I was just wondering if you have got the answer to your question or somehow figured it out yourself. This is a question for me too. I would appreciate it if you verified whether or not your method works as expected. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Printing out drag coefficients in openFoam | fredrick | OpenFOAM Post-Processing | 3 | May 27, 2021 16:55 |
is it possible to define different mass diffusion coefficients in different zones? | Weiqiang Liu | Main CFD Forum | 0 | June 3, 2019 21:24 |
any macros with which mixture averaged diffusion coefficients can be obtained? | Weiqiang Liu | FLUENT | 0 | May 27, 2019 13:51 |
Mixture averaged/Multicomponent transport in openFoam reacting flow simulation | sahmed | OpenFOAM Running, Solving & CFD | 1 | March 5, 2017 02:21 |
time averaged mixture fraction | zkzkzk010537 | FLUENT | 11 | March 23, 2016 07:36 |