📐 Grid Peclet Number Solver
Calculate local grid Peclet numbers to evaluate cell convection-diffusion numerical scheme limits.
⚡ Fortran 90 Engine
Double Precision (IEEE 754)
✓ ISO / ASME Validated
Cfd
📊 Solver Telemetry
● ACTIVE
📥 Flow Velocity & Cell Geometry
📖 Numerical Boundedness (Patankar, 1980):
$$\text{Pe} = \frac{\rho u L}{\Gamma} = \frac{u L}{\alpha}$$
$$\text{CDS Boundedness Limit: } |\text{Pe}| \le 2.0$$
$$\text{Exact 1D Solution: } \phi(x) = \frac{e^{\text{Pe}\, x/L} - 1}{e^{\text{Pe}} - 1}$$
15.000
Cell Péclet ($\text{Pe}_\Delta$)
Oscillatory ⚠️
CDS Stability (|Pe| ≤ 2)
6.67e-5 m
Max $\Delta x$ for CDS
Advection
Dominant Regime
💡
Recommended FVM Discretization: Upwind Differencing (UDS) / High-Resolution TVD
Transport is overwhelmingly advection-dominated. Upwind-based flux limiting or TVD limiter is required to enforce strict physical monotonicity.
📈 1D Scalar Advection-Diffusion Profile: $\phi(x)$
$\text{Pe} = 15.00$🔍 View Raw GNU Fortran Double-Precision Solver Output
MODE=1 MODE_NAME=Single Cell VELOCITY= 3.00000000E+00 CELL_SIZE= 5.00000000E-04 ALPHA= 1.00000000E-04 RHO= 1.22500000E+00 GAMMA= 0.00000000E+00 PE= 15.000000 PE_ABS= 15.000000 DIAGNOSIS=Strongly convective RECOMMENDATION=UDS or Power-Law; refine mesh to lower Pe
📘 Calculation Methodology: Grid Peclet Number & Convective Stability
Mathematical Model & Theory
The grid Peclet number $Pe_{grid}$ compares convective transport to diffusive transport across a single grid cell. If $Pe_{grid} > 2$, central differencing produces unphysical wiggles:
$$Pe_{grid} = \frac{\rho u \Delta x}{\Gamma} = \frac{u \Delta x}{\nu}$$
$$\text{Central Difference Boundedness Condition: } |Pe_{grid}| \le 2$$
Assumptions
- Finite volume advection-diffusion transport.
- Upwind differencing (UDS / QUICK / MUSCL) required when $Pe > 2$.
Academic References
- Patankar, S. V.: Numerical Heat Transfer and Fluid Flow.
- Versteeg, H. K., & Malalasekera, W.: An Introduction to CFD.
Worked Engineering Example
Problem Statement:
Air ($\nu = 1.5\times 10^{-5}\text{ m}^2/\text{s}$) flows at $u = 3.0\text{ m/s}$ with grid spacing $\Delta x = 20\text{ mm}$. Calculate $Pe_{grid}$.
Step-by-step Solution:
1. $Pe_{grid} = (3.0 \times 0.02) / (1.5 \times 10^{-5}) = 4000 \gg 2$.
Final Result:
$Pe_{grid} = \mathbf{4000}$ (upwind differencing required for stability).
Air ($\nu = 1.5\times 10^{-5}\text{ m}^2/\text{s}$) flows at $u = 3.0\text{ m/s}$ with grid spacing $\Delta x = 20\text{ mm}$. Calculate $Pe_{grid}$.
Step-by-step Solution:
1. $Pe_{grid} = (3.0 \times 0.02) / (1.5 \times 10^{-5}) = 4000 \gg 2$.
Final Result:
$Pe_{grid} = \mathbf{4000}$ (upwind differencing required for stability).