|
[Sponsors] |
June 29, 2010, 12:22 |
Crank Nicholson scheme for Laplacian failed
|
#1 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 |
Hi!
I modify the laplacian discretization scheme into the scalarTransportFoam solver as suggested in the programmer's guide P-43 to implement the crank Nicholson scheme : Code:
solve ( fvm::ddt(rho,T) + fvm::div(phi,T) - 0.5*(fvm::laplacian(rho*DT,T)+fvc::laplacian(rho*DT,T)) == ff ); In my fvSchemes I specify : Code:
laplacianSchemes { default none; laplacian ((rho*DT),T) Gauss linear corrected; } Regards, Cyp |
|
July 16, 2010, 12:12 |
|
#2 |
New Member
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16 |
Try moving
Code:
0.5*(fvm::laplacian(rho*DT,T)+fvc::laplacian(rho*DT,T)) |
|
July 19, 2010, 08:36 |
|
#3 |
Senior Member
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 299
Rep Power: 18 |
I tried your modification but I still have my problem...
How is your fvSchemes file ?? |
|
July 19, 2010, 08:58 |
|
#4 |
New Member
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16 |
The relevant part of my fvSchemes dictionary:
Code:
laplacianSchemes { default Gauss linear corrected; } Code:
ddtSchemes { default CrankNicholson 1; } Mike |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
About the Use of Crank Nicholson scheme | And | OpenFOAM | 1 | November 14, 2017 13:25 |
Transient boundary conditions | Jarrod Sinclair (Sinclair) | OpenFOAM Running, Solving & CFD | 133 | May 6, 2015 07:21 |
user subroutine error | CFDUSER | CFX | 2 | December 9, 2006 07:31 |
user defined function | cfduser | CFX | 0 | April 29, 2006 11:58 |
Crank Nicholson | noName | Main CFD Forum | 3 | September 15, 2004 10:10 |