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

Hyperbolic equasion solver or How OF works

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2013, 09:07
Question Hyperbolic equasion solver or How OF works?
  #1
Member
 
Avdeev Evgeniy
Join Date: Jan 2011
Location: Togliatty, Russia
Posts: 52
Blog Entries: 1
Rep Power: 12
j-avdeev will become famous soon enough
Send a message via Skype™ to j-avdeev
Hello,
I have working hyperbolic equasion solver (hyper1Foam in attachment - works in OF 2.1.0).

Main equation:
\frac{\partial\psi }{\partial t}+C{\frac{\partial \psi}{\partial x}}=0
accordingly
Code:
        solve
        (
	    fvm::ddt(psi)
	    +
	    fvm::div(phi,psi)
        );
and in fvShemes

Code:
ddtSchemes
{
    ddt(psi) Euler;
}
divSchemes
{
    div(phi,psi) Gauss linear;
}
Working geometry - coube 1x1x1 diveded on 5 equal parts along oZ axe.

Main equation in discretized form:

[LaTeX Error: Syntax error]
C_i=\frac{V\phi_i}{\Delta t}
where
\psi_i - old ~value ~of ~phi ~in ~cell ~i
\psi^n_i - new ~value

In my case
V = const = 1/5 = 0.2 ~~ volume ~of ~single ~cell
\Delta t = 0.1 ~time ~step

And when I insert result of phi from one of time directories...

for ~\phi_i = 0,434805
\frac{0.2 *0,749653 - 0.2*0,434805}{0.1 } + (0,2/0,1)*0,434805*{\frac{0,206803 - 1}{0.2}} =
=-2,81916421585

Not zero!
I khow that hyperbolic equation will not =0, because it is iterational method. Is it residual?

I want to understand - how OpenFOAM count values of phi from one iteration to another - how to get new phi from old values:
phi_new = ... ... ...

And what initial values of phi?

Results: Table of results on GoogleDocs
Attached Images
File Type: png hyperFoamMesh.png (6.6 KB, 18 views)
File Type: jpg Results.jpg (41.6 KB, 18 views)
File Type: png equation.png (20.2 KB, 16 views)
Attached Files
File Type: zip hyper1Foam.zip (8.3 KB, 5 views)
File Type: zip case.zip (7.9 KB, 4 views)

Last edited by j-avdeev; November 4, 2013 at 08:24.
j-avdeev is offline   Reply With Quote

Old   October 30, 2013, 09:12
Default
  #2
Member
 
Avdeev Evgeniy
Join Date: Jan 2011
Location: Togliatty, Russia
Posts: 52
Blog Entries: 1
Rep Power: 12
j-avdeev will become famous soon enough
Send a message via Skype™ to j-avdeev
second file of equatins
j-avdeev 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
smoothSolver diverges - solution in using PBiCG solver? makaveli_lcf OpenFOAM Running, Solving & CFD 3 September 11, 2013 12:44
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 08:45
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
Error during Solver cfd guy CFX 4 May 8, 2001 06:04


All times are GMT -4. The time now is 00:53.