diff --git a/validationCases/03_freeConvection/README.md b/validationCases/03_freeConvection/README.md
index c8b3596..e77f95f 100644
--- a/validationCases/03_freeConvection/README.md
+++ b/validationCases/03_freeConvection/README.md
@@ -3,4 +3,10 @@
 Two parallel, laminar gas flows split by a stainless steel plate. 
 Initially, hot gas (top) stream is at 1800 K, cold gas (bottom) is at 300 K, and plate is at 800 K.
 
+To demonstrate verification of the conservation equations, the integrated heat fluxes at the two fluid-solid interfaces converge to equal and opposite at steady-state.
+For validation, we examine the heat flux at the point location x/L = 0.98 near the trailing edge of the plate.
+The numerical convective heat flux from the solver is compared with the the theoretical heat flux calculated by the following method.
+
+1. Laminar correlations are used to calculate the local heat transfer coefficient (HTC) based on film properties. 
+2. The local heat flux is calculated as q = h(T\_{\infty} - T\_{wall}).
 
diff --git a/validationCases/03_freeConvection/results/integratedWallHeatFlux.pdf b/validationCases/03_freeConvection/results/integratedWallHeatFlux.pdf
new file mode 100644
index 0000000..3d374af
Binary files /dev/null and b/validationCases/03_freeConvection/results/integratedWallHeatFlux.pdf differ
diff --git a/validationCases/03_freeConvection/results/wallHeatFlux_xbyL0.98.pdf b/validationCases/03_freeConvection/results/wallHeatFlux_xbyL0.98.pdf
new file mode 100644
index 0000000..bd5f458
Binary files /dev/null and b/validationCases/03_freeConvection/results/wallHeatFlux_xbyL0.98.pdf differ