Here is an example of MATLAB code that was developed for 1D conduction heat transfer problems. This code was developed for MECH 302 – Finite Element Analysis at Bucknell. Specifics of this code include the following:

• User inputs for:
• Geometry, conduction coefficient, and mesh density
• Distributed internal heat generation as a function of location
• Essential, natural, and mixed boundary conditions
• A Gaussian Quadrature MATLAB function
• A plot of temperature and the derivative of temperature as a function of location
• The value of the maximum nodal temperature
• The temperature value of a specific location as specified by the user
• Boundary heat flux

The code is not published online because it’s an assignment for the course, but please feel free to email me (b.wheatley@bucknell.edu) and I would be happy to share.

Here is an example of the code at work: ```Input domain length in [m]: 8
Input number of elements: 8
Input conduction coefficient k in [W/m-C]: 200
Input initial area [m^2]: 0.1
Input final area [m^2]: 0.5
Input distributed internal heat generation [W/m^3] as a function of x: 10*x
Specify boundary conditions as follows: 1 - known temperature, 2 - known flux, 3 - convection.
Specify left boundary condition type: 3
Specify left convection coefficient in [W/m^2-C]: 50
Specify left fluid temperature in [deg C]: 10
Specify right boundary condition type: 2
Specify right flux in [W]: 0
The maximum temperature is 46.26 [deg C].
The maximum temperature location is 8 [m].
Specify location of desired approximation in [m]: 7.5
The temperature at x = 7.5 [m] is 46.13 [deg C].
Heat flux at left boundary is -210.4 [W].
Heat flux at right boundary is 0 [W].  ```