
[Sponsors] 
[Other] Multi species mass transport library [update] 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 10, 2018, 22:44 
fixedFlux boundary condition uploaded ?

#101  
Member
Haomin Yuan
Join Date: Jan 2012
Location: Madison, Wisconsin, USA
Posts: 59
Rep Power: 13 
Hi, all Foamers,
I just want to follow up this thread. Is the fixedFlux boundary condition available now ?? It is very important for the test cases in the paper. Thank you. Quote:


June 29, 2018, 23:27 

#102  
New Member
Gao Zhengwei
Join Date: Jan 2017
Location: HangZhou, P.R.China
Posts: 10
Rep Power: 8 
Quote:
I am new here and trying to modfity this code. Nikilas provided a good suggestion and I noticed that some codes are needed to be modified to ensure the consistence. In file multiSpeciesTransportModel.C and multiSpeciesTransportModel.H: I added a private member inertSpecie_ , which is initialized as inertSpecie_(thermo.lookup("inertSpecie")), In the function correct() in file Fick.C, a slight change is made: [I] volScalarField yt = 0.0*thermo_.composition().Y(0); surfaceScalarField nt = turbulence_.phi(); label inertIndex = 1; forAll(species(), i) { volScalarField& yi = thermo_.composition().Y(i); surfaceScalarField& ni = n_[i]; if (yi.name() != inertSpecie_) { tmp<fv::convectionScheme<scalar> > mvConvection ( fv::convectionScheme<scalar>::New ( mesh_, fields, turbulence_.phi(), mesh_.divScheme("div(phi,Yi_h)") ) ); if (mesh_.relaxField("Yi"))//Mohsen { yi.storePrevIter(); } tmp<fvScalarMatrix> yEqn ( fvm::ddt(thermo_.rho(), yi) // + fvm::div(turbulence_.phi(), yi, "div(phi,Yi_h)") + mvConvection>fvmDiv(turbulence_.phi(), yi)  fvm::laplacian(D_[i],yi, "laplacian(D,Yi)") == Sy_ ); eqnResidual = solve(yEqn() , mesh_.solver("Yi")).initialResidual(); maxResidual = max(eqnResidual, maxResidual); if (mesh_.relaxField("Yi"))//Mohsen { yi.relax(mesh_.fieldRelaxationFactor("Yi"));//Mohsen } yi.max(0.0); // yi.min(1.0); ni = yEqn().flux(); nt = ni; yt += yi; } else { inertIndex = i; } } // Calculate inert species volScalarField& yInert = thermo_.composition().Y()[inertIndex]; yInert = 1  yt; forAll(yInert.boundaryField(), boundaryI) { forAll(yInert.boundaryField()[boundaryI], faceI) { yInert.boundaryField()[boundaryI][faceI] = 1 yt.boundaryField()[boundaryI][faceI]; } } yInert.max(0.0); n_[inertIndex] = nt; updateMolarFractions(); return maxResidual; Note that, the inertIndex is not the same inertIndex_ in multiSpeciesTransportModel.H I hope this will be helpful. 

February 17, 2019, 18:43 
Update of library to OpenFOAM 6

#103 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 649
Rep Power: 20 
Hello,
I have created a noncomplete update of the library to OpenFOAM6. Main changes:  Handling of references, which are now const  The definition of turbulence.phi() seems to have changed: It is now a volume flux, but a mass flux is need, so that has been added as additional internal variable  Chemistry has changed so much, that I have (hopefully temporarily) disabled it. It compiles and runs (TM). I just tested it with the testCase given in the original library without chemistry. Any comments/improvements are welcome. 

February 18, 2019, 17:14 

#104  
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 649
Rep Power: 20 
Updates:
 To use transport const, thermo hConst change the commented lines in multiSpeciesTransportModels.C from Code:
// makeMultiSpeciesTransportModel(Fick, gasThermoPhysics); ... Code:
makeMultiSpeciesTransportModel(Fick, constGasHThermoPhysics); ... To start the attached test case, issue the following commands: Code:
blockMesh setFields modifiedReactingFoam Quote:


April 26, 2019, 09:39 

#105 
New Member
NKR
Join Date: Jan 2019
Posts: 2
Rep Power: 0 
hello mhsn, how did u solve that error /usr/bin/ld: cannot find llduSolvers.


May 14, 2019, 10:25 

#106 
Member
Join Date: Mar 2016
Posts: 73
Rep Power: 9 
Hello everyone,
thank you so much for sharing the library! This is very helpful! One thing I didn't understand is: There are models to calculate the diffusioncoefficient. What happens when I want to use Schmidtdiffusion for example? The diffusioncoefficient is calculated with mu, density and the schmidtnumber. Are all models for calculating the Diffusioncoefficient like champanenskog disabled for this model? Or are these different diffusioncoefficients? Thank you very much! 

May 19, 2019, 17:16 

#107 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 649
Rep Power: 20 
You set the diffusion model in the file constant/transportProperties, e.g.
Code:
multispeciesTransportModel SchmidtNumber<constGasHThermoPhysics>; Code:
Sc_(IOdictionary::lookupOrDefault<scalar>("Sc", 1)) Code:
this>D_[i] = this>turbulence_.muEff() / Sc_; If you use another mode, e.g. Fick, then you have to set Code:
multispeciesTransportModelFick<gasHThermoPhysics>; Have a look at the README.txt file in the archive and the pdf somewhere in this thread. 

August 23, 2019, 09:53 

#108 
Member
Thamali
Join Date: Jul 2013
Posts: 67
Rep Power: 11 
Hi Jherb and all,
Thanks a lot for all who contributed to this library. If you have updated the solver for OpenFOAM6 with chemistry, I would be grateful to have it. Otherwise I will try do it on my own. Thanks in advance. Thamali 

August 25, 2019, 15:06 

#109 
Senior Member
Joachim Herb
Join Date: Sep 2010
Posts: 649
Rep Power: 20 
You can find the current state of the port at https://github.com/jmozmoz/OFdevelopments
There are still bugs in treating the passive species. But based on a talk "Development and verification of a multispecies gas transport solver" by Vijaya Kumar G et al., at the OpenFOAM Workshop at Duisburg the problems with temperature differences due to diffusion should be solve. 

March 24, 2021, 01:44 
Tritium transport equation

#110 
New Member
Pratyush Kumar
Join Date: Jun 2019
Location: Mumbai
Posts: 19
Rep Power: 5 
Dear Valerio,
First of all, I would like to thank you for this updated library. I want to use this library for tritium transport using the ChampanEnskog model for diffusion coefficient and Fick’s model for flux modelling in a porous media domain. I am not interested in reaction mechanisms for now. As per my understanding to use it on some membrane domain, I should modify the solver and change the value of porosity and tortuosity in the respective libraries. Can I use this solver on the membrane without including the flux modelling? The reason why I am asking this question is because the solubility or entrapment effects have not been considered. If the answer is not, what modifications will be required? Also, I have a little confusion about the reaction mechanism of the solver. As I understand by going through solver, your solver can only solve the combustion reactions? Correct me if I'm wrong? So, if I want to model the copermeation of two molecules H2, D2 through membrane, a reaction like H2+D2> 2HD (where H is hydrogen, D is deuterium) can be modeled for the same, how should I begin to include it in the same solver? 

August 13, 2021, 12:43 

#111 
New Member
Ejaz Ahmed
Join Date: May 2021
Posts: 7
Rep Power: 4 
I am thanking you all for the contributions you have made. I'm currently working on icoReactingMultiphaseInterFoam solver on OpenFOAM  v2012 in order to model liquid evaporation which is using VOF (volume of fluid) based interface capturing method. I figured out that, the solver doesn't have any mass diffusivity in species transport equation. I need to modify the species transport equation for that. Can I install this "Multispecies Mass Transport Library" in OpenFOAM  v2012 & use the library in the icoReactingMultiphaseInterFoam solver? If not, what modifications do I need to make? I will really appreciate your help & suggestions. Thank you in advance!


November 9, 2021, 23:37 
Please Update the Library for New versions of OpenFOAM

#112 
New Member
Sk Hossen Ali
Join Date: Jul 2021
Location: India
Posts: 8
Rep Power: 3 
Dear Foamers,
Thanks a lot for all the efforts till now, this is to bring to your kind attention that the libraries that you provided is not compatible with new versions of openfoam. Please guide me on how to use the library if it is possible to run it in any currently supported versions of openfoam, or please try update it for the current versions of openfoam. I have tried to install openfoam2.3.0 in the ubuntu20.04LTS but did not compilled correctly. While I was trying to install Ubuntu13 it said it is no longer supported. Thanks in advance Sk Hossen Ali 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Species Mole and Mass Fraction Macro  combustion  FLUENT  17  February 9, 2023 10:04 
UDS and fluent internal species transport model  jinsong  FLUENT  0  May 3, 2018 11:37 
Inconsistencies in reading .dat file during run time in new injection model  Scram_1  OpenFOAM  0  March 23, 2018 22:29 
Surface instabilities with 2phase species transport  richard222  Fluent Multiphase  0  May 4, 2016 04:06 
Questions for a species transport problems (snapshots attached)  aleisia  FLUENT  2  October 9, 2011 04:40 