CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

How to modify the viscosity model

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

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

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, 15:14
Default Create your own viscositymodel
  #2
Senior Member
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26
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, 14:59
Default Hi, I've also started to
  #3
Member
 
sradl's Avatar
 
Stefan Radl
Join Date: Mar 2009
Location: Graz, Austria
Posts: 82
Rep Power: 18
sradl is on a distinguished road
Hi,


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?

br
Stefan
sradl is offline   Reply With Quote

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

have a look here:
http://www.cfd-online.com/OpenFOAM_D...es/1/2513.html

regards
markus
hartinger is offline   Reply With Quote

Old   October 13, 2010, 08:44
Post help please
  #5
Senior Member
 
Join Date: Sep 2010
Posts: 226
Rep Power: 16
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:

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

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




thanks a lot
T.D. is offline   Reply With Quote

Reply

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


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