search

MicroNanoFlows / OpenFOAM-2.4.0-MNF

OpenFOAM 2.4.0 plus the MicroNanoFlow Group Codes
http://www.micronanoflows.ac.uk/
Other
50 stars 37 forks source link

Boundary Conditions with calculated mole fractions gives FPEs #13

Open Previouslysane opened 2 years ago

Previouslysane commented 2 years ago

Hi MicroNanoFlows group,

This is regarding the pressure BCs with calculated molar masses. When initializing a simulation with no inlets and just 1 outlet, any BC with a calculated molar mass doesn't seem to work.

However, switching to a BC with calculated molar mass after starting the run with a specified molar mass condition for the first 10 steps seems to work. I am guessing this might have to do with 0 inserted particles in the first few steps, but didn't have time to investigate the code further.

Any help on the topic is greatly appreciated!

craigwhite87 commented 2 years ago

I have just tested this on the microScaleTestCase tutorial ( https://github.com/MicroNanoFlows/OpenFOAM-2.4.0-MNF/tree/devel-craig/tutorials/discreteMethods/dsmcFoamPlus/microScaleTestCase ) by turning the inlet patch into a diffuse wall and the changing the outlet to use dsmcLiouFangPressureOutletCalculatedMolarFraction and it works ok. Do you have a test case to try this on?