Simulation Laboratory M.E.Thermal Engineering] 5 | GCE, Salem - 11 Ex.No : 01 STEADY STATE CONDUCTION IN SOLIDS Date: Aim: To find the Temperature Distribution of the given solid by using ANSYS 10.1 Problem Description: 1) Length of the solid is 1cm. 2) Height of the solid is 1cm. Boundary conditions: Top is maintained Problem Description: Boundary conditions: 1) Top is maintained at 500C. 2) Other sides are maintained at 100C. 3) Thermal Conductivity (K) of the material is 10 W/mc Figure:
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 01STEADY STATE CONDUCTION IN SOLIDS
Date:
Aim: To find the Temperature Distribution of the given solid by using ANSYS 10.1
Problem Description:
1) Length of the solid is 1cm.2) Height of the solid is 1cm.
Boundary conditions:
Top is maintained Problem Description:
Boundary conditions:
1) Top is maintained at 500C.2) Other sides are maintained at 100C.3) Thermal Conductivity (K) of the material is 10 W/mc
Figure:
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
PROCEDURE:STARTING ANSYS: Click on ANSYS 10.1 in the programs menu
MODELLING THE STRUCTURE:1. Go to the ANSYS main menu2. Create geometry.
Now go to General Postproc > Plot Results > Vector Plot > Predefined and selectVelocity.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
RESULTS:
Thus the temperature distribution of the given solid cylinder is determined.
Ex.No : 03STEADY STATE RADIATION IN SOLIDS
Date:
Aim :
To find the temperature distribution of the given solid by using ANSYS 10.1
Problem Description:
1. Length of the strip is 2cm.2. Height of the strip is 1cm.
Boundary condition:
1. Left and Right side temperature is 900 C and Emissive is 0.7.2. Convection Co-efficient of the top is 100W/m2C and Bulk temperature is 50C.3. Stephen Boltze` man constant is 5.67e-8W/m2K4.
Material properties:
1. Thermal conductivity is 3W/m C.2. Density is 1600 kg/m3.3. Specific heat is 800 J/kg C.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Figure:
PROCEDURE:
STARTING ANSYS:
Click on ANSYS 10.1 in the programs menu
MODELLING THE STRUCTURE:
1. Go to the ANSYS Main Menu2. Create geometry
Pre-processor > Modeling > create > Areas > Rectangle > By 2 Corners > X =0, Y=0,Width=2, Height = 1.
ELEMENT PROPEERTIES
SELECTING ELEMENT TYPE:
Define the Type of element
Preprocessor > Element type > Add/Edit/Delete...> click `Add` > select ThermalMaterial properties
Preprocessor > Material Props > Material models > Thermal >Conductivity>Isotropic > kxx =3 (Thermal conductivity)
Specific heat =800Density =1600.
MESHING
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
DIVIDING THE SHANNEL INTO ELEMENTS:
Mesh size
Pre-processor > Meshing > Size cntrls > Manual size > Areas > All Areas > 0.05
Mesh
Pre processor > Meshing > Mesh > Areas > Free > Pick All
BOUNDARY CONDITIONS AND CONSTRAINTS
Analysis Type
Solution > Analysis Type > New Analysis > steady state
Constraints
Solution > Define Loads > Apply
Thermal > Temperature > On lines
Fill the window in as shown to constraint the side to constant temperature of 900
Solution > Define Loads > Apply
Thermal > Convection > Online
Fill the window in as shown to constrain the side to film coefficient of 100 and bulk Temperature of 50.
Solution > Define Loads > Apply
Thermal > Radiation > Online
Fill the window in as shown to constrain the side to a emissivity of 0.7 and enclosure no is2.
Preprocessor > Radiation opts > solution opt
Fill the window in as shown to constrain the side to a Stephen Boltze man constant of 5.67e-8.7 and enclosure no is 2.
Main menu > Radiation Opt > Matrix method > Other settings
Fill the window in as shown to constrain the side to a Stephen Boltze man constant of 5.67e-8.7 And type of geometry is 2D.
Thus the temperature distribution of the given Strip is determined.
Ex.No : 05CONDUCTION HEAT TRANSFER IN A PLATE
Date:
Aim :
To solve the 2-D heat conduction problem below, using ANSYS. And find the
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Temperature distribution.
Problem description:Rectangular plate based on the specified temperatures on the plate edges, and the plate
Dimensions. Note: Thermal conductivity of the plate, Kxx=401 W/(m-K)
Procedure:
Define element type:
Preprocessor -> Element Type -> Add/Edit/Delete> Thermal Solid> Quad 4node 55
Define Material Properties:
Preprocessor -> Material Properties -> Material Models> Thermal
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
ModelingPreprocessor -> Modeling-> Create -> Areas -> Rectangle -> By DimensionsFill in the fields, (X1, X2) as 0, 10 & (Y1, Y2) as 0, 20 and then click “OK”.
Specify mesh density controls.
We will specify numbers of element divisions along lines. Choose: Preprocessor -> -Meshing- Size Controls -> Manual Size -> Lines-> Picked Lines
The picking menu (below left) appears. On the graphics window, click on the bottomHorizontal line (this is one of the 10 meter lines), to highlight it. Then, click “OK” in thePicking menu. Then, the “Element Size” menu (below right) appears. Enter “10” for“NDIV”, as shown, then click “OK”. Now, repeat the above process to specify 20 divisionsalong either of the vertical lines.
Mesh the rectangle to create nodes and elements.Preprocessor -> Meshing -> Mesh -> Areas -> Mapped -> 3 or 4 Sided
A picking menu appears. Select “Pick All”. The rectangle will be meshed.
Solution:
Apply temperatures around the edges:Solution -> Define Loads-> Apply -> Thermal-> Temperature -> On Lines
A picking menu appears. Highlight the two vertical lines (the 20 meter lines), which have atemperature of 100 C, then click on “OK” in the picking menu. The box on the next page Appears. Highlight “TEMP” for “DOFs to be constrained”, and enter “100” for “VALUE”.Repeat the above process to apply the 100 C temperature to the bottom horizontal line, but inthis case, choose “Yes” for “Apply TEMP to endpoints?” Repeat the process once more, toapply the 200 C temperature to the top horizontal line, but in this case, choose “No” for“Apply TEMP to endpoints?”
Now, to address the fact that two corners do not have a specified temperature, as anApproximation, we will set the temperature at these to corners to 150 C.
Solution -> Define Loads-> Apply -> Thermal-> Temperature -> On Key points
A picking menu appears. Note that there are four “key points” in the model, one at each
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
corner of the large rectangular area. Click on the upper two corners, at the intersections of the100 C and 200 C lines. When these corner “key points” are highlighted, choose “OK” in thepicking menu, and the following box appears: Click on “TEMP” for “DOFs to beconstrained” and enter “150” for VALUE, then click “OK”.
Solve the problem:
Solution ->-Solve -> Current LS
Post processing:Plot the temperature distribution:
General Post proc -> Plot Results -> Contour Plot-> Nodal SolutionThe box below appears. Click on “DOF solution”, then “Temperature”, and then click “OK”.
List the temperatures at the middle of the plate:Utility Menu > Select > EntitiesSelect Nodes and By Location from the scroll down menus. Click X coordinates and type 5into the input box , then click OK. Then go to
General Post proc -> List Results -> Nodal SolutionIn the box that appears, click on “DOF Solution” and “Temperature”, as shown, then click“OK”.
Write down the results:
SI no Node number Temperature SI no Node number Temperature
Thus the given 2-d heat conduction problem is solved using ansys. and the temperature
distribution within the rectangular plate, based on the specified temperatures on the plate edges, and the
plate dimensions was found.
Ex.No : 06 THERMAL ANALYSIS USING MIXED BOUNDARY CONDITIONS
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Date:
Aim: To solve the 2-D heat conduction and convection problem below using ANSYS, and find the temperature
distribution.
Problem Description:
The Mixed Convection/Conduction/Insulated Boundary Conditions Example is constrained as shown in thefollowing figure (Note that the section is assumed to be infinitely long):
Procedure:Create geometry
Preprocessor > Modeling > Create > Areas > Rectangle > By 2 Corners > X=0, Y=0,Width=1, Height=1Define the Type of Element
Preprocessor > Element Type > Add/Edit/Delete... > click 'Add' > Select Thermal Mass Solid, Quad 4Node55ET, 1, PLANE55
Material Properties
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic >KXX = 10This will specify a thermal conductivity of 10 W/m*C.Element
Mesh SizePreprocessor > Meshing > Size Cntrls > Manual Size > Areas > All Areas > 0.05
MeshPreprocessor > Meshing > Mesh > Areas > Free > Pick All
Solution Phase: Assigning Loads and SolvingDefine Analysis Type
Solution > Analysis Type > New Analysis > Steady-State
Apply Conduction ConstraintsIn this Problem, all 2 sides of the block have fixed temperatures, while convection occurs on the other 2 sides.
Solution > Define Loads > Apply > Thermal > Temperature > On LinesSelect the top line of the block and constrain it to a constant value of 500 C Using the same method, constrain the left side of the block to a constant value of 100 C
Select the right side of the block.Fill in the window. This will specify a convection of 10 W/m2*C and an ambient temperature of 100 degrees Celsius. Note that VALJ and VAL2J have been left blank. This is because we have uniform convection across the line.Apply Insulated Boundary Conditions
Select the bottom of the block.Enter a constant Film coefficient (VALI) of 0. This will eliminate Convection through the side, thereby modeling an insulated wall. Note: you do not need to enter a Bulk (Or
ambient) temperature.
Solve the SystemSolution > Solve > Current LSSOLVE
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Post processing:Viewing the ResultsPlot Temperature
General Post proc > Plot Results > Contour Plot > Nodal Solu ... > DOF solution, Temperature
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
TEMP
List the temperatures at the middle of the plate:
Utility Menu > Select > Entities
Select Nodes and By Location from the scroll down menus. Click X coordinates and type0.5 into the input box , then click OK. Then go to
General Postproc -> List Results -> Nodal Solution
In the box that appears, click on “DOF Solution” and “Temperature”, as shown, then click“OK”.
Write down the results
SI no Node number Temperature SI no Node number Temperature
Thus the given 2-D heat conduction and convection problem is solved using ANSYS.
And the temperature distribution within the component based on the specified temperatures
on the edges, and the dimensions, was found.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 07LAMINAR FLOW OVER FLAT PLATE
Date:
Aim:
To find the Temperature Distribution of the Laminar flow over flat plate.
Problem Description:
Dimensions:
1) Length of the plate is 1 m long.2) The flow area is 1 m * .25 m. .
Boundary conditions:
1) The velocity of the air at infinite distance from the plate is 0.5 m/s.2) Atmospheric pressure is assumed on all faces except the face where velocity is
input into the system.
Figure:
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
PROCEDURESTARTING ANSYS
1) Click on ANSYS 10.1 in the programs menuMODELING THE STRUCTURE: Click
Click the right line and then OK. The following window will now appear.
Enter 0 for the constant pressure value for these faces and click OK. This sets the pressure to Atmospheric.
SOLUTION
Go to ANSYS Main Menu>Solution>Flotran Set Up>Execution Ctrl.
Go to Solution>Run FLOTRAN.
POST-PROCESSING
Plot the temperature distribution
General Postproc>Plot Results>Contour Plot>Nodal Solution.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
TEMPERATURE DISTRIBUTION
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
VECTOR PLOT
Now go to General Postproc>Plot Results>Vector Plot>Predefined andselect Velocity. Enter a scale factor (VRATIO) of 0.4. The vector plot looks like this:
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
VECTOR PLOT
RESULT :
Thus the temperature distribution of the laminar flow over the flat plate is determined.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 08COMPOSITE WALL
Date:
Aim:
To find the Temperature distribution of the given composite wall by using ANSYS 10.1
Problem Description:
Dimensions:
1) Length = 3 m
2) Width = 3 m
3) Thickness of each Layer = 1 m
Boundary Conditions:
1) The left side of the block has a constant temperature of 400 K.
2) The right side of the block has convection (h=20 W/m-K ; T= 300 K)
3) The Al section generates heat at a rate of 200 W/m3
4) The He section absorbs heat at a rate of 175 W/m3
Figure
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
PROCEDURE:
STARTING ANSYS:MODELING THE STRUCTURE:
1) Go to the ANSYS main menu2) Create geometry
Preprocessor > Modeling > Create > Areas > Arbitrary > Through KPs . Length = 3m,Width =3m.
ELEMENT PROPERTIES
Aluminum(1stlayer): KAl= 235 W/m*K
Helium(2ndlayer): KHe= 0.1513 W/m*K
Copper(3rdlayer): KCu= 400 W/m*
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
1) Now that we’ve defined what material ANSYS will be analyzing, we have to definehow ANSYS should analyze our block.
2) Click Preprocessor>Element Type>Add/Edit/Delete... In the 'Element Types'window that opens click on Add... The following window opens:
3) Type 1 in the Element Type reference number.
4) Click on Thermal Mass>Solid and select Quad 8node 77. Click OK. Close the'Element Types' window.
5) Now we have selected Element Type 1 to be a Thermal Solid 8node Element.
6) This finishes the section defining how the part is to be analyzed.
7) Now we have selected Element Type 1 to be a Thermal Solid 8node Element.
8) This finishes the section defining how the part is to be analyzed.
SELECTING ELEMENT TYPE:
Define the Type of Element
Preprocessor > Element Type > Add/Edit/Delete…>click Add > Thermal Mass >Solid > select Quad8node 77.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Element material properties
Preprocessor > Material props> Material models > Thermal > conductivity >Isotropic.
1) Fill in 235 for Thermal conductivity. Click OK. This is the Thermal Conductivity of Al.
2) Now repeat the steps of clicking Thermal>Conductivity>Isotropic and then Defining the Thermal Conductivityas 0.1513 for the Model 2.
3) You have now defined the k value of Helium.
4) Define the last section and this time use K = 400. This is the Thermal Conductivity of Copper.
And select the bottom line of the rectangle. Enter 0 inthe VY field. Leave VX and VZ blank.
Repeat the above and set the velocity to ZERO for the 2 lines that make up thehalf – circle. (VX =
VY=0,VZ can be left blank). This is due to the no –slip condition along the cylinder.
The last step is to apply atmospheric pressure to the outlet of the fluid region. Go to
Preprocessor > Loads > Define Loads > Apply >Fluid CFD > PressureDOF > On Lines and
select the right most vertical line and click OK. In the Window that appears, enter 0 for the constant pressure value
and click OK.
SOLUTION
Go to ANSYS main menu > Solution > Flotran set up > Execution Cntrl. The following window appears. Change the first input value to 1000, as shown. No other changes are
needed. Click OK.
Go to Solution > Run FLOTRAN.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
POST PROCESSING
1) Plotting the velocity distribution.
2) Go to General post proc > Read Results > Last set.
3) Then go to General Postproc > Plot Results > Contour plot > Nodal solution
select VSUM and click OK
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
TEMPERATURE DISTRIBUTION
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
VECTOR PLOT:
NOW GO TO General Post proc > Plot Results > vector Plot > Predefined and
select velocity. The vector plot looks like this.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
VECTOR PLOT
RESULTS:
Thus the velocity distribution of a flow of air over an infinite cylinder is computed andPlotted.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 10STEADY STATE RADIATION IN SOLID CYLINDER
Date:
Aim:
To find the Temperature Distribution of the given solid cylinder by using ANSYS 10.1
Problem Description:
1. Radius of the solid is 1m.2. Height of the solid is 1m.
Boundary conditions:
1. Top side is maintained at 700°C & emissivity 0.82. Bottom side is maintained at 500°C & emissivity 0.43. Circumference is maintained at 400°C & emissivity 1(Black body)
Now go to General PostProc > Plot Resuots > Vector Plot > Predefined and selectvelocity.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
RESULTS:
Thus the temperature distribution and thermal flux of the given fin is determined.
Ex.No : 12STEADY STATE CONVECTION IN COMPOSITE HOLLOW CYLINDER
Date:
Aim:
To find the temperature distribution of the given composite hollow cylinder by using ANSYS 10.1
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Problem Description:
1. Inner diameter of the cylinder 0.025m2. Outer diameter of the cylinder 0.038m3. Inner diameter of the cylinder (insulator) 0.038m4. Outer diameter of the cylinder (insulator) 0.058m5. Depth 2m
Boundary condition:
6. Hot gas flow inside the cylinder 330 C & h=400 W/m2 C7. Outer surface of the insulator is at 30C & h=60 W/m2 C
Now go to General PostProc > Plot Results > Vector Plot > Predefined and Thermal Flux
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
RESULTS:
Thus the temperature distribution of the given composite hollow cylinder is determined.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 13COMBINED CONDUCTION AND CONVECTION IN TRIANGULAR FIN
Date:
Aim:
To find the Temperature distribution of the given triangular fin by using ANSYS 10.1
Problem Description:
1. Length of the fin is 0.04 m.2. Side of the fin is 0.005m.3. Wall temperature is 400 C4. Convection Co-efficient of the wall is 90 W/m2C5. Bulk Temperature 50 C
Figure
PROCEDURE;
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
STARTING ANSYS:
Click on ANSYS 10.1 in the program menu.
MODELLING THE STRUCTURE:
1) Go to the ANSYS main menu.2) Create geometry.
Preprocessor > Modeling > create > Volume> Prism > By Side Lengths> Z=0.04, Side=0.005
ELEMENT PROPERTIES
SELECTING ELEMENT TYPE:
Define the type of element
Preprocessor > Element Type > Add/Edit/Delete >click `Add` > select
Thermal Mass Solid, 20 node90
Element material properties
Preprocessor> Material props > Material Models > Thermal conductivity>Isotropic > KXX = 54 (Thermal conductivity)
MESHING
DIVIDING THE CHANNEL INTO ELEMENTS:
Mesh Size
Preprocessor > Meshing > size cntrls > Manual size > Areas > All areas > 0.05
General PostProc > Plot Results > Vector Plot > Predefined and Thermal Flux
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
RESULT:
Thus the temperature distribution and thermal flux of the given fin is determined.
Simulation Laboratory
M.E.Thermal Engineering] 5 | GCE, Salem - 11
Ex.No : 14RADIATION HEAT TRANSFER IN TRIANGULAR FIN
Date:
Aim:
To find the Temperature distribution of the given triangular strip by using ANSYS 10.1
Problem Description:
1. Radius of the strip is 1 m.2. Angle of the strip is 60.3. Wall temperature at the sides are 900 C and 400 C respectively4. Emissivities of the sides are 0.8 for both sides.5. Heat transfer at the bottom side is zero6. K=30 W/m C