Question in definition of terms in solve

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

 March 19, 2009, 17:02 Question in definition of terms in solve #1 Senior Member   Antonio Martins Join Date: Mar 2009 Location: Porto, Porto, Portugal Posts: 112 Rep Power: 10 Sponsored Links Dear Foamers, I am trying to model a problem that involves solving an equation in terms of the stress tensor that includes the convected derivative of the stress tensor and the deformation rate tensor, in the form (convected derivative of the stress Tensor)+K*deformation rate & stress tensor)+K*(stress tensor & deformation rate)=viscosity*deformation rate where & means inner product. The convected derivative is a combination of four terms, in the form ddt(stress tensor)+gradient (velocity & stress tensor) - stress tensor & gradient (velocity) - gradient (transpose of the velocity) & stress tensor When I represent this equation in OpenFoam, i only consider implicitly two terms: - ddt (stress tensor): fvc::ddt(stresstensor) - gradient (velocity & stress tensor)= div(phi,stresstensor) The remaining terms are not discretized neither explicitly or implicitly. The final expression is in form solve ( fvc::ddt(stresstensor) + div(phi,stresstensor) = viscosity*Deformation rate + stresstensor & Deformation rate + transpose(deformation rate) & stresstensor-K*(deformation rate & stress tensor) - K*(stress tensor & deformation rate) ); In the right side none of the terms is written in the form fvc::.. Only the deformation rate involves the calculation of the gradient of the velocity field, in the form fvc::grad(velocity). Is this enough to ensure that the remaining terms are represented also as fvc::? The expression that I wrote compiles but is it O.K? Should I discretize, using perhaps fvc, the terms in the right side? I have tried to do that directly but with no sucess. I got errors of invalid operation or problems in the units. However, in all the solvers that I analysed when solving equations the terms in the right side involve allways fvc. Can anyone help me out? Titio A possible solution involves the definition of

 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 ivan_cozza OpenFOAM Running, Solving & CFD 7 June 30, 2011 10:53 sabitha FLUENT 0 April 6, 2006 04:58 Z FLUENT 1 March 31, 2005 15:08 Greg Perkins FLUENT 1 March 14, 2002 22:54 John C. Chien Main CFD Forum 36 October 5, 1999 12:58