Lie symmetry method for a nonlinear heat-diffusion equation

This paper applies the classical Lie symmetry method to a nonlinear heat-diffusion equation with variable coefficients to classify admitted symmetries, reduce the partial differential equation to ordinary differential equations, and construct invariant solutions for physically relevant cases such as Storm-type materials and power-law dependencies.

The Gist

Imagine you are trying to predict how heat spreads through a strange, magical material. In a normal metal, heat flows smoothly and predictably, like water running down a flat hill. But in this paper, the authors are studying a material where the rules of the game change depending on how hot it already is.

Think of it like a crowded dance floor. If the room is empty, people (heat) can move fast. But as the room gets crowded (hotter), the movement slows down, or maybe speeds up, depending on the specific "dance rules" of that material. Mathematically, this is a nonlinear heat-diffusion equation. It's a complex puzzle because the "friction" and "capacity" of the material change as the temperature changes.

The authors, Julieta, Ernesto, and Adriana, are using a powerful mathematical tool called the Lie Symmetry Method. Here is how they use it, explained simply:

1. The Detective's Magnifying Glass (Symmetry)

Imagine you are looking at a pattern on a piece of fabric. If you rotate the fabric 90 degrees and the pattern looks exactly the same, you say the fabric has "rotational symmetry."

In math, equations can have similar symmetries. The authors are asking: "If we stretch time, shrink space, or change the temperature in a specific way, does this heat equation still look the same?"

Finding these symmetries is like finding a secret "key" or a "shortcut" in the maze. If an equation has a symmetry, it means the problem is simpler than it looks. It's like realizing that a complicated 3D sculpture is actually just a 2D drawing rotated in space.

2. The Two Main Scenarios

The authors realized that the "dance rules" (the coefficients $C(u)$ and $K(u)$) determine how many symmetries exist. They split the problem into two main scenarios:

• Scenario A: The Changing Rules (Non-Constant Ratio)
  Here, the relationship between the material's heat capacity and its conductivity is constantly shifting.

  • The Result: They found that for most materials, there are 3 basic symmetries (like moving forward in time, moving sideways in space, or scaling the whole system).
  • The Special Cases: However, if the material follows very specific "recipes" (like the Storm-type materials or Power-law materials), the equation gains more symmetries (up to 5). It's like finding a hidden door in the maze that only opens if you wear a specific hat. When they find these special recipes, they can solve the equation much more easily.

• Scenario B: The Fixed Ratio (Constant Ratio)
  Here, the relationship between the two properties stays in a fixed proportion, like a perfectly balanced scale.

  • The Result: This is the "Golden Ticket." Because the rules are so balanced, the equation admits 6 symmetries. This is the maximum number of symmetries possible for this type of equation. It's like the maze suddenly turning into a straight hallway.

3. Folding the Map (Reduction)

Once they find these symmetries (the keys), they use them to "fold" the map.

• The original equation is a Partial Differential Equation (PDE). Think of this as a map with two dimensions: Space ($x$) and Time ($t$). It's a huge, complex landscape.
• By applying the symmetries, they can "fold" this 2D map into a 1D line (an Ordinary Differential Equation or ODE).
• The Analogy: Imagine you have a giant, wrinkled sheet of paper with a complex drawing on it. By folding the paper perfectly along the symmetry lines, you flatten it out until the drawing becomes a simple, straight line. Solving a line is much easier than solving a wrinkled sheet!

4. The Real-World Examples

The authors didn't just do abstract math; they tested their method on three real-world scenarios found in physics:

1. Power-Type Materials: Where heat capacity changes based on a power of the temperature (like $u^2$).
2. Storm-Type Materials: A specific condition named after a scientist (Storm) where the material properties change in a very specific, exponential way. This is common in certain metals.
3. Power-Law Mixtures: Where both properties change together in a specific polynomial way.

For each of these, they used their "folding" technique to find Exact Solutions. These are precise formulas that tell you exactly what the temperature will be at any point in space and time, without needing a computer to guess.

Why Does This Matter?

In the real world, we often use computers to simulate how heat moves (like in a car engine or a building). But computers can make mistakes, and they are slow.

This paper provides exact formulas (the "benchmarks").

• If you are an engineer designing a heat shield, you can use these formulas to check if your computer simulation is correct.
• If you are studying how ice melts (a "phase-change problem"), these formulas help you understand the exact behavior of the melting front.

Summary

Think of this paper as a master keymaker.

1. They looked at a locked door (a complex heat equation).
2. They figured out that the lock only opens if the key (the material's properties) has a specific shape.
3. They identified those shapes (Storm-type, Power-law).
4. Once the door is open, they showed us exactly how to walk through it to find the treasure (the exact solution).

They turned a messy, impossible-looking problem into a clean, solvable one by finding the hidden patterns (symmetries) inside the chaos.

Technical Summary

Here is a detailed technical summary of the paper "Lie symmetry method for a nonlinear heat-diffusion equation" by Bollati, Rodriguez, and Briozzo.

1. Problem Statement

The paper investigates the nonlinear heat-diffusion equation given by:
$$C(u)u_t = (K(u)u_x)_x$$
where $u(x,t)$ represents the dependent variable (e.g., temperature or concentration), and $C(u)$ and $K(u)$ are non-constant coefficients representing heat capacity and thermal conductivity (or diffusivity), respectively.

The core challenge addressed is the difficulty of finding exact analytical solutions for nonlinear partial differential equations (PDEs) with variable, solution-dependent coefficients. Standard linear methods fail due to the essential nonlinearity introduced by the functional dependence of $C(u)$ and $K(u)$ on $u$. The authors aim to identify specific functional relationships between $C(u)$ and $K(u)$ that allow the equation to admit Lie point symmetries, thereby enabling the reduction of the PDE to ordinary differential equations (ODEs) and the construction of exact invariant solutions.

2. Methodology

The authors employ the classical Lie symmetry method, a systematic group-theoretic approach developed by Sophus Lie. The methodology involves the following steps:

1. Infinitesimal Generator Formulation: They define a one-parameter Lie group of transformations acting on the space $(x, t, u)$ with infinitesimal generators $\xi_1(x,t,u)$, $\xi_2(x,t,u)$, and $\eta(x,t,u)$.
2. Invariance Condition: They apply the second prolongation of the infinitesimal generator, $X^{(2)}$, to the PDE. The condition $X^{(2)}(F) = 0$ (where $F$ is the PDE rewritten as $C(u)u_t - K'(u)u_x^2 - K(u)u_{xx} = 0$) yields a system of determining equations.
3. Classification of Cases: Solving the determining equations reveals that the structure of the admitted symmetries depends critically on the ratio $C(u)/K(u)$. The analysis is divided into two distinct cases:
   • Case 1: The ratio $C(u)/K(u)$ is non-constant.
   • Case 2: The ratio $C(u)/K(u)$ is a constant ($\alpha$).
4. Symmetry Reduction: For each admitted symmetry generator, the authors solve the characteristic equations to find similarity variables and invariant solutions, reducing the original PDE to ODEs or integral equations.
5. Application to Specific Models: The theoretical results are applied to three specific physical models found in the literature.

3. Key Contributions and Results

A. Classification of Symmetries (General Case)

The paper provides a complete classification of the Lie point symmetries based on the functional forms of $C(u)$ and $K(u)$:

• Case 1 ($C/K$ is non-constant):

  • Generic Case: For arbitrary non-constant $C/K$, the equation admits a 3-parameter Lie group generated by scaling and translation ($X_1, X_2, X_3$).
  • Specific Functional Dependencies:
    • If $C(u)$ and $K(u)$ satisfy a specific relation involving an integral of $K(u)$ (Eq. 38 or 40), a 4-parameter group is admitted (adding generator $X_4$).
    • If a more restrictive condition (Eq. 44) is met, a 5-parameter group is admitted (adding generator $X_5$).
  • The authors derive the explicit forms of $C(u)$ required for these higher-dimensional symmetries.

• Case 2 ($C/K = \alpha$ is constant):

  • This case corresponds to a generalized linear diffusion structure where coefficients vary proportionally.
  • The equation admits a 6-parameter Lie group.
  • The generators include complex transformations involving $x^2$, $t^2$, and mixed terms, leading to a rich algebraic structure (detailed in Table 2).

B. Construction of Invariant Solutions

For each generator identified, the authors derive the corresponding invariant solutions:

• Similarity Solutions: Solutions of the form $u = \phi(\xi)$ where $\xi$ is a similarity variable (e.g., $\xi = x/\sqrt{t}$ or $\xi = x/t$).
• Reduction to ODEs: The PDE is reduced to second-order ODEs or first-order ODEs. In many cases, the solutions are expressed as integral equations (e.g., Eq. 72, 108) or implicit relations involving integrals of $K(u)$.
• Explicit Solutions: For specific choices of $K(u)$ (e.g., exponential or power laws), the implicit relations are solved to yield explicit analytical expressions for $u(x,t)$.

C. Analysis of Particular Physical Cases

The authors apply their general framework to three specific models:

1. Power-type $C(u)$ with constant $K(u)$:

   • Model: $u_t/u^2 = k u_{xx}$.
   • Result: Admits a 4-parameter group. Explicit solutions are derived, including linear profiles and specific similarity solutions relevant to Stefan problems (phase-change).

2. Storm's Condition:

   • Model: $d/du \sqrt{C(u)/K(u)} = \lambda C(u)$.
   • Context: This condition is crucial for modeling heat conduction in metals where thermal properties vary exponentially with temperature.
   • Result: The authors show this condition corresponds to a specific case of their general classification ($B=-1/2$). They derive explicit solutions involving logarithmic functions for exponential $K(u)$ and $C(u)$.

3. Power-type $C(u)$ and $K(u)$:

   • Model: $C(u) = \rho c_0(1+\beta u^p)$, $K(u) = k_0(1+\beta u^p)$.
   • Context: Relevant for free boundary problems.
   • Result: Since $C/K$ is constant, this falls under Case 2. The authors derive solutions involving the error function (erf) for the linear case ($p=1$) and integral forms for general $p$.

4. Significance and Impact

• Systematic Framework: The paper provides a rigorous, unified classification of nonlinear heat-diffusion equations based on the ratio of their coefficients. It clarifies the conditions under which these complex equations become analytically tractable.
• Extension of Known Results: The work generalizes previous findings (e.g., those by Cherniha and others) by explicitly deriving the functional forms of $C(u)$ and $K(u)$ required for higher-dimensional symmetry groups.
• Physical Relevance: By analyzing Storm's condition and power-law dependencies, the paper offers exact solutions for real-world problems in heat transfer, phase-change materials, and porous media flow, where material properties are not constant.
• Benchmarking: The derived exact solutions serve as valuable benchmarks for validating numerical simulations of nonlinear diffusion processes, which are often difficult to solve analytically.

Conclusion

The authors successfully demonstrate that the Lie symmetry method is a powerful tool for analyzing nonlinear heat-diffusion equations. By categorizing the equations based on the relationship between $C(u)$ and $K(u)$, they identified specific classes of materials and conditions that admit rich symmetry structures, allowing for the derivation of exact similarity solutions that would otherwise be inaccessible. This work bridges the gap between abstract group theory and practical applications in thermal physics and transport phenomena.