🌀 nonlinear sciences

Exact self similar solutions for nonlinear coupled heat elastic wave systems with temperature dependent material properties

This paper employs Lie symmetry analysis to derive exact self-similar and traveling-wave solutions for a coupled nonlinear thermoelastic system with temperature-dependent material properties across various radial geometries, classifying the constitutive functions into exponential, power-law, and unified categories to reveal distinct thermal and displacement profiles.

Original authors: Manjit Singh, Radhika

Published 2026-07-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Manjit Singh, Radhika

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a solid object, like a metal rod or a rock, getting heated up. In the real world, when things get hot, they don't just get hot; they also expand, bend, and wiggle. This paper is about understanding the complex dance between heat and movement in these materials, but with a twist: the rules of the dance change depending on how hot it gets.

Here is a breakdown of what the authors did, using simple analogies:

1. The Problem: A Material That Changes Its Mind

Usually, scientists treat materials like a rigid robot: if you heat it, it expands by a fixed amount. But in reality, materials are more like a chameleon.

Heat Capacity: How much energy it takes to warm up.
Conductivity: How fast heat travels through it.
Stiffness: How hard it is to bend.

The authors realized that in extreme situations (like space rockets or deep geothermal rocks), these "chameleon" properties change drastically as the temperature rises. If you use the old, simple rules, your math breaks. This paper tries to solve the math for these "chameleon" materials.

2. The Method: Finding the "Secret Symmetry"

The equations describing this heat-and-wiggle dance are incredibly messy. They are like a tangled ball of yarn. To untangle them, the authors used a mathematical tool called Lie Symmetry Analysis.

Think of this like looking at a snowflake. Even though it's complex, you can rotate it or flip it, and it looks the same. The authors looked for these "rotations" and "flips" (symmetries) in the heat equations.

The Discovery: They found that the equations only have these secret symmetries if the material follows specific "laws" of how its properties change with heat.
The Three Laws: They identified three main families of materials that play nice with the math:
1. The Exponential Family: Properties that grow or shrink very fast (like a snowball rolling down a hill getting bigger).
2. The Power-Law Family: Properties that change in a steady, predictable rhythm (like a drumbeat).
3. The Unified Family: A master category that includes both of the above as special cases.

3. The Solution: Cracking the Code

Once they found these specific families, they could use the symmetries to simplify the messy equations. Imagine taking a 3D movie and flattening it into a 2D drawing without losing the story. They turned complex, multi-variable equations into simpler, single-variable equations.

They then solved these simpler equations to find Exact Solutions. These aren't approximations; they are perfect mathematical descriptions of how the heat and movement behave.

4. What the Solutions Look Like (The "Characters")

The paper describes two main types of solutions, which act like different characters in a story:

The Self-Similar Solutions (The "Zoom" Effect):
Imagine taking a photo of a ripple in a pond and zooming in or out. The shape of the ripple looks the same, just bigger or smaller.
- For the Exponential Family: The heat spreads out in a "logarithmic" pattern (a specific curve), while the physical movement of the material fades away as you get further from the source.
- For the Power-Law Family: The heat follows a smooth, power-law curve. The movement splits into three distinct parts:
  1. A steady stretch (the material is permanently bent).
  2. A growing wave (the wobble gets bigger over time).
  3. A coupling wave (a specific wiggle caused only by the heat pushing on the material). If there's no heat coupling, this third part disappears.
The Traveling-Wave Solutions (The "Surfing" Effect):
Imagine a wave moving down a beach. The shape stays the same, but it moves forward. The authors found exact formulas for these moving waves.
- The Cylinder Surprise: They found something unique about cylindrical shapes (like a pipe). In a cylinder, the math includes a "logarithmic" term (a specific type of curve) that doesn't appear in flat or spherical shapes. It's like the cylinder has a unique "signature" in the math that makes it slightly more complex to solve than a flat sheet or a ball.

5. Why This Matters (According to the Paper)

The authors aren't claiming this will cure diseases or build new engines right now. Instead, they are providing perfect benchmarks.

Think of it like a gold standard ruler. If a computer program tries to simulate a rocket engine or a geothermal drill, it needs to be tested against a known, perfect answer to see if it's working correctly.
These solutions act as that perfect ruler. They show exactly how heat and movement interact when the material's properties are changing, revealing the "dominant balance" (the main forces at play) in a way that computer simulations alone cannot easily show.

Summary

In short, this paper is a mathematical detective story. The authors looked for the hidden rules that allow complex heat-and-motion equations to be solved perfectly. They found three specific types of materials that follow these rules, solved the equations for them, and discovered that the shape of the object (flat, cylindrical, or spherical) changes the "music" of the solution, with cylinders having a unique, slightly more complex tune. These solutions serve as a perfect reference point for anyone trying to model these complex physical systems.

Technical Summary: Exact Self-Similar Solutions for Nonlinear Coupled Heat–Elastic Wave Systems

Problem Statement
This work addresses the mathematical modeling of nonlinear coupled thermoelastic systems where material properties—specifically heat capacity $C(u)$ , thermal conductivity $K(u)$ , and elastic modulus $E(u)$ —are dependent on the temperature field $u$ . The governing system consists of two coupled partial differential equations (PDEs) describing the thermal field $u(z, t)$ and the displacement field $w(z, t)$ in a generalized radial geometry parameterized by $\nu$ (encompassing planar, cylindrical, and spherical configurations). The authors note that while classical linear thermoelasticity assumes constant material parameters, many physically relevant scenarios (e.g., high-temperature aerospace applications, geothermal processes) require nonlinear models where these coefficients vary significantly with temperature. The primary challenge lies in the fact that for arbitrary constitutive functions, nontrivial similarity reductions are generally unavailable, necessitating a systematic group classification to identify specific functional forms that admit exact analytical solutions.

Methodology
The study employs Lie symmetry analysis to perform a complete group classification of the governing system. The methodology proceeds through the following steps:

Lie Symmetry Determination: The authors derive the determining equations for the infinitesimal generators of point symmetries acting on the space $(z, t, u, w)$ .
Constitutive Classification: By solving the overdetermined system of determining equations, the authors identify all admissible forms of the constitutive functions $C(u)$ , $K(u)$ , and $E(u)$ that lead to extensions of the principal Lie invariance algebra. This process branches into cases based on whether the elasticity is constant or temperature-dependent.
Optimal Systems: For each identified symmetry algebra, the authors construct optimal systems of one-dimensional subalgebras. This involves computing commutator tables and applying adjoint transformations to generate a minimal list of mutually inequivalent symmetry generators.
Similarity Reduction: Selected generators from the optimal systems are used to derive similarity variables. Substituting these transformations into the original PDEs reduces the system to coupled ordinary differential equations (ODEs).
Exact Solution Construction: The reduced ODE systems are solved using integrating factor methods, flux-variable substitutions, and specific ansätze (power-law and logarithmic) to obtain exact closed-form solutions.

Key Contributions and Results

Constitutive Classification: The analysis yields three principal classes of constitutive functions that admit extended symmetry algebras:
1. Exponential Class: $C(u) = C_0 e^{pu}$ , $K(u) = K_0 e^{qu}$ , $E(u) = E_0 e^{ru}$ , with specific constraints on the exponents $p, q, r$ .
2. Power-Law Class: $C(u) = C_0 u^p$ , $K(u) = K_0 u^{2p}$ , with constant elasticity $E(u) = E_0$ .
3. Unified Class: A general case where $E(u) = E_0$ and $K(u) = K_0 C(u)^2$ , which encompasses both the exponential and power-law branches as limiting cases depending on the form of $C(u)$ .
Symmetry Algebras and Optimal Systems: The paper explicitly determines the admitted Lie algebras for each class and constructs their optimal systems of subalgebras. The analysis distinguishes between cases where the geometry parameter $\nu \neq 0$ (cylindrical/spherical) and $\nu = 0$ (planar), noting that $\nu = 0$ admits additional translation symmetries.
Exact Solutions for the Exponential Class:
- For the scaling reduction, the thermal field exhibits a logarithmic similarity profile ( $u \sim \ln(z/t)$ ), while the displacement field decays algebraically ( $w \sim t^{3/2}/z$ ).
- The geometry parameter $\nu$ is absorbed into the leading-order transformation, meaning the functional form of the outer asymptotic solution remains identical for both cylindrical ( $\nu=1$ ) and spherical ( $\nu=2$ ) geometries.
- The solution is interpreted as a leading-order outer asymptotic solution, valid for $s \gg 1$ , requiring a separate inner expansion near the origin to ensure regularity.
Exact Solutions for the Power-Law Class:
- The thermal field follows a self-similar power-law profile $u \sim (z^2/t)^{1/p}$ .
- The displacement field decomposes into three distinct modes:
  1. A steady, time-independent elastic mode.
  2. A homogeneous growth mode reflecting the self-similar spread of the disturbance.
  3. A thermoelastically driven mode that grows linearly in time, proportional to the coupling constant $\beta$ . This term vanishes if $\beta = 0$ , highlighting its origin in thermoelastic coupling.
- Travelling-wave solutions are constructed, with the displacement field expressed via the Gauss hypergeometric function $_2F_1$ . Special cases include $p=0$ (constant conductivity), which yields elementary polynomial solutions, and $p \to \infty$ , which approaches an isothermal state with linear displacement.
Geometric Distinctions: The study identifies a structural distinction in the cylindrical case ( $\nu=1$ ). Unlike other geometries, the cylindrical Green's function introduces a logarithmic term ( $\ln z$ ) in the thermal profile when the flux constant is non-zero. This makes the $\nu=1$ solution structurally richer and distinct, often requiring quadrature rather than closed-form elementary expressions for the displacement field.

Significance and Scope
The paper claims to provide the first complete group classification of this specific nonlinear thermoelastic system with arbitrary temperature-dependent constitutive functions. The derived exact solutions serve as benchmark data for numerical solvers and offer transparent insight into the dominant physical balances between thermal diffusion, elastic stiffness, and thermoelastic coupling.

The authors explicitly limit the scope of their findings:

The similarity solutions presented are primarily leading-order outer asymptotic solutions valid in the far-field regime ( $s \gg 1$ ). A matched inner expansion for the source neighborhood is not constructed within this study.
The analysis assumes spatially and thermally uniform coupling constants ( $\beta, \gamma$ ) and density ( $\rho$ ), excluding functionally graded materials.
The stability of these solutions (whether they attract nearby trajectories) is not addressed, as it would require a separate spectral analysis.
Degenerate parameter cases (e.g., $p=0, q=0$ ) where adjoint weights vanish are excluded from the optimal system construction to maintain generality, though specific solutions for $p=0$ are provided as special cases.

In conclusion, the work establishes a systematic analytical framework for nonlinear thermoelasticity, demonstrating how specific constitutive laws enable the reduction of complex coupled PDEs to solvable ODEs, thereby revealing the interplay between material nonlinearity and geometric configuration.