# Momentum equation of Darcy's law

 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 17, 2009, 02:25 Momentum equation of Darcy's law #1 New Member   andrianomena Join Date: Mar 2009 Location: madagascar Posts: 23 Rep Power: 8 Hi there, I'm dealing with combustion modeling which takes into account Darcy's law for the momentum equation.I'm trying to solve thi momentum equation and it seems like I got trouble: the velocity equation is: U = -K*fvc:grad(p); where K is constant to be calculated the pressure equation is: div(-K1*grad(p2)) = source term to implement this into the code I tried this following syntax: -fvm:laplacian(-K1,pow(p,2)) == S where S is the source term to be calculated also and I came out with an error message like GeometricField:double....or something Can somebody out there help me out and tell me what I have to do to implement this pressure equation into the code? I have already posted this but nobody answered yet Please, Regards,

 April 17, 2009, 02:49 #2 Senior Member   Mark Olesen Join Date: Mar 2009 Location: http://olesenm.github.io/ Posts: 777 Rep Power: 18 You might wish to examine the rhoPorousSimpleFoam solver to gather some ideas.

 April 17, 2009, 03:27 #3 New Member   andrianomena Join Date: Mar 2009 Location: madagascar Posts: 23 Rep Power: 8 I have just taken a look at the rhoPorousSimpleFoam solver, but you know, I am new to Openfoam and there are some syntax in the code I dont know and I cant find in the programmers guide such as 'tmp<>', so I'm struggling a bit. My question is: is it possible to have this fvm:laplacian(k,pow(p,2)) in the lapalcian scheme? or do I have to cast 'p' in the coefficient 'k' and I get this fvm:laplacian(k1,p)? Regards,

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mihail CFX 7 September 7, 2014 06:27 kk Main CFD Forum 0 July 2, 2006 23:13 kk Main CFD Forum 0 July 2, 2006 23:10 George FLUENT 0 May 19, 2005 13:07 Paul Missel Main CFD Forum 6 July 20, 2001 05:26

All times are GMT -4. The time now is 04:59.