CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > 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 Search this Thread 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: 13
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 Search this Thread
Search this Thread:

Advanced Search
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 Off
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 13 October 24, 2019 17:07
How to solve energy eqn with enthalpy terms? sabitha FLUENT 0 April 6, 2006 05:58
UDS question Z FLUENT 1 March 31, 2005 16: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 13:58

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