|October 15, 2012, 11:11||
Most Appropriate Solver ??
Join Date: Sep 2009
Posts: 65Rep Power: 9
I am new to the CFD scene, although I have tried to dabble a little.
I'm now trying to model a real world problem. This problem is the opening of a valve and discharge of a compressible fluid into a piping system and through a restriction orifice. I do not intend to model the movement of the valve, just the ramped increase of the gas flow through the inlet boundary.
I have developed a 2D model of this domain which checkMesh seems happy with, showing no errors.
Reading the documentation that comes with OpenFOAM has lead me to believe that I should be looking to use rhoSimpleFoam or rhoPimpleFoam for this model although this is probably as far as I have got with this. I have tried to understand how to use these solvers. I want to model this problem as compressible and turbulent. This is my first problem:
The physical problem involves a very large pressure gradient across the orifice (or shock downstream). The valve inlet pressure will be 16.5 MPa, I have read that RANS K-e solvers would not be appropriate for such a large pressure gradient so I guess I should be utilising LES (but how?)
I have made the domain quite small, the pipe has an ID of 52mm (orifice 20.3 mm ID) and I have constructed the mesh to show about 1000mm upstream and downstream of the orifice. Can I leave OpenFOAM to calculate the outlet boundary Pressure since I do not know the pressure to specify? Essentially this would be a function of the inlet pressure and velocity through the orifice.
From a basic flow calculation across the orifice with 1kPa outlet pressure I would have a mass flow of approx. 10 kg/s. This creates some VERY large velocities downstream of the orifice, I have calculated the Reynolds number for this problem as 3.7e7. If I calculate the velocity in the low pressure side of the orifice I get about Mach 10!
I guess I have rambled on a little so basically I need to understand the solver to use to model a high velocity compressible fluid through a pipe/orifice system with a large pressure gradient. I would also like to model some of the thermophysical and transport properties but I'm happy to go inviscid and ideal as a first step if you advise this!
I hope someone can give me some pointers to get me moving with this!
|October 18, 2012, 06:46||
Join Date: Sep 2009
Posts: 65Rep Power: 9
I've done some more experimenting with various solvers and I have settled on rhoPimpleFoam using k-epsilon as my RAS model rather than LES. I'm currently just running with a low velocity while I experiment.
I would like to understand, though how I should set my boundary conditions? I have set my inlet BC for U as a fixedValue as I plan to implement a ramp and this seems the most obvious BC to use for this method. I have set my outlet U as zeroGradient as I would like OpenFOAM to simply report what it calculates rather than stipulating this flow (over specifying??) my outlet P is specified as a fixedValue.
When I run this, however, I notice that the velocity (Magnitude) begins to build from both inlet AND outlet boundaries which I know is wrong since I'm modelling a compressible fluid.
Have I mistakenly used zeroGradient? is there a more appropriate BC to use on the outlet boundary?
I have attached links to my files and would appreciate someones advice as to whether I am going about this in the right way?
btw I used the angledDuct example in the compressible/rhoPimpleFoam/ras directory as a template for this attempt.
|Thread||Thread Starter||Forum||Replies||Last Post|
|thobois class engineTopoChangerMesh error||Peter_600||OpenFOAM||4||August 2, 2014 09:52|
|A New Solver for Supersonic Combustion||nakul||OpenFOAM||6||March 30, 2014 03:00|
|A New Solver for Supersonic Combustion||nakul||OpenFOAM Announcements from Other Sources||18||February 19, 2013 08:48|
|Development of a Low mach PISO solver||nishant_hull||OpenFOAM Programming & Development||0||August 25, 2009 12:48|
|why the solver reject it? Anyone with experience?||bearcat||CFX||6||April 28, 2008 14:08|