Wednesday 1 December 2010

Thermal - Mixed Boundary Example (Conduction/Convection/Insulated)


Thermal - Mixed Boundary Example (Conduction/Convection/Insulated)


Introduction

This tutorial was created using ANSYS 7.0 to solve simple thermal examples. Analysis of a simple conduction as well a mixed conduction/convection/insulation problem will be demonstrated.
The Mixed Convection/Conduction/Insulated Boundary Conditions Example is constrained as shown in the following figure (Note that the section is assumed to be infinitely long):

Preprocessing: Defining the Problem

  1. Give example a Title
  2. Open preprocessor menuANSYS Main Menu > Preprocessor
    /PREP7
  3. Create geometryPreprocessor > Modeling > Create > Areas > Rectangle > By 2 Corners > X=0, Y=0, Width=1, Height=1
    BLC4,0,0,1,1
  4. Define the Type of Element
  5. Preprocessor > Element Type > Add/Edit/Delete... > click 'Add' > Select Thermal Mass Solid, Quad 4Node 55
    ET,1,PLANE55
    As in the conduction example, we will use PLANE55 (Thermal Solid, Quad 4node 55). This element has 4 nodes and a single DOF (temperature) at each node. PLANE55 can only be used for 2 dimensional steady-state or transient thermal analysis.

  6. Element Material Properties
  7. Preprocessor > Material Props > Material Models > Thermal > Conductivity > Isotropic > KXX = 10
    MP,KXX,1,10This will specify a thermal conductivity of 10 W/m*C.

  8. Mesh Size
  9. Preprocessor > Meshing > Size Cntrls > ManualSize > Areas > All Areas > 0.05
    AESIZE,ALL,0.05

  10. Mesh
  11. Preprocessor > Meshing > Mesh > Areas > Free > Pick All
    AMESH,ALL

Solution Phase: Assigning Loads and Solving


  1. Define Analysis Type
  2. Solution > Analysis Type > New Analysis > Steady-State
    ANTYPE,0

  3. Apply Conduction Constraints
  4. In this example, all 2 sides of the block have fixed temperatures, while convection occurs on the other 2 sides.

    • Solution > Define Loads > Apply > Thermal > Temperature > On Lines
    • Select 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

  5. Apply Convection Boundary Conditions

    • Solution > Define Loads > Apply > Thermal > Convection > On Lines
    • Select the right side of the block.The following window will appear:

    • Fill in the window as shown. This will specify a convection of 10 W/m2*C and an ambient temperature of 100 degrees Celcius. Note that VALJ and VAL2J have been left blank. This is because we have uniform convection across the line.

  6. Apply Insulated Boundary Conditions

    • Solution > Define Loads > Apply > Thermal > Convection > On Lines
    • 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
    You should obtain the following:


  7. Solve the System
  8. Solution > Solve > Current LS
    SOLVE

Postprocessing: Viewing the Results


  1. Results Using ANSYS
  2. Plot Temperature
    General Postproc > Plot Results > Contour Plot > Nodal Solu ... > DOF solution, Temperature TEMP

ANSYS Command Listing

/title, Simple Convection Example
/PREP7

! define geometry

length=1.0
height=1.0
blc4,0,0,length, height ! area - one corner, then width and height

! mesh 2D areas

ET,1, PLANE55 ! Thermal element only
MP,KXX,1,10 ! 10 W/mC
MAT,1
TYPE,1
ESIZE,length/20 ! number of element sub-divisions/side
AMESH,ALL

FINISH
/SOLU

ANTYPE,0 ! STEADY-STATE THERMAL ANALYSIS

! fixed temp BC's

NSEL,S,LOC,Y,height ! select nodes on top with y=height
D,ALL,TEMP,500 ! apply fixed temp of 500C
NSEL,ALL
NSEL,S,LOC,X,0 ! select nodes on three sides
D,ALL,TEMP,100 ! apply fixed temp of 100C
NSEL,ALL

! convection BC's

NSEL,S,LOC,X,length ! right edge
SF,ALL,CONV,10,100 ! apply fixed temp of 100C
NSEL,ALL

! Insulated BC's

NSEL,S,LOC,Y,0 ! bottom edge
SF,ALL,CONV,0 ! insulate edge
NSEL,ALL

SOLVE
FINISH

/POST1
PLNSOL,TEMP,,0, ! contour plot of temperatures

0 comments:

Post a Comment