CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Dirichlet BC downind strange behavior

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

LinkBack Thread Tools Display Modes
Old   August 1, 2006, 10:25
Default Dear All, I am trying to simu
Join Date: Mar 2009
Posts: 93
Rep Power: 10
sampaio is on a distinguished road
Dear All,
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.

sampaio is offline   Reply With Quote


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
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 Siemens 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

All times are GMT -4. The time now is 19:59.