diff --git a/.ci/start_fluent.py b/.ci/start_fluent.py
new file mode 100644
index 00000000000..26650d19ddf
--- /dev/null
+++ b/.ci/start_fluent.py
@@ -0,0 +1,42 @@
+import os
+import subprocess
+import sys
+import tempfile
+import time
+
+from ansys.fluent.core import EXAMPLES_PATH
+
+
+def start_fluent_container(args):
+    fd, sifile = tempfile.mkstemp(
+        suffix=".txt", prefix="serverinfo-", dir=EXAMPLES_PATH
+    )
+    os.close(fd)
+    timeout = 100
+    license_server = os.environ["ANSYSLMD_LICENSE_FILE"]
+    port = os.environ["PYFLUENT_FLUENT_PORT"]
+
+    subprocess.run(["docker", "run", "--name", "fluent_server", "-d", "--rm",
+                    "-p", f"{port}:{port}",
+                    "-v", f"{EXAMPLES_PATH}:{EXAMPLES_PATH}",
+                    "-e", f"ANSYSLMD_LICENSE_FILE={license_server}",
+                    "-e", f"REMOTING_PORTS={port}/portspan=2",
+                    "-e", "FLUENT_LAUNCHED_FROM_PYFLUENT=1",
+                    "ghcr.io/pyansys/pyfluent",
+                    "-g", f"-sifile={sifile}"] + args)
+
+    sifile_last_mtime = os.stat(sifile).st_mtime
+    while True:
+        if os.stat(sifile).st_mtime > sifile_last_mtime:
+            time.sleep(1)
+            break
+        if timeout == 0:
+            break
+        time.sleep(1)
+        timeout -= 1
+    if os.path.exists(sifile):
+        os.remove(sifile)
+
+
+if __name__ == "__main__":
+    start_fluent_container(sys.argv[1:])
diff --git a/.ci/stop_fluent.py b/.ci/stop_fluent.py
new file mode 100644
index 00000000000..ff6fa5f7435
--- /dev/null
+++ b/.ci/stop_fluent.py
@@ -0,0 +1,9 @@
+import subprocess
+
+
+def stop_fluent_container():
+    subprocess.run(["docker", "stop", "fluent_server"])
+
+
+if __name__ == "__main__":
+    stop_fluent_container()
diff --git a/.github/workflows/ci.yml b/.github/workflows/ci.yml
index 744346517ed..4a29eadd34c 100644
--- a/.github/workflows/ci.yml
+++ b/.github/workflows/ci.yml
@@ -10,6 +10,10 @@ on:
     branches:
       - main
 
+concurrency:
+  group: ${{ github.ref }}
+  cancel-in-progress: true
+
 jobs:
   stylecheck:
     name: Style Check
@@ -119,12 +123,25 @@ jobs:
       - name: Install pyfluent with post module
         run: make install-post
 
+      - name: Login to GitHub Container Registry
+        uses: docker/login-action@v1
+        with:
+          registry: ghcr.io
+          username: ${{ secrets.GH_USERNAME }}
+          password: ${{ secrets.REPO_DOWNLOAD_PAT }}
+
+      - name: Pull Fluent docker image
+        run: make docker-pull
+
       - name: Build Documentation
         run: |
           pip install -r requirements_docs.txt
           make -C doc html
           touch doc/_build/html/.nojekyll
           echo "fluentdocs.pyansys.com" >> doc/_build/html/CNAME
+        env:
+          ANSYSLMD_LICENSE_FILE: ${{ format('1055@{0}', secrets.LICENSE_SERVER) }}
+          PYFLUENT_FLUENT_PORT: 63084
 
       - name: Upload HTML Documentation
         uses: actions/upload-artifact@v2
diff --git a/.github/workflows/nightly-doc-build.yml b/.github/workflows/nightly-doc-build.yml
index 6e5bc42bd9b..4852e23202d 100644
--- a/.github/workflows/nightly-doc-build.yml
+++ b/.github/workflows/nightly-doc-build.yml
@@ -33,12 +33,25 @@ jobs:
       - name: Install pyfluent with post module
         run: make install-post
 
+      - name: Login to GitHub Container Registry
+        uses: docker/login-action@v1
+        with:
+          registry: ghcr.io
+          username: ${{ secrets.GH_USERNAME }}
+          password: ${{ secrets.REPO_DOWNLOAD_PAT }}
+
+      - name: Pull Fluent docker image
+        run: make docker-pull
+
       - name: Build Documentation
         run: |
           pip install -r requirements_docs.txt
           make -C doc html
           touch doc/_build/html/.nojekyll
           echo "dev.fluentdocs.pyansys.com" >> doc/_build/html/CNAME
+        env:
+          ANSYSLMD_LICENSE_FILE: ${{ format('1055@{0}', secrets.LICENSE_SERVER) }}
+          PYFLUENT_FLUENT_PORT: 63084
 
       - name: Deploy
         uses: JamesIves/github-pages-deploy-action@4.1.4
diff --git a/Makefile b/Makefile
index a7f0ebc1cc8..459a37f5c8b 100644
--- a/Makefile
+++ b/Makefile
@@ -18,6 +18,9 @@ install-pyvistaqt-requirements:
 	@sudo apt-get update
 	@sudo apt-get install libegl1 -y
 
+docker-pull:
+	@docker pull ghcr.io/pyansys/pyfluent:latest
+
 test-import:
 	@python -c "import ansys.fluent.core as pyfluent"
 
diff --git a/doc/source/conf.py b/doc/source/conf.py
index 4169b650e4b..b7db6a1aafc 100644
--- a/doc/source/conf.py
+++ b/doc/source/conf.py
@@ -1,5 +1,8 @@
 """Sphinx documentation configuration file."""
 from datetime import datetime
+import os
+import subprocess
+import sys
 
 from pyansys_sphinx_theme import pyansys_logo_black
 from sphinx_gallery.sorting import FileNameSortKey
@@ -93,6 +96,25 @@
 copybutton_prompt_is_regexp = True
 
 
+_THIS_DIR = os.path.dirname(__file__)
+_START_FLUENT_FILE = os.path.normpath(
+    os.path.join(_THIS_DIR, "..", "..", ".ci", "start_fluent.py")
+)
+_STOP_FLUENT_FILE = os.path.normpath(
+    os.path.join(_THIS_DIR, "..", "..", ".ci", "stop_fluent.py")
+)
+
+
+def _start_or_stop_fluent_container(gallery_conf, fname, when):
+    if when == "before":
+        if fname in ["mixing_elbow_settings_api.py",
+                     "mixing_elbow_tui_api.py"]:
+            args = ["3ddp", "-t4", "-meshing"]
+        subprocess.run([sys.executable, _START_FLUENT_FILE] + args)
+    elif when == "after":
+        subprocess.run([sys.executable, _STOP_FLUENT_FILE])
+
+
 # -- Sphinx Gallery Options ---------------------------------------------------
 sphinx_gallery_conf = {
     # convert rst to md for ipynb
@@ -113,6 +135,8 @@
     "doc_module": "ansys-fluent-core",
     "ignore_pattern": "flycheck*",
     "thumbnail_size": (350, 350),
+    'reset_modules_order': 'both',
+    'reset_modules': (_start_or_stop_fluent_container),
 }
 
 
diff --git a/examples/00-fluent/mixing_elbow_settings_api.py b/examples/00-fluent/mixing_elbow_settings_api.py
new file mode 100644
index 00000000000..dc6b9a6a230
--- /dev/null
+++ b/examples/00-fluent/mixing_elbow_settings_api.py
@@ -0,0 +1,529 @@
+"""
+.. _ref_mixing_elbow_settings_api:
+
+Fluid Flow and Heat Transfer in a Mixing Elbow
+----------------------------------------------
+This example illustrates the setup and solution of a three-dimensional
+turbulent fluid flow and heat transfer problem in a mixing elbow. The mixing
+elbow configuration is encountered in piping systems in power plants and
+process industries. It is often important to predict the flow field and
+temperature field in the area of the mixing region in order to properly design
+the junction.
+
+This example demonstrates how to do the following:
+
+- Use the Watertight Geometry guided workflow to:
+    - Import a CAD geometry
+    - Generate a surface mesh
+    - Decribe the geometry
+    - Generate a volume mesh
+- Launch Ansys Fluent.
+- Read an existing mesh file into Ansys Fluent.
+- Use mixed units to define the geometry and fluid properties.
+- Set material properties and boundary conditions for a turbulent
+  forced-convection problem.
+- Create a surface report definition and use it as a convergence criterion.
+- Calculate a solution using the pressure-based solver.
+- Visually examine the flow and temperature fields using the postprocessing
+  tools available in Ansys Fluent.
+
+Problem Description:
+A cold fluid at 20 deg C flows into the pipe through a large inlet, and mixes
+with a warmer fluid at 40 deg C that enters through a smaller inlet located at
+the elbow. The pipe dimensions are in inches and the fluid properties and
+boundary conditions are given in SI units. The Reynolds number for the flow at
+the larger inlet is 50, 800, so a turbulent flow model will be required.
+
+"""
+
+
+###############################################################################
+
+# First, connect with a Fluent server
+
+import ansys.fluent.core as pyfluent
+from ansys.fluent.core import examples
+
+
+import_filename = examples.download_file(
+    "mixing_elbow.pmdb", "pyfluent/mixing_elbow"
+)
+s = pyfluent.launch_fluent(start_instance=False)
+
+###############################################################################
+
+# Import the CAD geometry (mixing_elbow.scdoc). For Length Units, select "in".
+# Execute the Import Geometry task.
+s.workflow.InitializeWorkflow(WorkflowType="Watertight Geometry")
+
+###############################################################################
+
+# Import the CAD geometry. For Length Units, select "in".
+# Execute the Import Geometry task.
+
+s.workflow.TaskObject["Import Geometry"].Arguments = dict(
+    FileName=import_filename, LengthUnit="in"
+)
+
+s.workflow.TaskObject["Import Geometry"].Execute()
+
+###############################################################################
+
+# Add local sizing:
+# In the Add Local Sizing task, you are prompted as to whether or not you would
+# like to add local sizing controls to the faceted geometry. For the purposes
+# of this tutorial, you can keep the default setting. Execute to complete this
+# task and proceed to the next task in the workflow.
+
+s.workflow.TaskObject["Add Local Sizing"].AddChildToTask()
+s.workflow.TaskObject["Add Local Sizing"].Execute()
+
+###############################################################################
+
+# Generate the surface mesh:
+# In the Generate the Surface Mesh task, you can set various properties of the
+# surface mesh for the faceted geometry. Specify 0.3 for Maximum Size. Execute
+# the Surface Mesh to complete this task and proceed to the next task in the
+# workflow.
+
+s.workflow.TaskObject["Generate the Surface Mesh"].Arguments = {
+    "CFDSurfaceMeshControls": {"MaxSize": 0.3}
+}
+s.workflow.TaskObject["Generate the Surface Mesh"].Execute()
+
+###############################################################################
+
+# Describe the geometry:
+# When you select the Describe Geometry task, you are prompted with questions
+# relating to the nature of the imported geometry. Since the geometry defined
+# the fluid region. Select The geometry consists of only fluid regions with no
+# voids for Geometry Type. Execute Describe Geometry to complete this task and
+# proceed
+# to the next task in the workflow.
+s.workflow.TaskObject["Describe Geometry"].UpdateChildTasks(
+    SetupTypeChanged=False
+)
+s.workflow.TaskObject["Describe Geometry"].Arguments = dict(
+    SetupType="The geometry consists of only fluid regions with no voids"
+)
+s.workflow.TaskObject["Describe Geometry"].UpdateChildTasks(
+    SetupTypeChanged=True
+)
+s.workflow.TaskObject["Describe Geometry"].Execute()
+
+###############################################################################
+
+# Update Boundaries Task:
+# For the wall-inlet boundary, change the Boundary Type field to wall. Execute
+# Update Boundaries to complete this task and proceed to the next task in the
+# workflow.
+s.workflow.TaskObject["Update Boundaries"].Arguments = {
+    "BoundaryLabelList": ["wall-inlet"],
+    "BoundaryLabelTypeList": ["wall"],
+    "OldBoundaryLabelList": ["wall-inlet"],
+    "OldBoundaryLabelTypeList": ["velocity-inlet"],
+}
+s.workflow.TaskObject["Update Boundaries"].Execute()
+
+###############################################################################
+
+# Update your regions:
+# Select the Update Regions task, where you can review the names and types of
+# the various regions that have been generated from your imported geometry, and
+# change them as needed. Keep the default settings, and execute Update Regions.
+s.workflow.TaskObject["Update Regions"].Execute()
+
+###############################################################################
+
+# Add Boundary Layers:
+# Select the Add Boundary Layers task, where you can set properties of the
+# boundary layer mesh. Keep the default settings, and Add Boundary Layers.
+s.workflow.TaskObject["Add Boundary Layers"].AddChildToTask()
+s.workflow.TaskObject["Add Boundary Layers"].InsertCompoundChildTask()
+s.workflow.TaskObject["smooth-transition_1"].Arguments = {
+    "BLControlName": "smooth-transition_1",
+}
+s.workflow.TaskObject["Add Boundary Layers"].Arguments = {}
+s.workflow.TaskObject["smooth-transition_1"].Execute()
+
+###############################################################################
+
+# Generate the volume mesh:
+# Select the Generate the Volume Mesh task, where you can set properties of the
+# volume mesh. Select the poly-hexcore for Fill With. Execute Generate the
+# Volume Mesh.
+s.workflow.TaskObject["Generate the Volume Mesh"].Arguments = {
+    "VolumeFill": "poly-hexcore",
+    "VolumeFillControls": {
+        "HexMaxCellLength": 0.3,
+    },
+}
+s.workflow.TaskObject["Generate the Volume Mesh"].Execute()
+
+###############################################################################
+
+# Check the mesh in Meshing mode
+s.tui.meshing.mesh.check_mesh()
+
+###############################################################################
+
+# Save the mesh file (mixing_elbow.msh.h5).
+# s.tui.meshing.file.write_mesh('mixing_elbow.msh.h5')
+
+###############################################################################
+
+# Switch to Solution mode:
+# Now that a high-quality mesh has been generated using Ansys Fluent in meshing
+# mode, you can now switch to solver mode to complete the setup of the
+# simulation. We have just checked the mesh, so select Yes to switch to
+# solution mode.
+s.tui.meshing.switch_to_solution_mode("yes")
+
+###############################################################################
+
+# Check the mesh in Solver mode:
+# The mesh check will list the minimum and maximum x, y, and z values from the
+# mesh in the default SI unit of meters. It will also report a number of other
+# mesh features that are checked. Any errors in the mesh will be reported at
+# this time. Ensure that the minimum volume is not negative, since Ansys Fluent
+# cannot begin a calculation when this is the case.
+s.tui.solver.mesh.check()
+
+###############################################################################
+
+# The settings objects provide a natural way to access and modify settings.
+# The top-level settings object for a session can be accessed with the
+# get_settings_root() method of the session object.
+# Enabling the settings objects.
+root = s.get_settings_root()
+
+###############################################################################
+
+# Set the working units for the mesh:
+# select "in" to set inches as the working unit for length. Note:  Because the
+# default SI units will be used for everything except length, there is no need
+# to change any other units in this problem. If you want a different working
+# unit for length, other than inches (for example, millimeters), make the
+# appropriate change.
+s.tui.solver.define.units("length", "in")
+
+###############################################################################
+
+# Enable heat transfer by activating the energy equation.
+root.setup.models.energy.enabled = True
+
+###############################################################################
+
+# Create a new material called water-liquid.
+s.tui.solver.define.materials.copy("fluid", "water-liquid")
+
+###############################################################################
+
+# Set up the cell zone conditions for the fluid zone (elbow-fluid). Select
+# water-liquid from the Material list.
+s.tui.solver.define.boundary_conditions.fluid(
+    "elbow-fluid",
+    "yes",
+    "water-liquid",
+    "no",
+    "no",
+    "no",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "1",
+    "no",
+    "no",
+    "no",
+    "no",
+    "no",
+)
+
+###############################################################################
+
+# Set up the boundary conditions for the inlets, outlet, and walls for your CFD
+# analysis.
+
+# cold inlet (cold-inlet), Setting: Value:
+# Velocity Specification Method: Magnitude, Normal to Boundary
+
+# Velocity Magnitude: 0.4 [m/s]
+# Specification Method: Intensity and Hydraulic Diameter
+# Turbulent Intensity: 5 [%]
+# Hydraulic Diameter: 4 [inch]
+# Temperature: 293.15 [K]
+root.setup.boundary_conditions.velocity_inlet["cold-inlet"].vmag = {
+    "option": "constant or expression",
+    "constant": 0.4,
+}
+root.setup.boundary_conditions.velocity_inlet[
+    "cold-inlet"
+].ke_spec = "Intensity and Hydraulic Diameter"
+root.setup.boundary_conditions.velocity_inlet["cold-inlet"].turb_intensity = 5
+root.setup.boundary_conditions.velocity_inlet[
+    "cold-inlet"
+].turb_hydraulic_diam = "4 [in]"
+root.setup.boundary_conditions.velocity_inlet["cold-inlet"].t = {
+    "option": "constant or expression",
+    "constant": 293.15,
+}
+
+###############################################################################
+
+# hot inlet (hot-inlet), Setting: Value:
+# Velocity Specification Method: Magnitude, Normal to Boundary
+# Velocity Magnitude: 1.2 [m/s]
+# Specification Method: Intensity and Hydraulic Diameter
+# Turbulent Intensity: 5 [%]
+# Hydraulic Diameter: 1 [inch]
+# Temperature: 313.15 [K]
+
+root.setup.boundary_conditions.velocity_inlet["hot-inlet"].vmag = {
+    "option": "constant or expression",
+    "constant": 1.2,
+}
+root.setup.boundary_conditions.velocity_inlet[
+    "hot-inlet"
+].ke_spec = "Intensity and Hydraulic Diameter"
+root.setup.boundary_conditions.velocity_inlet[
+    "hot-inlet"
+].turb_hydraulic_diam = "1 [in]"
+root.setup.boundary_conditions.velocity_inlet["hot-inlet"].t = {
+    "option": "constant or expression",
+    "constant": 313.15,
+}
+
+###############################################################################
+
+# pressure outlet (outlet), Setting: Value:
+# Backflow Turbulent Intensity: 5 [%]
+# Backflow Turbulent Viscosity Ratio: 4
+root.setup.boundary_conditions.pressure_outlet[
+    "outlet"
+].turb_viscosity_ratio = 4
+
+###############################################################################
+
+# Enable the plotting of residuals during the calculation.
+s.tui.solver.solve.monitors.residual.plot("yes")
+
+###############################################################################
+
+# Create a surface report definition of average temperature at the outlet
+# (outlet) called outlet-temp-avg
+root.solution.report_definitions.surface["outlet-temp-avg"] = {}
+root.solution.report_definitions.surface[
+    "outlet-temp-avg"
+].report_type = "surface-massavg"
+root.solution.report_definitions.surface[
+    "outlet-temp-avg"
+].field = "temperature"
+root.solution.report_definitions.surface["outlet-temp-avg"].surface_names = [
+    "outlet"
+]
+root.solution.report_definitions.compute(report_defs=["outlet-temp-avg"])
+
+###############################################################################
+
+# Create a surface report file called outlet-temp-avg-rfile using
+# report-definition outlet-temp-avg
+# s.tui.solver.solve.report_files.add(
+#    "outlet-temp-avg-rfile",
+#    "report-defs",
+#    "outlet-temp-avg",
+#    "()",
+#    "file-name",
+#    "outlet-temp-avg-rfile.out",
+#    "print?",
+#    "yes",
+#    "file-name",
+#    "outlet-temp-avg-rfile.out",
+#    "frequency",
+#    "3",
+#    "frequency-of",
+#    "iteration",
+#    "itr-index",
+#    "1",
+#    "run-index",
+#    "0",
+#    "quit",
+# )
+
+###############################################################################
+
+# Create a convergence condition for outlet-temp-avg:
+# Provide con-outlet-temp-avg for Conditions. Select outlet-temp-avg Report
+# Definition. Provide 1e-5 for Stop Criterion. Provide 20 for Ignore Iterations
+# Before. Provide 15 for Use Iterations. Enable Print. Set Every Iteration to
+# 3.
+
+
+# These settings will cause Fluent to consider the solution converged when the
+# surface report definition value for each of the previous 15 iterations is
+# within 0.001% of the current value. Convergence of the values will be checked
+# every 3 iterations. The first 20 iterations will be ignored, allowing for any
+# initial solution dynamics to settle out. Note that the value printed to the
+# console is the deviation between the current and previous iteration values
+# only.
+# Change Convergence Conditions
+s.tui.solver.solve.convergence_conditions(
+    "conv-reports",
+    "add",
+    "con-outlet-temp-avg",
+    "initial-values-to-ignore",
+    "20",
+    "previous-values-to-consider",
+    "15",
+    "print?",
+    "yes",
+    "report-defs",
+    "outlet-temp-avg",
+    "stop-criterion",
+    "1e-05",
+    "quit",
+    "quit",
+    "condition",
+    "1",
+    "frequency",
+    "3",
+    "quit",
+)
+s.tui.solver.solve.convergence_conditions("frequency", "3", "quit").result()
+
+###############################################################################
+
+# Initialize the flow field using the Hybrid Initialization
+s.tui.solver.solve.initialize.hyb_initialization().result()
+
+###############################################################################
+
+# Solve for 150 Iterations.
+s.tui.solver.solve.iterate(150).result()
+
+###############################################################################
+
+# Save the case and data file (mixing_elbow1.cas.h5 and mixing_elbow1.dat.h5).
+# s.tui.solver.file.write_case_data('mixing_elbow1.cas.h5')
+
+###############################################################################
+
+# Examine the mass flux report for convergence: Select cold-inlet, hot-inlet,
+# and outlet from the Boundaries selection list.
+# Compute a Mass Flux Report for convergence
+root.solution.report_definitions.flux[
+    "report_mfr"
+] = {}  # Create a default report flux report
+root.solution.report_definitions.flux["report_mfr"].zone_names = [
+    "cold-inlet",
+    "hot-inlet",
+    "outlet",
+]
+root.solution.report_definitions.compute(report_defs=["report_mfr"])
+
+###############################################################################
+
+# Create and display a definition for velocity magnitude contours on the
+# symmetry plane:
+# Provide contour-vel for Contour Name. Select velocity magnitude. Select
+# symmetry-xyplane from the Surfaces list. Display contour-vel contour.
+
+root.results.graphics.contour["contour-vel"] = {}
+root.results.graphics.contour["contour-vel"].print_state()
+root.results.graphics.contour["contour-vel"].field = "velocity-magnitude"
+root.results.graphics.contour["contour-vel"].surfaces_list = [
+    "symmetry-xyplane"
+]
+# root.results.graphics.contour["contour-vel"].display()
+
+###############################################################################
+
+# Create and display a definition for temperature contours on the symmetry
+# plane:
+# Provide contour-temp for Contour Name. Select temperature. Select
+# symmetry-xyplane from the Surfaces list. Display contour-temp contour.
+
+root.results.graphics.contour["contour-temp"] = {}
+root.results.graphics.contour["contour-temp"].print_state()
+root.results.graphics.contour["contour-temp"].field = "temperature"
+root.results.graphics.contour["contour-temp"].surfaces_list = [
+    "symmetry-xyplane"
+]
+# root.results.graphics.contour["contour-temp"].display()
+
+###############################################################################
+
+# Create and display velocity vectors on the symmetry-xyplane plane:
+
+# Provide vector-vel for Vector Name. Select arrow for the Style. Select
+# symmetry-xyplane from the Surfaces selection list.
+root.results.graphics.vector["vector-vel"] = {}
+root.results.graphics.vector["vector-vel"].print_state()
+root.results.graphics.vector["vector-vel"].field = "temperature"
+root.results.graphics.vector["vector-vel"].surfaces_list = ["symmetry-xyplane"]
+root.results.graphics.vector["vector-vel"].scale.scale_f = 4
+root.results.graphics.vector["vector-vel"].style = "arrow"
+# root.results.graphics.vector["vector-vel"].display()
+
+###############################################################################
+
+# Create an iso-surface representing the intersection of the plane z=0 and the
+# surface outlet. Name: z=0_outlet
+s.tui.solver.surface.iso_surface(
+    "z-coordinate", "z=0_outlet", "outlet", "()", "()", "0", "()"
+).result()
+
+# Create Contour on the iso-surface
+root.results.graphics.contour["contour-z_0_outlet"] = {}
+root.results.graphics.contour["contour-z_0_outlet"].print_state()
+root.results.graphics.contour["contour-z_0_outlet"].field = "temperature"
+root.results.graphics.contour["contour-z_0_outlet"].surfaces_list = [
+    "z=0_outlet"
+]
+# root.results.graphics.contour["contour-z_0_outlet"].display()
+
+###############################################################################
+# s.tui.solver.file.write_case_data("mixing_elbow1_set.cas.h5").result()
+
+# Display and save an XY plot of the temperature profile across the centerline
+# of the outlet for the initial solution
+
+s.tui.solver.display.objects.create(
+    "xy",
+    "xy-outlet-temp",
+    "y-axis-function",
+    "temperature",
+    "surfaces-list",
+    "z=0_outlet",
+    "()",
+    "quit",
+).result()
+# s.tui.solver.display.objects.display("xy-outlet-temp").result()
+# s.tui.solver.plot.plot(
+#    "yes",
+#    "temp-1.xy",
+#    "no",
+#    "no",
+#    "no",
+#    "temperature",
+#    "yes",
+#    "1",
+#    "0",
+#    "0",
+#    "z=0_outlet",
+#    "()",
+#).result()
+
+###############################################################################
+
+# Write final case and data.
+# s.tui.solver.file.write_case_data('mixing_elbow2_set.cas.h5').result()
+
+###############################################################################
diff --git a/examples/00-fluent/mixing_elbow_tui_api.py b/examples/00-fluent/mixing_elbow_tui_api.py
new file mode 100644
index 00000000000..352699dbe73
--- /dev/null
+++ b/examples/00-fluent/mixing_elbow_tui_api.py
@@ -0,0 +1,534 @@
+"""
+.. _ref_mixing_elbow_tui_api:
+
+Fluid Flow and Heat Transfer in a Mixing Elbow
+----------------------------------------------
+This example illustrates the setup and solution of a three-dimensional
+turbulent fluid flow and heat transfer problem in a mixing elbow. The mixing
+elbow configuration is encountered in piping systems in power plants and
+processindustries. It is often important to predict the flow field and
+temperature field in the area of the mixing regionin order to properly design
+the junction.
+
+This example demonstrates how to do the following:
+
+- Use the Watertight Geometry guided workflow to:
+    - Import a CAD geometry
+    - Generate a surface mesh
+    - Decribe the geometry
+    - Generate a volume mesh
+- Launch Ansys Fluent.
+- Read an existing mesh file into Ansys Fluent.
+- Use mixed units to define the geometry and fluid properties.
+- Set material properties and boundary conditions for a turbulent
+  forced-convection problem.
+- Create a surface report definition and use it as a convergence criterion.
+- Calculate a solution using the pressure-based solver.
+- Visually examine the flow and temperature fields using the postprocessing
+  tools available in Ansys Fluent.
+
+Problem Description:
+A cold fluid at 20 deg C flows into the pipe through a large inlet, and mixes
+with a warmer fluid at 40 deg C that enters through a smaller inlet located at
+the elbow. The pipe dimensions are in inches and the fluid properties and
+boundary conditions are given in SI units. The Reynolds number for the flow at
+the larger inlet is 50, 800, so a turbulent flow model will be required.
+
+"""
+
+###############################################################################
+
+# First, connect with a Fluent server
+
+import ansys.fluent.core as pyfluent
+from ansys.fluent.core import examples
+
+
+import_filename = examples.download_file(
+    "mixing_elbow.pmdb", "pyfluent/mixing_elbow"
+)
+s = pyfluent.launch_fluent(start_instance=False)
+
+###############################################################################
+
+# Select the Watertight Geometry Meshing Workflow
+s.workflow.InitializeWorkflow(WorkflowType="Watertight Geometry")
+
+
+###############################################################################
+
+# Import the CAD geometry. For Length Units, select "in".
+# Execute the Import Geometry task.
+
+s.workflow.TaskObject["Import Geometry"].Arguments = dict(
+    FileName=import_filename, LengthUnit="in"
+)
+
+s.workflow.TaskObject["Import Geometry"].Execute()
+
+###############################################################################
+
+# Add local sizing:
+# In the Add Local Sizing task, you are prompted as to whether or not you would
+# like to add local sizing controls to the faceted geometry. For the purposes
+# of this example, you can keep the default setting. Execute to complete this
+# task
+# and proceed to the next task in the workflow.
+s.workflow.TaskObject["Add Local Sizing"].AddChildToTask()
+s.workflow.TaskObject["Add Local Sizing"].Execute()
+
+###############################################################################
+
+# Generate the surface mesh:
+# In the Generate the Surface Mesh task, you can set various properties of the
+# surface mesh for the faceted geometry. Specify 0.3 for Maximum Size. Execute
+# the Surface Mesh to complete this task and proceed to the next task in the
+# workflow.
+s.workflow.TaskObject["Generate the Surface Mesh"].Arguments = {
+    "CFDSurfaceMeshControls": {"MaxSize": 0.3}
+}
+s.workflow.TaskObject["Generate the Surface Mesh"].Execute()
+
+###############################################################################
+
+# Describe the geometry:
+# When you select the Describe Geometry task, you are prompted with questions
+# relating to the nature of the imported geometry. Since the geometry defined
+# the fluid region. Select The geometry consists of only fluid regions with no
+# voids for Geometry Type. Execute Describe Geometry to complete this task and
+# proceed
+# to the next task in the workflow.
+s.workflow.TaskObject["Describe Geometry"].UpdateChildTasks(
+    SetupTypeChanged=False
+)
+s.workflow.TaskObject["Describe Geometry"].Arguments = dict(
+    SetupType="The geometry consists of only fluid regions with no voids"
+)
+s.workflow.TaskObject["Describe Geometry"].UpdateChildTasks(
+    SetupTypeChanged=True
+)
+s.workflow.TaskObject["Describe Geometry"].Execute()
+
+###############################################################################
+
+# Update Boundaries Task:
+# For the wall-inlet boundary, change the Boundary Type field to wall. Execute
+# Update Boundaries to complete this task and proceed to the next task in the
+# workflow.
+s.workflow.TaskObject["Update Boundaries"].Arguments = {
+    "BoundaryLabelList": ["wall-inlet"],
+    "BoundaryLabelTypeList": ["wall"],
+    "OldBoundaryLabelList": ["wall-inlet"],
+    "OldBoundaryLabelTypeList": ["velocity-inlet"],
+}
+s.workflow.TaskObject["Update Boundaries"].Execute()
+
+###############################################################################
+
+# Update your regions:
+# Select the Update Regions task, where you can review the names and types of
+# the various regions that have been generated from your imported geometry, and
+# change them as needed. Keep the default settings, and execute Update Regions.
+s.workflow.TaskObject["Update Regions"].Execute()
+
+###############################################################################
+
+# Add Boundary Layers:
+# Select the Add Boundary Layers task, where you can set properties of the
+# boundary layer mesh. Keep the default settings, and Add Boundary Layers.
+
+s.workflow.TaskObject["Add Boundary Layers"].AddChildToTask()
+s.workflow.TaskObject["Add Boundary Layers"].InsertCompoundChildTask()
+s.workflow.TaskObject["smooth-transition_1"].Arguments = {
+    "BLControlName": "smooth-transition_1",
+}
+s.workflow.TaskObject["Add Boundary Layers"].Arguments = {}
+s.workflow.TaskObject["smooth-transition_1"].Execute()
+
+###############################################################################
+
+# Generate the volume mesh:
+# Select the Generate the Volume Mesh task, where you can set properties of the
+# volume mesh. Select the poly-hexcore for Fill With. Execute Generate the
+# Volume Mesh.
+s.workflow.TaskObject["Generate the Volume Mesh"].Arguments = {
+    "VolumeFill": "poly-hexcore",
+    "VolumeFillControls": {
+        "HexMaxCellLength": 0.3,
+    },
+}
+s.workflow.TaskObject["Generate the Volume Mesh"].Execute()
+
+###############################################################################
+
+# Check the mesh in Meshing mode
+s.tui.meshing.mesh.check_mesh().result()
+
+###############################################################################
+
+# Save the mesh file (mixing_elbow.msh.h5)
+# s.tui.meshing.file.write_mesh('mixing_elbow.msh.h5')
+
+###############################################################################
+
+# Switch to Solution mode:
+# Now that a high-quality mesh has been generated using Ansys Fluent in meshing
+# mode, you can now switch to solver mode to complete the setup of the
+# simulation. We have just checked the mesh, so select Yes to switch to
+# solution mode.
+s.tui.meshing.switch_to_solution_mode("yes").result()
+
+###############################################################################
+
+# Check the mesh in Solver mode:
+# The mesh check will list the minimum and maximum x, y, and z values from the
+# mesh in the default SI unit of meters. It will also report a number of other
+# mesh features that are checked. Any errors in the mesh will be reported at
+# this time. Ensure that the minimum volume is not negative, since Ansys Fluent
+# cannot begin a calculation when this is the case.
+s.tui.solver.mesh.check().result()
+
+###############################################################################
+
+# Set the working units for the mesh:
+# select "in" to set inches as the working unit for length. Note:  Because the
+# default SI units will be used for everything except length, there is no need
+# to change any other units in this problem. If you want a different working
+# unit for length, other than inches (for example, millimeters), make the
+# appropriate change.
+s.tui.solver.define.units("length", "in").result()
+
+###############################################################################
+
+# Enable heat transfer by activating the energy equation.
+s.tui.solver.define.models.energy("yes", ", ", ", ", ", ", ", ").result()
+
+###############################################################################
+
+# Create a new material called water-liquid.
+s.tui.solver.define.materials.copy("fluid", "water-liquid").result()
+
+###############################################################################
+
+# Set up the cell zone conditions for the fluid zone (elbow-fluid). Select
+# water-liquid from the Material list.
+s.tui.solver.define.boundary_conditions.fluid(
+    "elbow-fluid",
+    "yes",
+    "water-liquid",
+    "no",
+    "no",
+    "no",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "0",
+    "no",
+    "1",
+    "no",
+    "no",
+    "no",
+    "no",
+    "no",
+).result()
+
+###############################################################################
+
+# Set up the boundary conditions for the inlets, outlet, and walls for your CFD
+# analysis.
+
+# cold inlet (cold-inlet), Setting: Value:
+# Velocity Specification Method: Magnitude, Normal to Boundary
+
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "cold-inlet", [], "vmag", "no", 0.4, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "cold-inlet", [], "ke-spec", "no", "no", "no", "yes", "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "cold-inlet", [], "turb-intensity", 5, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "cold-inlet", [], "turb-hydraulic-diam", 4, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "cold-inlet", [], "temperature", "no", 293.15, "quit"
+).result()
+
+###############################################################################
+
+# hot inlet (hot-inlet), Setting: Value:
+# Velocity Specification Method: Magnitude, Normal to Boundary
+
+
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "hot-inlet", [], "vmag", "no", 1.2, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "hot-inlet", [], "ke-spec", "no", "no", "no", "yes", "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "hot-inlet", [], "turb-intensity", 5, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "hot-inlet", [], "turb-hydraulic-diam", 1, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.velocity_inlet(
+    "hot-inlet", [], "temperature", "no", 313.15, "quit"
+).result()
+
+###############################################################################
+
+# pressure outlet (outlet), Setting: Value:
+# Backflow Turbulent Intensity: 5 [%]
+# Backflow Turbulent Viscosity Ratio: 4
+
+s.tui.solver.define.boundary_conditions.set.pressure_outlet(
+    "outlet", [], "turb-intensity", 5, "quit"
+).result()
+s.tui.solver.define.boundary_conditions.set.pressure_outlet(
+    "outlet", [], "turb-viscosity-ratio", 4, "quit"
+).result()
+
+###############################################################################
+
+# Enable the plotting of residuals during the calculation.
+s.tui.solver.solve.monitors.residual.plot("yes")
+
+###############################################################################
+
+# Create a surface report definition of average temperature at the outlet
+# (outlet) called "outlet-temp-avg
+s.tui.solver.solve.report_definitions.add(
+    "outlet-temp-avg",
+    "surface-massavg",
+    "field",
+    "temperature",
+    "surface-names",
+    "outlet",
+    "()",
+    "quit",
+).result()
+
+###############################################################################
+
+# Create a surface report file called outlet-temp-avg-rfile using
+# report-definition outlet-temp-avg
+# s.tui.solver.solve.report_files.add(
+#    "outlet-temp-avg-rfile",
+#    "report-defs",
+#    "outlet-temp-avg",
+#    "()",
+#    "file-name",
+#    "outlet-temp-avg-rfile.out",
+#    "print?",
+#    "yes",
+#    "file-name",
+#    "outlet-temp-avg-rfile.out",
+#    "frequency",
+#    "3",
+#    "frequency-of",
+#    "iteration",
+#    "itr-index",
+#    "1",
+#    "run-index",
+#    "0",
+#    "quit",
+# ).result()
+
+###############################################################################
+
+# Create a convergence condition for outlet-temp-avg:
+# Provide con-outlet-temp-avg for Conditions. Select outlet-temp-avg Report
+# Definition. Provide 1e-5 for Stop Criterion. Provide 20 for Ignore Iterations
+# Before. Provide 15 for Use Iterations. Enable Print. Set Every Iteration to
+# 3.
+
+
+# These settings will cause Fluent to consider the solution converged when the
+# surface report definition value for each of the previous 15 iterations is
+# within 0.001% of the current value. Convergence of the values will be checked
+# every 3 iterations. The first 20 iterations will be ignored, allowing for any
+# initial solution dynamics to settle out. Note that the value printed to the
+# console is the deviation between the current and previous iteration values
+# only.
+s.tui.solver.solve.convergence_conditions(
+    "conv-reports",
+    "add",
+    "con-outlet-temp-avg",
+    "initial-values-to-ignore",
+    "20",
+    "previous-values-to-consider",
+    "15",
+    "print?",
+    "yes",
+    "report-defs",
+    "outlet-temp-avg",
+    "stop-criterion",
+    "1e-05",
+    "quit",
+    "quit",
+    "condition",
+    "1",
+    "frequency",
+    "3",
+    "quit",
+).result()
+s.tui.solver.solve.convergence_conditions("frequency", "3", "quit").result()
+
+###############################################################################
+
+# Initialize the flow field using the Hybrid Initialization
+s.tui.solver.solve.initialize.hyb_initialization().result()
+
+###############################################################################
+
+# Save the case file (mixing_elbow1.cas.h5).
+# s.tui.solver.file.write_case('mixing_elbow1.cas.h5').result()
+
+###############################################################################
+
+# Solve for 150 Iterations.
+s.tui.solver.solve.iterate(150).result()
+
+###############################################################################
+
+# Examine the mass flux report for convergence: Select cold-inlet, hot-inlet,
+# and outlet from the Boundaries selection list.
+# s.tui.solver.report.fluxes.mass_flow(
+#    "no", "cold-inlet", "hot-inlet", "outlet", "()", "yes", "mass-flux1.flp"
+# )
+
+###############################################################################
+
+# Save the data file (mixing_elbow1.dat.h5).
+# s.tui.solver.file.write_data('mixing_elbow1.dat.h5').result()
+
+###############################################################################
+
+# Create and display a definition for velocity magnitude contours on the
+# symmetry plane:
+# Provide contour-vel for Contour Name. Select velocity magnitude. Select
+# symmetry-xyplane from the Surfaces list. Display contour-vel contour.
+
+
+s.tui.solver.display.objects.create(
+    "contour",
+    "contour-vel",
+    "filled?",
+    "yes",
+    "node-values?",
+    "yes",
+    "field",
+    "velocity-magnitude",
+    "surfaces-list",
+    "symmetry-xyplane",
+    "()",
+    "coloring",
+    "banded",
+    "quit",
+).result()
+# s.tui.solver.display.objects.display("contour-vel").result()
+
+###############################################################################
+
+# Create and display a definition for temperature contours on the symmetry
+# plane:
+# Provide contour-temp for Contour Name. Select temperature. Select
+# symmetry-xyplane from the Surfaces list. Display contour-temp contour.
+
+s.tui.solver.display.objects.create(
+    "contour",
+    "contour-temp",
+    "filled?",
+    "yes",
+    "node-values?",
+    "yes",
+    "field",
+    "temperature",
+    "surfaces-list",
+    "symmetry-xyplane",
+    "()",
+    "coloring",
+    "smooth",
+    "quit",
+)
+# s.tui.solver.display.objects.display("contour-temp").result()
+
+###############################################################################
+
+# Create and display velocity vectors on the symmetry-xyplane plane:
+
+# Provide vector-vel for Vector Name. Select arrow for the Style. Select
+# symmetry-xyplane from the Surfaces selection list. Provide 4 for Scale. Set
+# Skip to 2.
+s.tui.solver.display.objects.create(
+    "vector",
+    "vector-vel",
+    "style",
+    "arrow",
+    "surface-list",
+    "symmetry-xyplane",
+    "()",
+    "scale",
+    "scale-f",
+    "4",
+    "quit",
+    "skip",
+    "2",
+    "quit",
+).result()
+# s.tui.solver.display.objects.display("vector-vel").result()
+
+###############################################################################
+
+# Create an iso-surface representing the intersection of the plane z=0 and the
+# surface outlet. Name: z=0_outlet
+s.tui.solver.surface.iso_surface(
+    "z-coordinate", "z=0_outlet", "outlet", "()", "()", "0", "()"
+).result()
+
+###############################################################################
+# s.tui.solver.file.write_case_data("mixing_elbow1_tui.cas.h5").result()
+
+# Display and save an XY plot of the temperature profile across the centerline
+# of the outlet for the initial solution
+s.tui.solver.display.objects.create(
+    "xy",
+    "xy-outlet-temp",
+    "y-axis-function",
+    "temperature",
+    "surfaces-list",
+    "z=0_outlet",
+    "()",
+    "quit",
+).result()
+# s.tui.solver.display.objects.display("xy-outlet-temp").result()
+# s.tui.solver.plot.plot(
+#    "yes",
+#    "temp-1.xy",
+#    "no",
+#    "no",
+#    "no",
+#    "temperature",
+#    "yes",
+#    "1",
+#    "0",
+#    "0",
+#    "z=0_outlet",
+#    "()",
+# ).result()
+
+###############################################################################
+
+# Write final case and data.
+# s.tui.solver.file.write_case_data("mixing_elbow2_tui.cas.h5").result()
+
+###############################################################################