CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Why do we need deferred correction when it can't help yield a valid result?

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 5, 2023, 20:36
Default
  #21
New Member
 
Jasmine Feng
Join Date: Feb 2023
Posts: 14
Rep Power: 3
JasmineFeng is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
One of the fundamental property of a convection-diffusion PDE is that it does not create new extrema in the solution. Therefore, it is clear that in a numerical solution the presence of new extrema is not physically meaningful. Monotone schemes have a property to guarantee that.

But, if a source term is added in the PDE, the action of such a term can actually produce new physical extrema.
I see. A positive (negative) source can create a local maximum (minimum), in which case the values of internal nodes will exceed those of the boundary nodes, namely unbounded results.
JasmineFeng is offline   Reply With Quote

Old   February 5, 2023, 20:45
Default
  #22
New Member
 
Jasmine Feng
Join Date: Feb 2023
Posts: 14
Rep Power: 3
JasmineFeng is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I don't understand what you are doing. Negative Peclet means an anti-diffusive problem.

Are you just trying to replicate the 1D stready convection-diffusion exercise shown in the textbook of Ferzigr, Peric and Street?
I don't think diffusion will change with the sign of the grid Peclect number.
Pe=ρ*u*Δx/Γ
So Pe > 0 denotes u>0 and Pe<0 denotes u<0, but the diffusive coefficient Γ keeps the same.

Yeah, this is a 1D steady convection-diffusion exercise without any source term. But from Versteeg and Malalasekera's textbook instead of Ferziger et al's.
JasmineFeng is offline   Reply With Quote

Old   February 6, 2023, 00:41
Default
  #23
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,286
Rep Power: 34
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by JasmineFeng View Post
I got it. Because schemes like SOU and MUSCL are unconditionally stable, so they will always produce bounded results without wiggles.


Yes. If you got stable and good results with CD on those peclet numbers you would have broken known mathematics. When we use unbounded CD all the solvers produce wiggle in solution or just diverge. This is why we have bounded central differencing schemes in solvers. This scheme blends enough upwind scheme to the make it stable.

Last edited by arjun; February 6, 2023 at 03:33.
arjun is offline   Reply With Quote

Old   February 6, 2023, 03:46
Default Ceribena Noche
  #24
New Member
 
ceribenanoche
Join Date: Feb 2023
Posts: 1
Rep Power: 0
ceribenanoche is on a distinguished road
Check out the latest Ceribena Noche results as https://ceribenanoche.com/ we provide the latest Ceribena Noche results daily.
ceribenanoche is offline   Reply With Quote

Reply

Tags
deferred correction, iterative method


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
error to check simulation result: could not find valid OpenFOAM mesh xiaoyoyo OpenFOAM 5 May 18, 2017 12:38
[snappyHexMesh] How to define to right point for locationInMesh Mirage12 OpenFOAM Meshing & Mesh Conversion 7 March 13, 2016 15:07
[snappyHexMesh] determining displacement for added points CFDnewbie147 OpenFOAM Meshing & Mesh Conversion 1 October 22, 2013 10:53
Timestep and Pressure Correction Relationship in SIMPLE rks171 Main CFD Forum 23 May 4, 2012 02:04
OF 1.6 | Ubuntu 9.10 (64bit) | GLIBCXX_3.4.11 not found piprus OpenFOAM Installation 22 February 25, 2010 14:43


All times are GMT -4. The time now is 05:30.