Finite Element Model of Transient Heat Conduction with Isothermal Phase Change (Two and Three Dimensional).

Abstract

The partial differential equation for transient heat conduction is solved by a finite element analog using a quadratic weighting function for the discretized spatial domain. The transient problem is solved by the Crank-Nicolson approximation. Two-dimensional and three-dimensional models incorporated in the same computer program are presented. The finite element method is reviewed, assumptions and limitations upon which the model is based are presented, and a complete derivation of the system analog is included. Certain problems can only be modeled as a three-dimensional system, e.g. thaw degradation around roadway culverts, embankment dams on permafrost where dam length is short relative to dam width, and thaw and freezeback under buildings. In most cases, however, the more economical two-dimensional model can be used. Numerical tests of both models have been accomplished but field verification has not been attempted. A user's manual and a FORTRAN IV computer listing of the program are presented. (Author)

Open PDF

Document Details

Document Type
Technical Report
Publication Date
Nov 01, 1977
Accession Number
ADA047369

Entities

People

  • G. L. Guymon
  • T. V. Hromadka Ii

Organizations

  • Cold Regions Research and Engineering Laboratory

Tags

Communities of Interest

  • Air Platforms
  • Energy and Power Technologies
  • Space

DTIC Thesaurus Topics

  • Cold Regions
  • Computational Fluid Dynamics
  • Computational Science
  • Computer Programs
  • Computer Simulations
  • Computers
  • Differential Equations
  • Engineers
  • Finite Element Analysis
  • Heat Energy
  • Latent Heat
  • Partial Differential Equations
  • Plastic Explosives
  • Self Assembly
  • Thermal Conductivity
  • Three Dimensional
  • Two Dimensional

Readers

  • Computational Modeling and Simulation
  • Finite Element Method (FEM) for solving Partial Differential Equations (PDEs)
  • Hydraulic Engineering.