diff --git a/doc/changelog.d/4830.fixed.md b/doc/changelog.d/4830.fixed.md
new file mode 100644
index 000000000000..70f643540f3b
--- /dev/null
+++ b/doc/changelog.d/4830.fixed.md
@@ -0,0 +1 @@
+Modeling solidification workflow and image
diff --git a/doc/source/_static/modeling_solidification_5.png b/doc/source/_static/modeling_solidification_5.png
index f5c825b956c5..fd9a25c50d24 100644
Binary files a/doc/source/_static/modeling_solidification_5.png and b/doc/source/_static/modeling_solidification_5.png differ
diff --git a/examples/00-fluent/Modeling_solidification_workflow.py b/examples/00-fluent/Modeling_solidification_workflow.py
index 4361bd1b6c98..b9d5b352b7c2 100644
--- a/examples/00-fluent/Modeling_solidification_workflow.py
+++ b/examples/00-fluent/Modeling_solidification_workflow.py
@@ -93,6 +93,7 @@
     Methods,
     Results,
     RunCalculation,
+    Scene,
     Setup,
     Solution,
     VelocityInlet,
@@ -344,7 +345,10 @@
     "option": "auto-range-off",
     "auto_range_off": {"clip_to_range": True, "minimum": 1100, "maximum": 1200},
 }
-mushy_temp.display()
+scene = Scene(solver, new_instance_name="scene-1")
+scene.graphics_objects.add(name="mesh-1")
+scene.graphics_objects.add(name="temperature-mushy")
+scene.display()
 
 graphics.views.restore_view(view_name="front")
 graphics.picture.save_picture(file_name="modeling_solidification_3.png")
@@ -362,6 +366,7 @@
 # %%
 # Enable transient flow and heat transfer
 # ---------------------------------------
+solver.settings.file.read_case_data(file_name="steady_state")
 solver_general_settings.solver.time = "unsteady-1st-order"  # First-order implicit
 
 controls.equations = {"flow": True, "w-swirl": True}
@@ -396,7 +401,9 @@
 solutions.run_calculation.calculate()
 
 # Liquid fraction at t = 5.0 s
-liquid_fraction_contour_t_5_sec = Contour(solver, new_instance_name="liquid-fraction")
+liquid_fraction_contour_t_5_sec = Contour(
+    solver, new_instance_name="liquid-fraction_t_5_sec"
+)
 
 liquid_fraction_contour_t_5_sec.field = "liquid-fraction"
 liquid_fraction_contour_t_5_sec.display()