一阶双曲型偏微分方程的数值解法——迎风格式、Lax-Friedrichs格式、Lax-Wendroff格式和Beam-Warming格式等

已于 2022-01-26 14:51:56 修改
阅读量1.5w 收藏 125
点赞数 14
分类专栏： 微分方程数值解 文章标签： 线性代数 偏微分方程 矩阵 python
于 2022-01-12 13:39:31 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/liuqihang11/article/details/122451568
版权

声明

本部分是一个学习笔记，主要内容来自于华冬英老师编写的《微分方程的数值解法与程序实践》。如果觉得内容不错，可自行购买价格良心的官方正版教材。http://www.hxedu.com.cn.上有配套的代码以PPT课件可供免费下载。另外，官方代码均用C语言编写，之后我也会陆陆续续上传自己编写的Python代码。

本部分内容公式太多实在懒得敲了，因此以图片形式呈现，基本能满足学习要求。不过也同时上传了文字可复制的PDF版本，内容与图片版本完全一致。感兴趣的话可在https://download.csdn.net/download/liuqihang11/77977916下载，不过需要付费1.99元，请按需下载。

syphomn
关注
  • 14
    点赞
  • 踩
  • 125
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
专栏目录
偏微分方程算法之一阶双曲差分法
欢迎各位大佬光临“寒舍”！请多多指点，共同进步！
04-18 681
主要介绍一阶双曲型偏微分方程（以一阶对流方程为例）的多种差分格式算法，并通过C++实现算例。
偏微分方程数值解迎风格式代码
04-29
差分格式是数值计算方法中微分以及偏微分导数的一种离散化方法，即用相邻两个或者多个数值点的差分取代偏微分方程中导数或者偏导数的一种算法。 选择差分格式是离散化偏微分方程的第一步。本代码是迎风格式的代码。具体参加李荣华的数值解书目。
UPW_utux0.rar_matlab 迎风格式_偏微分方程_迎风格式_迎风格式 matlab_迎风格式matlab
07-14
function [ue,un]=UPW_utux0(v,dt,t) 一个简单的双曲型偏微分方程： ut + ux = 0 初始条件为： u(x,0) = 1, x≤0 0, x＞0. 边界条件为： u(-1,t)=1,u(1,t)=0. 本题要求： 使用迎风格式，选择 v=0.5, 计算并画出当dt=0.01和0.0025时， 方程在t=0.5，x在(-1,1)时的数值解和精确解 输入： v--即a*dt/dx dt--数值格式的时间步 t--要求解的时间 输出： ue--在时间t时的1×N精确解矩阵 un--在时间t时的1×N数值解矩阵 输出图像： 精确解和数值解的图像
双曲型偏微分方程的数值解法（3）
01-26
4.1 二阶双曲型方程的紧差分法 4.2 二维双曲型方程的ADI格式 4.3 二维双曲型方程的紧ADI格式
双曲型偏微分方程的数值解法（1）
01-26
双曲型偏微分方程的初值依赖特性以及波传导特性 迎风格式 Leap-Frog Scheme格式 Lax-Friedrichs 格式 Lax-Wendroff 格式 Beam-Warming格式 隐格式设计
双曲型偏微分方程数值求解MATLAB程序
02-25
【内容介绍】本资源主要利用MATLAB的实时脚本编程实现了双曲型偏微分方程数值求解，以图-文-代码三者互相嵌套的形式介绍实现过程，一目了然。对数值解和解析解进行作图对比分析。 【适用对象】工科生、数学专业等。 【算法涵盖】迎风格式、Lax-Friedrichs 格式、Lax-Wendroff 格式 感谢支持！
间断初值双曲守恒问题的Lax-Friedrichs和后向欧拉数值解法
LSEC小陆的博客
12-04 3555
间断初值双曲守恒问题的Lax-Friedrichs和后向欧拉数值解法我们用求解PDE常用的Backwand Euler和Lax-Friedrichs数值方法来求解如下双曲问题:首先回顾一下双曲问题中一些经典的有限差分格式：- 该黎曼问题的精确解为：计算精确解Matlab代码：function u = ExSolu( x, t ) %计算精确解，这里x支持向量输入，输出多个精确解 a = 1;%设置a
常用守恒格式求解双曲问题（Lax-Friedrichts、local Lax-F、Roe、Engquist-Osher、Godunov、Lax-Wendroff）
LSEC小陆的博客
12-09 1万+
Burgers方程的守恒格式求解（Proj No.1）陆嵩 计算数学问题描述我们对下述Burgers方程的初边值问题 ⎧⎩⎨⎪⎪⎪⎪ut+(u2/2)x=0u(x,0)=u0(x)=sin(πx+π)periodicB.C.forx−1≤x≤1,t>0\label{q} \left \{ \begin{aligned} & u_t + (u^2/2)_x = 0 & -1\leq x
CFD有限体积5-对流项离散
qq_41817886的博客
02-22 1530
FVM，对流项离散
3.2.1 对流项离散格式的基本介绍 | 3.2.2 中心离散格式与一阶迎风格式（OpenFOAM理论笔记系列）
CloudBird07的博客
08-28 5138
3.2 对流项的离散格式 3.2.1 对流项离散格式的基本介绍 在开始本节的讨论开始，笔者首先要说明，本节所介绍的内容相对于整个对流项离散格式的开发的历史可以说是极其简略的。可以这么说，整个计算流体动力学发展的历史有一半以上是由对流项离散格式写成的。本节的内容以Versteeg和Malalasekera的著作《An Introduction to Computational Fluid Dynamic: The Finite Volume Method Second edition》中的对流项离散格式介绍内
一维激波管问题——计算流体力学（Lax-Friedrichs格式）
05-12
本文用python3.6对一维激波管（lax激波管）问题进行实现，用到的格式为Lax-Friedrichs格式，激波管问题的初值采用了lax激波管的初值（具体见代码），推进时间为0.14s，网格数为1000
Lax Wendoff 1D Burgers：从 Taylor-Galerkin 过程派生的 La​​x Wendroff 方案-matlab开发
05-29
一维 Burgers 方程的 Taylor-Galerkin 离散化的简单实现，当单元尺寸恒定时，它简化为 Lax-Wendroff 方案。 包括派生说明。 对于实际使用，运行粗网格/时间步长组合，并基于 max(abs(u))，使用更细的网格和 CFL = 0.8 重新估计 dt。 尽管更高的值在技术上更准确，但实际上虚假振荡（双曲线守恒的 LW 方案的典型特征）会对稳定性施加额外的、难以预测的惩罚。 “真实”粘性项的实现有助于控制振荡，但对网格/时间步长大小比率（B = u_max*((dt^2)/(dx^2)) ）施加了另一个限制。 请注意，即使在不连续点附近使用振荡解决方案，它仍然比一阶逆风方法更能捕获冲击。
一维激波管问题——计算流体力学（Lax-Wendroff格式）
05-12
本文用python3.6对一维激波管lax激波管问题进行了计算，使用了Lax-Wendoff格式，激波管问题的初值采用了lax激波管的初值（具体见代码），推进时间为0.14s，网格数为500，相较于Lax-Friedrichs格式在稀疏波与激波处更加符合解析解, 接触间断面处会发生耗散。
Lax-wendroff_Lax-Wendroff_
09-30
使用Lax-wendroff格式解方程 Ut+3Ux=0
双曲型方程LaxWendroff的差分格式程序(Matlab).rar_LaxWendroff_Matla代码_Wendroff
07-14
双曲型方程Lax-Wendroff的差分格式程序(Matlab)，非常好的例子，具体算例，有理论分析，解题过程分析，Matlab代码，注释清晰明了，有计算结果输出和图形显示，按照Word文档说明，很快就能学会，操作简单明了，非常...
论文研究 - 一阶双曲型方程的一维显式Tolesa数值格式及其在宏观交通流模型中的应用
05-22
提出了一种求解一阶双曲型偏微分方程的数值方案，并将其应用于具有恒定速度和线速度-密度关系的宏观交通流模型的仿真研究中。 宏观交通流模型由Lighthill Whitham and Richards（LWR）首先开发，用于通过诸如流量，...
一维激波管（Lax shock tube）问题的数值求解
热门推荐
LSEC小陆的博客
05-28 1万+
一维激波管（Lax shock tube）问题的数值求解 文章目录一维激波管（Lax shock tube）问题的数值求解@[toc] 问题描述一般的守恒算法格式精确解常用守恒格式（先介绍标量）Lax-Friedrichs格式Roe格式Lax-Wendroff格式ENO格式简介（先介绍标量）5阶WENO格式空间离散Lax-Friedrichs通量分裂原函数方法重构左值右值空间算子和时间三阶Rung...
一阶双曲型偏微分方程matlab代码
最新发布
05-13
一阶双曲型偏微分方程一般的形式为： $$ \frac{\partial u}{\partial t}+a\frac{\partial u}{\partial x}=0 $$ 其中 $a$ 为常数。 可以用有限差分法来数值求解这个方程，其中 $u_{i,j}$ 表示在位置 $x_i$ 和时间 $t_j$ 处的解。 我们可以选择用向前差分、向后差分或中心差分来离散化偏微分方程。下面以中心差分法为例： $$ \frac{u_{i,j+1}-u_{i,j}}{\Delta t}+a\frac{u_{i+1,j}-u_{i-1,j}}{2\Delta x}=0 $$ 整理得到： $$ u_{i,j+1}=u_{i,j}-\frac{a\Delta t}{2\Delta x}(u_{i+1,j}-u_{i-1,j}) $$ 根据时空离散化的方法，可以用以下 MATLAB 代码实现一阶双曲型偏微分方程的求解： ```matlab % 离散化参数 Nx = 100; % 空间离散化步数 Nt = 200; % 时间离散化步数 a = 1; % 常数 a % 区间参数 x_start = 0; x_end = 1; t_start = 0; t_end = 1; % 离散化步长 dx = (x_end - x_start) / Nx; dt = (t_end - t_start) / Nt; % 初始条件 u0 = sin(pi * linspace(x_start, x_end, Nx+1)); % 数值求解 u = u0; for j = 1:Nt u_new = u; for i = 2:Nx u_new(i) = u(i) - a * dt / (2 * dx) * (u(i+1) - u(i-1)); end u = u_new; end % 可视化 figure(); plot(linspace(x_start, x_end, Nx+1), u0, 'r--', 'LineWidth', 1.5); hold on; plot(linspace(x_start, x_end, Nx+1), u, 'b-', 'LineWidth', 1.5); legend('t=0', 't=1'); xlabel('x'); ylabel('u'); title(sprintf("1-order hyperbolic PDE, dx=%.2f, dt=%.2f, a=%.2f", dx, dt, a)); ``` 这个代码用的是中心差分法，实现了一阶双曲型偏微分方程的离散化求解，并把结果可视化出来。其中，离散化步数和常数 $a$ 都可以根据具体问题进行调整。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

syphomn

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值