Applied Partial Differential Equations Solutions

Working Through PDEs by Hand Is Still Worth Doing

I spent years running finite element models in commercial software, but honestly the things that broke most often were the ones nobody thought to check first. Separation of variables still comes up in real engineering work, even if it feels like something you only see in textbooks. The trick is knowing when it applies and when it quietly stops working. I once had a student who got completely stuck on a heat equation problem because the boundary conditions were non-homogeneous on two sides instead of just one, and they kept trying to force separation anyway. That is not a trick question. That is just a normal thing that happens. The core idea is straightforward enough. You take a PDE, assume the solution can be written as a product of functions each depending on only one variable, and split the problem into ordinary differential equations. It works beautifully for constant coefficient linear equations on domains with simple geometry. Rectangle. Circle. Cylinder. Sphere. Those are the shapes where everything lines up. Once you cross into irregular boundaries or nonlinear terms, the whole approach starts falling apart and you need something else. Take the one-dimensional heat equation, which is probably the first example you encounter. u sub t equals alpha times u sub xx. You assume u of x, t equals X of x times T of t, plug it in, divide through by XT, and suddenly you have one side depending only on x and the other only on t. Both have to equal a constant. That constant becomes your eigenvalue, and the sign of that constant determines whether you get exponential decay, oscillatory behavior, or both. The boundary conditions decide which eigenvalues are actually allowed. Dirichlet conditions give you a discrete spectrum. Neumann conditions change the allowed modes. Mixed conditions are where people usually make mistakes because they forget that zero is sometimes an eigenvalue.

I ran into a case a few years back involving a composite wall problem where two different materials met at an interface inside the domain. Standard separation of variables gives you eigenfunctions for each region, but the coupling conditions at the interface turned the whole thing into a transcendental eigenvalue problem. There was no closed form for the eigenvalues. I ended up writing a small Newton iteration to find them numerically and then just used those values in the series solution. It took about ten minutes to code and saved me from trying to force an analytic approach that would never converge. Fourier series come next because they are how you actually satisfy arbitrary initial conditions. Once you have your eigenfunctions, you expand the initial temperature distribution or displacement profile in terms of them. The coefficients come from orthogonality relations. That part is mechanical. What people miss is that convergence can be really slow if your initial condition has discontinuities or sharp gradients. Gibbs phenomenon is not a theoretical curiosity. It shows up in real calculations and makes truncated series look wrong even when they are technically correct. I usually recommend checking convergence by computing the first twenty coefficients and seeing whether they are decaying like one over n or one over n squared. If they decay slowly, your solution will need many more terms than you expect. Transform methods are another tool worth knowing about, especially when your domain is infinite or semi-infinite. Laplace transforms handle time domains well. Fourier transforms work for spatial domains that extend to infinity. The trick is matching the transform to whichever variable has the simpler boundary behavior. I once solved a wave equation on a semi-infinite string using a Laplace transform in time and found that the reflection at the boundary appeared naturally as a time-shifted term in the inverse transform. No need to invent image methods or second-guess yourself. The transform does the reflection algebra for you.

Green functions deserve more practical attention than they usually get. They are essentially the impulse response of your differential operator with specified boundary conditions. Once you have the Green function, any source term and any boundary data can be folded into an integral representation of the solution. The downside is that Green functions are rarely available in closed form outside the simplest geometries. For a rectangular domain with Dirichlet conditions you can build one from eigenfunction expansions. For more complex setups you typically compute it numerically and then convolve. This approach is also sensitive to boundary condition type. Mixing Dirichlet and Neumann on different parts of the boundary requires a Green function that satisfies the correct condition on each part, and that constraint alone can make construction very tedious. One thing that almost nobody warns beginners about is what happens when your coefficients are not constant. Separation of variables breaks down immediately. You might think of a change of variables to fix it, but that only works for a narrow set of coefficient structures. When I encountered a diffusion equation with position-dependent diffusivity on a graduation project, the standard approach gave nothing useful. We ended up using a shooting method combined with a fourth-order Runge-Kutta integrator for the radial ODE that resulted from assuming a power-law diffusivity. The eigenvalues had to be found by scanning for zeros of the solution at the boundary. It was slower than a clean analytic solution but it actually converged and matched experimental data within five percent. Numerical methods fill the gap where analytic approaches cannot reach. Finite difference schemes are the easiest to implement but they require structured grids, which means they struggle with curved boundaries. Finite element methods handle complex geometry much better but the assembly process is computationally heavier and the resulting matrix systems can become ill-conditioned if your mesh quality drops. I found that for most structural vibration problems on irregular domains, a moderate quality triangular mesh with linear elements gave acceptable results with a setup time of roughly fifteen minutes per model, while a higher order quadrilateral mesh took about twice as long to set up but cut the solution error in half. The tradeoff is real and it depends entirely on how much time you have versus how much precision you need.

Boundary layer problems represent another category where analytic intuition fails. When a small parameter multiplies the highest derivative, like in convection-diffusion equations at high Peclet numbers, the solution develops thin regions near boundaries where gradients become enormous. Standard numerical methods either smear these out or oscillate wildly unless the mesh is refined specifically in those layers. I learned this the hard way when simulating heat transfer in a duct with a high Reynolds number. My uniform grid produced a physically impossible temperature profile near the wall. Switching to an adaptive grid that concentrated points near the boundary reduced the error from forty percent to under three percent and took about eight minutes to re-run. Nonlinear PDEs do not respond to any of the linear techniques described here. Method of characteristics still works for first-order equations, but second-order nonlinear problems usually require perturbation methods or complete numerical treatment. The Burgers equation is one of the few nonlinear PDEs that can be linearized through the Cole-Hopf transformation, and knowing that trick can save you a lot of headaches on exam problems. Most nonlinear PDEs in practice need a numerical solver, and even then you should validate your results against known solutions whenever possible. I once had a colleague who trusted a finite volume code for a nonlinear reaction-diffusion system without checking against an analytical limit case. The code produced smooth results that looked reasonable until we realized the reaction term had been implemented with the wrong sign. Smooth wrong results are worse than oscillatory wrong results because they do not trigger any alarm bells. If you need reference materials, the classic texts by Haberman and Strauss cover the standard separation of variables and Fourier method material very well for undergraduate level work. For graduate level treatment with more emphasis on modern techniques, Haberman again or Evans if you want a rigorous PDE foundation. Numerical approaches are well covered in LeVeque for finite differences and finite volumes and in Verfurth for finite elements. Online lecture notes from MIT OpenCourseWare and Stanford are freely available and often include worked examples that match homework problems closely.

When to Use Which Approach

Constant coefficient linear PDE on a rectangle with homogeneous boundary conditions. Use separation of variables and Fourier series. This is the bread and butter case and it should take you under thirty minutes from start to finish if you know your eigenfunction expansions. Constant coefficient linear PDE on an infinite or semi-infinite domain. Use transform methods. Laplace for time, Fourier for space. This avoids the eigenfunction matching hassle entirely and handles discontinuous initial data more gracefully. Variable coefficients or irregular geometry. Numerical methods are your only realistic option. Choose finite elements for complex boundaries and finite differences for simple structured domains. Validate against an analytic solution in a limiting case before trusting the numbers.

Nonlinear PDE. Check whether a transformation exists first. If not, use numerical methods and validate aggressively. The space of valid numerical solutions for nonlinear problems is much narrower than for linear ones, and convergence does not guarantee correctness. The biggest mistake I see students make is treating every PDE as if it belongs in the same toolbox. The method you pick should be driven by the structure of the equation and the geometry of the domain, not by what you practiced most recently. Applied Partial Differential Equations Solutions are not a single method. They are a collection of techniques that overlap in coverage and diverge sharply at the edges. Knowing where each one breaks down matters more than knowing how to apply it when everything goes smoothly.