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

How to modify the viscosity model

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

LinkBack Thread Tools Display Modes
Old   July 25, 2005, 20:05
Default Hi, I am working with simpl
New Member
Jing Wang
Join Date: Mar 2009
Location: Toronto, Ontario, Canada
Posts: 7
Rep Power: 8
mpml is on a distinguished road

I am working with simpleFoam on polymer extrusion problem, using Bird-Carreau viscosity model. As is seen from source code, you used the velocity field (U) and the surface integration field (phi) to initialise the viscosity model:

autoPtr<transportmodel> laminarTransport
transportModel::New(U, phi)

My question is: how to add new factors into the viscosity model? For example, viscosity is also a function of pressure and temperature (called the WLF-Arrhenius equation).

If I want to implement this effect, I think I have to modify the sourcecode of transportModel and add pressure and/or temperature into the member functions. Can I have some guideline on the implementation? In addition, is the solution speed and convergence going to be affected seriously?

Thanks very much.
mpml is offline   Reply With Quote

Old   July 26, 2005, 14:14
Default Create your own viscositymodel
Super Moderator
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Create your own viscositymodel and when you need your T look it up from the mesh (the mesh is also a database and holds all the fields on it)

Something like:

const volScalarField& T = U_.mesh().lookupObject<volscalarfield>("T");

The solution speed should not be affected. Convergence only if there is a strong temperature dependence.
mattijs is offline   Reply With Quote

Old   July 1, 2007, 13:59
Default Hi, I've also started to
sradl's Avatar
Stefan Radl
Join Date: Mar 2009
Location: Graz, Austria
Posts: 82
Rep Power: 9
sradl is on a distinguished road

I've also started to implement a solver for polymerflow. Dissipation in that case is important - has anybody experience in adding this to the energy equation?

sradl is offline   Reply With Quote

Old   July 2, 2007, 06:30
Default Hi Stefan, have a look here
Senior Member
Markus Hartinger
Join Date: Mar 2009
Posts: 102
Rep Power: 8
hartinger is on a distinguished road
Hi Stefan,

have a look here:

hartinger is offline   Reply With Quote

Old   October 13, 2010, 07:44
Post help please
Senior Member
Join Date: Sep 2010
Location: France
Posts: 194
Rep Power: 7
T.D. is on a distinguished road
Hi Foamers,

I was wondering, how one can directly implement a new visosityLaw even as function of a Field like "T" inside the solver directly, and without going to call the predefined viscosityModels, without even to use them.
The problem is just to re-update each time the new viscosity at the new time step due to the new values of the Field "T" .
How this can be done ?

tell me if this example is true or not, where should be the update for varables?

the Example for the problem inside the solver:

if i am solving:

//My new viscosity law as function of field T
volScalaField new_nu = 2*T-5;

fvm::ddt(T) + fvm::div(phi, T) == fvm::laplacian(new_nu, T)

thanks a lot
T.D. 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
Modify kEpsilon model harly OpenFOAM 5 December 3, 2008 18:43
modify viscosity tian CFX 0 April 15, 2002 23:39
how to modify model wang FLUENT 1 April 14, 2002 03:40
how to modify the surface reaction rate model luo FLUENT 2 April 11, 2002 07:03
How to modify the model? cfxbeginer CFX 4 November 30, 2001 23:20

All times are GMT -4. The time now is 17:36.