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

block solver already available in OpenFOAM?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 11, 2010, 06:39
Default block solver already available in OpenFOAM?
  #1
Senior Member
 
Stefan Herbert
Join Date: Dec 2009
Location: Darmstadt, Germany
Posts: 128
Rep Power: 7
herbert is on a distinguished road
Hi FOAMers,

I've read a lot of times that a block matix solver was planned to be implemented. It should be something comparable to coupledMatrix but capable for multi-variable-coupling instead of multi-region-coupling.

Does anyone know about the proceeding for that. I'd like to use it for finite-element stress analysis, because when I'm solving the equations segregated, the solution converges far too slow (or even not at all).

Thanks,
Stefan
herbert is offline   Reply With Quote

Old   May 11, 2010, 11:37
Default
  #2
Member
 
Patricio Bohorquez
Join Date: Mar 2009
Location: Jaén, Spain
Posts: 84
Rep Power: 7
pbohorquez is on a distinguished road
Let's see what is going on the "New Implementations in OpenFOAM" to be introduced at the "5th OpenFOAM Workshop". It seems that there are several talks planned on that topic.

Exciting!

Last edited by pbohorquez; May 12, 2010 at 16:32.
pbohorquez is offline   Reply With Quote

Old   May 12, 2010, 16:17
Default
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,718
Rep Power: 19
hjasak will become famous soon enough
Yes, it is already in use in a few projects (neutronics simulations, coupled multi-phase VOF) and to come out at the Workshop. There will be a presentation by Julia to show what Kathrin and Julia have done at the Summer School in Zagreb last year + I am seeing further work by Ivor Clifford in neutronics at Penn State.

In short: it is implemented, parts of it are due for release and we need to apply it to more cases.

Hrv
__________________
Hrvoje Jasak
hjasak is offline   Reply With Quote

Old   July 9, 2010, 05:00
Default
  #4
New Member
 
Eric.wang
Join Date: Jun 2009
Location: beijing china
Posts: 10
Rep Power: 6
wang.zhy is on a distinguished road
Hi Hrv
I downloaded your presentation "OpenFOAM: Year in Review"。and I found these big good news
New Features in Upcoming Release:
Block matrix implementation with parallelisation support
Re-meshing with tetrahedral edge swapping: Sandeep Menon
......

It may be rude ,but I can't help myself.
I really need it !
where and when could I find it ?
wang.zhy is offline   Reply With Quote

Old   July 9, 2010, 05:11
Default
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,718
Rep Power: 19
hjasak will become famous soon enough
Well, I am finishing off the merge and 1.7-ext will follow up immediately afterwards. However, the block matrix is not so easy (have a look at the work Julia and Kathrin have done with it to have an idea of what to do).

As for Sandeep's contribution, this will slot right in, with tutorial examples etc. If we ask him nicely, Sandeep will give us some documentation as well

I will keep you posted,

Hrv
__________________
Hrvoje Jasak
hjasak is offline   Reply With Quote

Old   July 9, 2010, 05:34
Default
  #6
New Member
 
Eric.wang
Join Date: Jun 2009
Location: beijing china
Posts: 10
Rep Power: 6
wang.zhy is on a distinguished road
Thanks Hrv!
wang.zhy is offline   Reply With Quote

Old   July 9, 2010, 18:01
Default
  #7
Senior Member
 
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 378
Rep Power: 9
marupio is on a distinguished road
Quote:
Originally Posted by wang.zhy View Post
Hi Hrv
I downloaded your presentation "OpenFOAM: Year in Review"。and ...
Really? I've been looking for it. Where can I find it? I was unable to go to the workshop, and was looking forward to reading that one.

-Dave
marupio is offline   Reply With Quote

Old   July 10, 2010, 02:15
Default
  #8
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 80
Rep Power: 7
juho is on a distinguished road
Quote:
Originally Posted by marupio View Post
Really? I've been looking for it. Where can I find it? I was unable to go to the workshop, and was looking forward to reading that one.

-Dave
You can find most of the presentations here:
http://web.student.chalmers.se/groups/ofw5/Program.htm

-Juho
juho is offline   Reply With Quote

Old   July 10, 2010, 07:11
Default
  #9
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,880
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by hjasak View Post
Well, I am finishing off the merge and 1.7-ext will follow up immediately afterwards. However, the block matrix is not so easy (have a look at the work Julia and Kathrin have done with it to have an idea of what to do).
I saw the slides. Very interesting work! If I understood it correctly, they
  • Estimate U as usual, with U = H/A.
  • Solve the coupled p, alpha1, alpha2 equations
  • Correct p and fluxes
  • "Re-evaluate" U
I have a couple of questions:
  1. Is the re-evaluation of U done reconstructing the flux (explicitly or implicitly?)?
  2. Would it be possible to deal with coupled U equations (like in FLUENT phase-coupled SIMPLE for Euler-Euler multi-fluid model)? In that algorithm, n equations for U are coupled to predict U, then a pressure equation based on total phase continuity is solved and phase velocities are corrected using a (point-implicit in FLUENT, but I'm not that worried about this point. I'm more interested in the coupling of U equations :-)) flux reconstruction.
Thanks
__________________
Alberto

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
GeekoCFD 32bit - The 32bit edition of GeekoCFD.
GeekoCFD text mode - A smaller version of GeekoCFD, text-mode only, with only OpenFOAM. Available in a variety of virtual formats.
alberto is offline   Reply With Quote

Old   July 11, 2010, 01:58
Default
  #10
New Member
 
Eric.wang
Join Date: Jun 2009
Location: beijing china
Posts: 10
Rep Power: 6
wang.zhy is on a distinguished road
Be greedy ,It would be great if OpenFOAM-1.7-dev have a full coupled solver.There is paper "A coupled finite volume solver for the solution of incompressible flows on unstructured grids" in JCP .and I saw L. Mangani and C. Bianchini have accomplish this algrithm with a third party linear solver.In my view of their reslult ,both stablish and convergence have a huge enhancement compare to segregated solver.I guess computation cost will also be saving .

In my test,for a turbulence Flat plate case (like the case in http://www.lerc.nasa.gov/WWW/wind/valid/archive.html),

In general CFX 11 need about 60 iteration to achieve a 1e-6 convergence ,while OpenFOAM1.6 need 3000 or more iteration to reach the same convergence .

CFX11 use a full coupled solver,and OpenFOAM 1.6 use segragated solver ,maybe that is the main cause of the difference .
wang.zhy is offline   Reply With Quote

Old   September 29, 2010, 11:43
Default
  #11
Senior Member
 
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 378
Rep Power: 9
marupio is on a distinguished road
Has any progress been made on the block matrix solvers? I have a biochemistry model with 30 scalars that need to be coupled with the transport equations. I'm currently using an ugly mix of an ODE solver with standard transport.
marupio is offline   Reply With Quote

Old   December 30, 2010, 14:18
Default
  #12
Senior Member
 
Nilesh Rane
Join Date: Apr 2010
Posts: 119
Rep Power: 6
nileshjrane is on a distinguished road
Hello FOAMers..

First of all,
Wishing you all a very Happy new year, 2011...May we all as the open source CFD community prosper in the coming year like non ever before...

My question to Hrv,

Where can i find the literature regarding the block matrix solver?? I am currently looking for the various options for coupled solvers...

Thank you,
Nilesh
__________________
Imagination is more important than knowledge..
nileshjrane is offline   Reply With Quote

Old   December 30, 2010, 14:21
Default
  #13
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,718
Rep Power: 19
hjasak will become famous soon enough
M&C2009 paper by Ivor Clifford
ECCOMAS 2010 by Kathrin Kissling and Julia Springer

Enjoy,

Hrv
__________________
Hrvoje Jasak
hjasak is offline   Reply With Quote

Old   December 30, 2010, 15:13
Default
  #14
Senior Member
 
Nilesh Rane
Join Date: Apr 2010
Posts: 119
Rep Power: 6
nileshjrane is on a distinguished road
That was quick..

I have one more question. I have worked on a research code. The code is like, all the governing equations (including turbulence, species conservation equations, momentum and energy) written in vector form in transformed co-ordinates. Then a Jacobian is formed and inverted using a specialised (for parallel computation) algorithm. This is fully coupled solver system. Looking at the present state of OpenFoam, is it possible to model such solver in it?? i mean does it have the required machinery (building blocks) for that?? i am concerned with this because this would be the ideal way (may not be for industry but for research at least) of solving complex flows like hypersonic reactive flow etc. Please correct me if i am wrong.

And if i am not guessing incorrectly, the block matrix solvers have clear advantage over such implicit methods in terms of computational costs.


Regards,
Nilesh
__________________
Imagination is more important than knowledge..
nileshjrane 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
Integrated conjugate heat transfer solver in OpenFOAM hjasak OpenFOAM Running, Solving & CFD 168 June 27, 2013 01:56
Creating New Solver: For particle-laden compressible jets sankarv OpenFOAM 0 April 4, 2010 18:06
The OpenFOAM extensions project mbeaudoin OpenFOAM 16 October 9, 2007 09:33
Block Tridiagonal Solver Abdulhafid M. Elfaghi Main CFD Forum 2 December 23, 2006 12:20
Block tridiagonal matrix solver Grace Main CFD Forum 0 February 23, 2002 14:18


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