663Project

Created: 2012-03-20 17:31
Updated: 2013-01-30 21:18
c++

README.md

Kenton's MecE 663 Project

This is a custom deal.ii application that has been developed to solve a heat transfer problem. While the problem would be extremely difficult to solve analytically due to the "complex" boundary conditions, the solved problem is quite simple for a finite-element analysis.

Problem Domain

The problem is a small rectangle with dimensions 2.5cm x 1cm, and is composed of 3 distinct boundary conditions. On y = 1, from x = 1 to x = 2.5, a constant temperature of 400K is applied. On x = 0, from y = 0 to y = 1, a heat flux in the x direction only that is applied according to the equation: B(1-y) where B is a constant.

There is a heat generation term in the given PDE, equal to Ax^2 where A is a constant.

Default Constants

Kx = 15 W/m^2
Ky = 25 W/m^2
A is in (0, 5) x 1000 W/m^5
B is 50 W/m^3
T (temperature at the Dirichlet boundary) = 400K

Implementation

The program is written with a static geometry and boundary locations / types, however the above constants may be changed through the command-line (as well as several other options pertaining to the finite-element formulation itself.

Command-line

The program may be invoked and produce a solution using the following syntax: ./project [OPTIONS] [SOLUTION_FILE.vtk] Where both [OPTIONS] and [SOLUTION_FILE.vtk] are optional parameters. [OPTIONS] follows unix-style command-line options, with the following values being accepted:

  • -? shows a help menu
  • -r refines the mesh times (defaults to 0)
  • -f write the generated mesh to file for visualization
  • -s actually solve the problem
  • -q <num|auto> allows you to manually specify the number of quadrature points used (defaults to auto)
  • -i set the maximum number of CG iterations performed for the solution (defaults to 1000)
  • -t set the tolerance for the residual in the CG solution (defaults to 1e-12)
  • -o set the Lagrange interpolation polynomial order (defaults to 1)
  • -x set the Kx value (defaults to 15)
  • -y set the Ky value (defaults to 25)
  • -A set the A value (defaults to 5000)
  • -B set the B value (defaults to 50)
  • -T set the boundary temperature at the Dirichlet boundary (defaults to 400)

If [SOLUTION_FILE.vtk] is specified, the solution will be graphically written to the specified file (in VTK format). If it is not specified and a solution is performed, the solution will be written to "solution.vtk", overwriting that file.

Cookies help us deliver our services. By using our services, you agree to our use of cookies Learn more