# Heat equation calculator differential equation

A partial differential diffusion equation of the form. Physically, the equation commonly arises in situations where is the thermal diffusivity and the temperature. Dividing both sides by gives. Anticipating the exponential solution inwe have picked a negative separation constant so that the solution remains finite at all times and has units of length.

The solution is. Since the general solution can have any. Now, if stanley wiki disney are given an initial conditionwe have. Multiplying both sides by and integrating from 0 to gives. Using the orthogonality of and. Following the same procedure as before, a similar answer is found, but with sine replaced by cosine:. Weisstein, Eric W. Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own. Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet. Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

MathWorld Book.

Terms of Use. Contact the MathWorld Team. Nonsteady-State Heat Conduction in a Cylinder.But its applicability is very limited. It is simply the rate equation in this heat transfer mode, where the temperature gradient is known.

### Separable Differential Equations Calculator

But a major problem in most conduction analyses is to determine the temperature field in a medium resulting from conditions imposed on its boundaries. In engineering, we have to solve heat transfer problems involving different geometries and different conditions such as a cylindrical nuclear fuel element, which involves internal heat source or the wall of a spherical containment. These problems are more complex than the planar analyses we did in previous sections.

Therefore these problems will be the subject of this section, in which the heat conduction equation will be introduced and solved. Note that heat flux may vary with time as well as position on a surface. In nuclear reactorslimitations of the local heat flux is of the highest importance for reactor safety. The heat conduction equation is a partial differential equation that describes the distribution of heat or the temperature field in a given body over time.

Detailed knowledge of the temperature field is very important in thermal conduction through materials. That is:. This equation is also known as the Fourier-Biot equationand provides the basic tool for heat conduction analysis. From its solution, we can obtain the temperature field as a function of time. At any point in the medium the net rate of energy transfer by conduction into a unit volume plus the volumetric rate of thermal energy generation must equal the rate of change of thermal energy stored within the volume.

The thermal conductivity of most liquids and solids varies with temperature. For vapors, it also depends upon pressure. In general:. From the foregoing equation, it follows that the conduction heat flux increases with increasing thermal conductivity and increases with increasing temperature difference. In general, the thermal conductivity of a solid is larger than that of a liquid, which is larger than that of a gas.

This trend is due largely to differences in intermolecular spacing for the two states of matter. In particular, diamond has the highest hardness and thermal conductivity of any bulk material. See also: Thermal Conductivity. These pellets are then loaded and encapsulated within a fuel rod or fuel pinwhich is made of zirconium alloys due to its very low absorption cross-section unlike the stainless steel.

The surface of the tube, which covers the pellets, is called fuel cladding. Fuel rods are base element of a fuel assembly. The thermal conductivity of uranium dioxide is very low when compared with metal uranium, uranium nitride, uranium carbide and zirconium cladding material. The thermal conductivity is one of parameters, which determine the fuel centerline temperature. This low thermal conductivity can result in localised overheating in the fuel centerline and therefore this overheating must be avoided.

Overheating of the fuel is prevented by maintaining the steady state peak linear heat rate LHR or the Heat Flux Hot Channel Factor â€” F Q z below the level at which fuel centerline melting occurs. Expansion of the fuel pellet upon centerline melting may cause the pellet to stress the cladding to the point of failure. Klimenko and V. MEI Press, ISBN â€”92â€”0â€”â€”7.

In other words, it is the measure of thermal inertia of given material. Additional simplifications of the general form of the heat equation are often possible. One of most powerful assumptions is that the special case of one-dimensional heat transfer in the x-direction.

In this case the derivatives with respect to y and z drop out and the equations above reduce to Cartesian coordinates :. In engineering, there are plenty of problems, that cannot be solved in cartesian coordinates.Wolfram Language Revolutionary knowledge-based programming language.

Wolfram Science Technology-enabling science of the computational universe. Wolfram Notebooks The preeminent environment for any technical workflows. Wolfram Engine Software engine implementing the Wolfram Language.

### Select a Web Site

Wolfram Data Framework Semantic framework for real-world data. Wolfram Universal Deployment System Instant deployment across cloud, desktop, mobile, and more. Wolfram Knowledgebase Curated computable knowledge powering Wolfram Alpha. Solve a Dirichlet Problem for the Laplace Equation. Solve a Dirichlet Problem for the Helmholtz Equation. Solve a Basic Sturm â€” Liouville Problem. Study the Vibrations of a Stretched String. Model the Flow of Heat in an Insulated Bar. Study Dispersion in Quantum Mechanics.

Observe a Quantum Particle in a Box. Generate Oscillations in a Circular Membrane. Study the Formation of a Shock Wave. Find the Value of a European Call Option. Construct a Complex Analytic Function. All rights reserved.

Enable JavaScript to interact with content and submit forms on Wolfram websites.In mathematics and physicsthe heat equation is a certain partial differential equation. Solutions of the heat equation are sometimes known as caloric functions. The theory of the heat equation was first developed by Joseph Fourier in for the purpose of modeling how a quantity such as heat diffuses through a given region.

As the prototypical parabolic partial differential equationthe heat equation is among the most widely studied topics in pure mathematicsand its analysis is regarded as fundamental to the broader field of partial differential equations. The heat equation can also be considered on Riemannian manifoldsleading to many geometric applications.

Certain solutions of the heat equation known as heat kernels provide subtle information about the region on which they are defined, as exemplified through their application to the Atiyahâ€”Singer index theorem. The heat equation, along with variants thereof, is also important in many fields of science and applied mathematics.

In probability theorythe heat equation is connected with the study of random walks and Brownian motion via the Fokkerâ€”Planck equation.

In image analysisthe heat equation is sometimes used to resolve pixelation and to identify edges. Following Robert Richtmyer and John von Neumann 's introduction of "artificial viscosity" methods, solutions of heat equations have been useful in the mathematical formulation of hydrodynamical shocks.

Solutions of the heat equation have also been given much attention in the numerical analysis literature, beginning in the s with work of Jim Douglas, D. Peaceman, and Henry Rachford Jr. It is typical to refer to t as "time" and x 1The collection of spatial variables is often referred to simply as x.

As such, the heat equation is often written more compactly as. In physics and engineering contexts, especially in the context of diffusion through a medium, it is more common to fix a Cartesian coordinate system and then to consider the specific case of a function u xyzt of three spatial variables xyz and time variable t. One then says that u is a solution of the heat equation if. In addition to other physical phenomena, this equation describes the flow of heat in a homogeneous and isotropic medium, with u xyzt being the temperature at the point xyz and time t.

In mathematical terms, one would say that the Laplacian is "translationally and rotationally invariant.

## Heat Conduction Equation

This can be taken as a significant and purely mathematical justification of the use of the Laplacian and of the heat equation in modeling any physical phenomena which are homogeneous and isotropic, of which heat diffusion is a principal example. This is not a major difference, for the following reason. Let u be a function with. Then, according to the chain ruleone has.We are going to give several forms of the heat equation for reference purposes, but we will only be really solving one of them.

Note that with this assumption the actual shape of the cross section i. Note that the 1-D assumption is actually not all that bad of an assumption as it might seem at first glance. If we assume that the lateral surface of the bar is perfectly insulated i. This means that heat can only flow from left to right or right to left and thus creating a 1-D temperature distribution. The assumption of the lateral surfaces being perfectly insulated is of course impossible, but it is possible to put enough insulation on the lateral surfaces that there will be very little heat flow through them and so, at least for a time, we can consider the lateral surfaces to be perfectly insulated.

As indicated we are going to assume, at least initially, that the specific heat may not be uniform throughout the bar. Note as well that in practice the specific heat depends upon the temperature.

Differential Equations: The Heat Equation

While this is a nice form of the heat equation it is not actually something we can solve. As noted the thermal conductivity can vary with the location in the bar. Also, much like the specific heat the thermal conductivity can vary with temperature, but we will assume that the total temperature change is not so great that this will be an issue and so we will assume for the purposes here that the thermal conductivity will not vary with temperature.

First, we know that if the temperature in a region is constant, i. Next, we know that if there is a temperature difference in a region we know the heat will flow from the hot portion to the cold portion of the region.

For example, if it is hotter to the right then we know that the heat should flow to the left. Finally, the greater the temperature difference in a region i. Note that we factored the minus sign out of the derivative to cancel against the minus sign that was already there.

In this case we generally say that the material in the bar is uniform.

Under these assumptions the heat equation becomes. There are four of them that are fairly common boundary conditions. The first type of boundary conditions that we can have would be the prescribed temperature boundary conditions, also called Dirichlet conditions. The prescribed temperature boundary conditions are.Documentation Help Center.

In a partial differential equation PDEthe function being solved for depends on several variables, and the differential equation can include partial derivatives taken with respect to each of the variables.

Partial differential equations are useful for modelling waves, heat flow, fluid dispersion, and other phenomena with spatial behavior that changes over time. You can think of these as ODEs of one variable that also change with respect to time.

Equations with a time derivative are parabolic. Equations without a time derivative are elliptic.

In other words, at least one equation in the system must include a time derivative. A 1-D PDE includes a function u xt that depends on time t and one spatial variable x. The spatial interval [ ab ] must be finite.

The diagonal elements of this matrix are either zero or positive. An element that is zero corresponds to an elliptic equation, and any other element corresponds to a parabolic equation. There must be at least one parabolic equation. An element of c that corresponds to a parabolic equation can vanish at isolated values of x if they are mesh points points where the solution is evaluated.

Discontinuities in c and s due to material interfaces are permitted provided that a mesh point is placed at each interface. To solve PDEs with pdepeyou must define the equation coefficients for cfand sthe initial conditions, the behavior of the solution at the boundaries, and a mesh of points to evaluate the solution on.

Together, the xmesh and tspan vectors form a 2-D grid that pdepe evaluates the solution on. You must express the PDEs in the standard form expected by pdepe. Written in this form, you can read off the values of the coefficients cfand s.

If there are multiple equations, then cfand s are vectors with each element corresponding to one equation. If there are multiple equations, then u0 is a vector with each element defining the initial condition of one equation. Note that the boundary conditions are expressed in terms of the flux frather than the partial derivative of u with respect to x. Also, of the two coefficients p xtu and q xtonly p can depend on u. In this case bcfun defines the boundary conditions.

If there are multiple equations, then the outputs pLqLpRand qR are vectors with each element defining the boundary condition of one equation. In some cases, you can improve solver performance by overriding these default values. To do this, use odeset to create an options structure. Then, pass the structure to pdepe as the last input argument:. Of the options for the underlying ODE solver ode15sonly those shown in the following table are available for pdepe.

InitialStepMaxStep. After you solve an equation with pdepeMATLAB returns the solution as a 3-D array solwhere sol i,j,k contains the k th component of the solution evaluated at t i and x j. The time mesh you specify is used purely for output purposes, and does not affect the internal time steps taken by the solver.

However, the spatial mesh you specify can affect the quality and speed of the solution. After solving an equation, you can use pdeval to evaluate the solution structure returned by pdepe with a different spatial mesh. The goal is to solve for the temperature u xt. The temperature is initially a nonzero constant, so the initial condition is.This website uses cookies to ensure you get the best experience. By using this website, you agree to our Cookie Policy.

Learn more Accept. Conic Sections Trigonometry. Conic Sections. Matrices Vectors. Chemical Reactions Chemical Properties.