Author: Benjamin Mullen, Cornell University

Problem Specification

1. Pre-Analysis & Start-Up

2. Geometry

3. Mesh

4. Physics Setup

5. Numerical Solution

6. Numerical Results

7. Verification & Validation

Exercises

Comments

# Numerical Results

### Open the Post Processor

In the *Project Schematic* double click ** Results** to open the post processor. When the

*A6: Fluid Flow (FLUENT) - CFD - Post*Window opens, look at the geometry by clicking the +Z axis on the compass

### Velocity Vectors

In the Post Processing window, click the Vector icon to create a vector result. When prompted, name the result `Velocity Vector`

. In the *Details of Velocity Vector* window, begin on the ** Geometry** tab. Under

**, select**

*Locations***. This will show the velocity along the entire geometry surface periodically. Next, click on the**

*Periodic 1***tab. Change the**

*Symbol***to**

*Symbol Size*`0.1`

. Finally, move to the **tab. We want to see the entire geometry of the pipe: not just half of it like we currently see. To see the whole pipe, check the box next to**

*View***, and change the**

*Apply Reflection/Mirroring***to**

*Method***. Because the pipe is long and skinny, it will be difficult to see the results. This post processor allows us to stretch the results to make the results easier to see. To apply a scaling, check the box next to**

*ZX Plane***, and change the**

*Apply Scale***to**

*Scale*`1,10,1`

(this will scale the y-direction by 10). When finished, press **to see the result. If you wish to see the result without the wireframe of the pipe, uncheck the box next to**

*Apply***under**

*Wireframe**User Location and Plots*.

In ANSYS version 14.5 and later, only one half of the pipe cross-section is displayed after using the mirroring option. You can work around this by applying the mirroring condition in the "Default transform" setting and not in the "View" Tab. To do this select "Default Transform" in the left-hand menu, uncheck "Instancing Info from Domain", check "Apply Reflection" and select to mirror about the ZX Plane.

### Velocity Contour

In the Post Processing window, click the Contour icon to create a Contour result. When prompted, name the result `Velocity Contour`

. In the *Details of Velocity Contour* window, begin on the ** Geometry** tab. Under

**, again select**

*Locations***. Also, change the**

*Periodic 1***to**

*Variable***. Next, move to the**

*Velocity***tab. Check the box next to**

*View***, and change the**

*Apply Reflection/Mirroring***to**

*Method***and again,check the box next to**

*ZX Plane***, and change the**

*Apply Scale***to**

*Scale*`1,10,1`

. When finished, press **to see the result. Finally, we need to remove the Velocity Vectors from the Graphic Window. Do this by unchecking the box next to**

*Apply***in the**

*Velocity Vector**Outline*window under

*User Location and Plots*.

### Temperature Contour

In the Post Processing window, click the Contour icon to create another Contour result. When prompted, name the result `Temperature Contour`

. In the *Details of Temperature Contour* window, begin on the ** Geometry** tab. Under

**, select**

*Locations***. This time, change the**

*Periodic 1***to**

*Variable***. Next, move to the**

*Temperature***tab. Apply the same mirroring and scaling as we did for the Velocity Contours. When finished, press**

*View***. Uncheck the box next to**

*Apply***to only see the Temperature Contours.**

*Velocity Contour*### Pressure Contour

Create another contour result, and name `Pressure Contour`

. Use all of the same settings as the previous results but this time choosing ** Variable > Pressure** in the

**tab.**

*Geometry*### Graph of Temperature along Centerline

To graph the temperature along the centerline, we first need to create the centerline as a path. To accomplish this, click on the Location icon , select ** Line**, and name the line

`Centerline`

. In the *Details of Centerline*window, set the

**to two points. Point 1 is (0,0,0), and Point 2 is (8.64,0,0). Enter these values into the**

*Method**details*window. Next, change the number of

**to 100. Press**

*Samples***once finished.**

*Apply*To create a chart, press the chart icon . When prompted, name the page `Temperature Along Centerline`

. In the *Details of Temperature Along Centerline* window, begin on the ** General** tab. In the

**, enter**

*Title*`Temperature Along Centerline`

. Next, click on the **tab. Under**

*Data Series**Data Source*, in the drop down menu next to

**, select**

*Location***. Now move to the**

*Centerline***tab. In the drop down menu next to**

*X Axis***, scroll all the way down and select**

*Variable***. In the**

*X***tab, change the**

*Y Axis***to**

*Variable***. when finished, press**

*Temperature***to see the chart.**

*Apply*### Graph of Temperature along Outlet

To graph the temperature along the outlet, we need to create the outlet as a path much like we did with the centerline. Click on the Location icon , select ** Line**, and name the line

`Outlet`

. In the *Details of Centerline*window, set the

**to two points. Point 1 is (8.64,0,0), and Point 2 is (8.64,0.06,0). Enter these values into the**

*Method**details*window. Next, change the number of

**to 100. Press**

*Samples***once finished.**

*Apply*Next, press the chart icon . When prompted, name the page `Temperature Along Outlet`

. In the *Details of Temperature Along Outlet* window, begin on the ** General** tab. In the

**, enter**

*Title*`Temperature Along Outlet`

. Next, click on the **tab. Under**

*Data Series**Data Source*, in the drop down menu next to

**, select**

*Location***. Now move to the**

*Outlet***tab. In the drop down menu next to**

*X Axis***, and select**

*Variable***. In the**

*Temperature***tab, change the**

*Y Axis***to**

*Variable***. when finished, press**

*Y***to see the chart.**

*Apply*### Graph of Nusselt Number along the heated section of the pipe

The ** Nusselt number** is a non-dimensional parameter that provides a measure of the convection heat transfer at a surface. It is the ratio of convection to pure conduction heat transfer. We will now derive the Nusselt number as a function of the given parameters and temperature.

The convection heat transfer at the pipe wall is:

We can rearrange terms to find an expression for h, the

**:**

*convection coefficient*Substitute the

**expression into the**

*convection coefficient***expression:**

*Nusselt Number*where

h is the convection coefficient.

k is the thermal conductivity.

L is the length scale. Similar to the Reynold's Number, the length scale is the diameter of the pipe for an internal pipe flow.

q''_w is the heat flux at the heated surface, 37.5 W/m^2.

Tw is the pipe wall temperature at a given location along the pipe.

Tm is the mean temperature in the pipe at the location where Tw is defined.

*Wall Temperature*

To find the temperature at the wall, click on ** insert** >>

**>>**

*location***, and name it**

*point**Tw exit*. In the Details of

*Tw exit*window, set Method to XYZ and enter (8.64, 0.06, 0) in Point.Click Apply to create a point at the upper right corner of the pipe.

Click on Expression right below and right click in the window to create a new expression named *Tw*. Under ** Details of Tw** panel, enter

**maxVal(Temperature)@Tw exit**in the

**tab.**

*Definition**Tw* now gives the temperature at the location (8.64, 0.06, 0), which is on the exit pipe wall.

*Mixed Mean Temperature*

To find the mean temperature at a given location in the pipe, click on ** insert** >>

**>>**

*location***, and name it exit. In the Details of exit window, set Method to Two Points and enter (8.64, 0, 0) for Point 1 and (8.64, 0.06, 0) for Point 2. Click Apply to create a line at the exit of the pipe. The mean temperature is the area weighted average temperature and we can use integral to find the appropriate mean Temperature:**

*line*Click on the ** Calculators** tab and double click on

**. Select**

*Function Calculator***for the Function,**

*lengthInt***for the Location, and**

*exit***for the Variable. Check**

*Velocity u***and click Calculate. The expression "lengthInt(Velocity u)@exit is essentially the integral of u*dr and can be conveniently used to calculate the mean temperature.**

*show equivalent expression*Under

**, right click in the window to create a new expression and name it Tm. In the Details of Tm window, enter the following:**

*Expressions*lengthInt(Velocity u*Y*Temperature)@exit/lengthInt(Velocity u*Y)@exit

This expression will now give the mean temperature at the location in which we called "exit". Recall the pipe radius r is defined in the Y direction in FLUENT. Hence we will use Y to define the radial position in the pipe, as shown in the expression above.

*Nusselt Number*

We are now ready to find the Nusselt Number. Create another expression and name it ** Nu exp**. Under the Definition tab, enter the Nusselt Number expression shown in the equation above. The units are entered in square brackets, this is done to ensure the expression for the Nusselt Number is dimensionless.

You may get a slightly higher or lower value for the Nusselt Number here.

We would like to compare the Nusselt Number along the heated section of the pipe. We can generate the Nusselt Number at a different location by simply changing the x-coordinate of ** exit** and

**, which we defined earlier. Once the new coordinates defined in**

*Tw exit***and**

*exit***are updated, the associated expression**

*Tw exit***,**

*Tw***, and**

*Tm***will be updated automatically.**

*Nu exp*We can expect a maximum and dominant convection heat transfer at the entrance of the heated section of the pipe. The convection heat transfer raises the temperature inside the pipe, as well as mean temperature, along the downstream direction. The mean temperature near the exit is higher relative to the entrance and therefore a lower convection heat transfer is expected at the exit. Again, the Nusselt Number is a measure of convection heat transfer relative to conduction heat transfer. Thus we should expect the Nusselt Number to decrease along the length of the pipe.

To export the data, click on the "export" button. Comma Seperated Value (.csv) is able to be read by matlab and Excel, so it should be fine.

We are now ready to validate and verify our results.

**Go to Step 7: Verification & Validation**

Go to all FLUENT Learning Modules