Moving boundary problems of Stefan-type for a novel third order nonlinear evolution equation with temporal modulation are here shown to be amenable to exact Airy-type solution via a classical Ermakov equation with its admitted nonlinear superposition principle. Application of the latter together with a class of involutory transformations sets the original moving boundary problem in a wide class with temporal modulation. As an appendix, reciprocally associated exactly solvable moving boundary problems are derived.
论文ID : 2508.17421标题 : Moving boundary problems with Ermakov symmetry reduction: nonlinear superposition principle and reciprocal transformation applications作者 : Colin Rogers (University of New South Wales), Adriana C. Briozzo (Universidad Austral)分类 : math.AP (数学分析)发表期刊 : Open Communications in Nonlinear Mathematical Physics, Vol.5 (2025) pp 90–100接收时间 : 2025年9月28日论文链接 : https://arxiv.org/abs/2508.17421 本文研究了一类新颖的三阶非线性演化方程的Stefan型移动边界问题,该方程具有时间调制特性。通过经典Ermakov方程及其非线性叠加原理,证明了这类问题可以获得精确的Airy型解。结合对合变换类,将原始移动边界问题嵌入到具有时间调制的广泛问题类中。作为附录,还导出了互反关联的精确可解移动边界问题。
Ermakov方程的重要性 : Ermakov方程作为多分量Ermakov型系统的经典基础成员,在非线性物理和连续介质力学中有广泛应用,包括超弹性管的非线性弹性动力学、海洋涡旋脉冲演化、磁气体动力学等。移动边界问题的挑战 : Stefan型移动边界问题在物理学和工程学中普遍存在,如相变问题、热传导、多孔介质渗流等,但解析求解通常极其困难。现有方法局限性 :传统数值方法计算复杂度高 缺乏统一的解析框架处理时间调制的非线性演化方程 对于三阶非线性方程的移动边界问题研究相对较少 本文旨在通过Ermakov对称约化和非线性叠加原理,为一类新的三阶非线性演化方程建立精确的解析解框架,并探索互反变换在扩展解空间中的应用。
引入新的三阶非线性演化方程 : 提出了具有时间调制的三阶非线性演化方程 u t + u x x x + λ ( t + a ) μ u − 4 u x = 0 u_t + u_{xxx} + \lambda(t+a)^{\mu}u^{-4}u_x = 0 u t + u xxx + λ ( t + a ) μ u − 4 u x = 0 建立Ermakov对称约化 : 通过ansatz u = ( t + a ) m Ψ ( x ( t + a ) n ) u = (t+a)^m\Psi(\frac{x}{(t+a)^n}) u = ( t + a ) m Ψ ( ( t + a ) n x ) 实现对称约化,将原方程约化为经典Ermakov型方程获得Airy型精确解 : 利用非线性叠加原理,构造了基于Airy函数的精确解析解构建广义问题类 : 通过对合变换 T ∗ T^* T ∗ 将单一问题嵌入到具有时间调制的广泛问题类中导出互反关联问题 : 通过互反变换 R ∗ R^* R ∗ 建立了与原问题关联的新的可解移动边界问题类研究Stefan型移动边界问题:
输入 : 三阶非线性演化方程及其边界条件输出 : 精确解析解约束 : 移动边界 x = S ( t ) x = S(t) x = S ( t ) 和相应的边界条件从原始方程:
u t + u x x x + λ ( t + a ) μ u − 4 u x = 0 u_t + u_{xxx} + \lambda(t+a)^{\mu}u^{-4}u_x = 0 u t + u xxx + λ ( t + a ) μ u − 4 u x = 0
采用对称ansatz:
u = ( t + a ) m Ψ ( x ( t + a ) n ) u = (t+a)^m\Psi\left(\frac{x}{(t+a)^n}\right) u = ( t + a ) m Ψ ( ( t + a ) n x )
通过参数匹配 m = − 1 / 3 m = -1/3 m = − 1/3 , n = 1 / 3 n = 1/3 n = 1/3 , μ = − 2 \mu = -2 μ = − 2 ,得到约化方程:
Ψ ′ ′ ′ − 1 3 ( ξ Ψ ) ′ + λ Ψ − 4 Ψ ′ = 0 \Psi''' - \frac{1}{3}(\xi\Psi)' + \lambda\Psi^{-4}\Psi' = 0 Ψ ′′′ − 3 1 ( ξ Ψ ) ′ + λ Ψ − 4 Ψ ′ = 0
积分后得到:
Ψ ′ ′ − 1 3 ξ Ψ − λ 3 Ψ − 3 = ζ \Psi'' - \frac{1}{3}\xi\Psi - \frac{\lambda}{3}\Psi^{-3} = \zeta Ψ ′′ − 3 1 ξ Ψ − 3 λ Ψ − 3 = ζ
通过尺度变换 ξ = ϵ z \xi = \epsilon z ξ = ϵz (其中 ϵ 3 = − 3 2 \epsilon^3 = -\frac{3}{2} ϵ 3 = − 2 3 ),转换为标准Ermakov方程:
Ψ z z + z 2 Ψ = k ∗ Ψ − 3 \Psi_{zz} + \frac{z}{2}\Psi = k^*\Psi^{-3} Ψ zz + 2 z Ψ = k ∗ Ψ − 3
Ermakov方程的解具有非线性叠加形式:
Ψ = c 1 Ω 1 2 + 2 c 2 Ω 1 Ω 2 + c 3 Ω 2 2 \Psi = \sqrt{c_1\Omega_1^2 + 2c_2\Omega_1\Omega_2 + c_3\Omega_2^2} Ψ = c 1 Ω 1 2 + 2 c 2 Ω 1 Ω 2 + c 3 Ω 2 2
其中 Ω 1 , Ω 2 \Omega_1, \Omega_2 Ω 1 , Ω 2 是辅助线性方程的解:
Ω z z + z 2 Ω = 0 \Omega_{zz} + \frac{z}{2}\Omega = 0 Ω zz + 2 z Ω = 0
约束条件:c 1 c 3 − c 2 2 = k ∗ W 2 c_1c_3 - c_2^2 = \frac{k^*}{W^2} c 1 c 3 − c 2 2 = W 2 k ∗ (W W W 是Wronskian行列式)
最终解表示为:
Ψ = c 1 Ai 2 ( − 2 1 / 3 ξ / ϵ ) + 2 c 2 Ai ( − 2 1 / 3 ξ / ϵ ) Bi ( − 2 1 / 3 ξ / ϵ ) + c 3 Bi 2 ( − 2 1 / 3 ξ / ϵ ) \Psi = \sqrt{c_1\text{Ai}^2(-2^{1/3}\xi/\epsilon) + 2c_2\text{Ai}(-2^{1/3}\xi/\epsilon)\text{Bi}(-2^{1/3}\xi/\epsilon) + c_3\text{Bi}^2(-2^{1/3}\xi/\epsilon)} Ψ = c 1 Ai 2 ( − 2 1/3 ξ / ϵ ) + 2 c 2 Ai ( − 2 1/3 ξ / ϵ ) Bi ( − 2 1/3 ξ / ϵ ) + c 3 Bi 2 ( − 2 1/3 ξ / ϵ )
统一框架 : 首次将Ermakov对称约化应用于三阶非线性演化方程的移动边界问题参数确定方法 : 通过系统的参数匹配确定对称约化的关键参数对合变换扩展 : 利用对合性质 T ∗ ∗ = I T^{**} = I T ∗∗ = I 的变换类扩展问题族互反变换应用 : 构造互反变换生成新的可解问题类考虑移动边界 S ( t ) = γ ( t + a ) 1 / 3 S(t) = \gamma(t+a)^{1/3} S ( t ) = γ ( t + a ) 1/3 上的Stefan型边界条件:
边界条件I : u x x ( S ( t ) , t ) − λ 3 ( t + a ) − 2 u − 3 ( S ( t ) , t ) = L m S i ( t ) S ˙ ( t ) u_{xx}(S(t),t) - \frac{\lambda}{3}(t+a)^{-2}u^{-3}(S(t),t) = L_mS^i(t)\dot{S}(t) u xx ( S ( t ) , t ) − 3 λ ( t + a ) − 2 u − 3 ( S ( t ) , t ) = L m S i ( t ) S ˙ ( t ) 边界条件II : u ( S ( t ) , t ) = P m S j ( t ) u(S(t),t) = P_mS^j(t) u ( S ( t ) , t ) = P m S j ( t ) 边界条件III : u x x ( 0 , t ) − λ 3 ( t + a ) − 2 u − 3 ( 0 , t ) = H 0 ( t + a ) k u_{xx}(0,t) - \frac{\lambda}{3}(t+a)^{-2}u^{-3}(0,t) = H_0(t+a)^k u xx ( 0 , t ) − 3 λ ( t + a ) − 2 u − 3 ( 0 , t ) = H 0 ( t + a ) k 通过对称性要求确定:
i = − 1 i = -1 i = − 1 , j = − 1 j = -1 j = − 1 , k = − 1 k = -1 k = − 1 边界条件简化为代数约束 精确解构造 : 成功构造了基于Airy函数的精确解析解,包含可调参数 c 1 , c 2 , c 3 , k ∗ c_1, c_2, c_3, k^* c 1 , c 2 , c 3 , k ∗ 边界条件满足 : 验证了所构造的解满足所有Stefan型边界条件参数关系 : 建立了边界参数与Ermakov方程参数之间的明确关系通过变换:
d t ∗ = ρ − 2 ( t ) d t , x ∗ = x , u ∗ = ρ − 1 ( t ) u dt^* = \rho^{-2}(t)dt, \quad x^* = x, \quad u^* = \rho^{-1}(t)u d t ∗ = ρ − 2 ( t ) d t , x ∗ = x , u ∗ = ρ − 1 ( t ) u
成功将单一问题嵌入到广泛的时间调制问题族中。
通过互反变换:
d x ∗ = u d x + [ − u x x + λ 3 u − 3 ( t + a ) − 2 ] d t dx^* = u dx + \left[-u_{xx} + \frac{\lambda}{3}u^{-3}(t+a)^{-2}\right]dt d x ∗ = u d x + [ − u xx + 3 λ u − 3 ( t + a ) − 2 ] d t
导出了新的移动边界问题,其中移动边界具有对数形式:
S ∗ ( t ) = γ 3 [ Ψ ( γ ) − L m γ − 1 ] ln ( t + a ) S^*(t) = \frac{\gamma}{3}[\Psi(\gamma) - L_m\gamma^{-1}]\ln(t+a) S ∗ ( t ) = 3 γ [ Ψ ( γ ) − L m γ − 1 ] ln ( t + a )
Ermakov方程 : 由Ermakov在1880年引入,后来发展为多分量系统移动边界问题 : Stefan问题的现代发展,在相变理论中的应用互反变换 : 由Bateman引入,后来发展为Bäcklund变换的特殊情况孤子理论中的应用 海洋学中的应用 等离子体物理中的应用 多孔介质渗流问题 解析可解性 : 证明了一类新的三阶非线性演化方程的Stefan型移动边界问题具有精确的Airy型解系统性方法 : 建立了基于Ermakov对称约化的系统性求解框架问题族构造 : 通过变换群理论构造了广泛的可解问题族特定形式限制 : 方法仅适用于具有特定对称性的方程边界条件约束 : 需要特定形式的Stefan型边界条件参数关系 : 边界参数之间存在严格的代数约束扩展到高维 : 将方法推广到多维移动边界问题其他对称性 : 探索其他类型的对称约化数值验证 : 通过数值方法验证解析解的稳定性理论创新 : 首次将Ermakov理论系统应用于三阶非线性方程的移动边界问题方法完整 : 提供了从对称约化到精确解构造的完整框架数学严谨 : 所有变换和约化过程都有严格的数学证明应用潜力 : 为相关物理问题提供了新的解析工具适用范围 : 方法的适用性受到对称性要求的严格限制物理解释 : 缺乏对所构造解的物理意义的深入讨论数值比较 : 未提供与数值方法的比较验证学术价值 : 为移动边界问题理论提供了新的数学工具实用价值 : 可应用于相变、热传导等实际物理问题方法论贡献 : 展示了群论方法在非线性偏微分方程中的强大应用相变问题 : Stefan型相变过程的精确分析热传导 : 非线性热传导中的移动边界问题流体力学 : 自由边界流动问题材料科学 : 材料生长和相变过程建模论文引用了56篇重要文献,涵盖:
Ermakov方程理论的经典文献 孤子理论和可积系统 移动边界问题的数学理论 互反变换和Bäcklund变换 相关的物理应用 主要参考文献包括Ablowitz等人的孤子理论工作、Rogers在移动边界问题方面的系列研究,以及Ermakov方程在各个物理领域的应用文献。