1 引言
本文研究定常非齐次不可压 Euler 流通过三维柱对称无穷管道的解的适定性和无穷远渐近速率, 管道中可包含障碍物. 三维定常非齐次不可压 Euler 方程为
(1.1) $\begin{equation} \left\{\begin{array}{ll} (\rho u_{1})_{x_{1}}+(\rho u_{2})_{x_{2}}+ (\rho u_{3})_{x_{3}}=0,\\ (\rho u_{1}u_{1})_{x_{1}}+(\rho u_{1}u_{2})_{x_{2}}+ (\rho u_{1}u_{3})_{x_{3}}+p_{x_{1}}=0,\\ (\rho u_{1}u_{2})_{x_{1}}+(\rho u_{2}u_{2})_{x_{2}}+ (\rho u_{2}u_{3})_{x_{3}}+p_{x_{2}}=0,\\ (\rho u_{1}u_{3})_{x_{1}}+(\rho u_{2}u_{3})_{x_{2}}+ (\rho u_{3}u_{3})_{x_{3}}+p_{x_{3}}=0,\\ (u_{1})_{x_{1}}+( u_{2})_{x_{2}}+ ( u_{3})_{x_{3}}=0, \end{array}\right.\label{1} \end{equation}$
其中$\rho$ , $(u_{1}, u_{2},u_{3})$ 和 $p$ 分别表示密度, 速度和压力.
如图1 所示, 流体区域是一个内部包含障碍物 $\Omega'$ 的无穷管道区域 $\widetilde \Omega$ .
图 1
$\widetilde \Omega = \Big\{(x_{1},x_{2},x_{3})\in\Bbb R^{3}|\sqrt{x_{2}^{2}+x_{3}^{2}}\in[0,f_1(x_{1})), \ x_{1}\in (-\infty,+\infty)\Big\} $
$ \Omega'=\Big\{(x_{1},x_{2},x_{3})\in\Bbb R^{3}|\sqrt{x_{2}^{2}+x_{3}^{2}}\in[L_{2}]\Big\}, $
其中 $L_{1}$ 和 $L_{2}$ 为常数. 令
${f_2}(x) = \left\{ \begin{array}{ll} {\tilde f_2}(x),\quad& x \in \left[ {L_1,L_2} \right],\\ 0,\quad& x \in \Bbb R \setminus [L_{1},L_{2}]. \end{array} \right.$
函数 $f_1(x_1)$ 和 $f_2(x_1)$ 满足
(1.2) $\begin{equation}\label{1.9} 0\leq f_2(x_1)\leq C\quad \mbox{和}\quad\frac{1}{C}\leq \big(f_1-f_2\big)(x_1)\leq C,\quad x_1\in \Bbb R, \quad\mbox{其中 $C$ 是常数}, \end{equation}$
(1.3) $\begin{equation}\label{1.11} \lim\limits_{x_{1}\rightarrow-\infty}f_1(x_{1})=1,\ \lim\limits_{x_{1}\rightarrow+\infty}f_1(x_{1})=a,\ \quad \mbox{其中 }\ a>0, \end{equation}$
(1.4) $\begin{equation}\label{1.10} \|f_1\|_{C^{2,\alpha}(\Bbb R\setminus \Pi)}\leq C \quad \mbox{和}\quad\|f_2\|_{C^{2,\alpha}(\Bbb R\setminus \Pi)}\leq C, \quad\mbox{其中}\ \alpha \in\left(0,1\right), \end{equation}$
这里$\Pi=\{P_{ i}\}_{i=1}^m$ 是边界上所有角点的集合, 其中 $\theta_i$ 和 $\beta_i$ 都是严格小于 $\pi$ . 记 $\Omega=\tilde\Omega\setminus\Omega'$ . 流体在边界满足条件
(1.5) $\begin{equation} (u_{1},u_{2},u_{3})\cdot\overrightarrow{n}=0, \quad(x_1,x_2,x_3)\in \partial \Omega\backslash\Pi, \label{1.12} \end{equation}$
其中 $\overrightarrow{n}$ 是 $\partial \Omega\backslash\Pi$ 上的单位外法向量. 在本文中的结论都包含 $f_0\equiv0$ , 即没有障碍物的情况.
定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程
(1.6) $\begin{equation} \Delta\psi=f(\psi), \end{equation}$
这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法.
本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性.
定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性.
在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件.
对于三维柱对称定常非齐次不可压 Euler 方程管道问题解的适定性和无穷远渐近速率的研究主要有以下三个困难: 第一, 方程的解可能在对称轴附近出现奇性, 对应于流线在对称轴附近可能会分叉, 特别是在管道中放置障碍物的情况; 第二, 当涡旋速度 $W\neq 0$ 时, 漩涡速度会产生离心力, 可能导致流场出现静止或回流; 第三, 边界的角点会影响解的正则性. 针对困难一, 首先我们在对称轴附近引入比较函数, 证明方程在对称轴附近不会出现奇性; 再证明水平速度 $U>0$ , 其中: 先证明无旋流的流线非退化性, 再通过估计有旋流和无旋流的差由连续性方法得到有旋流的流线不回流, 最后结合正负无穷远处的渐近行为和极大值原理得到流线的非退化性. 针对困难二, 本文引入了对于涡旋速度 $W$ 的流线守恒量 $rW$ . 针对困难三, 我们在角点处构造了合适的闸函数, 利用极大值原理得到一致估计.
本文主要的创新点如下: 第一, 对于含有障碍物的三维柱对称管道问题, 首次证明了涡旋速度 $W \neq 0$ 的非齐次不可压 Euler 流解的适定性; 第二, 在两种边界条件下, 得到了无穷远的渐近收敛速率: 若无穷管道在有限长度以外是平边界, 则方程的解以指数速率收敛到渐近状态; 若无穷管道以多项式速率收敛到平边界, 则方程的解以相同的多项式速率收敛到渐近状态.
本文的章节安排如下: 第 2 节将三维柱对称齐次不可压 Euler 方程的无穷长管道问题转化为一个无界区域中非线性椭圆方程的边值问题并介绍本文的主要结果; 第 3 节证明解的适定性; 第 4 节讨论解在无穷远处的渐近收敛速率.
2 问题设定与主要定理
2.1 柱对称变换
$ x=x_{1},\qquad r=\sqrt{x_{2}^{2}+x_{3}^{2}}, \qquad x_{2}=r\sin{\theta},\qquad x_{3}=r\cos{\theta}. $
$\rho(x,r,\theta)=\rho(x,r),\ (u_{1},u_{2},u_{3})(x,r,\theta)=(u_{1},u_{2},u_{3})(x,r),\ p(x,r,\theta)=p(x,r).$
(2.1) $\begin{equation} \label{symmcondtion} \lim_{(x_2^2+x_3^2)\rightarrow0}\frac{x_2}{\sqrt{x_2^2+x_3^2}}u_2(x)+\frac{x_3}{\sqrt{x_2^2+x_3^2}}u_3(x)=0. \end{equation}$
$ u_{1}=U(x,r),\ u_{2}=V(x,r)\frac{x_{2}}{r}+W(x,r)\frac{x_{3}}{r},\ u_{3}=V(x,r)\frac{x_{3}}{r}-W(x,r)\frac{x_{2}}{r}, $
其中 $U$ , $V$ 和 $W$ 分别是轴向速度, 径向速度和漩涡速度. 直接计算得到
$u_{1}^{2}+u_{2}^{2}+u_{3}^{2}=U^{2}+V^{2}+W^{2}.$
不可压 Euler 方程的 Bernoulli 函数为
$B=\frac{U^{2}+V^{2}+W^{2}}{2}+\frac{p}{\rho}.$
在柱对称坐标系下, 不可压 Euler 方程(1.1)变为
(2.2) $\begin{equation} \left\{\begin{array}{ll} (r\rho U)_{x}+(r\rho V)_{r}=0,\\ (r\rho U^{2})_{x}+(r\rho UV)_{r}+rp_{x}=0,\\ (r\rho UV)_{x}+(r\rho V^{2})_{r}+r p_{r}=\rho W^{2},\\ (r\rho U(rW))_{x}+(r\rho V(rW))_{r}=0,\\ (rU)_{x}+(rV)_{r}=0. \end{array}\right.\label{1.19} \end{equation}$
$\Gamma_{1}=\{(x,r)|x\in(-\infty,+\infty),r=f_2(x)\}$
$\ \Gamma_{2}=\{(x,r)|x\in(-\infty,+\infty),r=f_1(x)\}.$
(2.3) $\begin{equation} (U,V,W)\cdot\overrightarrow{n}=0,\qquad (x,r)\in \Gamma_{1}\backslash\Pi,\label{1.14b} \end{equation}$
其中 $\overrightarrow{n}$ 是 $\Gamma_{1}$ 上的单位外法向量. 滑动边界条件(1.5)变为
(2.4) $\begin{equation} (U,V,W)\cdot \overrightarrow{n}=0,\qquad (x,r)\in \Gamma_{2}\backslash\Pi,\label{1.15a} \end{equation}$
其中 $\overrightarrow{n}$ 是 $\Gamma_{2}\backslash\Pi$ 上的单位外法向量.
2.2 主要定理
本文的主要结果包括适定性和无穷远渐近速率两部分. 适定性部分的定理陈述如下.
定理2.1 如果管道满足(1.2)-(1.5)式, 给定 $\left(W^{2}_-(r),U_-(r),\rho_-(r)\right)\in \left(C^{1,1}\left[1\right]\right)^{3}$ 满足如下条件
(2.5) $\begin{equation}\label{AC2.14} U'_-(0)=\rho'_-(0)=W_-(0)=0, \end{equation}$
(2.6) $\begin{equation}\label{AC2.15} \rho_-(r)>0,\quad U_-(r)>0,\quad r\in\left(0,1\right), \end{equation}$
(2.7) $\begin{equation}\label{MT2.15} \left(\left(U^2_-+W^2_-\right)\rho \right)'(r)>0,\quad \delta r^2\left(\left(U^2_-+W^2_-\right)\rho_- \right)'(r)\geq \left(r^2W^2_-\rho_- \right)'(r), \quad r\in\left(0,1\right), \end{equation}$
(2.8) $\begin{equation}\label{MT2.16} \bigg(\frac{(U_-^2+W_-^2)'+\frac{2W_-^2}{r}}{rU_-}\bigg)'(r)\geq0,\qquad r\in[0,1], \end{equation}$
(2.9) $\begin{equation}\label{MT2.17} \bigg(\frac{(r^2W_-^2)'}{rU_-}\bigg)'(r)\leq0, \quad r\in[0,1], \end{equation}$
(2.10) $\begin{equation}\label{T2.17} \|\rho'_{-}\|_{C^{1,1}}\leq \delta, \end{equation}$
其中 $\delta$ 是依赖于 $U_-$ 和 $W_-^2$ 的正常数, 那么存在唯一的解 $(\rho, U, V, W,p)$ , 其中 $(\rho, U, V, W^2,p)\in (C^{\alpha}(\overline{\Omega}) \cap C^{1,\alpha}(\Omega))^5$ 满足方程(2.2)和边界条件(2.3)以及(2.4)并且满足
(2.11) $\begin{equation}\label{non-degeneracy} U>0,\qquad (x,r)\in \overline{\Omega}\setminus \Pi. \end{equation}$
(ii) 当 $x\rightarrow \pm \infty$ ,
(2.12) $\begin{equation}\label{AC1} \rho\rightarrow \rho_{\pm}, \qquad U\rightarrow U_\pm, \qquad V\rightarrow 0, \qquad W\rightarrow W_\pm, \qquad p\rightarrow p_\pm, \end{equation}$
这里的 $\rho_{+},U_{+},W_{+},p_{+}$ 由 $\rho_{-},U_{-},W_{-},p_{-}$ 和 $m$ 唯一确定.
注2.1 此处的条件(2.5)为对称轴处的对称条件; 条件(2.6)为管道入口处密度和水平速度的正性条件; 条件(2.7)可推出流体轴向速度具有非负性, 从而保证流线在对称轴附近不出现分叉或者回流的条件; 条件(2.8)-(2.10)是为了保证解的唯一性以及渐近状态存在性唯一的条件. 对于不具有涡旋的齐次不可压Euler流, 即: $W_-=0$ , $\rho=1$ , 条件(2.6)-(2.10)可简化为轴向速度的凸性条件, 即
$r U_-''(r)\geq U_-'(r)\geq 0. $
注2.2 此处的条件为: $W_-(r)^{2}\in C^{1,1}\left[1\right]$ , 而不需要 $W_-(r)$ 是连续函数, 即允许存在间断点 $r_0\in(0, 1)$ , 使得
$ \lim_{r\rightarrow r_0^+}W_-(r)=-\lim_{r\rightarrow r_0^-}W_-(r). $
对应于在 $\Omega$ 中存在一条 $C^{2, \alpha}$ 曲线 $\Gamma$ , 其切线平行于 $(rU, rV)$ . $\Gamma$ 将 $\Omega$ 划分为 $\Omega_+$ 和 $\Omega_-$ , 使得
$ \lim_{(x,r)\rightarrow\Gamma\atop (x,r)\subset \Omega_+}W(x,r)=-\lim_{(x,r)\rightarrow\Gamma\atop (x,r)\subset \Omega_-}W(x,r). $
此类关于 $W(x, r)$ 的间断为一类特殊的接触间断, 可参考文献[18 ].
在适定性的基础上, 我们进一步的研究了解的无穷远渐近速率.
定理2.2 如果$W_-(r)$ , $U_-(r)$ 和 $\rho_-(r)$ 满足定理2.1 中的条件,有如下两种情况
情况1 如果边界是在有限长度以外为平直, 即存在 $K>0$ 使得
(2.13) $\begin{equation}\label{ping} f_1(x)=1, \quad x\le-K, \end{equation}$
那么存在 $C>0$ , $\epsilon>0$ 和 $K^{'}> K$ 满足
(2.14) $\begin{equation} \Vert{U-U_{-}}\Vert_{C^{\alpha}(\Omega_{x-1,x+1})}+\Vert{V}\Vert_{C^{\alpha}(\Omega_{x-1,x+1})}\le C{\rm e}^{\epsilon x},\quad x \le -K^{'}. \end{equation}$
情况2 如果边界以代数速率收敛到平直边界, 即 $f_1(x)$ 满足
(2.15) $\begin{equation}\label{DS} \Vert f_1(x)-1\Vert_{C^\alpha([x-1,x+1])}\le C|x|^{-l}, \quad x\le-K^{'},\ \mbox{其中 }\ l>0. \end{equation}$
那么存在 $C>0$ 和 $K^{'}>K$ 满足
(2.16) $\begin{equation} \Vert{U-U_-}\Vert_{C^{\alpha}(\Omega_{x-1,x+1})}+\Vert{V}\Vert_{C^{\alpha}(\Omega_{x-1,x+1})}\le C|x|^{-l},\quad x \le -K^{'}, \end{equation}$
其中 $\Omega_{x-1,x+1}=\Omega\cap\left[x-1,x+1\right]\times\left[1\right]$ .
注2.3 值得注意的是, 相比于柱对称可压缩Euler方程已有的结果, 此定理可以包括漩涡速度 $W$ 不为 0 的情况, 并且在解的适定性条件基础上不需要添加额外的条件.
2.3 负无穷远处的状态
我们先假设 $V_-\equiv0$ , 这将在命题3.2中证明. 当 $x\rightarrow-\infty$
(2.17) $\begin{equation} \label{1.24} \rho\rightarrow \rho_-(r), \end{equation}$
(2.18) $\begin{equation}\label{1.25} U\rightarrow U_-(r), \end{equation}$
(2.19) $\begin{equation} \frac{1}{2}(U^2+V^2+W^2)+\frac{p}{\rho}\rightarrow B_-(r) \end{equation}$
(2.20) $\begin{equation}\label{1.26x} W\rightarrow W_-(r). \end{equation}$
给定的 $U_-(r)$ , $\rho_-(r)$ , $W_{-}(r)$ , 那么 $B_{-}(r)$ 在入口处满足
(2.21) $\begin{equation}\label{mass-infty} B_-(r)= \frac{(U_{-}(r))^2}{2}+\frac{(W_{-}(r))^2}{2}+\frac{ p_{-}(r)}{\rho_-(r)}. \end{equation}$
由式$(2.2)_3$ 和 $V_-\equiv0$ 得
(2.22) $\begin{equation} \frac{\partial {p_-(r)}}{\partial{r}}=\frac{\rho_{-}(r)W_-^{2}(r)}{r}, \end{equation}$
(2.23) $\begin{equation}\label{equ:3.6} p_{-}(r)=p_{-c}+\int_0^rH_{-}(s){\rm d}s, \end{equation}$
其中 $H_{-}(r)$ =$\frac{\rho_{-}(r)W_-^{2}(r)}{r}$ 和 $p_{-c}$ 是任意常数表示在 $r=0$ 处的压力. 由 $(2.2)_4$ 式可知, 存在守恒量 $m$ 满足
(2.24) $\begin{equation} m=\int_0^1s U_{-}(s){\rm d}s. \end{equation}$
2.4 流线守恒量和流函数的方程
流线守恒量来自于 Euler 方程组的线性退化部分
(2.25) $\begin{equation} \left\{\begin{array}{ll} (rU)_{x}+(rV)_{r}=0,\\ (r\rho U)_{x}+(r\rho V)_{r}=0,\\ (rU(rW))_{x}+(rV(rW))_{r}=0,\\ (rUB)_{x}+(rVB)_{r}=0. \label{LDF} \end{array}\right. \end{equation}$
根据方程 $(2.25)_1$ 和 $\Omega$ 是单连通区域得, 存在一个流函数 $\psi$ : 满足
(2.26) $\begin{equation} \psi_{x}=-rV,\ \psi_{r}=rU.\label{2.1} \end{equation}$
对于 $r\in (0, 1)$ , 定义一个从 $\left[1\right]$ 到 $\left[m\right]$ 一一映射如下
(2.27) $\begin{equation}\label{equ:3.8xx} \psi_-(r)=\int_0^r sU_{-}(s) ds. \end{equation}$
令$ \kappa(m):=\psi^{-1}_{-}(m) $ 为 $\left[m\right]$ 到 $\left[1\right]$ 的一一映射. 直接计算得到
(2.28) $\begin{equation} \kappa'(\psi_-)=\frac{1}{\kappa U_{-}(\kappa)} \quad\mbox{和}\quad \kappa''(\psi_-)=\left(-\frac{1}{\kappa(\psi_-)}-\frac{U_{-}'(\kappa)}{U_{-}(\kappa)}\right)(\kappa')^2(\psi_-). \end{equation}$
由于$(2.25)_2$ - $(2.25)_4$ 式以及 $rW$ , $B$ 和 $\rho$ 是流线守恒量, 于是, 当 $0\leq\psi\leq m$ 时,定义
$ {\cal G}(\psi):=\frac{1}{\rho_-(\kappa(\psi))},\ {\cal B}(\psi):=B_-(\kappa(\psi))\ \mbox{和}\ {\cal W}(\psi):=\kappa(\psi)W_-(\kappa(\psi)). $
(2.29) $\begin{equation}\label{equ:3.4} \left\{\begin{array}{ll} \psi_{x}=-rV,\\ \psi_{r}=rU,\\ {\cal B}(\psi)=\frac{1}{2}(U^{2}+V^{2}+W^2)+\frac{p}{\rho},\\ {\cal G}(\psi)=\frac{1}{\rho},\\ {\cal W}(\psi)=rW. \end{array}\right. \end{equation}$
(2.30) $\begin{equation} \rho\partial_{x}B=\rho UU_{x}+\rho VV_{x}+\rho WW_{x}+p_{x}-\frac{p\rho_{x}}{\rho}. \end{equation}$
$ U_r-V_x={\cal K}\left(\left|\frac{\nabla_{(x, r)}\psi}{r}\right|^2,\psi,r\right):=r\left( {\cal B}'-p{\cal G}'\right)-\frac{{\cal W}{\cal W}'}{r}.$
(2.31) $\begin{equation} {\cal B}=\frac{1}{2}\left|\frac{\nabla_{(x, r)}\psi}{r}\right|^2+\frac{{\cal W}^2}{2r^2}+p{\cal G}. \end{equation}$
(2.32) $\begin{equation} p=\frac{1}{{\cal G}} \bigg({\cal B}-\frac{1}{2}\left|\frac{\nabla_{(x, r)}\psi}{r}\right|^2-\frac{{\cal W}^2}{2r^2}\bigg). \end{equation}$
$ U_r-V_x=\left(\frac{\psi_r}{r}\right)_{r}+\left(\frac{\psi_x}{r}\right)_x.$
因此, 我们得到关于 $\psi$ 的如下二阶方程
(2.33) $\begin{equation} \mbox{div}_{(x, r)}\left(\frac{\nabla_{(x, r)} \psi}{r}\right)={\cal K}\left(\left|\frac{\nabla_{(x, r)}\psi}{r}\right|^2,\psi,r\right). \label{2.8} \end{equation}$
由边界条件(2.3)和(2.4), 得 $\psi$ 在 $\Gamma_1$ 和 $\Gamma_2$ 上分别都是常数, 即
$ \psi=0,\quad (x,r)\in \Gamma_1\qquad \mbox{和}\qquad \psi=m \quad (x,r)\in \Gamma_2. $
因此, 将不可压 Euler 方程定值问题转化为关于 $\psi$ 的二阶方程的边值问题.
(2.34) $\begin{equation} \left\{\begin{array}{ll} \mbox{div}_{(x, r)}\left(\frac{\nabla_{(x, r)} \psi}{r}\right)={\cal K}\left(\left|\frac{\nabla_{(x, r)}\psi}{r}\right|^2,\psi,r\right),\quad& (x,r)\in\Omega,\\ \psi|_{\Gamma_1}=0, \psi|_{\Gamma_2}=m. \end{array}\right.\label{2.9} \end{equation}$
(2.35) $\begin{equation}\label{equ:3.43} \psi_r>0,\qquad (x,r)\in\overline{\Omega}\setminus\pi. \end{equation}$
令$\partial_1:=\partial_x$ 和$\partial_2:=\partial_r$ . 将方程(2.33)写成非散度形式.
(2.36) $\begin{equation}\label{E3.13} \Delta\psi+b_{i}\partial_{i}\psi={\cal F}(\psi,r), \end{equation}$
$\begin{eqnarray*} &&b_1(\psi,r)=-\frac{\partial_1\psi{\cal G}'}{2{\cal G}},\\ &&b_2(\psi,r)=-\frac{\partial_2\psi{\cal G}'}{2{\cal G}}-\frac{1}{r},\\ &&{\cal F}(\psi,r)=r^2{\cal G}\left(\frac{{\cal B}}{{\cal G}}\right)'-\frac{{\cal G}}{2}\left(\frac{{\cal W}^2}{{\cal G}}\right)'. \end{eqnarray*}$
3 定理1.1的证明
本节中, 我们将求解问题1 即证明定理 2.1. 首先考虑没有边界没有角点的情况, 即 $\Gamma_1$ 和 $\Gamma_2$ 为 $C^{2, \alpha}$ 边界.
将 ${\cal G}(\psi)$ , ${\cal B}(\psi)$ 和 ${\cal W}(\psi)$ 从 $[m]$ 延拓到 $\Bbb R$ . 令
(3.1) $\begin{equation}\dot{{\cal G}}(s)=\left\{ \begin{array}{ll} {\cal G}'(s),& \mbox{当}\ 0\leq s\leq m,\\ {\cal G}'(m)\frac{2m-s}{m},&\mbox{当}\ m\leq s\leq 2m,\\ 0,& \mbox{当}\ s\geq2m\ \mbox{或}\ s\leq 0, \end{array}\right. \end{equation}$
(3.2) $\begin{equation} \label{4.3} \dot{{\cal B}}(s)=\left\{ \begin{array}{ll} \left(\frac{{\cal B}}{{\cal G}}\right)'(s),& \mbox{当}\ 0\leq s\leq m,\\ \left(\frac{{\cal B}}{{\cal G}}\right)'(m)\frac{2m-s}{m},&\mbox{当}\ m\leq s\leq 2m,\\ 0, & \mbox{当}\ s\geq2m\ \mbox{或}\ s\leq-0 \end{array}\right. \end{equation}$
(3.3) $\begin{equation}\dot{{\cal W}^2}(s)=\left\{ \begin{array}{ll} \left(\frac{{\cal W}^2}{{\cal G}}\right)'(s),& \mbox{当}\ 0\leq s\leq m,\\ \left(\frac{{\cal W}^2}{{\cal G}}\right)'(m)\frac{2m-s}{m}, & \mbox{当}\ m\leq s\leq 2m,\\ 0,& \mbox{当}\ s\geq2m\ \mbox{或}\ s\leq 0. \end{array}\right. \end{equation}$
(3.4) $\begin{matrix} \label{3.16} \tilde{{\cal G}}(s)&=&{\cal G}(0)+\int^{s}_{0}\dot{{\cal G}}(t){\rm d}t, \\ \tilde{{\cal B}}(s)&=&\tilde{{\cal G}}(s)\left(\frac{{\cal B}}{{\cal G}}(0)+\int^{s}_{0}\dot{{\cal B}}(t){\rm d}t\right),\\ \tilde{{\cal W}^2}(s)&=&\tilde{{\cal G}}(s)\left(\frac{{\cal W}^2}{{\cal G}}(0)+\int^{s}_{0}\dot{{\cal W}}^{2}(t){\rm d}t\right). \end{matrix}$
显然, $\tilde{{\cal G}}(s),\ \tilde{{\cal B}}(s),\ \tilde{{\cal W}^2}(s)\in C^{1,1}(\Bbb R)$ . 当 $0\leq s\leq m$ 时, 直接计算得到
$\begin{eqnarray*} {\cal K}(0,s,r)&=&\left(\frac{\tilde{{\cal B}}}{\tilde{{\cal G}}}\right)'(s) -\frac{1}{2r^2}\left(\frac{\tilde{{\cal W}}^{2}}{\tilde{{\cal G}}}\right)'(s) \\ &=&\kappa'\left(\left(\frac{\left(U^2_-+W^2_-\right)\rho_-}{2}\right)'(\kappa)+\frac{\rho_{-}W_{-}^{2}}{\kappa}(\kappa) -\frac{1}{2r^2}\left(\rho_{-}W_-^{2}\right)'(\kappa)\right). \end{eqnarray*}$
$\begin{eqnarray*} &&\left(\frac{\left(U^2_-+W^2_-\right)}{2}\rho_-\right)'(\kappa)+\frac{\rho_{-}W_{-}^{2}}{\kappa}(\kappa) -\frac{1}{r^2}\left(\rho_{-}\frac{W_-^{2}}{2}\right)'(\kappa)\\ &\geq& \left(1-\delta\right)\left(\frac{\left(U^2_-+W^2_-\right)}{2}\rho_-\right)'(\kappa)+\frac{\rho_{-}W_{-}^{2}}{\kappa}(\kappa)\\ &\geq& \delta'>0. \end{eqnarray*}$
根据 $\tilde{{\cal G}}$ , $\tilde{{\cal B}}$ 和 $\tilde{{\cal W}}^{2}$ 的定义得
(3.5) $\begin{equation}\label{max} \left(\frac{\tilde{{\cal B}}}{\tilde{{\cal G}}}\right)'(s) -\frac{1}{2r^2}\left(\frac{\tilde{{\cal W}}^{2}}{\tilde{{\cal G}}}\right)'(s)\geq 0, \quad s\in\Bbb R, \end{equation}$
当 $0<s<m$ 时, ${\cal K}(0,s,r)$ 严格大于$0$ .
3.1 解的存在性
(3.6) $\begin{equation}\label{approximated} \left\{\begin{array}{ll} \Delta\psi^{(k)}_L+b_{i}^{(k)}\partial_i\psi^{(k)}_L={\cal F}\left(\psi^{(k)}_L,r+k\right), & (x,r)\in \Omega_L,\\ \psi^{(k)}_L=\frac{r^{2}m}{f_{1}^{2}(x)},\ \ \ & (x,r)\in \partial \Omega_L. \end{array}\right. \end{equation}$
$ \{(x,r)|(x,r)\in\Omega,\ |x|<L\}\subset\Omega_{L}\subset\{(x,r)|(x,r)\in\Omega,\ |x|<4L\}, $
$0<k\ll1$ , $L \geq\max\{|L_1|,|L_2|\}$ 和
$ b_i^{(k)}=b_i\left(\psi^{(k)}_{L},r+k\right). $
这里$\partial\Omega_{L}$ 是 $C^{2,\alpha_1}$ 边界, $0<\alpha_1<\alpha<1$ , 并满足一致外球条件. 对于固定的 $k>0$ 和 $L>0$ , 由文献[33 ]知,问题(3.6)存在唯一解 $\psi^{(k)}_L\in C^{2,H}(\Omega_L)\cap C^0(\bar{\Omega}_L)$ .
步骤2 证明解在 $r=0$ 处不会有奇性. 通过式(3.5)得到
$\begin{eqnarray*} {\cal F}\left(\psi^{(k)}_L,r+k\right) =\left(r^2\tilde{{\cal G}}\left(\frac{\tilde{{\cal B}}}{\tilde{{\cal G}}}\right)'-\frac{1}{2}\tilde{{\cal G}}\left(\frac{\tilde{{\cal W}}^{2}}{\tilde{{\cal G}}}\right)'\right) \geq 0. \end{eqnarray*}$
由强极大值原理得 $\psi^{(k)}_L$ 不能在内部取得极大极小值, 即在内部满足$0<\psi^{(k)}_L<m$ . 这意味着 $\tilde{{\cal G}}={\cal G}$ , $\tilde{{\cal B}}={\cal B}$ 和 $\tilde{{\cal W}}={\cal W}$ . 令 $\phi(r)=\frac{m}{b^2}(r+k)^2$ , 其中 $b=\inf\limits_{x\in\Bbb R}f(x)$ . 如果存在一个临界点 $P_{\mbox{max}}=(x_0,r_0)$ 使得
$\psi^{(k)}_L(P_{\mbox{max}})\geq \phi(P_{\mbox{max}}),\quad \nabla\psi^{(k)}_{L}=\nabla\phi= \left(0,\frac{2m}{b^2}(r+k)\right),$
$\begin{eqnarray*} \Delta\left(\psi^{(k)}_{L}-\phi\right)+b_{i}^{(k)}\partial_i\left(\psi^{(k)}_{L}-\phi\right) =r^2{\cal G}\left(\frac{{\cal B}}{{\cal G}}\right)'-\frac{{\cal G}}{2}\left(\frac{{\cal W}^{2}} {{\cal G}}\right)' -\frac{{\cal G}'}{2{\cal G}}\frac{4m^2}{b^4}(r+k)^2. \end{eqnarray*}$
(3.7) $\begin{equation}\label{C3x} \Delta(\psi^{(k)}_{L}-\phi)+b_{i}^{(k)}\left(\psi^{(k)}_{L},r+k\right)\partial_i(\psi^{(k)}_{L}-\phi)\geq0. \end{equation}$
显然, 在边界上 $\psi^{(k)}_{L}\leq \frac{m}{b^2}(r+k)^2$ . 因此, 通过比较原理得到
$ 0\leq \psi^{(k)}_{L}\leq \frac{m}{b^2}(r+k)^2.\label{C3} $
于是, 存在 $b>0$ , 当 $k\rightarrow0$ 时, $\psi^{(k)}_{L}\rightarrow\psi_L$ 并满足
$\left|\frac{\nabla\psi_{L}}{r}\right|\leq Cm,\qquad 0<r<\frac{b}{2}.$
首先考虑无旋流的非退化性. 对于无旋流, $\underline{B}$ 和 $\underline{G}$ 是分别取值于 ${\cal B}$ 和 ${\cal G}$ 值域的正常数, 以及漩涡速度 $W=0$ . 当流体是无旋流, Euler 方程(2.2)变为
(3.8) $\begin{equation}\label{equ:5.2} \left\{\begin{array}{ll} (r U)_{x}+(r V)_{r}=0,\\ U_r-V_x=0, \end{array}\right. \end{equation}$
对应的 $\bar{\psi}_{L}$ 满足下列边值问题
(3.9) $\begin{equation}\label{irrotational} \left\{\begin{array}{ll} \mbox{div}_{(x, r)}\left(\frac{\nabla\bar{\psi}}{r}\right)=0, \quad &(x, r)\in \Omega_{L},\\ \bar{\psi}=\frac{mr^2}{f^{2}_{1}(x)},\ \quad &(x, r)\in \partial\Omega_{L}. \end{array}\right. \end{equation}$
由步骤2可知在 $\Omega^o$ 内 $0<\bar{\psi}_{L}<m$ . 因为 $\bar{\psi}_{L}$ 在 $\Gamma_1\cap\partial\Omega_{L}$ 上有极小值, 利用 Hopf 引理得到
(3.10) $\begin{equation} \frac{\partial\bar{\psi}_{L}}{\partial\eta}=-\partial_{1}\bar{\psi}_{L} f'_{2}(x)+\partial_{2}\bar{\psi}_{L}>0, \end{equation}$
其中 $\eta$ 是 $\Omega_{L}$ 上的单位外法向量. 又因为 $\bar{\psi}_{L}$ 在 $\Gamma_1\cap\partial\Omega_{L}$ 等于0, 于是得到
(3.11) $\begin{equation} \partial_{1}\bar{\psi}_{L}+\partial_{2}\bar{\psi}_{L} f'_{2}(x)=0, \end{equation}$
(3.12) $\begin{equation} \frac{\partial\bar{\psi}_{L}}{\partial\eta}=\partial_{2}\bar{\psi}_{L}(1+f'_{2}(x)^{2})>0. \end{equation}$
这表明在 $\Gamma_1\cap\partial\Omega_{L}$ 上 $\partial_{2}\bar{\psi}_{L}>0$ . 同理, 得到在其他边界上 $\partial_{2}\psi_{L}>0$ , 即
(3.13) $\begin{equation}\label{L-positve} \mbox{在 }\ \partial \Omega_L\ \mbox{上, }\ U>0. \end{equation}$
根据 $U_r-V_x=0$ , 得到一个位势函数 $\phi$ 满足
(3.14) $\begin{equation} \phi_x=U, \quad \phi_x=V. \end{equation}$
令 $\Phi\left(x_{1},x_{2},x_{3}\right)=\phi\left(x,r\right)$ 得
(3.15) $\begin{equation} \Phi_{x_{1}}=\phi_x, \quad \Phi_{x_{2}}=\phi_r\frac{x_{2}}{r},\quad \Phi_{x_{3}}=\phi_r\frac{x_{3}}{r}. \end{equation}$
由不可压条件 $\mbox{div }u=0$ 得
(3.16) $\begin{equation} \Delta\Phi=0. \end{equation}$
(3.17) $\begin{equation} \Delta U=0. \end{equation}$
由(3.13)式, 利用强极大值原理得到在 $\overline{\Omega}_{L}$ 上 $U>0$ .
$\begin{eqnarray*} \phi=\psi_L-\bar{\psi}_L\qquad \mbox{和}\qquad\psi_{\tau}=\bar{\psi}_L+\tau(\psi_L-\bar{\psi}_L),\quad \mbox{其中 }\ \tau\in(0,1). \end{eqnarray*} $
(3.18) $\begin{equation} \mbox{div}\left(\frac{\nabla\psi_L-\nabla\bar{\psi}_L}{r}\right)=r\left({\cal B}'-p\left({\cal G}\right)'\right)-\frac{{\cal W}{\cal W}'}{r}. \end{equation} $
(3.19) $\begin{equation} \Delta \phi-\frac{1}{r}\partial_2\phi=r^2\left({\cal B}'(\psi)-p\left({\cal G}\right)'(\psi)\right)-{\cal W}{\cal W}'. \end{equation}$
因此, 对于任意的紧子集 $K\subset \Omega_{L}$ , 如果 $\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])}\leq \delta$ , 那么
(3.20) $\begin{matrix} \|\phi\|_{C^{2,\alpha}(K)}\leq C\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])} \leq C\delta.\label{7.16} \end{matrix} $
$ \|\psi_L-\bar{\psi}_L\|_{C^{2, \alpha}(K)}\leq C\delta. $
由于 $\partial_2\bar{\psi}_L>0$ , 选择足够小的 $\delta$ , 得到 $\partial_2\psi_L>0$ . 假设存在最大的有界正数 $\delta_0$ , 使得当 $\delta\leq\delta_0$ 时, $\partial_2\psi_L>0$ ; 当 $\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])}= \delta_0$ 时, 在某些点上 $\partial_2\psi_L\leq 0$ . 那么通过选取一个序列 $({\cal G}_n, {\cal B}_n, {\cal W}_n)$ , 其中对于每一个 $n\in \mathbb{N}$ , $\|({\cal G}'_n,{\cal B}'_n,{\cal W}'_n)\|_{C^{0,1}([m])}\leq \delta_0$ , 当 $n\rightarrow \infty$ 时, $({\cal G}_n, {\cal B}_n, {\cal W}_n)\rightarrow ({\cal G}, {\cal B}, {\cal W})$ , 其中 $\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])}=\delta_0$ . 记 $\psi_{L}^{(n)}$ 和 $\psi_L$ 分别是 $({\cal G}_n, {\cal B}_n, {\cal W}_n)$ 和 $({\cal G}, {\cal B}, {\cal W})$ 对应问题(3.6)的唯一解. 由有界区域解的唯一性和 $C^{2,\alpha}$ 估计, 当 $n\rightarrow\infty$ 得到 $\psi_{L}^{(n)}\rightarrow\psi_L$ . 由 $\partial_2 \psi_{L}^{(n)}>0$ , 可知 $\partial_2 \psi_{L}\geq0$ . 由文献[34 ]可知如果 $\partial_2 \psi_{L}\geq0$ , 那么 $\partial_2 \psi_{L}>0$ , 与 $\delta_0$ 的定义矛盾. 因此, 不存在有界的 $\delta_0$ . 从而, 对于任意的 $\delta>0$ 都有 $\partial_r \psi_{L}>0$ .
步骤4 证明 $\psi_{L}$ 关于 $L$ 的一致估计.
(3.21) $\begin{equation} \left|\psi_L\right|\leq m+C\sup\limits_{\Omega_L}\tilde{{\cal F}}, \end{equation}$
其中 $C={\rm e}^{\bar{f}}-1$ . 因此
(3.22) $\begin{equation} -C|\sup\limits_{\Omega_L}\tilde{{\cal F}}|\leq\psi_L\leq m+C|\sup\limits_{\Omega_L}\tilde{{\cal F}}|. \end{equation}$
利用文献[33 ,引理 6.20], 存在 $C$ 满足
(3.23) $\begin{equation} [\psi_L]_{1,\mu;\Omega_L}\leq C\left(1+m+\|\tilde{{\cal F}}\|_{0;\Omega_L}\right) \end{equation}$
(3.24) $\begin{equation} \|\psi_L\|^{(-1-\alpha)}_{2,\alpha,\Omega_L}\leq C\left(\|\tilde{{\cal F}}\|_{0;\Omega_L}+\|\psi_L\|_{0;\Omega_L}\right), 0<\alpha<\mu. \end{equation}$
对于任意的紧子集 $K\in\bar{\Omega}$ , 利用 Arzela-Ascoli 引理, $\{\psi_{L}\}$ 存在一个子列 $\psi_{L_n}$ 在 $C^{1,\alpha/2}(K)$ 意义下收敛到 $\psi$ . 又因为 $\psi_{L}\in C^{2, \alpha}(\Omega)$ . 显然, $\psi$ 是(4.22)的一个解且$\psi_r \geq 0$ . 证毕.
3. 解的唯一性
命题3.1 给定$U_-(r)$ , $\rho_-(r)$ ,$W_-(r)$ 满足定理 2.1 的假设, 那么问题 1的解是唯一的.
证 如果 $\psi_1$ 和 $\psi_2$ 都上述存在性结论的解, 即满足
(3.25) $\begin{equation} \left\{\begin{array}{ll} \mbox{div}_{(x, r)}\left(\frac{\nabla \psi_1-\nabla \psi_2}{r}\right)={\cal K}\left(\left|\frac{\nabla\psi_1}{r}\right|^2,\psi_1,r\right)-{\cal K}\left(\left|\frac{\nabla\psi_2}{r}\right|^2,\psi_2,r\right),\quad &(x,r)\in \Omega,\\ \psi_1-\psi_2=0,\quad& (x,r)\in\partial\Omega. \end{array}\right. \end{equation}$
$\begin{eqnarray*} \phi=\psi_1-\psi_2\qquad \mbox{和}\qquad\psi_{\tau}=\psi_1+\tau(\psi_1-\psi_2). \end{eqnarray*}$
(3.26) $\begin{equation}\label{UN} \Delta \phi-\frac{1}{r}\phi_{r}-r\widehat c_i\partial_i\phi-r\widehat d\phi=0, \end{equation}$
$\begin{eqnarray*} \widehat c_{i}\left(\frac{\nabla\psi_\tau}{r},\psi_\tau,r\right) =\int_{0}^{1}\frac{{\cal G}'\partial_{i}\psi_\tau}{2{\cal G}r}{\rm d}\tau,\qquad \widehat d\left(\frac{\nabla\psi_\tau}{r},\psi_\tau\right)=\int_{0}^{1}\partial_{\psi_\tau}{\cal K}\left(\left|\frac{\nabla\psi_\tau}{r}\right|^{2},\psi_\tau,r\right){\rm d}\tau. \end{eqnarray*}$
取 $\eta\in C^{\infty}_0(\Bbb R)$ 满足: 当 $|s|<L$ 时, $\eta=1$ ; 当 $|s|>L+1$ 时, $\eta=0$ . 将式(3.26)乘 $\eta^2(x)\phi$ , 再在 $\Omega$ 上积分得
$ \iint_{\Omega}\Delta \phi\eta^2(x)\phi {\rm d}x{\rm d}r=\iint_{\Omega}\left(\frac{1}{r}\phi \phi_{r}+r\widehat c_i\partial_{i}\phi\eta^2(x)\phi+r\widehat d\eta^2(x)\phi^2\right){\rm d}x{\rm d}r. $
$\begin{eqnarray*} &&\iint_{\Omega}\eta^2|\nabla \phi|^2{\rm d}x{\rm d}r\nonumber\\ &=&-2\iint_{\Omega}\phi\eta\nabla \eta \cdot\nabla\phi {\rm d}x{\rm d}r-\iint_{\Omega}\left(\frac{1}{r}\phi_{r}+r\widehat c_i\partial_{i}\phi\right)\eta^2(x)\phi {\rm d}x{\rm d}r-\iint_{\Omega}r\widehat d \eta^2(x)\phi^2{\rm d}x{\rm d}r.\nonumber \end{eqnarray*}$
对于右端最后一项, 我们考虑 $ \widehat d$ 中积分内的项.
$\begin{eqnarray*} d&=&r{\cal B}''-rp{\cal G}''-r\frac{{\cal B}'{\cal G}'-r^{-2}{\cal W}{\cal W}'{\cal G}'-p({\cal G}')^2}{{\cal G}}-\frac{({\cal W}')^2+{\cal W}{\cal W}''}{r}\nonumber\\ &=&r(B_{-}''(\kappa')^2+B_{-}'(\kappa''))\nonumber\\ &&-(\kappa^2 W_{-}^2)''(\kappa')^2+ r (\kappa^2 W_{-}^2)'(\kappa)''\nonumber\\ &&-(\kappa')^2\left(-r\rho_{-}'B_{-}'+r^{-1}(\kappa^2W_{-})'{\rho_{-}}'+\rho_{-}^{-3}p_-\left(\rho_{-}'\right)^2\right)\nonumber\\ &&-rp\left(\left(\frac{1}{\rho_{-}}\right)''(\kappa')^2+\left(\frac{1}{\rho_{-}}\right)'\kappa''\right)\nonumber\\ &:=&\sum\limits_{i=1}^{4}I_{i}.\nonumber \end{eqnarray*}$
$\begin{eqnarray*} &&{\cal B}'(\psi)=\kappa'B_-',\qquad {\cal B}''(\psi)=(\kappa')^2B_-''+\kappa''B_-';\nonumber\\ &&{\cal G}'(\psi)=\kappa'\rho_-',\qquad {\cal G}''(\psi)=(\kappa')^2\rho_-''+\kappa''\rho_-';\nonumber\\ & &{\cal W}'(\psi)=\kappa'W_-+\kappa\kappa'W_-',\qquad {\cal W}''(\psi)=\kappa(\kappa')^2W_-''+\kappa\kappa''W_-'+2(\kappa')^2W_-'+\kappa''W_-.\nonumber \end{eqnarray*}$
$\begin{eqnarray*} &&B_{-}''(\kappa)(\kappa')^2+B_{-}'\kappa''\nonumber\\ &=&(\kappa')^2\left(\left(-\frac{1}{\kappa}-\frac{U_{-}'}{U_-}\right)\left(\left(\frac{U_{-}^2+W_{-}^2}{2}\right)'+\frac{W_{-}^2}{\kappa}\right)(\kappa)+\frac{W_{-}^2}{\kappa}'+\left(\frac{U_{-}^2+W_{-}^2}{2}\right)''(\kappa)\right)\nonumber\\ \nonumber\\ &&+(\kappa')^2\left(p_{-}\left(\frac{1}{\rho_{-}}\right)'\right)'(\kappa)+\kappa''p_{-}\left(\frac{1}{\rho_{-}}\right)' \nonumber\\ &=&(\kappa')^2\kappa U_{-}(\kappa)\left(\frac{\left(\frac{U_{-}^2+W_{-}^2}{2}\right)'(\kappa)+\frac{W_{-}^2}{\kappa}}{\kappa U_{-}}\right)'+(\kappa')^2\left(p_-\left(\frac{1}{\rho_{-}}\right)'\right)'+\kappa(\psi_-)''p_-\left(\frac{1}{\rho_{-}}\right)'.\nonumber \end{eqnarray*}$
(3.27) $\begin{equation} B_{-}''(\kappa)(\kappa')^2+B_{-}'\kappa''\geq -\frac{\delta'}{4}. \end{equation}$
由条件(2.9)得 $I_{2}\geq -\frac{\delta'}{4}$ . 对于$I_{3}$ 和 $I_{4}$ , 可由(2.10)条件控制. 综上可知: $ d\geq -\delta'$ . 于是
$-\iint_{\Omega}r\widehat d \eta^2(x)\phi^2{\rm d}x{\rm d}r\leq \delta'\iint_{\Omega}r\eta^2(x)\phi^2{\rm d}x{\rm d}r.$
$ \iint_{\Omega}\frac{\eta^2 \phi^2}{r^{2}}dxd\leq C\iint_{\Omega}|\nabla(\eta \phi)|^2{\rm d}x{\rm d}r, $
$\iint_{\Omega}\phi^2{\rm d}x{\rm d}r\leq \iint_{\Omega}|\nabla \phi|^2{\rm d}x{\rm d}r,$
以及 $\delta'$ 的小性, 我们可以得到: 对于任意 $L>0$ ,
(3.28) $\begin{equation} \iint_{\Omega_{-L,L}}|\nabla \phi|^2{\rm d}x{\rm d}r\leq 0. \end{equation}$
因为 $\phi$ 在边界上等于0, 于是得, 在 $\Omega$ 上, $\phi= 0$ , 即 $\psi_1=\psi_2$ . 证毕.
命题3.2 对于$a_{B}>0$ , 如果 $\bar{\psi}$ 是下列方程的解
(3.29) $\begin{equation} \left\{\begin{array}{ll}\label{FlatProblem} {\rm div}_{(x, r)}\left(\frac{\nabla_{(x, r)}\bar{\psi}}{r}\right)={\cal K}\left(\left|\frac{\nabla\bar{\psi}}{r}\right|^2,\bar{\psi},r\right) \quad (x,r)\in D,\\ \bar{\psi}=0,\quad\ r=0,\quad \bar{\psi}=m, \quad r=a_B, \end{array}\right. \end{equation}$
其中 $D=\left\{\left(x,r\right)|x\in \left(-\infty,+\infty\right),r \in\left(0, a_B\right) \right\}$ , 则在 $D$ 上 $\partial_x\bar{\psi}=0$ .
(3.30) $\begin{equation} \partial_x\bar{\psi}=\lim_{\Delta x\rightarrow0}\frac{\bar{\psi}\left(x+\Delta x\right)-\bar{\psi}\left(x\right)}{\Delta x}. \end{equation}$
如果 $\bar{\psi}$ 是方程(3.2)的解, 那么 $\bar{\psi}\left(x+\Delta x\right)$ 也是方程(3.2)的解. 根据命题3.1得到
(3.31) $\begin{equation} \bar{\psi}\left(x+\Delta x\right)=\bar{\psi}\left(x\right). \end{equation}$
因此, 在 $D$ 上 $\partial_x\bar{\psi}=0$ . 证毕.
令 $\psi$ 是问题 1 的解, 那么存在一个序列 $\psi^{(n)}=\psi\chi_{\{0<r<f(x-n)\}}$ 在 $C^{2,\vartheta}(K)$ 意义下收敛到 $\bar{\psi}_{-}$ , 其中 $\bar{\psi}_{-}$ 满足: 对应于$a_B=1$ 的方程(3.29). 由(3.29)解的唯一性可知 $\bar{\psi}_{-}=\psi_{-}$ . 结合命题 3.2 可知: $\partial_x\bar{\psi}_{-}=V_-\equiv0$ , 并且, 当 $x\rightarrow - \infty$ ,
(3.32) $\begin{equation} \rho\rightarrow \rho_{-}, \qquad U\rightarrow U_-, \qquad V\rightarrow 0, \qquad W\rightarrow W_-, \qquad p\rightarrow p_-. \end{equation}$
相似地, 存在一个序列 $\psi^{(n)}=\psi\chi_{\{0<r<f(x+n)\}}$ 收敛到 $\bar{\psi}_{+}$ . $\bar{\psi}_{+}$ 满足: 对应于$a_B=a$ 的(3.29). 由(3.29)解的唯一性可知
$ \lim_{x\rightarrow+\infty}\psi(x,r)=\bar{\psi}_{+}, $
由此可以确定, 当 $x\rightarrow +\infty$ ,
(3.33) $\begin{equation} \rho\rightarrow \rho_{+}, \qquad U\rightarrow U_+, \qquad V\rightarrow 0, \qquad W\rightarrow W_+, \qquad p\rightarrow p_+. \end{equation}$
下面我们将证明 $U_{+}(r)>0$ . 由已有的分析可以得到 $U_{+}(r)=\partial_r\bar{\psi}_{+}(r)\geq0$ . 与式(3.13)类似, 利用 Hopf 引理, 得到 $U_+(0)>0$ 和 $U_+(a)>0$ . 假设存在一点 $r_0\in\left(0,1\right)$ 使得 $ U_{+}(r_0)=0$ , 那么$r_0$ 是 $ U_{+}(r)$ 的内部极小值点,由此可知$ \partial _r U_{+}(r_0)\leq 0$ . $ \partial _r U_+={\cal K}\left(0,\bar{\psi}_{+},r\right). $ 由(3.5)得 $\partial _r U_{+}(r_0)>0$ , 矛盾. 故 $U_{+}(r)>0$ .
进一步的由极大值原理, 类似于文献[23 ], 我们可以得到在 $\Omega$ 上 $U(x, r)>0$ . 由此对于 $f_1$ 和 $f_2$ 不具有角点的情况证明了定理1.
3.3 具有角点的情况
我们用一组光滑的边界去逼近角点附近的区域. 首先考虑下边界 $\Gamma_{2}$ . 令 $\Omega^{\varepsilon}$ 是 $C^{2,\alpha}$ 边界, $r=f^{\varepsilon}_{2}(x)$ 的正则性取决于 $\varepsilon$ . $\Omega^{\varepsilon}\subset\Omega$ 和 $f^{\varepsilon}_{2}(x)\rightarrow f_{2}(x)$ , 并且存在足够大的 $X>0$ 满足对于任意的 $\varepsilon>0$ 和 $|x|\geq X$ , 都有 $f^{\varepsilon}_{2}=f_{2}$ . 不失一般性, 我们考虑角度为 $\theta_0$ 的角点 $P$ , 这个区域为 $\Omega(R,\varepsilon):=\Omega^{\varepsilon}\cap B_R(P)$ . 令 $\Gamma_{+}$ 和 $\Gamma_-$ 是 $P$ 点引出来的左右两条切线. 令 $(s,\theta)$ 是以 $P$ 为原点, $\Gamma_+$ 上 $\theta=0$ 的极坐标系. $\Gamma_-$ 上的角度为 $\theta_0\in(0,\pi)$ . 记 $\theta_*=\frac{\pi-\theta_0}{2}$ . 考虑如下闸函数
(3.34) $\begin{equation} w=\frac{C_1m}{R^{1+\beta}}s^{1+\beta}\sin(A\theta+\theta_*), \end{equation}$
其中 $A\in(1,\frac{\pi+\theta_0}{2\theta_0})$ , $C_1=\max(\csc\theta_*,\csc(A\theta_0+\theta_*))$ , 和 $\beta\in(0, 1)$ .
在 $\Omega^{\varepsilon}\cap \partial B_R(x_0)$ , $\psi^{\varepsilon}\leq m\leq w$ 和在 $ B_R(x_0)\cap\partial\Omega_{\varepsilon}$ , $\psi^{\varepsilon}=0\leq w$ . 在 $\Omega^{\varepsilon}\cap B_R(x_0)$ 内, 当 $\beta$ 足够小时
$\begin{eqnarray*} \Delta w&=&\frac{C_1m}{R^{1+\beta}}(-A^2+1+\beta)s^{\beta-1}\sin(A\theta+\theta_*)+\frac{C_1m}{R^{1+\beta}}\beta(\beta+1)s^{\beta-1}\sin(A\theta+\theta_*)\\ & \leq&-Cs^{\beta-1} \leq\Delta\psi^{\varepsilon}. \end{eqnarray*}$
通过极大值原理得到, 在 $\Omega(R,\varepsilon)$ 上
(3.35) $\begin{equation} 0 \leq\psi^{\varepsilon}\leq Cs^{1+\beta}, \end{equation}$
其中$C$ 是不依赖于 $\varepsilon$ 的常数. 利用标准的伸缩变换和 Schauder 估计, 得到 $\psi^{\varepsilon}$ 在 $C^{1,\beta}$ 意义下收敛 $\psi$ . 显然, $\psi\in C^{1,\beta}$ 和 $|D\psi|\leq Cs^{\beta}$ . 当角点出现在上边界时, 选择闸函数
(3.36) $\begin{equation} w=m-\frac{C_1m}{R^{1+\beta}}s^{1+\beta}\sin(A\theta+\theta_*), \end{equation}$
4 定理1.2的证明
下面我们将证明解在无穷远处的渐进速率. 这里$\psi$ 是问题1的解,
$\psi_-$ 是问题(3.29)中$a_B$ 的解, 记
$ \Psi=\psi-\psi_-,\quad \Psi_\tau=\psi+\tau\Psi, \quad \tau\in[0,1]. $
直接计算得到, 存在 $K^{'}>K$ , 当 $x\le-K^{'}$ 时, $\Psi$ 满足
(4.1) $\begin{equation} L_0\Psi=\Delta\Psi-\frac{1}{r}\partial_2\Psi-r{\cal C}_{i}\partial_i\Psi-r{\cal D}\Psi=0, \end{equation}$
(4.2) $\begin{equation} {\cal C}_{i}=\widehat c_{i} \left(\frac{\nabla\Psi_\tau}{r},\Psi_\tau,r\right), \end{equation}$
(4.3) $\begin{equation} {\cal D}=\widehat d\left(\frac{\nabla\Psi_\tau}{r},\Psi_\tau,r\right). \end{equation}$
4.1 边界是在有限长度以外为平直的情况
令 $\eta(x,r)={\rm e}^{\mu r+\epsilon x}$ , 其中 $\mu$ 和 $\epsilon$ 是正常数, 对 $\eta$ 作用 $L_0$ 得
(4.4) $\begin{matrix} L_{0}\eta&=&\partial_{11}\eta+\partial_{22}\eta-r{\cal C}_{1}\partial_{1}\eta-\left(\frac{1}{r}+r{\cal C}_{2}\right)\partial_{2}\eta-r{\cal D}\eta\nonumber\\ &=&{\rm e}^{\mu r+\epsilon x}\left(\epsilon^{2}+\mu\left(\mu-\frac{1}{r}\right)-\epsilon c_1-\mu r{\cal C}_{2}-r{\cal D}\right) \end{matrix}$
令 $\mu$ 足够小, 有 $\mu-\frac{1}{r}\leq -\frac{1}{2}$ . 令 $\epsilon$ 足够小, 有 $\epsilon^{2}+\mu\left(\mu-\frac{1}{r}\right)\leq -\frac{\mu}{4}$ . 由 ${\cal D}>0$ , 令 $\mu$ 足够小, 有 $-\mu r{\cal C}_{2}-r{\cal D}\leq \frac{\mu}{8}$ . 于是,
(4.5) $\begin{equation} L_0\eta\leq -\frac{\mu}{8}{\rm e}^{\mu r+\epsilon x}. \end{equation}$
令 $\nu={\rm e}^{\mu r}$ , 同理得
(4.6) $\begin{equation} L_0\nu\leq -\frac{\mu}{8}{\rm e}^{\mu r}. \end{equation}$
令 $K_{2}$ 是一个常数且 $K_{2}>K_{1}$ , 存在一个常数 $ C_0$ 满足
(4.7) $\begin{equation} |\Psi(K_2,r)|\le C_0 {\rm e}^{\epsilon K_2} {\rm e}^{\mu r}. \end{equation}$
对于任意固定的 $\beta>0$ , 存在 $s_0<K_1$ 满足
(4.8) $\begin{equation} |\Psi(x,r)|\le \beta {\rm e}^{\mu r}, \quad x<s_0. \end{equation}$
(4.9) $\begin{equation} \Psi(x,0)=\Psi(x,1)=0, \quad x<K_2. \end{equation}$
(4.10) $\begin{equation} \left\{\begin{array}{ll} L_0\Psi=0>L_0\left(C_0\eta+\beta {\rm e}^{\mu r}\right), \quad& (x,r)\in\Omega_{s,-K_2},\\ \Psi\le C_0\eta+\beta {\rm e}^{\mu r},\quad &(x,r)\in\partial\Omega_{s,-K_2}. \end{array}\right. \end{equation}$
(4.11) $\begin{equation} \Psi\le C_0\eta+\beta {\rm e}^{\mu r}, \ \ (x,r)\in\Omega_{s,-K_2}. \end{equation}$
(4.12) $\begin{equation} \Psi\le C_0\eta+\beta {\rm e}^{\mu r}, \ \ (x,r)\in\Omega_{-\infty,-K_2}. \end{equation}$
(4.13) $\begin{equation} -\Psi\le C_1\eta+\beta {\rm e}^{\mu r}, \ \ (x,r)\in\Omega_{-\infty,-K_2}. \end{equation}$
(4.14) $\begin{equation} |\Psi|\le C_2\eta+\beta {\rm e}^{\mu r}, \ \ (x,r)\in\Omega_{-\infty,-K_2}. \end{equation}$
这里的 $C_2=\max \{ {C_0},{C_1}\} $ . 由 $\beta$ 的任意性得
(4.15) $\begin{equation} |\Psi|\le C_2 {\rm e}^{\epsilon x},\ \ (x,r)\in\Omega_{-\infty,-K_2}. \end{equation}$
(4.16) $\begin{equation} \Vert\Psi\Vert_{C^{2,\alpha}(\Omega_{x-1,x+1})}\le C |\Psi|_{C^{0}(\Omega_{x-1,x+1})}\le C {\rm e}^{\epsilon x}, \ \ \ x\le-K_2, \end{equation}$
(4.17) $\begin{equation} \Vert U-U_0\Vert_{C^{1,\alpha}(\Omega_{x-1,x+1})}+\Vert V\Vert_{C^{1,\alpha}(\Omega_{x-1,x+})}\le C {\rm e}^{\epsilon x},\ \ \ x\le-K_2. \end{equation}$
4.2 边界以代数速率收敛到平直边界的情况
令 $\Upsilon(x,r)=|-x|^{-l}{\rm e}^{\mu r}$ . 其中 $ \mu$ 是正常数. 直接计算得到
$\begin{eqnarray*} \bar{L}\Upsilon&=&\partial_{11}\Upsilon+\partial_{22}\Upsilon-r{\cal C}_{1}\partial_{1}\Upsilon-\left(\frac{1}{r}+r{\cal C}_{2}\right)\partial_{2}\Upsilon-r{\cal D}\Upsilon\nonumber\\ &=&|x|^{-l}{\rm e}^{\mu r}\left(\frac{rl\left(l+1\right)}{x^{2}}+\mu\left(\mu-\frac{1}{r}\right)-\epsilon c_1-\mu r{\cal C}_{2}-r{\cal D}\right). \end{eqnarray*}$
(4.18) $\begin{equation} \bar{L}\Upsilon(x,r)\leq -\frac{\mu}{8}|x|^{-l}{\rm e}^{\mu r}. \end{equation}$
记 $\hat{ f} =\min \{ {1},{f_2(x)}\} $ 和 $\hat{\Omega}=\{(x,r)| 0\le r \le \hat{f}\}$ . 于是, 存在足够大的 $J_2>J_1$ , 当 $x<-J_2$ , 得
(4.19) $\begin{equation} |\Psi|=|\psi-\psi_0|\le|\psi-\frac{m_0}{2\pi}|+|\psi_0-\frac{m_0}{2\pi}|\le(\nabla\psi+\nabla\psi_0)|\hat{f}-1|\le C|x|^{-l}, \quad x\in\partial \hat{\Omega}. \end{equation}$
令 $J_3$ 是一个大于$J_2$ 的常数, 存在常数 $C$ 满足
(4.20) $\begin{equation} |\Psi(J_3, r)|\le C|J_3|^{-l}{\rm e}^{\mu r}. \end{equation}$
对于任意的 $\beta>0$ , 存在 $s_1<-J_3$ 满足
(4.21) $\begin{equation} |\Psi|\le \beta {\rm e}^{\mu r},\quad x<s_1. \end{equation}$
(4.22) $\begin{equation} \left\{\begin{array}{ll} L_0\Psi=0>L_0\left(C\upsilon+\beta {\rm e}^{\mu r}\right), \ \ &(x,r)\in \Omega_{s,-J_3},\\ \Psi\le C\upsilon+\beta {\rm e}^{\mu r}, \ \ &(x,r)\in \partial\Omega_{s,-J_3}. \end{array}\right. \end{equation}$
(4.23) $\begin{equation} |\Psi|\le C\upsilon+\beta {\rm e}^{\mu r}, \ \ (x,r)\in \Omega_{s,-J_3}. \end{equation}$
(4.24) $\begin{equation} |\Psi|\le C|x|^{-l}, \ \ (x,r)\in \Omega_{s,-J_3}. \end{equation}$
(4.25) $\begin{equation} \Vert\Psi\Vert_{C^{2,\alpha}(\Omega_{x-1,x+1})}\le C |\Psi|_{C^{0}(\Omega_{x-1,x+1})}\le C|x|^{-l},\ \ \ x\le-J_3. \end{equation}$
(4.26) $\begin{equation} \Vert U-U_-\Vert_{C^{1,\alpha}(\Omega_{x-1,x+1})}+\Vert V\Vert_{C^{1,\alpha}(\Omega_{x-1,x+1})}\le C|x|^{-l},\ \ \ \ x\le-J_3. \end{equation}$
参考文献
View Option
[1]
Badiani T V . Existence of steady symmetric vortex pairs on a planar domain with an obstacle
Mathematical Proceedings of the Cambridge Philosophical Society , 2000 , 132 (2 ): 365 -384
[本文引用: 1]
[2]
Burton G R . Steady symmetric vortex pairs and rearrangements
Proceedings of the Royal Society of Edinburgh Section a Mathematics , 1988 , 108 (3/4 ): 269 -290
[本文引用: 1]
[3]
Burton G R . Variational problems on classes of rearrangements and multiple configurations for steady vortices
Annales De Linstitut Henri Poincare , 1989 , 6 (4 ): 295 -319
[本文引用: 1]
[4]
Elcrat A R , Miller K G . Steady vortex flows with circulation past asymmetric obstacles
Communications in Partial Differential Equations , 1987 , 12 (10 ): 1095 -1115
DOI:10.1080/03605308708820520
URL
[本文引用: 1]
[5]
Miller E . Rearrangements in steady vortex flows with circulation
Proceedings of the American Mathematical Society , 1991 , 111 (4 ): 1051 -1055
[本文引用: 1]
[6]
Ambrosetti A , Struwe M . Existence of steady vortex rings in an ideal fluid
Archive for Rational Mechanics and Analysis , 1989 , 108 (2 ): 97 -109
DOI:10.1007/BF01053458
URL
[本文引用: 1]
[7]
Cao D M , Liu Z , Wei J . Regularization of point vortices for the Euler equation in dimension two
Archive for Rational Mechanics and Analysis , 2012 , 212 (1 ): 179 -217
DOI:10.1007/s00205-013-0692-y
URL
[本文引用: 1]
[9]
Berger M S , Fraenkel L E . Nonlinear desingularization in certain free-boundary problems
Communications in Mathematical Physics , 1980 , 77 (2 ): 149 -172
DOI:10.1007/BF01982715
URL
[本文引用: 1]
[12]
Bers L . Mathematical Aspects of Subsonic and Transonic Gas Dynamics . New York : John Wiley and Sons , 1958
[本文引用: 1]
[13]
Xie C J , Xin Z P . Global subsonic and subsonic-sonic flows through infinitely long nozzles
Indiana University Mathematics Journal , 2007 , 56 (6 ): 2991 -3023
DOI:10.1512/iumj.2007.56.3108
URL
[本文引用: 1]
[14]
Xie C J , Xin Z P . Global subsonic and subsonic-sonic flows through infinitely long axially symmetric nozzles
Journal of Differential Equations , 2010 , 248 (11 ): 2657 -2683
DOI:10.1016/j.jde.2010.02.007
URL
[本文引用: 1]
[16]
Xie C J , Xin Z P . Existence of global steady subsonic Euler flows through infinitely long nozzles
SIAM Journal on Mathematical Analysis , 2010 , 42 (2 ): 751 -784
DOI:10.1137/09076667X
URL
[本文引用: 1]
[17]
Du L L , Xie C J , Xin Z P . Steady subsonic ideal flows through an infinitely long nozzle with large vorticity
Communications in Mathematical Physics , 2014 , 328 (1 ): 327 -354
DOI:10.1007/s00220-014-1951-y
URL
[本文引用: 2]
[18]
Chen G Q , Huang F M , Wang T Y , et al. Steady Euler flows with large vorticity and characteristic discontinuities in arbitrary infinitely long nozzles
Advances in Mathematics , 2019 , 346 : 946 -1008
DOI:10.1016/j.aim.2019.02.002
URL
[本文引用: 3]
[19]
Du L L , Duan B . Global subsonic Euler flows in an infinitely long axisymmetric nozzle
Journal of Differential Equations , 2011 , 250 (2 ): 813 -847
DOI:10.1016/j.jde.2010.06.005
URL
[本文引用: 1]
[21]
Duan B , Luo Z . Subsonic non-isentropic Euler flows with large vorticity in axisymmetric nozzles
Journal of Mathematical Analysis and Applications , 2015 , 430 (2 ): 1037 -1057
DOI:10.1016/j.jmaa.2015.05.023
URL
[本文引用: 1]
[22]
Duan B , Weng S K . Global smooth axisymmetric subsonic flows with nonzero swirl in an infinitely long axisymmetric nozzle
Ztschrift für angewandte Mathematik und Physik , 2018 , 69 (5 ): 1 -17
[本文引用: 2]
[23]
Deng X M , Wang T Y , Xiang W . Three-dimensional full Euler flows with nontrivial swirl in axisymmetric nozzles
SIAM Journal on Mathematical Analysis , 2018 , 50 (3 ): 2740 -2772
DOI:10.1137/16M1107991
URL
[本文引用: 4]
[24]
Shiffman M . On the existence of subsonic flows of a compressible fluid
Journal Rational Mechanics Analysis , 1952 , 1 : 605 -652
[本文引用: 1]
[25]
Bers L . Existence and uniqueness of a subsonic flow past a given profile
Communications on Pure and Applied Mathematics , 1954 , 7 (3 ): 441 -504
DOI:10.1002/cpa.3160070303
URL
[本文引用: 1]
[26]
Bers L . Results and conjectures in the mathematical theory of subsonic and transonic gas flows
Communications on Pure and Applied Mathematics , 1954 , 7 (1 ): 79 -104
DOI:10.1002/cpa.3160070107
URL
[本文引用: 1]
[27]
Finn R , Gilbarg D . Three-dimensional subsonic flows and asymptotic estimates for elliptic partial differential equations
Acta Mathematica , 1957 , 98 (1 ): 265 -296
DOI:10.1007/BF02404476
URL
[本文引用: 1]
[28]
董光昌 . 空间亚音速流及此边值问题在更高维情况的推广
浙江大学学报 , 1979 , 1 : 33 -63
[本文引用: 1]
Dong G C . Three-dimensinal subsonic ows and their boundary value problems extended to higher dimensions
Journal of Zhejiang University , 1979 , 1 : 33 -63
[本文引用: 1]
[29]
Chen C , Du L L , Xie C J , et al. Two dimensional subsonic Euler flows past a wall or a symmetric body
Archive for Rational Mechanics and Analysis , 2016 , 221 (2 ): 559 -602
DOI:10.1007/s00205-016-0968-0
URL
[本文引用: 1]
[30]
Ou B . An irrotational and incompressible flow around a body in space
Journal of Partial Differential Equations , 1994 , 7 (2 ): 160 -170
[本文引用: 1]
[31]
Ma L , Xie C . Existence and optimal convergence rates of multi-dimensional subsonic potential flows through an infinitely long nozzle with an obstacle inside
Journal of Mathematica Physics , 2020 , 61 (7 ): 1 -23
[本文引用: 1]
[32]
Ma L . The optimal convergence rates of non-isentropic subsonic Euler flows through the infinitely long three-dimensional axisymmetric nozzles
Mathematical Methods in the Applied Sciences , 2020 , 48 (10 ): 6553 -6565
[本文引用: 2]
[33]
Gilbarg D , Trudinger N . Elliptic Partial Differential Equations of Second Order
Berlin: Springer , 1983
[本文引用: 3]
[34]
Chen G Q , Deng X M , Xiang W . Global steady subsonic flows through infinitely long nozzles for the full Euler equations
SIAM Journal on Mathematical Analysis , 2012 , 44 (4 ): 2888 -2919
DOI:10.1137/11085325X
URL
[本文引用: 1]
Existence of steady symmetric vortex pairs on a planar domain with an obstacle
1
2000
... 定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程 ...
Steady symmetric vortex pairs and rearrangements
1
1988
... 定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程 ...
Variational problems on classes of rearrangements and multiple configurations for steady vortices
1
1989
... 定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程 ...
Steady vortex flows with circulation past asymmetric obstacles
1
1987
... 定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程 ...
Rearrangements in steady vortex flows with circulation
1
1991
... 定常不可压 Euler 方程是椭圆-双曲特征耦合的系统, 已有一系列的研究工作. 当密度 $\rho$ 为常数时,方程(1.1)为齐次定常 Euler 方程. 对于齐次定常 Euler 方程的研究主要有两种方法: 涡方法和流函数方法. 涡方法是通过涡/旋函数研究速度场的存在唯一性和相关性质. 利用涡方法研究二维齐次不可压定常 Euler 方程的结果, 可参见文献[1 ⇓ ⇓ ⇓ -5 ]. 流函数方法是由不可压方程引入的流函数, 将齐次定常 Euler 方程转化为关于流函数的椭圆方程. 对于二维的情况, 流函数 $\psi$ 满足椭圆方程 ...
Existence of steady vortex rings in an ideal fluid
1
1989
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
Regularization of point vortices for the Euler equation in dimension two
1
2012
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
Planar vortex patch problem in incompressible steady flow
1
2015
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
Nonlinear desingularization in certain free-boundary problems
1
1980
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
An elliptic problem related to planar vortex pairs
1
2005
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
On the existence of global vortex rings
1
1980
... 这里的非线性项 $f$ 是 $\psi$ 的函数, 表示流场的旋度, 相关结果见文献[6 ⇓ ⇓ ⇓ ⇓ -11 ]. 对于三维柱对称非齐次 Euler方程, 经过柱坐标变换可知其旋度由 $\partial_{r}U-\partial_{x}V$ 和 $W$ 两个量表示, 而二维的旋度由单量表示, 因此在旋度方面, 三维柱对称和二维有本质的区别. 对于三维柱对称非齐次的 Euler 方程, 流函数 $\psi$ 满足椭圆方程的非线性项 $f=f\left(\nabla \psi, \psi, r\right)$ 比二维齐次 Euler 方程的非线性项 $f=f\left(\psi\right)$ 更加复杂, 具体推导请见第 2 节. 因此, 三维柱对称非齐次 Euler 方程不能简单的套用二维的方法. ...
1
1958
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Global subsonic and subsonic-sonic flows through infinitely long nozzles
1
2007
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Global subsonic and subsonic-sonic flows through infinitely long axially symmetric nozzles
1
2010
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Subsonic flows in a multi-dimensional nozzle
1
2011
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Existence of global steady subsonic Euler flows through infinitely long nozzles
1
2010
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Steady subsonic ideal flows through an infinitely long nozzle with large vorticity
2
2014
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... 在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
Steady Euler flows with large vorticity and characteristic discontinuities in arbitrary infinitely long nozzles
3
2019
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... [18 ]和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... 此类关于 $W(x, r)$ 的间断为一类特殊的接触间断, 可参考文献[18 ]. ...
Global subsonic Euler flows in an infinitely long axisymmetric nozzle
1
2011
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Three-dimensional full Euler flows in axisymmetric nozzles
1
2013
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Subsonic non-isentropic Euler flows with large vorticity in axisymmetric nozzles
1
2015
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
Global smooth axisymmetric subsonic flows with nonzero swirl in an infinitely long axisymmetric nozzle
2
2018
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... 在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
Three-dimensional full Euler flows with nontrivial swirl in axisymmetric nozzles
4
2018
... 本文研究的是由 Bers[12 ] 提出的管道问题. 此类问题对应于物理中的风洞实验、实际工程中的油气传输管道等实际场景. 对于可压缩无旋流的研究已经有了大量结果, 可参见文献[13 ⇓ -15 ]. 相比于无旋流的情况, 有旋流的研究更加复杂. 当旋度满足小性条件时, Xie 和 Xin[16 ] 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于大旋度情况, Du, Xie 和 Xin[17 ] 在水平速度的凸性条件下, 证明了二维有旋可压 Euler 流亚音速解的适定性. Chen, Huang, Wang 和 Xiang[18 ] 去除凸性条件和小性要求, 证明了二维有旋可压 Euler 流亚音速解的适定性. 对于三维柱对称流, Du 和 Duan[19 ] 在凸性条件下, 证明了等熵可压 Euler 流亚音速解的适定性, 其后 Du 和 Luo [20 -21 ] 在凸性假设下, 证明了非等熵可压 Euler 流亚音速解的适定性. 当漩涡速度 $W$ 不等于 0 时, Duan 和 Weng[22 ] 证明了三维柱对称可压 Euler 流具有小漩涡的亚音速解的适定性. Deng, Wang 和 Xiang[23 ] 证明了三维柱对称可压 Euler 流具有大漩涡解的适定性. 通过补偿紧性框架, 由绝热指数极限, 即当绝热指数趋向于无穷大时, 定常可压缩 Euler 流体收敛到定常不可压缩 Euler 流体, Chen, Huang, Wang 和 Xiang[18 ] 和 Deng, Wang 和 Xiang[23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... [23 ] 分别得到二维和三维柱对称不可压缩 Euler 方程极限解的存在唯一性. ...
... 在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
... 进一步的由极大值原理, 类似于文献[23 ], 我们可以得到在 $\Omega$ 上 $U(x, r)>0$ . 由此对于 $f_1$ 和 $f_2$ 不具有角点的情况证明了定理1. ...
On the existence of subsonic flows of a compressible fluid
1
1952
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
Existence and uniqueness of a subsonic flow past a given profile
1
1954
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
Results and conjectures in the mathematical theory of subsonic and transonic gas flows
1
1954
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
Three-dimensional subsonic flows and asymptotic estimates for elliptic partial differential equations
1
1957
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
空间亚音速流及此边值问题在更高维情况的推广
1
1979
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
空间亚音速流及此边值问题在更高维情况的推广
1
1979
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
Two dimensional subsonic Euler flows past a wall or a symmetric body
1
2016
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
An irrotational and incompressible flow around a body in space
1
1994
... 定常 Euler 流研究中的另一个重要问题是研究流体流经障碍物的问题. 无旋流的情况, Shiffman 和 Bers[24 ⇓ -26 ] 证明二维可压 Euler 流亚音速解的适定性. 对于三维(或更高维)的情况, Finn 和 Gilabarg[27 ] 证明了当流速低于某一速度时, 三维无旋 Euler 流亚音速解的适定性和渐近性行为. 董光昌[28 ] 将 Finn 等人的结果扩展到空间维数大于等于三的一般情况. 对于有旋光滑解适定性问题的研究成果较少. Chen, Du, Chen 和 Xin[29 ] 对于二维对称翼形面问题, 证明了有旋亚音速解的存在性和唯一性. 对于无旋不可压缩 Euler 方程, 二维的情况可由复变函数的经典理论讨论, 空间维数大于等于三的情况由 Ou[30 ] 通过变分法证明了解的存在性和唯一性. ...
Existence and optimal convergence rates of multi-dimensional subsonic potential flows through an infinitely long nozzle with an obstacle inside
1
2020
... 在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
The optimal convergence rates of non-isentropic subsonic Euler flows through the infinitely long three-dimensional axisymmetric nozzles
2
2020
... 在适定性的基础上, 对于无界区域, 一个自然的问题就是研究解在无穷远处的渐近收敛速率. Du, Xie 和 Xin[17 ] 研究了二维可压 Euler 方程管道问题有旋解的无穷远的渐近收敛速率. Ma 和 Xie[31 ] 考虑了高维 ($n\geq 3$ ) 可压 Euler 方程管道问题无旋解的无穷远的渐近收敛速率. Ma[32 ] 在漩涡速度 $W=0$ 的情况下得到了三维柱对称可压 Euler 方程有旋解的无穷远的渐近收敛速率. 值得指出的是: Ma[32 ] 的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
... [32 ]的结果在文献[22 -23 ]的条件的基础上需要两个额外的条件. ...
Elliptic Partial Differential Equations of Second Order
3
1983
... 这里$\partial\Omega_{L}$ 是 $C^{2,\alpha_1}$ 边界, $0<\alpha_1<\alpha<1$ , 并满足一致外球条件. 对于固定的 $k>0$ 和 $L>0$ , 由文献[33 ]知,问题(3.6)存在唯一解 $\psi^{(k)}_L\in C^{2,H}(\Omega_L)\cap C^0(\bar{\Omega}_L)$ . ...
... 利用文献[33 ,定理 3.7] 得 ...
... 利用文献[33 ,引理 6.20], 存在 $C$ 满足 ...
Global steady subsonic flows through infinitely long nozzles for the full Euler equations
1
2012
... 由于 $\partial_2\bar{\psi}_L>0$ , 选择足够小的 $\delta$ , 得到 $\partial_2\psi_L>0$ . 假设存在最大的有界正数 $\delta_0$ , 使得当 $\delta\leq\delta_0$ 时, $\partial_2\psi_L>0$ ; 当 $\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])}= \delta_0$ 时, 在某些点上 $\partial_2\psi_L\leq 0$ . 那么通过选取一个序列 $({\cal G}_n, {\cal B}_n, {\cal W}_n)$ , 其中对于每一个 $n\in \mathbb{N}$ , $\|({\cal G}'_n,{\cal B}'_n,{\cal W}'_n)\|_{C^{0,1}([m])}\leq \delta_0$ , 当 $n\rightarrow \infty$ 时, $({\cal G}_n, {\cal B}_n, {\cal W}_n)\rightarrow ({\cal G}, {\cal B}, {\cal W})$ , 其中 $\|({\cal G}',{\cal B}',{\cal W}')\|_{C^{0,1}([m])}=\delta_0$ . 记 $\psi_{L}^{(n)}$ 和 $\psi_L$ 分别是 $({\cal G}_n, {\cal B}_n, {\cal W}_n)$ 和 $({\cal G}, {\cal B}, {\cal W})$ 对应问题(3.6)的唯一解. 由有界区域解的唯一性和 $C^{2,\alpha}$ 估计, 当 $n\rightarrow\infty$ 得到 $\psi_{L}^{(n)}\rightarrow\psi_L$ . 由 $\partial_2 \psi_{L}^{(n)}>0$ , 可知 $\partial_2 \psi_{L}\geq0$ . 由文献[34 ]可知如果 $\partial_2 \psi_{L}\geq0$ , 那么 $\partial_2 \psi_{L}>0$ , 与 $\delta_0$ 的定义矛盾. 因此, 不存在有界的 $\delta_0$ . 从而, 对于任意的 $\delta>0$ 都有 $\partial_r \psi_{L}>0$ . ...