|August 1, 2006, 10:25||
Dear All, I am trying to simu
Join Date: Mar 2009
Posts: 93Rep Power: 9
I am trying to simulate a "pathologic" (cause it has no physical meaning, and its purpose is just to test numerical squemes) scalarTransportFoam case, in which a unitary square box, 20 divisions in each edge, has prescribed (fixed) temperature, in all 4 edges: 0 in the bottom and left edge; -1 in the top edge; 1 in the right edge.
The velocity is uniform (1, 1, 0).
The initial temperature is also 0 everywhere.
termal diffusivity DT is set to 10^-10.
I am solving it as a transient problem with timestep = 0.005, with backward differences squeme.
One should expect that the downwind BC would not influence the interior domain, at least when Upwind is used. When Central Differences is used, we get the expected "checkered board" problem.
However, with QUICk and even with Upwind, I get an unbounded solution that grows indefinitely. In the upwind solution, the first interior cell-center, just upwindward of the downwind BC (top and right edges) grows beyond the bounding limits (-1, 1) as the time is advanced.
Is this something expected, that I am missing because of my ignorance?
I can post the case here if anyone wants to look at it.
|Thread||Thread Starter||Forum||Replies||Last Post|
|strange solver behavior---Treatment of Rough Wall||Yi||CFX||6||April 14, 2010 18:31|
|Strange incorrect inlet behavior for reactingFoam||mrangitschdowcom||OpenFOAM Running, Solving & CFD||4||August 29, 2007 15:55|
|very strange behavior of residuals||Mohammad El-Alti||CD-adapco||0||March 30, 2007 04:28|
|Dirichlet Pressure condition||Saud Khashan||FLUENT||0||June 25, 2001 09:58|
|Dirichlet B.C. for k-epsilon||Jens Wyrwa||FLUENT||0||December 18, 2000 06:18|