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

lduMatrix problem in High order scheme

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

Like Tree1Likes
  • 1 Post By santiagomarquezd

Reply
 
LinkBack Thread Tools Display Modes
Old   March 27, 2013, 22:41
Default lduMatrix problem in High order scheme
  #1
New Member
 
ren XG
Join Date: Sep 2012
Posts: 10
Rep Power: 4
rxgrch is on a distinguished road
Hi foamers:

what I am doing is to implement a high order scheme for the FVM's face interplation, but now I face a problem with the lduMatrix, as we know, the matrix values just store for the face shared by two cells, and indexed by face, which defined by fvMeshLduAddressing. to implement my high order face interplation, the value not shared by two adjacent cells also had to be stored in the lduMtrix, which result in the modification of fvMeshLduAddressing. however, I found this modify of fvMeshLduAddressing may trigger a chain reaction in the code, which I had to modify the existing low order schemes, a heavy task. So can anyone help me with this problem, and point out if there are something wrong with my implement, or some advices? Thank you very much!
rxgrch is offline   Reply With Quote

Old   April 1, 2013, 11:07
Default
  #2
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Hi, I would suggest to read this material:

https://docs.google.com/file/d/1Y0cX...it?usp=sharing

http://www.sciencedirect.com/science...17931002003307

in order to check if is it possible to write your scheme in the FOAM's TVD framework. In addition take into account that schemes with double upwind cell information cannot be implemented in non-structured meshes in general since this cell cannot be always univocally determined.

Regards.
fumiya likes this.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd is offline   Reply With Quote

Old   April 2, 2013, 03:22
Default
  #3
New Member
 
ren XG
Join Date: Sep 2012
Posts: 10
Rep Power: 4
rxgrch is on a distinguished road
Quote:
Originally Posted by santiagomarquezd View Post
Hi, I would suggest to read this material:

https://docs.google.com/file/d/1Y0cX...it?usp=sharing

http://www.sciencedirect.com/science...17931002003307

in order to check if is it possible to write your scheme in the FOAM's TVD framework. In addition take into account that schemes with double upwind cell information cannot be implemented in non-structured meshes in general since this cell cannot be always univocally determined.

Regards.
Thanks for your reply, and all my work is now under the structured mesh, your document and paper seemed to be of great value for my work. it seemed that the four TVD schemes you mentioned in the paper still used just two cells C and D in Fig 2(b), the cell U actually isn't used, right? Did all the four sheme had been implement in the current OpenFOAM version? Or just your in-house code? Thank you!
rxgrch is offline   Reply With Quote

Old   April 2, 2013, 10:30
Default
  #4
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 14
santiagomarquezd will become famous soon enough
Hi, even if you are using structured meshes it is possible to remove the dependency of the U cell (take for example the QUICK limiter) and write the method in terms of C and P (or P and N in FOAM nomenclature), using (16) in the work of Darwish. This work explains the way TVD is implemented in FOAM. I think you would have to try this way first before changing the lduAddressing.

Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D.
Post-doctoral Fellow
Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL
T.E.: 54-342-4511594 Ext. 1005
Güemes 3450 - (3000) Santa Fe
Santa Fe - Argentina
http://www.cimec.org.ar
santiagomarquezd 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
Solution scheme: 1st, 2nd order, etc in momentum eqs nottern Main CFD Forum 1 May 11, 2012 00:30
Using Runge kutta 4th order scheme in FSI using CFX vmlxb6 CFX 0 April 26, 2011 14:55
Fifth order anti diffusive scheme for 1D hyperbolic heat transfer areffallah Main CFD Forum 0 November 3, 2010 05:48
High order FVM scheme .... cfd_newbie Main CFD Forum 2 February 6, 2008 07:46
Discretisation scheme in CFX-1st or 2nd order? Pete CFX 10 January 12, 2005 13:48


All times are GMT -4. The time now is 13:07.