Include Page |
---|
...
|
...
|
Step 6: Results
Please make sure your project is saved in Workbench. Double click on Results in the Project Schematic window. This will open CFD-Post (the program used to analyze results from FLUENT computation.)
Overview
You may have noticed in previous sections, that the pipe looks extremely long and thin on the screen. In fact, due to the symmetry condition of the experiment, we have only modeled half the pipe in our analysis. To be able to make full use of the results, we must:
1) Generate the results for the parameter investigated (e.g. temperature, pressure, velocity)
2) Mirror the result to reflect the result of the full pipe
3) Adjust the scaling of the resulting graphs so that the variation of the pipe in y direction is more pronounced.
Temperature Contour
Our first challenge is the temperature contour. On the top menu, click on contour . We will be calling this contour "Temperature Contour", OK when done. On the left hand side, Details of Temperature Contour will allow you to select parameters relevant to the results we're looking for. In this example, the Locations is periodic 1, the Variable is Temperature. The number of contours is a personal preference, in this example, we have selected 100. This step tells CFD-Post we are looking to plot contours of temperature.
...
The next step is to mirror the image, this will make the results more intuitive and easier to understand. From the previous screen, select the View tab. This tab will allow us to adjust the appearance of the contour plot we have just generated. Check Apply Reflection/Mirroring. Select ZX Plane for Method. Choosing this option reflects the current model in the ZX plane and allows us to view the "full" pipe.
Finally, we want to view the variation of parameters along the y axis of the pipe. This means we will need to "stretch out" the plot. Select Apply Scale. Enter 30 for y-axis. This will stretch our model in the y direction. Click Apply.
After you click Apply, you will see that under Outline > User Locations and Plots, Temperature Contour is created. You will also see that the Temperature Contour is plotted in the Graphics window on the right. Under Outline > User Locations and Plots, uncheck Wireframe to see just the Temperature Contour in the Graphics window.
...
Include Page | ||||
---|---|---|---|---|
|
Numerical Results
Note |
---|
Cornell MAE 4272 Fall 2020: For the canvas quiz, you should use the FLUENT inputs from this tutorial. Later, you will need to repeat the FLUENT simulation with inputs from YOUR MEASUREMENTS in the lab and compare the FLUENT results with the experiment. |
Temperature Contour
HTML |
---|
<iframe width="560" height="315" src="https://www.youtube.com/embed/7_P7arvK-4Q?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the Above Video:
- Create a Contour Plot and name it Temperature Contours
- Under Geometry > Locations, choose periodic 1 and set the variable of interest as Temperature.
- Under View > Apply Scale, enter (1, 30, 1)
- Default Transform
- Apply Reflection about ZX Plane
You can save the image to a file using the camera icon indicated in the image below or using the Snipping Tool in Windows (you can search for it under Start > Programs).
In developing the experiment, it was assumed that by the end of the adiabatic mixing stage, the flow will be well-mixed. Do the results from the numerical solution simulation support this assumption?
Velocity Vectors
Our next challenge is to produce velocity vectors. This is a very similar process to creating the temperature contours above. On the top menu, click on vector . Please name it "Velocity Vector" and click OK. Under Details of Velocity Vector, select periodic 1 for Locations. Select Velocity for Variable. This tells CFD-post we are looking for vector plots of velocity.
In the next step, we will specify the appearance of vector arrows. Select the Symbol tab. Enter 0.05 for Symbol Size. This again is dependent on personal preference.
Finally click Apply. You will see that under Outline > User Locations and Plots, Velocity Vector is created. Uncheck Temperature Contour so that Graphics window shows just the Velocity Vector plot.
It would be beneficial to repeat the previous steps involved with mirroring and stretching the plot:
...
https://confluence.cornell.edu/download/attachments/111221576/velocity%20vector.png
Please focus on the first section of the pipe as it shows flow development. Can you see at which point the flow becomes fully developed?
Centerline Temperature Plot
Now let's look at the temperature variation along the center-line of the pipe. In the experiment, we are only able to measure the temperature at two points. First, at the inlet of the pipe and second, after the adiabatic mixing stage. The simulation can show us the variation of temperate in between these two points.
To create the desired plot:
Insert > Chart
Please name this chart "Centerline Temperature". You will see Details of Centerline Temperature appear on the lower left panel. Select the General tab and name the chart "Temperature Variation along Pipe Axis".
Moving on, please select the Data Series tab. This tab will help us specify the source of the chart data. Create a new data series . Change the name from Series 1 to FLUENT. Under Data Source, specify Centerline as Location. Click Apply. This should result in a plot like:
Wiki Markup \[insert the picture here\]
We would like to verify our simulation results, to do this we are comparing to experimental data. Experimental data is can be downloaded here. Download it to a directory of your choice. Now, click a new data series . Name it "Experiment". Under Data Source, select File and browse for the downloaded experimental data.
Now that we have our data sources, we will proceed by specifying the axes. Starting with the X-axis:
Click on X Axis tab. Next to Variable, choose X.
Now we will specify the Y Axis parameter. Click on Y Axis tab. Next to Variable, choose Temperature.
Now we will specify how we want to the chart to display. The default setting is to display the data series in lines. Since we only have 3 experimental points, we want them to be displayed in data points. Click on Line Display. Then click on experimental tab. Next to Line Style, change Automatic to None. Next to Symbols, change None to Diamond. Change the color to red. Click Apply. You will see Centerline Temperature created under Report in the Outline tab.
This is what you should see in the Graphics window.
...
https://confluence.cornell.edu/download/attachments/111221576/Centerline%20Temperature.png
...
Pressure Plot
Now let's us look at the pressure variation at the centerline. First, we will create a line call centerline.
Insert > Location > Line
Name it "Centerline" and click OK. On the lower left panel, you will see Details of Centerline. Enter the following coordinates.
Point 1 (0,0,0)
Point 2 (6.096,0,0)
Enter 50 for Samples. (This will be the number of sample points used when plotting data)
Click Apply.
You will see centerline created under User Locations and Plots.
Next, we will create a chart using this Location data.
Insert > Chart
Enter "Axial Pressure" as Name. You will see Details of Axial Pressure appear on the lower left panel. Under General, give the chart Title as "Pressure Variation along Pipe Axis".
Now click on Data Series tap to specify the location of the chart data. Create a new data series . Change the name from Series 1 to FLUENT. Under Data Source, specify Centerline as Location.
We would also like to compare our simulation result with experimental data. Experimental data is can be downloaded here. Download it to the directory that you like. Now, click a new data series . Name it Experiment. Under Data Source, select File and browse to the downloaded experimental data.
Now we will specify the X Axis parameter. Click on X Axis tab. Next to Variable, choose X.
Now we will specify the Y Axis parameter. Click on Y Axis tab. Next to Variable, choose Pressure.
Now we will specify how we want to the chart to display. The default setting is to display the data series in lines. Since we only have 3 experimental points, we want them to be displayed in data points. Click on Line Display. Then click on experimental tab. Next to Line Style, change Automatic to None. Next to Symbols, change None to Diamond. Change the color to red. Click Apply. You will see Axial Pressure created under Report in the Outline tab.
This is what you should see in the Graphics window.
...
https://confluence.cornell.edu/download/attachments/111221576/centerline%20pressure.png
Does the simulation result compares well with the experimental data?
Wall Temperature Plot
Now let's us look at the temperature variation along the wall. First, we will create a line call wall.
Insert > Location > Line
Name it "Wall" and click OK. On the lower left panel, you will see Details of Wall. Enter the following coordinates.
Point 1 (0,0.0294,0)
Point 2 (6.096,0.0294,0)
Enter 50 for Samples. (This will be the number of sample points used when plotting data)
Click Apply.
You will see wall created under User Locations and Plots.
Next, we will create a chart using this Location data.
Insert > Chart
Enter "Wall Temperature" as Name. You will see Details of Wall Temperature appear on the lower left panel. Under General, give the chart Title as "Wall Temperature".
Now click on Data Series tap to specify the location of the chart data. Create a new data series . Change the name from Series 1 to FLUENT. Under Data Source, specify Wall as Location.
...
HTML |
---|
<iframe width="560" height="315" src="https://www.youtube.com/embed/HX27T78mQfg?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the Above Video:
- Click Vectors > Name Velocity Vectors
- In Details plot Velocity
- In Variable > Color by Variable
- Color by temperature
- Color by temperature
We see that the flow speeds up as the density decreases in order to keep the same mass flow rate.
Does the flow become fully developed at the end of the first section?
Wall Temperature Plot
HTML |
---|
<iframe width="560" height="315" src="https://www.youtube.com/embed/-Cvn7HPp9eY?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the Above Video:
- Create location > Line
- Point 1: (0,.0294,0) Point 2: (6.045,.0294,0)
- Create Chart
- Under Data Series > Location Select line you just created
- X Tab
- Plot X location
- Y Tab
- Plot Temperature
The wall temperature values can be exported to a csv file (suitable for Excel) using the Export button highlighted below.
Importing Experimental Data and Chart Formatting
We can import experimental data from a csv file using the procedure shown below. The csv file should contain two columns, corresponding to the x-axis and y-axis values, respectively
HTML |
---|
<iframe width="560" height="315" src="https://www.youtube.com/embed/oy3xJrmWNLo?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the Above Video:
- In details of Chart > New Series > Import Data from file
- Find file and click okay
- Change how experimental data looks in CFD Post
- In Line Display under Series 2 > No Line > Symbols > Ellipse
- Chart Details > General
- Give Title
- Chart Details > Chart Display
- Font Sizes
- Grid Sizes
Note that the centerline temperature and pressure variations can be plotted by duplicating this plot as mentioned in the video.
Mixed Mean Temperature at an Axial Location
From energy conservation, we can show that the mixed mean temperature is constant in the flow development and mixing sections and varies linearly in the heated section. This is shown schematically in the following figure from the MAE 4272 lab manual.
The slope of Tm in the heated section can be obtained from the following equation which is derived from energy balance in the heated section:
Using the above equation, calculate the mixed mean temperature Tm at x=2.67 m. Remember to add the inlet temperature, otherwise you will just end up with the temperature difference between the mixed mean temperature and the inlet temperature (where we assumed the flow was fully mixed). An alternate procedure to calculate Tm involves integrating the temperature profile. This procedure is covered in the Verification & Validation section in the video entitled Check Energy Conservation via Mixed Mean Temperature Variation. If energy is conserved in the FLUENT simulation, the values calculated using the two procedures should match.
Wall Temperature at an Axial Location
For calculating the Nusselt no. at an axial location, we need the wall temperature at that location. The wall temperature at an axial location can be calculated in two ways:
- By interpolating the Tw vs. x values exported to Excel from the wall temperature plot obtained above.
- By directly extracting the wall temperature at the desired location using the Probe function in the post-processor.
The following video shows you the procedure for extracting the wall temperature at x=2.67 m using the Probe function; this value can then be used to calculate the Nusselt number. To repeat the calculation at a different axial location, you can right-click on appropriate items in the tree, duplicate and modify as necessary. You need to double-click on an item in the tree to modify it; this is easy to overlook.
HTML |
---|
<iframe width="640" height="360" src="//www.youtube.com/embed/mv3uDy7ZuCY?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the above video:
- Create a line at x = 2.67
- Go to Expressions, right click and click New
- Name Tw267
- Right click in Definition box > Functions > CFD-Post > maxVal
- Right click within the parenthese of maxVal > Variables > Temperature
- Right click after the @ sign > Locations > x267
Nusselt No.
Cornell MAE 4272 students, Fall 2020: Please use the following procedure for calculating the Nusselt number at two axial locations in the heated section where the flow is thermally fully-developed.
The convective heat transfer coefficient, h can be determined from Newton's law of cooling:
The wall heat flux on the left hand side is known from the boundary condition. We have shown you how to get the wall temperature T_w at an axial location. Calculate the mixed mean temperature at the same axial location by evaluating the following integral using the procedure shown in the video in the Verification and Validation section. The title of the video is Check Energy Conservation via Mixed Mean Temperature Variation.
Once you determine h, you can non-dimensionalize it to get the Nusselt number.
The Nusselt number in the thermally fully-developed region should be the same at different axial locations. Why? Do your Nu values at the two axial locations compare reasonably well? Why or why not?
Wall Shear
We plot the wall shear using the procedure shown in the video below.
HTML |
---|
<iframe width="640" height="360" src="//www.youtube.com/embed/WiK1uBTdK-M?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the above video:
- Click on the Chart Viewer tab
- Click chart in the top toolbar
- Name it Wall Shear
- Click on Data Series Tab
- In Location dropdown menu, choose Wall
- Click on X Axis tab
- under Variable, choose X
- Click on Y Axis tab
- under Variable, choose Wall Shear X or Wall Shear
- Go to Location in the tree
- Double click on Wall
- Increase Sampling from 50 to 200
We then consider the trends in the wall shear in the heated, mixing and flow development sections and try to justify them through physical reasoning.
HTML |
---|
<iframe width="640" height="360" src="//www.youtube.com/embed/NbdGrJTesZg?rel=0" frameborder="0" allowfullscreen></iframe> |
You can spiff up your plot using the tips discussed below. This video also shows you how you could read in experimental results for comparing the wall shear between simulation and experiment.
HTML |
---|
<iframe width="640" height="360" src="//www.youtube.com/embed/6RNykoM86xA?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the above video:
- To edit how the Wall shear graph is displayed
- Select Wall shear in tree
- Click on General, check Title and enter title in Title blank
- in Data Series, enter series name in Name blank
- in X Axis, enter x axis label in Custom Label blank
- same for Y Axis
- in Line Display, uncheck "Use series..." and type within Legend Name blank
- in Chart Display, under Sizes, toggle with the line sizes and font
- To add another data series
- go to Data series tab of Wall Shear
- click on New button
- Scroll down, click File as Data Source > browse for your file
- To export the chart
- Click the button next to undo
- Select location to save
- Rename
When the simulation was repeated for conditions for which experimental data are available, we got the comparison shown below. The difference in the average wall shear in the heated section between the simulation and experiment is a respectable 4%. Note that the wall shear in turbulent flows is difficult to predict accurately due to the steep velocity gradients at the wall.
Fanning Friction Factor
The Fanning friction factor, also called the skin friction coefficient, is obtained by non-dimensionalizing the wall shear. It can be calculated and plotted using the procedure outlined below.
HTML |
---|
<iframe width="640" height="360" src="//www.youtube.com/embed/Jcht7CAPzAc?rel=0" frameborder="0" allowfullscreen></iframe> |
Summary of the above video:
- Create Location > Point
- (2.67, .0294, 0)
- Probe wall shear at point:
- Create Expression > Right Click > New Expression
- probe(Wall Shear)@w267 / 1.483 (Pa/K) /Tm267
- Can duplicate for different locations for a plot
Final Plots
HTML |
---|
<iframe width="560" height="315" src="https://www.youtube.com/embed/x-flaXXU7xg?rel=0" frameborder="0" allowfullscreen></iframe> |
Input Summary
In older versions of FLUENT, you can view the input summary (model, material properties, boundary conditions, etc) by clicking on Report in the menu bar of FLUENT. A small window will pop up and you can print the selected input summary directly in FLUENT.
Go to Step 7: Verification & Validation
Now we will specify the X Axis parameter. Click on X Axis tab. Next to Variable, choose X.
...
https://confluence.cornell.edu/download/attachments/111221576/Wall%20Temperature.png
Does the simulation result compares well with the experimental data?
Axial Velocity Profile
Now let's us look at the axial velocity profile at various location in the pipe. We are interested in the flow development before the heated section and would like to observe how heat addition in the heated section will affect the flow development. To do this, let's start with investigation the flow before the heated section.
...
Axial Velocity Profile before Heated Section
...
The heated section is from 1.83m to 4.27m. Let's create 4 lines before 1.83m in the pipe.
Insert > Location > Line
Name it "Inlet" and click OK. On the lower left panel, you will see Details of Inlet. Enter the following coordinates.
Point 1 (0,0,0)
Point 2 (0,0.0294,0)
Enter 50 for Samples. (This will be the number of sample points used when plotting data)
Click Apply.
Create Preheat 1.
Insert > Location > Line
Name it "Preheat 1" and click OK. On the lower left panel, you will see Details of Preheat 1. Enter the following coordinates.
Point 1 (0.6,0,0)
Point 2 (0.6,0.0294,0)
Enter 50 for Samples. (This will be the number of sample points used when plotting data)
Click Apply.
Continue the same step for creating line Preheat 2 (x=1.2m), Preheat 3 (x=1.8m).
Check that you have the following under Outline.
Now we will have enough interval to look at the flow development before the heating. Let's create a chart to investigate this.
...
Now we will specify the X Axis parameter. Click on X Axis tab. Next to Variable, choose Velocity u. Next we will specify the Y Axis parameter. Click on Y Axis tab. Next to Variable, choose Y. Click Apply. You will see First Section Axial Velocity Profile created under Report in the Outline tab.
This is what you should see in the Graphics window.
...
https://confluence.cornell.edu/download/attachments/111221576/Axial%20Velocity%20Profile.png
...
Axial Velocity Profile before and after Heated Section
...
Let's create lines after heated section
Insert > Location > Line
Name it "Postheat 1" and click OK. On the lower left panel, you will see Details of Postheat 1. Enter the following coordinates.
Point 1 (4.27,0,0)
Point 2 (4.27,0.0294,0)
Enter 50 for Samples. (This will be the number of sample points used when plotting data) Click Apply.
Create Postheat 2.
Insert > Location > Line
Name it "Postheat 2" and click OK. On the lower left panel, you will see Details of Postheat 2. Enter the following coordinates.
Point 1 (5,0,0)
Point 2 (5,0.0294,0)
Enter 50 for Samples. Click Apply.
Continue the same step for creating line Outlet (x=6.096m)
Check that you have the following under Outline.
Now we will have enough interval to look at the flow development before and after the heating. Let's create a chart to investigate this.
Insert > Chart
Enter "Second Section Axial Velocity Profile" as Name. You will see Details of Second Section Axial Velocity Profile appear on the lower left panel. Under General, give the chart Title as "Axial Velocity Profile".
Now click on Data Series tap to specify the location of the chart data. Create a new data series. Under Data Source, specify Preheat 3 as Location. Change the name to x=1.8m. Continue adding Data Source until we added all Preheat 3, Postheat 1, Postheat 2, and Outlet. Name them according to the figure shown below.
Now we will specify the X Axis parameter. Click on X Axis tab. Next to Variable, choose Velocity u. Next we will specify the Y Axis parameter. Click on Y Axis tab. Next to Variable, choose Y. Click Apply. You will see First Section Axial Velocity Profile created under Report in the Outline tab.
This is what you should see in the Graphics window.
...
https://confluence.cornell.edu/download/attachments/111221576/Axial%20Velocity%20Profile%20after.png
The flow is accelerated in the heated section.
Temperature Profile
Now let's us look at the temperature profile before and after the heating section.
Insert > Chart
Enter "Temperature Profile" as Name. You will see Details of Temperature Profile appears on the lower left panel. Under General, give the chart Title as "Temperature Profile".
Now click on Data Series tab to specify the location of the chart data. Create a new data series. Under Data Source, specify Preheat 3 as Location. Change the name to x=1.8m. Similarly, add the locations: Preheat 3, Postheat 1, Postheat 2, and Outlet. Name them according to the figure shown below.
Now we will specify the X Axis parameter. Click on X Axis tab. Next to Variable, choose Temperature. Next we will specify the Y Axis parameter. Click on Y Axis tab. Next to Variable, choose Y. Click Apply. You will see Temperature Profile created under Report in the Outline tab.This is what you should see in the Graphics window.
...
https://confluence.cornell.edu/download/attachments/111221576/Temperature%20Profile%202.png
The plot shows temperature is nearly uniform at the outlet (end of mixing section).
Go to Step 7: Verification & Validation