|
[Sponsors] |
March 12, 2016, 00:40 |
Sysyem of linear equations
|
#1 |
Member
Soyol
Join Date: Jun 2014
Location: Berkeley, CA
Posts: 39
Rep Power: 12 |
Hello everybody,
I am using finite difference to solve linear internal wave 2d equations. I discretize the equation but at the end the linear system of equations do not come to a well format to be able to code it and solve it. Last edited by soes; March 12, 2016 at 13:58. |
|
March 12, 2016, 04:34 |
|
#2 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
Sorry, but it is very hard to read from your pages what you did...
It seems a second order PDE, is it hyperbolic? You can then reduce to a system of two equations of first order and solve with classical methods you can find in many textbook. For example http://www.cambridge.org/se/academic...on-2nd-edition |
|
March 12, 2016, 14:06 |
|
#3 | |
Member
Soyol
Join Date: Jun 2014
Location: Berkeley, CA
Posts: 39
Rep Power: 12 |
Quote:
Thanks in advance |
||
March 12, 2016, 14:37 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,897
Rep Power: 73 |
but Wxx and Wzz are second derivatives in space? Therefore, you wrote a fourth-order PDE ....
I suggest to decompose in a lower order system of equations to analyse the nature of the PDE. Then you can simply apply any of the method in the book I linked. |
|
March 12, 2016, 14:43 |
|
#5 | |
Member
Soyol
Join Date: Jun 2014
Location: Berkeley, CA
Posts: 39
Rep Power: 12 |
Quote:
Yes they are spatial derivatives which go into time derivative (i.e. laplacian of time and space both). But what do you mean exactly by decompose to lower order system? Is my discretization wrong? I cannot proceed with discretization of the equation without modifying it? |
||
March 14, 2016, 14:47 |
|
#6 |
Senior Member
duri
Join Date: May 2010
Posts: 245
Rep Power: 17 |
The LHS in the final equation has Wxx. I think it is typographical error, it should be W. Then the given equation is correct. By the way, what equation is this?. It looks very complicated.
|
|
March 18, 2016, 01:56 |
|
#7 |
Member
Soyol
Join Date: Jun 2014
Location: Berkeley, CA
Posts: 39
Rep Power: 12 |
No it is unfortunately correct
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Heat flux BC | lfrigeri3 | OpenFOAM | 9 | February 22, 2014 12:59 |
Oscillating behaviour when mesh gets too small | james91 | OpenFOAM Running, Solving & CFD | 0 | January 9, 2014 06:04 |
Solving Linear System of Equations in UDF | mohsen zendehbad | FLUENT | 5 | May 8, 2010 11:39 |
Solving Linear System of Equations in UDF | mohsen zendehbad | FLUENT | 0 | April 25, 2010 02:23 |
IdeasUnvToFoam Bug amp Fix | benru | OpenFOAM Bugs | 42 | November 13, 2009 08:59 |