I'm interested in solving a heat equation with a discontinuous source term in one dimension on the real line:
where is the Heaviside function. The boundary conditions are that The initial condition is arbitrary, but even for a symmetric Gaussian or a symmetric exponential an analytical solution would be interesting.
Because of the discontinuity in the right-hand side induced by , my intuition would be to solve in the two domains x>0 and x<0 and patch them together at x=0. However, the source terms have the effect of causing exponential decay with rate x at the point x for x>0 and corresponding growth at the point −x. This introduces an asymmetry, which should cause a diffusive flux back to the side x>0. Solving for the solution in the two domains does not seem to capture this feedback, and I'm not sure how to solve globally. What are some avenues to proceed?