I'm looking for the analytical solution for the 1D convection diffusion equation with a constant heat flux.(adsbygoogle = window.adsbygoogle || []).push({});

Boundary conditions:

The domain I'm looking at is x from 0 meters to 1 meter. The temperature at x=0 is T=0 degrees Celsius. At x=1, T=100 C.

I'm given the equation:

d/dx(rho*Cp*u*T)=d/dx(k dT/dx) + q'''

rho=density

Cp=specific heat

u=velocity in x-dir

T=temperature

k=thermal conductivity

q=generated heat (constant)

alpha=k/(rho*Cp)

P=(rho*u*L)/k (Peclet number)

To=0 C=273 K

T_L=100 C=373 K

I get it in the form:

u/alpha (dT/dx) - d^2T/dx^2 = q'''/k (equation 2)

I need to solve for T as a function of x (space).

I found the no generation solution to be

T(x)=(T_L-To)*(exp(P*x/L)-1)/(exp(P)-1)+To;

I tried plugging my equation 2 into a DE solver and I got

T(x)=[a*c1*e^(bx/a)]/(b) + cx/b + c2

c1 and c2 are constants

a=-1

b=u/alpha

c=q'''/k

I tried plugging in my boundary conditions to solve for c1 and c2 multiple times, but the solutions are not coming close to my numerical solutions at all.

I have also searched numerous sites trying to find this solution.

If anyone knows what T(x) is for this type of problem with constant heat flux, please respond to this post. Thank you.

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Convection diffusion equation 1D exact solution

**Physics Forums | Science Articles, Homework Help, Discussion**