I am modelling gas combustion in porous media. Most contemporary models assume that the pressure drop from the porous media is small enough to disregard, but I want to include that in my investigation. Most studies alos assume that the gas is filtered through the media at a constant speed $u$.

I use Darcy's law to couple the velocity and pressure feilds:

$u=-\frac{k}{\mu}\nabla p$.

Plus we have the ideal gas euqation of state and the continuity equation:

$p=\rho RT/M_r$

$\frac{\partial \rho}{\partial t}+\nabla \cdot (\rho u)=0$.

Since we are talking about combustion the case is not isothermal, two energy euqations are also present. All the quations are solved with finite differences.

I first of all want to qualitatively reproduce the results obtained in previous studies to validate my model. The question is how do I simulate a constant velocity in the domain. It seems like the Darcy equation implues that I should apply a proper pressure drop, i.e first order boundary conditions on the pressure. This implies a big pressure drop if I want to obtain something around $u=0.5$ m/s in the camera. Moreover the velocity appears to be increasing rapidly towards the far end, where atmosperic pressure is prescribed.

I was thinking perhaps that I ran into this problem limitation as a payment for not solving a real Navier-Stockes. But I am not sure, so here I am :).

P.S. I do realize that I can through away the gasodynamical part of the model and just prescribe the speed in the energy equations as far as validation goes, I have indeed done that and got good results.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.