Physics benchmarks

Diffusion in Poiseuille Flow

Validation: 2D Finite Difference, Advection-Diffusion, Laminar flow

Description

Advection-diffusion of a scalar inside a 2D channel with a prescribed parabolic velocity profile, run to steady state. Half of the channel is modelled, with the centreline at y = 0.

The wall is held at scalar value T = 1 and the inlet at T = 0. The centreline and the outlet are zero gradient.

T= 1 T= 0 ux(y) ∂T ∂x = 0 ∂T ∂y = 0 x y
Figure 1: Flow domain: half channel

Simulation Inputs

InputValue
Material propertiesDiffusivity D = 7.143 × 10⁻² m²/s. Péclet number Pe = u0H/D = 27.65.
GeometryL = 2.0 m, H = 0.4 m.
Boundary ConditionsPrescribed velocity profile at inlet ux(y)=u0[1−(y/H)2] with centre-line velocity u0 = 4.937 m/s.

Assumptions and Notes

Flow is steady and laminar, with a fixed velocity field (no momentum solve). Profile is compared at x = 1.0 m (mid-length).

Results

Scalar field in the half channel
Figure 2: Scalar field

Scalar profile at x = 1.0 m vs Graetz solution

Hover or tap the chart to read the values at any point. Click a legend item to show or hide it. FlowDeck values are at cell centres; the lines are completed with the wall boundary values.

Figure 3: Scalar profile at x = 1.0 m vs Graetz solution

Scalar at x = 1.0 m: Graetz solution vs FlowDeck

y (m)Graetz TFlowDeck TAbs. error
0.001250.0864380.0830640.00337
0.003750.0865200.0831430.00338
0.003750.0865200.0836000.00292
0.006250.0866850.0833000.00339
0.006250.0866850.0837570.00293
0.008750.0869320.0835350.00340
0.008750.0869320.0839930.00294
0.011250.0872610.0838490.00341
0.011250.0872610.0843070.00295
0.013750.0876730.0842410.00343
0.013750.0876730.0847000.00297
0.016250.0881670.0847130.00345
0.016250.0881670.0851720.00300
0.018750.0887440.0852620.00348
0.018750.0887440.0857220.00302
0.021250.0894030.0858910.00351
0.021250.0894030.0863520.00305
0.023750.0901440.0865980.00355
0.023750.0901440.0870600.00308
0.026250.0909690.0873850.00358
0.026250.0909690.0878470.00312
0.028750.0918760.0882500.00363
0.028750.0918760.0887140.00316
0.031250.0928650.0891950.00367
0.031250.0928650.0896600.00321
0.033750.0939380.0902190.00372
0.033750.0939380.0906850.00325
0.036250.0950930.0913220.00377
0.036250.0950930.0917900.00330
0.038750.0963320.0925050.00383
0.038750.0963320.0929740.00336
0.041250.0976540.0937680.00389
0.041250.0976540.0942390.00341
0.043750.0990590.0951110.00395
0.043750.0990590.0955840.00348
0.046250.1005470.0965350.00401
0.046250.1005470.0970080.00354
0.048750.1021190.0980380.00408
0.048750.1021190.0985140.00361
0.051250.1037740.0996220.00415
0.051250.1037740.1000990.00367
0.053750.1055130.1012870.00423
0.053750.1055130.1017660.00375
0.056250.1073350.1030330.00430
0.056250.1073350.1035130.00382
0.058750.1092420.1048590.00438
0.058750.1092420.1053420.00390
0.061250.1112320.1067680.00446
0.061250.1112320.1072520.00398
0.063750.1133060.1087570.00455
0.063750.1133060.1092440.00406
0.066250.1154650.1108290.00464
0.066250.1154650.1113180.00415
0.068750.1177070.1129820.00473
0.068750.1177070.1134730.00423
0.071250.1200340.1152170.00482
0.071250.1200340.1157110.00432
0.073750.1224460.1175350.00491
0.073750.1224460.1180310.00442
0.076250.1249420.1199360.00501
0.076250.1249420.1204330.00451
0.078750.1275220.1224190.00510
0.078750.1275220.1229190.00460
0.081250.1301880.1249850.00520
0.081250.1301880.1254870.00470
0.083750.1329380.1276340.00530
0.083750.1329380.1281390.00480
0.086250.1357720.1303670.00541
0.086250.1357720.1308730.00490
0.088750.1386920.1331830.00551
0.088750.1386920.1336920.00500
0.091250.1416960.1360820.00561
0.091250.1416960.1365940.00510
0.093750.1447850.1390660.00572
0.093750.1447850.1395800.00521
0.096250.1479590.1421330.00583
0.096250.1479590.1426490.00531
0.098750.1512180.1452850.00593
0.098750.1512180.1458030.00541
0.101250.1545620.1485210.00604
0.101250.1545620.1490410.00552
0.103750.1579900.1518410.00615
0.103750.1579900.1523640.00563
0.106250.1615040.1552460.00626
0.106250.1615040.1557710.00573
0.108750.1651020.1587350.00637
0.108750.1651020.1592620.00584
0.111250.1687850.1623090.00648
0.111250.1687850.1628380.00595
0.113750.1725520.1659670.00658
0.113750.1725520.1664990.00605
0.116250.1764040.1697110.00669
0.116250.1764040.1702440.00616
0.118750.1803400.1735390.00680
0.118750.1803400.1740740.00627
0.121250.1843600.1774510.00691
0.121250.1843600.1779880.00637
0.123750.1884650.1814490.00702
0.123750.1884650.1819880.00648
0.126250.1926530.1855310.00712
0.126250.1926530.1860720.00658
0.128750.1969250.1896980.00723
0.128750.1969250.1902400.00668
0.131250.2012810.1939490.00733
0.131250.2012810.1944930.00679
0.133750.2057190.1982850.00743
0.133750.2057190.1988300.00689
0.136250.2102410.2027050.00754
0.136250.2102410.2032520.00699
0.138750.2148450.2072090.00764
0.138750.2148450.2077580.00709
0.141250.2195320.2117980.00773
0.141250.2195320.2123470.00718
0.143750.2243000.2164700.00783
0.143750.2243000.2170210.00728
0.146250.2291510.2212260.00792
0.146250.2291510.2217770.00737
0.148750.2340820.2260650.00802
0.148750.2340820.2266170.00746
0.151250.2390950.2309870.00811
0.151250.2390950.2315410.00755
0.153750.2441880.2359920.00820
0.153750.2441880.2365460.00764
0.156250.2493610.2410800.00828
0.156250.2493610.2416350.00773
0.158750.2546130.2462490.00836
0.158750.2546130.2468050.00781
0.161250.2599450.2515010.00844
0.161250.2599450.2520570.00789
0.163750.2653550.2568340.00852
0.163750.2653550.2573900.00797
0.166250.2708430.2622470.00860
0.166250.2708430.2628030.00804
0.168750.2764090.2677420.00867
0.168750.2764090.2682980.00811
0.171250.2820520.2733160.00874
0.171250.2820520.2738720.00818
0.173750.2877700.2789700.00880
0.173750.2877700.2795250.00825
0.176250.2935650.2847020.00886
0.176250.2935650.2852570.00831
0.178750.2994340.2905130.00892
0.178750.2994340.2910670.00837
0.181250.3053770.2964020.00898
0.181250.3053770.2969550.00842
0.183750.3113940.3023680.00903
0.183750.3113940.3029200.00847
0.186250.3174840.3084100.00907
0.186250.3174840.3089620.00852
0.188750.3236460.3145280.00912
0.188750.3236460.3150790.00857
0.191250.3298790.3207220.00916
0.191250.3298790.3212700.00861
0.193750.3361820.3269890.00919
0.193750.3361820.3275370.00865
0.196250.3425560.3333310.00922
0.196250.3425560.3338760.00868
0.198750.3489980.3397450.00925
0.198750.3489980.3402880.00871
0.201250.3555070.3462310.00928
0.201250.3555070.3467720.00873
0.203750.3620840.3527880.00930
0.203750.3620840.3533280.00876
0.206250.3687270.3594160.00931
0.206250.3687270.3599530.00877
0.208750.3754350.3661130.00932
0.208750.3754350.3666470.00879
0.211250.3822080.3728790.00933
0.211250.3822080.3734100.00880
0.213750.3890430.3797120.00933
0.213750.3890430.3802410.00880
0.216250.3959410.3866120.00933
0.216250.3959410.3871380.00880
0.218750.4029000.3935770.00932
0.218750.4029000.3941000.00880
0.221250.4099200.4006070.00931
0.221250.4099200.4011270.00879
0.223750.4169990.4077010.00930
0.223750.4169990.4082170.00878
0.226250.4241360.4148580.00928
0.226250.4241360.4153700.00877
0.228750.4313300.4220760.00925
0.228750.4313300.4225840.00875
0.231250.4385800.4293540.00923
0.231250.4385800.4298590.00872
0.233750.4458850.4366920.00919
0.233750.4458850.4371930.00869
0.236250.4532440.4440890.00916
0.236250.4532440.4445850.00866
0.238750.4606560.4515420.00911
0.238750.4606560.4520340.00862
0.241250.4681200.4590520.00907
0.241250.4681200.4595390.00858
0.243750.4756340.4666160.00902
0.243750.4756340.4670990.00854
0.246250.4831980.4742350.00896
0.246250.4831980.4747120.00849
0.248750.4908100.4819060.00890
0.248750.4908100.4823780.00843
0.251250.4984690.4896280.00884
0.251250.4984690.4900950.00837
0.253750.5061740.4974000.00877
0.253750.5061740.4978620.00831
0.256250.5139240.5052220.00870
0.256250.5139240.5056780.00825
0.258750.5217170.5130910.00863
0.258750.5217170.5135420.00818
0.261250.5295530.5210070.00855
0.261250.5295530.5214520.00810
0.263750.5374310.5289690.00846
0.263750.5374310.5294080.00802
0.266250.5453480.5369740.00837
0.266250.5453480.5374070.00794
0.268750.5533050.5450230.00828
0.268750.5533050.5454500.00786
0.271250.5613000.5531130.00819
0.271250.5613000.5535340.00777
0.273750.5693310.5612440.00809
0.273750.5693310.5616580.00767
0.276250.5773980.5694140.00798
0.276250.5773980.5698220.00758
0.278750.5855000.5776230.00788
0.278750.5855000.5780230.00748
0.281250.5936340.5858680.00777
0.281250.5936340.5862620.00737
0.283750.6018010.5941490.00765
0.283750.6018010.5945360.00727
0.286250.6099990.6024640.00754
0.286250.6099990.6028440.00715
0.288750.6182270.6108120.00741
0.288750.6182270.6111850.00704
0.291250.6264830.6191930.00729
0.291250.6264830.6195590.00692
0.293750.6347670.6276040.00716
0.293750.6347670.6279630.00680
0.296250.6430780.6360460.00703
0.296250.6430780.6363970.00668
0.298750.6514150.6445150.00690
0.298750.6514150.6448590.00656
0.301250.6597750.6530120.00676
0.301250.6597750.6533480.00643
0.303750.6681600.6615360.00662
0.303750.6681600.6618640.00630
0.306250.6765660.6700840.00648
0.306250.6765660.6704050.00616
0.308750.6849940.6786560.00634
0.308750.6849940.6789690.00602
0.311250.6934420.6872510.00619
0.311250.6934420.6875570.00589
0.313750.7019100.6958680.00604
0.313750.7019100.6961660.00574
0.316250.7103950.7045060.00589
0.316250.7103950.7047950.00560
0.318750.7188990.7131640.00574
0.318750.7188990.7134450.00545
0.321250.7274180.7218400.00558
0.321250.7274180.7221130.00531
0.323750.7359540.7305340.00542
0.323750.7359540.7307980.00516
0.326250.7445030.7392440.00526
0.326250.7445030.7395000.00500
0.328750.7530670.7479700.00510
0.328750.7530670.7482180.00485
0.331250.7616430.7567110.00493
0.331250.7616430.7569510.00469
0.333750.7702320.7654650.00477
0.333750.7702320.7656970.00454
0.336250.7788320.7742330.00460
0.336250.7788320.7744560.00438
0.338750.7874420.7830120.00443
0.338750.7874420.7832270.00422
0.341250.7960620.7918030.00426
0.341250.7960620.7920090.00405
0.343750.8046910.8006040.00409
0.343750.8046910.8008020.00389
0.346250.8133280.8094150.00391
0.346250.8133280.8096040.00372
0.348750.8219730.8182350.00374
0.348750.8219730.8184150.00356
0.351250.8306250.8270620.00356
0.351250.8306250.8272340.00339
0.353750.8392820.8358970.00339
0.353750.8392820.8360600.00322
0.356250.8479460.8447390.00321
0.356250.8479460.8448930.00305
0.358750.8566150.8535860.00303
0.358750.8566150.8537320.00288
0.361250.8652880.8624390.00285
0.361250.8652880.8625760.00271
0.363750.8739650.8712970.00267
0.363750.8739650.8714250.00254
0.366250.8826460.8801590.00249
0.366250.8826460.8802780.00237
0.368750.8913300.8890250.00230
0.368750.8913300.8891350.00219
0.371250.9000160.8978940.00212
0.371250.9000160.8979960.00202
0.373750.9087050.9067660.00194
0.373750.9087050.9068590.00185
0.376250.9173960.9156400.00176
0.376250.9173960.9157240.00167
0.378750.9260880.9245160.00157
0.378750.9260880.9245910.00150
0.381250.9347820.9333940.00139
0.381250.9347820.9334600.00132
0.383750.9434760.9422730.00120
0.383750.9434760.9423300.00115
0.386250.9521710.9511520.00102
0.386250.9521710.9512010.00097
0.388750.9608670.9600330.00083
0.388750.9608670.9600730.00079
0.391250.9695630.9689140.00065
0.391250.9695630.9689450.00062
0.393750.9782590.9777960.00046
0.393750.9782590.9778180.00044
0.396250.9869550.9866770.00028
0.396250.9869550.9866910.00026
0.398750.9956520.9955590.00009
0.398750.9956520.9955640.00009

Reference Data

The reference is the classical Graetz series [1], taking the first 10 eigenfunctions, with coefficients Tm found through integrals of the hypergeometric functions ₁F₁ [2]:

T(x,y)=1−∑m=010Tme−m4x/(HPe)e−(my/2H)2F11(−m24+14,12,m2y2H2)

References

  1. L. Graetz, “Ueber die Wärmeleitungsfähigkeit von Flüssigkeiten,” Annalen der Physik 254(1), 79–94, 1883. doi:10.1002/andp.18832540106
  2. R. K. Shah and A. L. London, Laminar Flow Forced Convection in Ducts, Advances in Heat Transfer, Supplement 1, Academic Press, New York, 1978.

Data licence. The model results on this page and in the CSV downloads are published by FlowDeck under the Creative Commons Attribution 4.0 licence (CC BY 4.0). The reference values are from the sources listed under References.