CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Question in definition of terms in solve

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

LinkBack Thread Tools Display Modes
Old   March 19, 2009, 17:02
Question Question in definition of terms in solve
Senior Member
Antonio Martins
Join Date: Mar 2009
Location: Porto, Porto, Portugal
Posts: 112
Rep Power: 8
titio is on a distinguished road
Send a message via MSN to titio Send a message via Skype™ to titio
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?


A possible solution involves the definition of
titio is offline   Reply With Quote


Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
CellZones to solve different eq in different zones ivan_cozza OpenFOAM Running, Solving & CFD 7 June 30, 2011 10:53
How to solve energy eqn with enthalpy terms? sabitha FLUENT 0 April 6, 2006 04:58
UDS question Z FLUENT 1 March 31, 2005 15:08
UDS Flux Question - HARD! Greg Perkins FLUENT 1 March 14, 2002 22:54
Is CFD Science or Art ? John C. Chien Main CFD Forum 36 October 5, 1999 12:58

All times are GMT -4. The time now is 11:58.