Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

axial scaling of navier-stokes and brinkman equations

Please login with a confirmed email address before reporting spam

Hi,

im using the free and porous media flow physics system in Comsol 4.0. The thing is im considering hollow fibres where the length of the domain is around 1000 times bigger than the radius (2d axisymmetric coordinates). Trying to mesh and solve something like this is troublesome. I can use convection-diffusion equations in a scaled form but im wondering if anyone knows how to do it with the NS and brinkman equations?

Thanks,
Adam

2 Replies Last Post 2010/06/23 7:30 GMT-4

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010/06/17 14:09 GMT-4
I'm using a mapped mesh for much the same thing: brinkman in porous and RANS NS in fluid. My aspect ratio also about 1000. Results are surprisingly stable using a quad mesh, mapped so as to have lots of elements near the interface. If I mesh too finely at the interface then it becomes unstable again, presumably reaching some limit of aspect ratio.

I have a different question though: what is the difference between using the free/porous interface versus having two separate physics models. Are the equations on the interface different? I force continuity of pressre and velocity and basically allow a discontinuty in viscous stress. I dont know if that is correct!!

Thanks, John
I'm using a mapped mesh for much the same thing: brinkman in porous and RANS NS in fluid. My aspect ratio also about 1000. Results are surprisingly stable using a quad mesh, mapped so as to have lots of elements near the interface. If I mesh too finely at the interface then it becomes unstable again, presumably reaching some limit of aspect ratio. I have a different question though: what is the difference between using the free/porous interface versus having two separate physics models. Are the equations on the interface different? I force continuity of pressre and velocity and basically allow a discontinuty in viscous stress. I dont know if that is correct!! Thanks, John

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 2010/06/23 7:30 GMT-4
Hi John,

thanks for the reply. ive started using the mapped mesh now too and it seems ok although its impossible to show a 2d surface plot with non-equal axes, even though you can view the geometry this way. i was getting strange results for the free and porous media flow physics so I thought id try it your way. can i ask how do you set up the continuity of flow and pressure at the boundary between the free medium and porous matrix?

Thanks,
Adam
Hi John, thanks for the reply. ive started using the mapped mesh now too and it seems ok although its impossible to show a 2d surface plot with non-equal axes, even though you can view the geometry this way. i was getting strange results for the free and porous media flow physics so I thought id try it your way. can i ask how do you set up the continuity of flow and pressure at the boundary between the free medium and porous matrix? Thanks, Adam

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.