0% found this document useful (0 votes)
17 views17 pages

CFD FVM Techniques for Convection-Diffusion

Uploaded by

Muthamil Arasu
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
17 views17 pages

CFD FVM Techniques for Convection-Diffusion

Uploaded by

Muthamil Arasu
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

MCDM608L: COMPUTATIONAL

FLUID DYNAMICS (CFD)


[Link] (WS24-25)

Sivakumar, R
SMEC, VIT Chennai
Module-4

FVM for Convection-Diffusion Problems:


FVM for 1D steady state convection-diffusion,
Central differencing scheme, Conservativeness,
Boundedness, Transportiveness, Upward
differencing scheme, Hybrid differencing scheme for
2D convection-diffusion, Power-law scheme, QUICK
scheme. (6 hours)
Improved Formulation Schemes

 Upwind scheme
 Exponential scheme
 Hybrid scheme
 Power law scheme
 QUICK Scheme
Upwind Scheme
 Inadequacies of the central differencing scheme is its inability to
identify flow direction
 The value of property  at a west cell face is always influenced by
both P & W  Unsuitable
 Upwind differencing scheme/Donor cell differencing scheme takes
into account the flow direction when determining the value at a cell
face
 Formulation of diffusion term is unchanged
 The convected value of  at a cell face is taken to be equal to the
value at the upstream node
Upwind Scheme

 When the flow is in the positive direction

 When the flow is in the negative direction


Upwind Scheme

 No negative coefficient would arise. Solutions will always be


physically realistic
Upwind Scheme
Example:
 De = 1 ; Dw = 1; Fe = 4; Fw = 4
 E = 200; W = 100, P = ? (117)
 E = 100; W = 200, P = ? (183)
 Realistic values for P
 Satisfies the rules
 First order accurate since the scheme is based on the backward
differencing formula
 Produces erroneous results when the flow is not aligned with
the grid lines  FALSE DIFFUSION
Upwind Scheme
False DIffusion
 False diffusion is numerically introduced diffusion and arises in
convection dominated flows, i.e. high Pe number flows
 False diffusion will occur due to the oblique flow direction and
non-zero gradient of temperature in the direction normal to
the flow

Diffusion set to zero


k=0
Hot fluid
T = 100ºC

Cold fluid
T = 0ºC
Upwind Scheme
False DIffusion
 Grid refinement will minimize the false diffusion

First-order Upwind Second-order Upwind

8x8

64 x 64
Hybrid Scheme
 Drawback of upwind scheme: overestimate diffusion at large values
of Peclet number
 Combination of central and upwind differencing schemes
 Central differencing scheme: -2  Pe  2
 Upwind differencing scheme: Pe > 2
 Uses piecewise formula based on the local Peclet number to
evaluate the net flux through each control volume face
 The Peclet number is evaluated at the face of the control volume
Hybrid Scheme
Hybrid Scheme
 Scheme is conservative
 Coefficients are always positive  unconditionally bounded
 Satisfies the transportiveness requirement by using upwind
formulation for large value of Peclet number
 Produces physically realistic solutions
 First order accurate
Hybrid Scheme
 Conservativeness: Global conservation of the fluid
property  must be ensured.
 Boundedness: Values predicted by the scheme should be
within realistic bounds. For linear problems without
sources, those would be the maximum and minimum
boundary values. Fluid flow is non-linear and values in the
domain may be outside the range of boundary values.
 Transportiveness: Diffusion works in all directions but
convection only in the flow direction. The numerical
scheme should recognize the direction of the flow as it
affects the strength of convection versus diffusion.
Power law Scheme
 It is premature to set the diffusion effects equal to zero
when Pe > 2
 Power law scheme is more accurate approximation to the
one dimensional exact solution
 Produces better results than hybrid scheme
 Diffusion is set to zero when cell Pe > 10
 For 0 < Pe < 10, flux is evaluated by using a polynomial
expression
Power law Scheme
 The face value is determined from an exponential profile through
the cell values. The exponential profile is approximated by the
following power law equation:

1  0.1Pe 5 
e   P  E  P 
Pe

 First order accurate


QUICK Scheme
 Quadratic Upwind interpolation Scheme
 Interpolation through a parabola: three points necessary P, E and
point in upstream side

• g coefficients in terms of
nodal coordinates
• third order accurate
End of Module-4

335

Common questions

Powered by AI

The Hybrid Scheme combines central and upwind differencing schemes, utilizing a piecewise formula based on the local Peclet number evaluation. It is conservative, ensuring global conservation of the fluid property, and hits transportiveness by recognizing flow direction . This scheme is useful when Peclet numbers are not exceedingly high, granting physical realism and robustness . On the other hand, the Power Law Scheme excels in scenarios with Peclet numbers greater than 2, as it provides a more accurate approximation to the one-dimensional exact solution, particularly by setting diffusion effects to zero when cell Pe > 10 . It uses a polynomial expression for flux evaluation in the range 0 < Pe < 10, offering improved accuracy over the Hybrid Scheme for such conditions .

Transportiveness in numerical schemes refers to a scheme's ability to accurately model transport phenomena like convection and diffusion, particularly reflecting the dominant transport mechanisms in the direction of flow. It is important because diffusion acts in all directions in a fluid, while convection acts predominantly in the flow direction. In convection-dominated scenarios, adhering to this requirement ensures that numerical models accurately represent the correct transport properties without introducing non-physical artifacts like false diffusion and other inaccuracies . This helps in achieving realistic and reliable simulations of fluid behavior, aligning numerical predictions with actual physical behaviors .

Boundedness in computational fluid dynamics schemes refers to the principle that predicted values should remain within realistic physical bounds. This is essential for ensuring numerical solutions that are consistent with physical reality, avoiding non-physical behaviors in the solution domain . For linear problems without source terms, boundedness implies that values should not exceed the specified boundary values. In convection-dominated schemes, such as the Hybrid Scheme, boundedness ensures coefficients are positive and solutions remain stable and realistic, particularly when dealing with non-linear fluid flow where solution values might otherwise stray outside the expected range due to numerical inaccuracies .

The Peclet number is crucial in computational fluid dynamics to guide the choice or modification of differencing schemes based on flow conditions. For example, the Hybrid Scheme adapts its approach based on the Peclet number: it uses central differencing when -2 ≤ Pe ≤ 2, while resorting to the Upwind scheme when Pe > 2 to ensure realistic solutions by properly representing directional convective effects . Similarly, the Power Law Scheme alters diffusion modeling, setting diffusion to zero when the Peclet number exceeds 10 due to domination by convection . Thus, the Peclet number helps tailor the numerical approaches for convection-diffusion problems, ensuring stability and physical realism across varied flow regimes .

The QUICK Scheme is advantageous over the Upwind Scheme in scenarios demanding higher accuracy due to its third-order accuracy and ability to better capture flow behavior through quadratic upwind interpolation . It is particularly suitable for flows where higher-order precision is crucial due to complex geometries or when capturing small-scale phenomena, offering greater accuracy than the first-order Upwind Scheme would provide . Its polynomial approach uses three points (P, E, and upstream point), which enhances solution accuracy in capturing the curvature of the flow path over the simple direct approach of the Upwind Scheme .

The principal drawback of the central differencing scheme is its inability to accurately identify and account for flow direction, which makes it unsuitable for convection-dominant problems . It may provide oscillatory and non-physical solutions due to numerical instability when applied in high Peclet number scenarios . In contrast, schemes like the Upwind or Hybrid are more reliable in these situations as they incorporate directional flow information, leading to stable and physically realistic solutions . Whereas, the central differencing scheme could be suited for low Peclet number flows where diffusion dominates, but its application in convection-dominated flows is limited by potential inaccuracies and non-physical solutions .

The Power Law Scheme enhances accuracy in convection-dominated problems by accurately mimicking one-dimensional solutions across a range of Peclet numbers. It implements an exponential profile approximation to adjust the face value between nodes when 0 < Pe < 10, and sets diffusion to zero for Pe > 10, thus reducing the overestimation of diffusion inherent in other schemes like the Hybrid Scheme . This refined approach results in improved accuracy when dealing with varying degrees of convective dominance, ensuring the predictions closely align with physical expectations under high Peclet number conditions .

The Upwind differencing scheme accounts for flow direction by using values from the upstream node to influence the face value of a component in a convection-diffusion scenario, thus ensuring that solutions reflect the flow's true direction . This approach prevents negative coefficients and maintains physically realistic solutions. However, it is only first-order accurate, leading to potentially erroneous results when the flow is misaligned with grid lines due to false diffusion, especially in high Peclet number flows .

Grid refinement plays a crucial role in reducing false diffusion by aligning the computational grid more precisely with the flow streamlines, which minimizes numerical dispersion errors inherent in discretization methods like the Upwind Scheme . By increasing the grid's resolution, the effect of misalignment between grid lines and the flow direction is lessened, thereby reducing the introduction of artificial diffusion. This adjustment significantly enhances the fidelity and accuracy of the simulation results, ensuring that they are more reflective of the true physical phenomena . Grid refinement is thus an effective strategy to address numerical inaccuracies, providing clearer insights into the behavior of convection-dominated flows .

False diffusion is a numerical anomaly that introduces artificial or excess diffusion in computational models of convection-dominated flows, typically in cases of high Peclet number . It occurs when the numerical method improperly disperses the quantity being solved for, due to alignment issues between the flow direction and computational grid. This can lead to significant errors in predicted flow characteristics and malfunctioning of the numerical model, as the scheme might underestimate the convective transport effects . Refining the computational grid can minimize false diffusion, ensuring that the directions of grid lines align well with flow streamlines, thus reducing numerical diffusion artifacts and maintaining result accuracy .

You might also like