一类变分波动方程经典解的爆破
Blow-Up of Classical Solution to a Class of Variational Wave Equations
摘要: 针对一类一维非线性变分波动方程初值问题,应用特征线法和能量估计方法证明了即使初始能量任意小,经典解也会在有限时间内发生爆破。该模型源于受电场影响的向列相液晶体的研究。
Abstract: We consider the initial value problem for a class of one-dimensional nonlinear variational wave equations. The method of characteristics and energy estimate is used to prove that the classical solutions will blow up in finite time even if the initial energy is arbitrarily small. This model is derived from the study of the nematic liquid crystals affected by electric fields.
文章引用:景晓萍. 一类变分波动方程经典解的爆破[J]. 理论数学, 2025, 15(3): 112-119. https://doi.org/10.12677/pm.2025.153082

1. 引言

液晶作为液体和晶体的中间状态而存在,它像液体一样流动,同时又具备晶体的某些有向结构特性。根据分子排列的形式和有序性的不同,液晶分为三种类型:向列相、近晶相和胆甾相。其中,向列相是最简单的液晶,它的棒状分子质心位置分布杂乱无章,但是分子整体取向一致,其长轴大致平行。在外力作用下很容易流动,被广泛应用于液晶显示器、生物芯片和柔性电子材料等领域。因此,对刻画向列型液晶分子运动的数学模型的研究具有重要理论意义和现实应用价值。

液晶材料的分子取向对磁场、电场和表面结构等微弱的外部刺激高度敏感,从而导致材料性质发生改变,随着分子局部方向的变化,材料中经常出现缺陷或者奇点。孙新涛等[1]指出,当向列相液晶受外加电场作用时,如果液晶材料含有杂质或带电粒子,或者外加电场很强,液晶材料会发生强烈的电光效应,此时讨论向列相液晶棒状分子方向向量是很困难的。实际应用中采用高纯度绝缘材料和较弱电场,此时液晶发生简单的展曲、扭曲或弯曲形变,棒状分子的方向向量由电场对液晶分子施加的转矩和形变转矩的平衡来确定。本文正是基于这个基本原理建立数学模型来讨论一定弱电场作用下向列相液晶的分子运动情况,从而建立一类变分波动方程经典解的爆破条件,为之后的实际应用和数值模拟做准备。

Hunter等[2]指出,向列液晶体的Frank-Oseen势能密度为

W( n,n )= 1 2 α ( n ) 2 + 1 2 β ( n×n ) 2 + 1 2 γ | n×( ×n ) | 2 , (1)

其中 α,β,γ 是液晶体的正弹性常数,分别代表展曲,扭曲和弯曲。当向列相液晶处于一定弱电场中时,引入电能密度为 W elec =λ g 2 ( n ) ,其中 λ 为介电常数, n 是单位向量,表示向列相液晶棒状分子的长轴取向(后面简称“方向 n ”), g( n ) 与方向 n 有关,见[3]。由最小作用原理

δ ( 1 2 t n t nW( n,n )λ g 2 ( n ) )dxdt =0,     nn=1 (2)

可知,当方向 n=( cosu( x,t ),sinu( x,t ),0 ) 时(这里,因变量 u 1 表示方向 n x轴的夹角, ( t,x ) 分别是时间自变量和空间自变量),Euler-Lagrange方程为

u tt c( u ) [ c( u ) u x ] x + λ 2 ( u+ u 3 )=0, (3)

其中波速 c 2 ( u )=α sin 2 u+γ cos 2 u 。这就是本文研究的一维非线性变分波动方程。

λ=0 时,方程(3)是最简单的一维变分波动方程

u tt c( u ) [ c( u ) u x ] x =0, (4)

[4]。假设波速c是严格正函数,Glassey等[5]最早使用特征线证明了方程(4)对于一般的光滑初值不存在整体光滑解;Duan等[6][5]中的结果拓展到多维。这种正则性的缺失迫使研究者放弃经典解而去考虑弱解。Zhang等[7] [8]和Bressan等[9]系统地研究了方程(4)耗散型弱解的整体存在性。Bressan等[10]-[12]用能量坐标法证明了方程(4)守恒型弱解的存在唯一性以及其他相关的性质。近年来,Song [13]利用[5]的方法证明了一维变分Sine-Gordon方程光滑解的爆破;Qin等[14]进一步证明了多维Sine-Gordon方程在特定初值下的光滑解在有限时间内爆破。最近,在不假设双曲方程的波速c是一个正函数的情况下,Hu [15]证明了一维非线性退化双曲抛物耦合系统对于任意小的初始能量,光滑解在有限时间内发生爆破。更多爆破结果见文献[16] [17]。本文将研究受电场影响的向列相液晶的变分波动方程并证明其初值问题经典解发生爆破的充要条件,为后续数值计算和弱解的研究奠定理论基础。

本文假设 c( u ) C 2 满足

0< c 0 c( u ) c 1 | c ( u ) | c 1 , (5)

这里 c 0 c 1 是正常数。主要结论如下:

定理1 考虑一维非线性变分波动方程初值问题

{ u tt c( u ) [ c( u ) u x ] x + λ 2 ( u+ u 3 )=0, u( 0,x )= u 0 +εφ( x x 0 ε ), u t ( 0,x )=( c( u( 0,x ) )+ε ) u x ( 0,x ), (6)

其中 u 0 , x 0 ,ε 均为正常数。假设 c ( u 0 ) 为正常数且光滑函数 φ( z ) 满足

φ( z ) C c 1 ( ( 1,1 ) ), max{ 3 ε , 64 c 1 2 c ( u 0 )( x 0 ε )ε }< φ ( 0 )< 2 c 0 λ x 0 ε . (7)

如果 u( t,x ) C 2 ( [ 0,T )× ) 是初值问题(6)的光滑解,那么对于任给的 ε>0 ,存在足够小的正数 ε 0 ,当 ε< ε 0 时, T<

本文使用[5] [6]的特征线法证明定理1,具体而言,为方便假设特征线,需引入新变量,把二阶方程转化为一阶方程组,将原经典解 u( t,x ) C 2 ( [ 0,T )× ) 在有限时间内爆破的问题转换为新变量沿特征线在一定时间内爆破和有界的问题。其中,证明难点是 c ( u ) 符号的控制。根据一阶方程组假设特征线得到特征三角区域,在该区域内对能量方程积分得到新变量的估计,从而控制 c ( u ) 沿某一条特征线恒正。

2. 定理1的证明

为了证明定理1,首先引入新变量,将二阶方程转化成一阶方程组,将原初值问题(6)转换为新变量的初值问题;其次,推导新变量满足的能量方程并建立能量估计;进一步,在特征三角区域对新变量进

行估计从而控制 c ( u ) 的符号;最后通过引理1证明初值问题(6)的光滑解在 t= x 0 ε c 1 之前发生爆破。

2.1. 变量变换

引入新变量

R 1 := u t +c( u ) u x ,   S 1 := u t c( u ) u x , (8)

则有

u x = R 1 S 1 2c( u ) u t = R 1 + S 1 2 , (9)

且初值问题(3) (6)等价于以下初值问题

{ t R 1 c( u ) x R 1 = c ( u ) 4c( u ) ( R 1 2 S 1 2 ) λ 2 ( u+ u 3 ), t S 1 +c( u ) x S 1 = c ( u ) 4c( u ) ( S 1 2 R 1 2 ) λ 2 ( u+ u 3 ), (10)

{ R 1 ( 0,x )=ε u x ( 0,x ), S 1 ( 0,x )=( 2c( u( 0,x ) )+ε ) u x ( 0,x ) (11)

2.2. 能量估计

考虑足够小的 ε=min{ c 0 , x 0 3 } ,由于 φ( ) C c 1 ( ( 1,1 ) ) ,当 x\[ x 0 ε, x 0 +ε ] 时,

R 1 ( 0,x )= S 1 ( 0,x )=0. (12)

将方程组(10)中第一个方程乘以 2 R 1 ,第二个方程乘以 2 S 1 ,然后相加,得到能量方程

t [ R 1 2 + S 1 2 +λ( u 2 + u 4 2 ) ] x [ c( u )( R 1 2 S 1 2 ) ]=0, (13)

定义能量函数

E( t ):= [ R 1 2 + S 1 2 +λ( u 2 + u 4 2 ) ] ( t,x )dx, (14)

则由(5) (6) (12)得初始能量

E( 0 )= + [ R 1 2 ( 0,x )+ S 1 2 ( 0,x )+λ( u 2 ( 0,x )+ u 4 2 ( 0,x ) ) ]dx = x 0 ε x 0 +ε { [ ε 2 + ( 2c( u( 0,x ) )+ε ) 2 ] φ 2 ( x x 0 ε ) +λ[ ( u 0 +εφ( x x 0 ε ) ) 2 + 1 2 ( u 0 +εφ( x x 0 ε ) ) 4 ] }dx ε 1 1 [ ( ε 2 + ( 2 c 1 +ε ) 2 )( φ 2 ( z ) ) ]dz +2λ M 1 ε { 2 φ ( z ) L 2 [ c 0 2 + ( 2 c 1 + c 0 ) 2 ]+2λ M 1 }ε:= M 2 ε, (15)

这里 M 1 = max z[ 1,1 ] [ ( u 0 +εφ( z ) ) 2 + 1 2 ( u 0 +εφ( z ) ) 4 ] M 2 为正常数。

易得 E( t ) 与时间t无关,则有

E( t )=E( 0 ) M 2 ε,

下面对方程(13)在特征三角区域积分。

x轴上取两点 ( x 1 ,0 ),( x 2 ,0 ) 且满足 0< x 1 < x 2 。过该两点定义两组特征线

{ d x + ( t ) dt =c( u( t, x + ( t ) ) ),   x + ( 0 )= x 1 , d x ( t ) dt =c( u( t, x ( t ) ) ),   x ( 0 )= x 2 . (16)

由(5)可知,如果 x 2 x 1 2 c 0 ( x 0 ε ) c 1 ,那么正特征线 x + ( t ) 和负特征线 x ( t ) 会相交于点 ( x A , t A ) 且有 x 1 < x A < x 2 t A < x 0 ε c 1 。设特征线 x + ( t ) x ( t ) x轴所围区域为 Ω ,如图1所示。

Figure 1. Characteristic plane

图1. 特征平面

在区域 Ω 上对方程(13)应用格林公式,可得:

Ω [ R 1 2 + S 1 2 +λ( u 2 + u 4 2 ) ] t [ c( u )( R 1 2 S 1 2 ) ] x ( t,x )dxdt = Ω [ c( u( t,x ) )( R 1 2 ( t,x ) S 1 2 ( t,x ) ) ]dt+[ R 1 2 ( t,x )+ S 1 2 ( t,x )+λ( u 2 + u 4 2 )( t,x ) ]dx = x 2 x A [ 2 S 1 2 +λ( u 2 + u 4 2 ) ] ( t ( x ),x )dx+ x A x 1 [ 2 R 1 2 +λ( u 2 + u 4 2 ) ] ( t + ( x ),x )dx + x 1 x 2 [ R 1 2 + S 1 2 +λ( u 2 + u 4 2 ) ] ( 0,x )dx=0.

则有

x A x 2 [ S 1 2 + λ 2 ( u 2 + u 4 2 ) ] ( t ( x ),x )dx+ x 1 x A [ R 1 2 + λ 2 ( u 2 + u 4 2 ) ] ( t + ( x ),x )dx = 1 2 x 1 x 2 [ R 1 2 + S 1 2 +λ( u 2 + u 4 2 ) ] ( 0,x )dx M 2 ε. (17)

x轴上取点 ( x 0 ,0 )( x 1 < x 0 < x 2 ) ,定义特征线 x ^ ( t )( t< x 0 ε c 1 )

d x ^ ( t ) dt =c( u( t, x ^ ( t ) ) ),  x ^ ( 0 )= x 0 , (18)

图1所示。下面证明当 ε 足够小时, c ( u( t, x ^ ( t ) ) )>0

由定义(8)可知

du( t, x ^ ( t ) ) dt = S 1 ( t, x ^ ( t ) ). (19)

对(19)从 t=0 t< x 0 ε c 1 积分并应用(17)得

| u( t, x ^ ( t ) )u( 0, x 0 ) | 0 t | S 1 ( t, x ^ ( t ) ) |dτ t ( 0 t S 1 2 ( t, x ^ ( t ) )dτ ) 1 2 x 0 ε c 0 c 1 ( x ^ ( t ) x 0 S 1 2 ( t ^ ( x ),x )dx ) 1 2 M 2 ( x 0 ε ) c 0 c 1 ε .

由于函数 c( ) C 2 ,因此,对于足够小的 δ>0 ,当 ε( 0,δ ) 时,有

c ( u( t, x ^ ( t ) ) ) c ( u( 0, x 0 ) ) 2 c ( u 0 ) 4 >0. (20)

2.3. 光滑解的爆破

引理1 在定理1条件下 ε 0 =min{ c 0 ,δ, x 0 3 , 1 3 M 2 M 3 , c ( u 0 ) x 0 48 c 1 2 M 2 M 3 } ,其中 M 3 = 1 c 0 + c 1 c 0 2 ,则当 ε< ε 0 时,初值问题(10) (11)的光滑解 R 1 ( t, x ^ ( t ) ) t< x 0 ε c 1 趋于正无穷,而光滑解 S 1 ( t, x ^ ( t ) ) ,在区间 [ 0, T ) ( T < x 0 ε c 1 ) 一致有界。

证明:由(7) (11)及引理1条件可知

max{ 3, 64 c 1 2 c ( u 0 )( x 0 ε ) }< R 1 ( 0, x 0 )< 2 c 0 λ x 0 R 2 ( 0, x 0 )>1,  S 2 ( 0, x 0 )<0.

第1步,证明 R 1 ( t, x ^ ( t ) ) t[ 0, x 0 ε c 1 ) 恒成立。

否则,由 R 1 ( 0, x 0 )>3 且函数 R 1 ( t, x ^ ( t ) ) t[ 0, x 0 ε c 1 ) 光滑,设 t 1 ( 0, x 0 ε c 1 ) 使得 R 1 ( t, x ^ ( t ) ) C 1 ( [ 0, t 1 ] ) 且满足 R 1 ( t, x ^ ( t ) )>1 t[ 0, t 1 ) R 1 ( t 1 , x ^ ( t 1 ) )=1 ,则由(5) (10) (20)有

d dt ( 1 R 1 ( t, x ^ ( t ) ) )= 1 R 1 2 [ c ( u ) 4c( u ) ( R 1 2 S 1 2 ) λ 2 ( u+ u 3 ) ] c ( u 0 ) 16 c 1 + 1 R 1 2 [ c 1 4 c 0 S 1 2 + λ 2 | u+ u 3 | ] c ( u 0 ) 16 c 1 + 1 R 1 2 [ c 1 4 c 0 S 1 2 + λ 2 | u+ u 3 | ]. (21)

对(21)从 t=0 t 1 积分并应用引理1条件得

1= 1 R 1 ( t 1 , x ^ ( t 1 ) ) 1 R 1 ( 0, x 0 ) 0 t 1 c ( u 0 ) 16 c 1 dt + 0 t 1 [ c 1 c 0 S 1 2 + λ 2 | u+ u 3 | ] ( t, x ^ ( t ) )dt 1 R 1 ( 0, x 0 ) + λ 2 c 0 x ^ ( t 1 ) x 0 1+3 u 2 + u 4 2 ( t ^ ( x ),x )dx + c 1 c 0 2 x ^ ( t 1 ) x 0 S 1 2 ( t ^ ( x ),x )dx 1 R 1 ( 0, x 0 ) + λ 2 c 0 c 1 t 1 + 1 2 c 0 x ^ ( t 1 ) x 0 2λ( u 2 + u 4 2 ) ( t ^ ( x ),x )dx+ c 1 c 0 2 M 2 ε < 1 R 1 ( 0, x 0 ) + λ x 0 2 c 0 +( 1 c 0 + c 1 c 0 2 ) M 2 ε< 2 R 1 ( 0, x 0 ) + M 2 M 3 ε 0 < 2 3 + M 2 M 3 1 3 M 2 M 3 =1,

这里 M 3 := 1 c 0 + 3 c 1 c 0 2 ,矛盾。因此 R 1 ( t, x ^ ( t ) )>1 t[ 0, x 0 ε c 1 ) 恒成立。

第2步,证明 R 1 ( t, x ^ ( t ) )+ t[ 0, x 0 ε c 1 ) 成立。

对(21)从 t=0 t< x 0 ε c 1 积分并应用引理1条件得

1 R 1 ( t, x ^ ( t ) ) 1 R 1 ( 0, x 0 ) c ( u 0 ) 16 c 1 t+ 0 t [ c 1 c 0 S 1 2 + λ 2 | u+ u 3 | ] ( t, x ^ ( t ) )dτ 1 R 1 ( 0, x 0 ) c ( u 0 ) 16 c 1 t+ λ 2 c 0 x ^ ( t ) x 0 1+3 u 2 + u 4 2 ( t ^ ( x ),x )dx + c 1 c 0 2 x ^ ( t ) x 0 S 1 2 ( t ^ ( x ),x )dx 1 R 1 ( 0, x 0 ) c ( u 0 ) 16 c 1 t+ λ x 0 2 c 0 +[ 1 c 0 + c 1 c 0 2 ] M 2 ε< 2 R 1 ( 0, x 0 ) c ( u 0 ) 16 c 1 t+ M 2 M 3 ε < c ( u 0 )( x 0 ε ) 32 c 1 2 c ( u 0 ) 16 c 1 t+ M 2 M 3 c ( u 0 ) x 0 48 c 1 2 M 2 M 3 = c ( u 0 ) 16 c 1 ( x 0 ε 2 c 1 t+ x 0 3 c 1 ).

x 0 ε 2 c 1 t+ x 0 3 c 1 =0 ,得 t= 5 x 0 3ε 6 c 1 < 5 x 0 3ε+ x 0 3ε 6 c 1 = x 0 ε c 1 。因此, T = 5 x 0 3ε 6 c 1 < x 0 ε c 1 ,使得当 ε< ε 0 t T R 1 ( t, x ^ ( t ) )+

第3步,证明 S 1 ( t, x ^ ( t ) ) [ 0, T ) 一致有界。

对于任意 t [ 0, T ) ,过曲线 x= x ^ ( t ) 上的点 ( t , x ) 作正特征线:

d x ( t ) dt =c( u( t, x ( t ) ) ),  x ( 0 )= x 0 . (22)

对光滑解 S 1 ( t,x ) 沿特征线 x= x ( t ) t=0 t= t 积分得

S 1 ( t , x ^ ( t ) )= S 1 ( 0, x 0 )+ 0 t [ c ( u ) 4c( u ) ( S 1 2 R 1 2 ) λ 2 ( u+ u 3 ) ] ( t, x ( t ) )dt,

| S 1 ( t , x ^ ( t ) ) || S 1 ( 0, x 0 ) |+ 0 t [ | c ( u ) | 4c( u ) ( S 1 2 + R 1 2 )+ λ 2 | u+ u 3 | ] ( t, x ( t ) )dt ( 2 c 1 +ε ) max z[ 1,1 ] | φ ( z ) |+ c 1 c 0 2 x 0 x ( S 1 2 + R 1 2 ) ( t ( x ),x )dx + 1 c 0 x 0 x λ( u 2 + u 4 2 ) ( t ( x ),x )dx+ λ x 0 2 c 0 ( 2 c 1 + ε 0 ) max z[ 1,1 ] | φ ( z ) |+ c 1 + c 0 c 0 2 M 2 ε 0 + λ x 0 2 c 0 .

( t , x ) 在曲线 x= x ^ ( t ) 的任意性知,当 ε< ε 0 时, S 1 ( t, x ^ ( t ) ) [ 0, T ) 一致有界。

综上,引理1证毕。

由引理1以及(5) (9)可知,对于任给的 ε>0 ,存在正数 T = 5 x 0 3ε 6 c 1 < x 0 ε c 1 和足够小的正数 ε 0 =min{ c 0 ,δ, x 0 3 , 1 3 M 2 M 3 , c ( u 0 ) x 0 48 c 1 2 M 2 M 3 } ,使得当 ε< ε 0 t T - 时,有 u+ 。即当初始能量 E( 0 ) 任意小时,初值问题(6)的经典解会在有限时间内发生爆破,定理1证毕。

致 谢

感谢审稿专家和编辑老师对本文提出宝贵的修改意见。

参考文献

[1] 孙新涛, 钱祥忠. 电场对向列相液晶指向矢取向的影响[J]. 淮南工业学院学报, 1999, 19(2): 71-74.
[2] Hunter, J.K. and Zheng, Y. (1995) On a Nonlinear Hyperbolic Variational Equation: I. Global Existence of Weak Solutions. Archive for Rational Mechanics and Analysis, 129, 305-353.
https://doi.org/10.1007/bf00379259
[3] Stewart, I.W. (2004) The Static and Dynamic Continuum Theory of Liquid Crystals: A Mathematical Introduction. CRC Press, 26-28.
[4] Alì, G. and Hunter, J.K. (2009) Orientation Waves in a Director Field with Rotational Inertia. Kinetic & Related Models, 2, 1-37.
https://doi.org/10.3934/krm.2009.2.1
[5] Glassey, R.T., Hunter, J.K. and Zheng, Y. (1996) Singularities of a Variational Wave Equation. Journal of Differential Equations, 129, 49-78.
https://doi.org/10.1006/jdeq.1996.0111
[6] Duan, W., Hu, Y. and Wang, G. (2020) Singularity and Existence for a Multidimensional Variational Wave Equation Arising from Nematic Liquid Crystals. Journal of Mathematical Analysis and Applications, 487, Article 124026.
https://doi.org/10.1016/j.jmaa.2020.124026
[7] Zhang, P. and Zheng, Y. (2003) Weak Solutions to a Nonlinear Variational Wave Equation. Archive for Rational Mechanics and Analysis, 166, 303-319.
https://doi.org/10.1007/s00205-002-0232-7
[8] Zhang, P. and Zheng, Y. (2005) Weak Solutions to a Nonlinear Variational Wave Equation with General Data. Annales de l’Institut Henri Poincaré C, Analyse non linéaire, 22, 207-226.
https://doi.org/10.1016/j.anihpc.2004.04.001
[9] Bressan, A. and Huang, T. (2016) Representation of Dissipative Solutions to a Nonlinear Variational Wave Equation. Communications in Mathematical Sciences, 14, 31-53.
https://doi.org/10.4310/cms.2016.v14.n1.a2
[10] Bressan, A. and Zheng, Y. (2006) Conservative Solutions to a Nonlinear Variational Wave Equation. Communications in Mathematical Physics, 266, 471-497.
https://doi.org/10.1007/s00220-006-0047-8
[11] Bressan, A. (2016) Uniqueness of Conservative Solutions for Nonlinear Wave Equations via Characteristics. Bulletin of the Brazilian Mathematical Society, New Series, 47, 157-169.
https://doi.org/10.1007/s00574-016-0129-y
[12] Bressan, A. and Chen, G. (2017) Lipschitz Metrics for a Class of Nonlinear Wave Equations. Archive for Rational Mechanics and Analysis, 226, 1303-1343.
https://doi.org/10.1007/s00205-017-1155-7
[13] Song, K. (2003) On Singularity of a Nonlinear Variational Sine-Gordon Equation. Journal of Differential Equations, 189, 183-198.
https://doi.org/10.1016/s0022-0396(02)00150-x
[14] Qin, F.M., Song, K. and Wang, Q. (2023) Singularity Formation for Nonlinear Variational Sine-Gordon Equation in a Multidimensional Space. Bulletin of the Korean Mathematical Society, 60, 1697-1704.
https://doi.org/10.4134/BKMS.b220859
[15] Hu, Y. (2022) Singularity for a Nonlinear Degenerate Hyperbolic-Parabolic Coupled System Arising from Nematic Liquid Crystals. Advances in Nonlinear Analysis, 12, Article 20220268.
https://doi.org/10.1515/anona-2022-0268
[16] Chen, G. and Zheng, Y. (2013) Singularity and Existence for a Wave System of Nematic Liquid Crystals. Journal of Mathematical Analysis and Applications, 398, 170-188.
https://doi.org/10.1016/j.jmaa.2012.08.048
[17] Chen, G. and Sofiani, M. (2023) Singularity Formation for the General Poiseuille Flow of Nematic Liquid Crystals. Communications on Applied Mathematics and Computation, 5, 1130-1147.
https://doi.org/10.1007/s42967-022-00190-5