<p>We establish optimal error bounds on an exponential wave integrator (EWI) for the space fractional nonlinear Schrödinger equation (SFNLSE) with low regularity potential and/or nonlinearity. For the semi-discretization in time, under the assumption of <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(L^\infty \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>L</mi> <mi>∞</mi> </msup> </math></EquationSource> </InlineEquation>-potential, <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(C^1\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>C</mi> <mn>1</mn> </msup> </math></EquationSource> </InlineEquation>-nonlinearity, and <InlineEquation ID="IEq3"> <EquationSource Format="TEX">\(H^\alpha \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>H</mi> <mi>α</mi> </msup> </math></EquationSource> </InlineEquation>-solution with <InlineEquation ID="IEq4"> <EquationSource Format="TEX">\(1&lt;\alpha \le 2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>1</mn> <mo>&lt;</mo> <mi>α</mi> <mo>≤</mo> <mn>2</mn> </mrow> </math></EquationSource> </InlineEquation> being the fractional index of <InlineEquation ID="IEq5"> <EquationSource Format="TEX">\((-\Delta )^\frac{\alpha }{2}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mrow> <mo stretchy="false">(</mo> <mo>-</mo> <mi mathvariant="normal">Δ</mi> <mo stretchy="false">)</mo> </mrow> <mfrac> <mi>α</mi> <mn>2</mn> </mfrac> </msup> </math></EquationSource> </InlineEquation>, we prove an optimal first-order <InlineEquation ID="IEq6"> <EquationSource Format="TEX">\(L^2\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>L</mi> <mn>2</mn> </msup> </math></EquationSource> </InlineEquation>-norm error bound <InlineEquation ID="IEq7"> <EquationSource Format="TEX">\(O(\tau )\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mi>τ</mi> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> and a uniform <InlineEquation ID="IEq8"> <EquationSource Format="TEX">\(H^\alpha \)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>H</mi> <mi>α</mi> </msup> </math></EquationSource> </InlineEquation>-norm bound of the semi-discrete numerical solution, where <InlineEquation ID="IEq9"> <EquationSource Format="TEX">\(\tau \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>τ</mi> </math></EquationSource> </InlineEquation> is the time step size. We further discretize the EWI in space by the Fourier spectral method and obtain an optimal error bound in <InlineEquation ID="IEq10"> <EquationSource Format="TEX">\(L^{2}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>L</mi> <mn>2</mn> </msup> </math></EquationSource> </InlineEquation>-norm <InlineEquation ID="IEq11"> <EquationSource Format="TEX">\(O(\tau +h^{m})\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mi>τ</mi> <mo>+</mo> <msup> <mi>h</mi> <mi>m</mi> </msup> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> without introducing any CFL-type time step size restrictions, where <i>h</i> is the spatial step size, <i>m</i> is the regularity of the exact solution. Moreover, under slightly stronger regularity assumptions, we obtain optimal error bounds <InlineEquation ID="IEq12"> <EquationSource Format="TEX">\(O(\tau )\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mi>τ</mi> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq13"> <EquationSource Format="TEX">\(O(\tau +h^{m-{\frac{\alpha }{2}}})\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mi>τ</mi> <mo>+</mo> <msup> <mi>h</mi> <mrow> <mi>m</mi> <mo>-</mo> <mfrac> <mi>α</mi> <mn>2</mn> </mfrac> </mrow> </msup> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> in <InlineEquation ID="IEq14"> <EquationSource Format="TEX">\(H^\frac{\alpha }{2}\)</EquationSource> <EquationSource Format="MATHML"><math> <msup> <mi>H</mi> <mfrac> <mi>α</mi> <mn>2</mn> </mfrac> </msup> </math></EquationSource> </InlineEquation>-norm, which is the norm associated to the energy. Extensive numerical examples are provided to validate the optimal error bounds and show their sharpness. We also find distinct evolving patterns between the SFNLSE and the classical nonlinear Schrödinger equation.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Optimal error bounds on an exponential wave integrator Fourier spectral method for fractional nonlinear Schrödinger equations with low regularity potential and nonlinearity

  • Junqing Jia,
  • Xiaoyun Jiang

摘要

We establish optimal error bounds on an exponential wave integrator (EWI) for the space fractional nonlinear Schrödinger equation (SFNLSE) with low regularity potential and/or nonlinearity. For the semi-discretization in time, under the assumption of \(L^\infty \) L -potential, \(C^1\) C 1 -nonlinearity, and \(H^\alpha \) H α -solution with \(1<\alpha \le 2\) 1 < α 2 being the fractional index of \((-\Delta )^\frac{\alpha }{2}\) ( - Δ ) α 2 , we prove an optimal first-order \(L^2\) L 2 -norm error bound \(O(\tau )\) O ( τ ) and a uniform \(H^\alpha \) H α -norm bound of the semi-discrete numerical solution, where \(\tau \) τ is the time step size. We further discretize the EWI in space by the Fourier spectral method and obtain an optimal error bound in \(L^{2}\) L 2 -norm \(O(\tau +h^{m})\) O ( τ + h m ) without introducing any CFL-type time step size restrictions, where h is the spatial step size, m is the regularity of the exact solution. Moreover, under slightly stronger regularity assumptions, we obtain optimal error bounds \(O(\tau )\) O ( τ ) and \(O(\tau +h^{m-{\frac{\alpha }{2}}})\) O ( τ + h m - α 2 ) in \(H^\frac{\alpha }{2}\) H α 2 -norm, which is the norm associated to the energy. Extensive numerical examples are provided to validate the optimal error bounds and show their sharpness. We also find distinct evolving patterns between the SFNLSE and the classical nonlinear Schrödinger equation.