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

Fiber Suspension Solver laminar AdvaniFolgarTucker model

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

Like Tree2Likes
  • 1 Post By kerstin
  • 1 Post By kerstin

Reply
 
LinkBack Thread Tools Display Modes
Old   January 9, 2008, 16:15
Default Hi OpenFOAM-Users, In my th
  #1
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
Hi OpenFOAM-Users,

In my thesis I have programmed a solver for laminar fiber suspension flow in OpenFOAM and before everything gets lost, I would like to contribute it to the community.
Unfortunately I have never written an english paper and didn't clean the code up. And I have not written a documentation. So everything is very messed up and done with a lot of workarounds and clumsy programming. It is possible, that there are some files already dispensable in the archive...But I give the short facts here:

-Basis: icoFoam (version OF-1.3)
-last version I compiled was with OF-1.3. I get some errors with OF-1.4...didn't try to resolve it so far

- name of the solver: nonNewton
-Fiber suspension Model:
- equations are documented in my (german) thesis
http://hdl.handle.net/2003/24559 (chapter 6)
- closed model, closure for 4th order configuration tensor, differential equation for second order configuration tensor
- hybrid closure from Advani and Tucker
- Invariant based closure IBOF-5 from Chung and Kwon
- Fiber-Fiber interactions with approach from Folgar and Tucker, interaction coefficient with approach from Bay and Phan-Thien

- models are selected with switches in constant/fiberProperties and constant/RheologyControl and constant/MatrixProperties

-Phan-Thien/Tanner (PTT) model for viscoelastic matrix. (I never really tested the implementation and didn't use it!!!)

I don't read here regularly, so if someone is really interested, I will do my best to help. Unfortunately the solver is not as ready to use, as I would prefer to have it...
So you may write an email to
thesis at kerstinheinen.de

- solver packed in nonNewton.tar.gz
download at www.kerstinheinen.de/nonNewton.tar.gz
- example 2D channel with a cylinder www.kerstinheinen.de/nonNewton_example.tar.gz

(sry, actually I am to stupid, how to upload the files in this post)

Thanks to the OpenFOAM developers. Keep it running! It's great work!

Kerstin
awacs likes this.
kerstin is offline   Reply With Quote

Old   January 29, 2009, 04:40
Default Hi all, I updated and renew
  #2
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
Hi all,

I updated and renewed the implementation of the viscoelastic Phan-Thien/Tanner Model for laminar flows in OF-1.5.

This is a multi-modular implementation for use of superposition of different parameter sets.

Attached a description of equations and implementation and also related shortcomings. This is more or less an unsophisticated hack, but still better than nothing for all interested in PTT actually.




kerstin is offline   Reply With Quote

Old   January 29, 2009, 04:43
Default http://www.cfd-online.com/Ope
  #3
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
laminarPTTMultimod.tar.gz
case_laminarPTTMultiMod.tar.gz
Description_of_PTT_Model_in_OpenFOAM-1.5.zip
awacs likes this.
kerstin is offline   Reply With Quote

Old   December 22, 2009, 09:28
Default laminar fiber Suspension code for OF-1.6
  #4
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
Hi all,
I cleaned up my old implementation of a laminar fiber suspension code and updated to OF-1.6.

Now its hopefully more straight forward to understand and include in other codes..
This code did not have extensive testing so far. It still needs more validation...So handle with care and don't blame me, if you get strange results

Actually I just copied the fiber files into icoFoam and nonNewtonianIcoFoam. Long term plan would be to have it in the same general way like the viscoelastic solver...

The code writes out also the eigenvalues of A2. lambdamax is the largest eigenvalue... Still missing is a post-processing utility to view the fiber orientation as vectors or orientation ellipsoids...

- The test case attached is a channel with cylinder.
- just run the script runAll.
- general advice: the fiber suspension solver should always have as a first step running for a newtonian velocity field with icoFoam...otherwise it can be quite unstable...

- the example case is not perfect. You will get some warnings for negative eigenvalues, which dissapear after some time...its maybe due to bad initialisation...

One nice extension would be a model for turbulent flows...where orientation of the fibers is coupled to turbulence properties... Is there actually anyone out there working on fiber suspension models (based on the orientation tensor approach) for turbulent flows?

Kerstin
Attached Files
File Type: gz fiberFoam.tar.gz (16.4 KB, 64 views)
File Type: gz icoFiberFoam_channel.tar.gz (5.3 KB, 55 views)
kerstin is offline   Reply With Quote

Old   January 10, 2010, 12:34
Default
  #5
New Member
 
Kerstin Heinen
Join Date: Mar 2009
Location: Ludwigshafen, Germany
Posts: 27
Rep Power: 8
kerstin is on a distinguished road
Next step done: attached a pdf doc with the equations used in the solver uploaded above...



Kerstin
Attached Files
File Type: pdf laminarFiberFoam.pdf (80.7 KB, 108 views)
kerstin is offline   Reply With Quote

Old   November 21, 2014, 23:49
Default
  #6
New Member
 
pavan
Join Date: Nov 2014
Posts: 1
Rep Power: 0
pawanhv is on a distinguished road
Hello, I am a beginner here.. I am interested to know the difference in results in the example here if I change the duration from 2mins to 4mins? , It crashes after 4.1mins I guess
pawanhv is offline   Reply With Quote

Reply

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
Continuous Fiber Model sunyafeng FLUENT 0 June 6, 2008 09:24
SimpleFoam as laminar solver anja OpenFOAM Pre-Processing 3 April 5, 2007 09:28
algebraic slip model with stirred suspension? Manuel Krause CFX 0 April 15, 2005 09:10
Laminar solver? LES? or DNS? Ray Main CFD Forum 5 March 19, 2003 11:19
i wanna model draw glass fiber Stone D. CFX 0 May 27, 2002 02:20


All times are GMT -4. The time now is 19:25.