Applying Fourier cosine transfer method with this Neumann boundary condition $\frac{dh(0,t)}{dx}=\frac{1}{a}(g(0,t)-h(x,t)) $ to solve heat equation!

398 Views Asked by At

Consider heat equation $$\frac{\partial h}{\partial t}= \text{D}\frac{\partial^{2} h}{\partial x^{2}}$$ I want to solve heat\diffusion equation with this Neumann boundary condition $$\frac{dh(0,t)}{dx}=\frac{1}{a}{g(0,t)-h(x,t)}$$ Fourier cosine transfer for second derivative of x is given as $$F_{c}(\frac{d^{2}h(x,t)}{dx^{2}})= U_{c}(w,t)= -w^{2}U_{c}(w,t)-\sqrt{\frac{2}{\pi}}\frac{dh(0,t)}{dx} $$ Consider Fc to be fourier cosine operator which converts h(x,t) into Uc(w,t). Fourier cosine transfer of time derivative is given as $$F_{c}(\frac{dh}{dt})=\frac{d}{dt}U_{c}{(w,t)}$$ The problem is, by solving heat equation using Fourier cosine transform with above mentioned neumann boundary condition we will end up having h(x,t) at both sides of resulting equation. $$\frac{\partial}{\partial t}U_{c}{(w,t)}=\text{D}(-w^{2}U_{c}(w,t)-\sqrt{\frac{2}{\pi}}\frac{\partial h(0,t)}{\partial x})$$ by substituting values of above mentioned boundary condition we get following ODE $$\frac{d}{dt}U_{c}{(w,t)}=\text{D}(-w^{2}U_{c}(w,t)-\sqrt{\frac{2}{\pi}}\frac{1}{a}{g(0,t)-h(x,t)})$$ Because Uc is a transformed form of h(x,t) which we just assumed and it will be transformed back when we apply inverse fourier transform How can I solve this ODE so that h(x,t) can be later taken to LHS. We have initial condition $$h(0,t)=0 $$