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

Edit linear solver

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

Like Tree4Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 23, 2013, 12:38
Default Matrix in OpenFoam
  #21
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Hello Mr. Bruno

According to the reference that it's address link is below, in page 53;
http://www.diva-portal.org/smash/get...T01.pdf‎
I can,t undrestand that what is the placeholder matrix in openFoam? and how get this equations:

u(0)=[2,0] u(1)=[3,1] l(0)=[0,2] l(1)=[1,3]

UpperAddr=(2,3) LowerAddr=(0,1)

If it is possible for you, please explain a little about this concept.
Best regards,
Hassan
luckycfd is offline   Reply With Quote

Old   August 24, 2013, 20:14
Default
  #22
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,691
Blog Entries: 33
Rep Power: 75
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Hi Hassan,

Sorry, but this is beyond my experience level.
My suggestion is that you'll have to search for more information on this topic, possibly among the reference thesis made around OpenFOAM and the OpenFOAM Programmers Guide: http://foam.sourceforge.net/docs/Gui...mmersGuide.pdf

Good luck! Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   August 25, 2013, 02:13
Default
  #23
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
thank you Mr. Bruno
luckycfd is offline   Reply With Quote

Old   August 26, 2013, 04:12
Default Ueqn
  #24
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Hi,
In the icoFoam solver, where is the top level of Ueqn where Navier-stoke equation carries out?

best regards,
Hassan
luckycfd is offline   Reply With Quote

Old   August 26, 2013, 18:12
Default
  #25
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,691
Blog Entries: 33
Rep Power: 75
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Hi Hassan,

The best I can do on this question is to point you to this wiki page: http://openfoamwiki.net/index.php/IcoFoam - it basically explains how the icoFoam solver works.

Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   September 9, 2013, 15:25
Default
  #26
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Hi friends,
I can't understand the PCG linear solver how work in the icoFoam solver? In other words, which part of the icoFoam.C related to the PCG linear solver and use it? please show me.

Regards,
Hassan
luckycfd is offline   Reply With Quote

Old   September 9, 2013, 15:38
Default
  #27
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,691
Blog Entries: 33
Rep Power: 75
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Hi Hassan,

The PCG solver is not available directly to icoFoam. It's used very deep into the solve() methods and equation solving algorithms.

Some more reading material:
Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   September 9, 2013, 16:06
Default
  #28
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Thanks.
Regards,
Hassan.

Last edited by luckycfd; September 20, 2013 at 15:08.
luckycfd is offline   Reply With Quote

Old   October 10, 2013, 14:34
Default
  #29
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Dear Foamers,
How can I extract the linear system (AX=B) from the cavity case in tutorials using OpenFOAM?


Regards,
Hassan
luckycfd is offline   Reply With Quote

Old   December 9, 2013, 02:29
Default extract the linear system
  #30
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Hi dear Bruno,
Could I use the PCG linear solver for my linear system (Ax=b) separately? In the other words, I want to send two input files A and b, then get the x from this solver? yes or no please?

regards,
Hassan
luckycfd is offline   Reply With Quote

Old   December 28, 2013, 15:22
Default
  #31
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 7,691
Blog Entries: 33
Rep Power: 75
wyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the roughwyldckat is a jewel in the rough
Hi Hassan,

I'm sorry to say that I don't know the answer. In theory, it should be possible.
AFAIK, do keep in mind that PCG will not be able to give you an exact solution. At most, it's able to give you a partial/possible near solution, but not an exact one: http://www.cfd-online.com/Wiki/Conju...adient_methods

If I had to start investigating into this, I would have a look into the source code experiments that try to use one or more GPUs to help solve simulations faster, such as:
I say this because these libraries will sort-of hack directly into the LDU matrix management/solving system.


Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   January 15, 2014, 08:42
Default
  #32
New Member
 
Meher
Join Date: Jan 2014
Posts: 3
Rep Power: 3
ssmeherk is on a distinguished road
Quote:
Originally Posted by luckycfd View Post
Hi dear Bruno,
Could I use the PCG linear solver for my linear system (Ax=b) separately? In the other words, I want to send two input files A and b, then get the x from this solver? yes or no please?

regards,
Hassan
Hi Hassan,

I am also working on the same problem of solving sparse linear matrix (AX=B) using PCG solver.
Did you find any solution ? Any useful inputs from your side apart from what is available here?

Thanks,
Meher.
ssmeherk is offline   Reply With Quote

Old   January 24, 2014, 14:28
Default
  #33
Member
 
Hassan
Join Date: Oct 2012
Location: Iran
Posts: 60
Rep Power: 4
luckycfd is on a distinguished road
Hi Meher,
thanks for your attention, but no.
luckycfd is offline   Reply With Quote

Reply

Tags
compile, linear solver, openfoam

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
How to avoid checkerboard? dohnie OpenFOAM Programming & Development 6 January 23, 2015 12:18
Fatal overflow in linear solver (RSM) Luis CFX 9 May 24, 2013 07:19
ERROR #004100018; Fatal overflow in linear solver Attila CFX 1 April 13, 2012 23:22
ParaView Compilation jakaranda OpenFOAM Installation 3 October 27, 2008 13:46
free C code for large sparse matrix linear solver ztdep Main CFD Forum 7 May 24, 2007 15:14


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