对流方程的差分法求解-间断面的展示以及拟线性方程的求解

已于 2022-11-28 19:14:09 修改
阅读量633 收藏 1
点赞数
分类专栏： 传统数值方法 文章标签： 矩阵 线性代数
于 2020-12-30 22:05:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/111408657
版权
传统数值方法 专栏收录该内容
35 篇文章 20 订阅 ¥199.90 ¥299.90
订阅专栏
本文探讨了线性对流方程的差分法求解，包括一阶和二阶精度的数值格式，如迎风格式、LF格式、LW格式和BW格式，并通过具体算例展示了数值解的计算。同时，文章详细阐述了边界条件的处理方法，特别是针对不同初值条件的处理策略。此外，还讨论了拟线性方程的求解及其解的间断点展示和数值误差分析。
摘要由CSDN通过智能技术生成

可以通过本人知乎主页了解博客内容添加链接描述

三个线性对流方程

{ u t + u x = 0 , x

了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
专栏目录
订阅专栏
使用jasmin求解线性对流方程
syh704364538的博客
09-22 455
线性对流方程： ∂u(x⃗,t)∂t+∇⋅(a⃗)⋅u(x⃗,t)=0,x⃗∈ω\frac{\partial u(\vec{x} ,t)}{\partial t} + \nabla \cdot (\vec{a}) \cdot u(\vec{x}, t) = 0, \quad \vec{x} \in \omega∂t∂u(x,t)​+∇⋅(a)⋅u(x,t)=0,x∈ω ...
一维稳态对流扩散问题，无源项，一阶迎风差分格式的python程序
Fo*(Bi)的博客
11-15 2337
一维稳态对流扩散无源项一阶迎风差分格式 参考书籍：陶文铨的数值传热学+李人宪的有限体积法基础 例题 一个长度L为1的规则物体，左边界温度恒为1，右边界温度恒为0，流速u=2，扩散系数Γ为0.1，密度 ρ为1。 本题划分5个网格。 代码： #一维稳态对流扩散问题 #无源项，一阶迎风差分格式 from numpy import * from numpy.linalg import * import matplotlib.pyplot as plt L = 1 #物体总长度 n = 5 #划分网格数 I
【有限差分法】（三）一维和二维抛物方程CN格式以及长时间稳定性分析（附算例与Python代码）
qq_39482771的博客
11-28 8758
前言 我们知道差分里的CN格式是无条件稳定的。但是最近在学习有限元结合CN格式算长时间抛物问题的时候（在时间方向用差分空间方向用有限元）发现稳定性却不能保证，其数值解和真解误差会随着时间越来越大。并且实际能算的时间只有几秒，这几秒钟几乎是没有实际意义的。所以就想到一个问题，纯CN差分格式在长时间计算的时候，是否也是个理论看还行，实践臭弟弟的“花瓶”呢？ 在以前做差分课本上的算例大多也是只有几秒，还真没注意到长时间稳定性的问题(¬_¬ )。所以本文找了最简单的一维抛物问题，把时间层增加至1000S，来看看C.
关于CN离散格式
armstrongd的博客
10-01 2492
CN离散之后是五对角线矩阵 要解五对角线矩阵，就要把右边的b写进来，写进来之后，要把 初始值矩阵是已知的，所以直接代入进行迭代求解，期间，把每次的边界条件都进行修正 难点是如何把矩阵转换为一个n-1平方和一列的一个向量。...
CN格式matlab程序,征集MATLAB程序开发说明较好的格式
weixin_33583401的博客
03-18 188
%SIZE Size of array.% D = SIZE(X), for M-by-N matrix X, returns the two-element row vector% D = [M,N] containing the number of rows and columns in the matrix.% For N-D arrays, SIZE(X) returns ...
convection-diffusion2_对流扩散方程有限差分法求解_
10-03
在计算机科学和数值分析领域，对流扩散方程是描述许多...总之，通过对流扩散方程的有限差分法求解，我们可以模拟和预测各种物理现象，而"convection-diffusion2.f90"文件则提供了实现这一过程的具体算法和代码框架。
一维线性对流有限差分：使用五种有限差分方案的选择来求解一维线性对流方程。-matlab开发
06-01
一维线性平流方程使用五种有限差分方案（全部显式）来求解。 First Order Upwind、Lax-Friedrichs、Lax-Wendroff、Adams Average (Lax-Friedrichs) 和 Adams Average (Lax-Wendroff)。 使用启发式时间步长。 使用...
MATLAB求解扩散方程有限差分法 源程序代码.zip_一维 扩散_一维扩散_交叉扩散matlab_扩散方程_有限差分
07-14
MATLAB是进行数值计算的强大工具，其内置的矩阵运算功能非常适合求解线性方程组，这是有限差分法求解偏微分方程的关键步骤。源程序代码中可能包含以下部分： 1. 定义网格：创建空间和时间的离散网格。 2. 初始化...
一维对流方程编程求解_一维对流方程编程求解_源码
09-30
本项目提供了中心差分、向前差分和向后差分三种数值方法来近似求解该方程，下面将详细介绍这些方法及其应用。 一维对流方程通常写作： \[ \frac{\partial u}{\partial t} + c \frac{\partial u}{\partial x} = 0 \]...
【有限差分法】（一）有限差分法的基本流程与常用格式
热门推荐
qq_39482771的博客
11-25 2万+
2.有限差分法的引入与介绍 有限差分法：通过将求解区域离散化得到离散的微分算子，将微分方程定解问题转化为代数方程组求解问题。 有限差分法的原理：我们在用数值方法求解微分方程时，由于计算机不能处理连续的数学方程，所以在用计算机处理问题时，必须先将一个无限的连续问题通过一些特殊的方法离散化，使之成有限形式的线性代数方程组，再将这个线性代数方程组交给计算机求解以达到我们的目的。在偏微分方程中，通过离散化区域获得导数（或偏导数）的手段就构成了我们的有限差分法。 3.有限差分法的流程 以一维情况为例来说明一个有限
抛物线差分格式求解matlab程序
01-26
抛物线差分格式求解，包括一维古典显格式，DFF格式，CN格式，局部一维方法和预测校正格式的求解；文件中有具体的题目及解法说明并还有matlab程序供大家参考；易读易懂。
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数值解矩阵 输出图像： 精确解和数值解的图像
matlab解对流方程
04-09
用拉克斯-温德洛夫多步格式解对流方程 matlab源程序
非线性发展方程定解问题
图灵的猫的博客
09-30 847
非线性发展方程定解问题 文章目录非线性发展方程定解问题前言Burgers方程正则长波方程Korteweg-de Vries方程Camassa-Holm (C-H)方程Schrödinger 方程Kuramoto-Tsuzuki 方程ZakharOv方程Ginzburg-Landau 方程Cahn-Hilliard 方程外延增长模型方程相场晶体模型方程 前言 非线性现象的研究是自然科学领域甚至社会科学领域也十分关心的问题. 由于自然界中许多的现象本质上是非线性的, 所以非线性现象引起了工程师、物理 学家、数学
数学物理方程学习笔记5
qq_51557623的博客
11-11 2066
数学物理方程学习笔记5 2.热传导方程 2.2Cauchy问题的解 Cauchy问题： n=1 {ut−a2uxx=0,x∈R,t>0u(x,0)=φ(x),x∈R \begin{cases} u_t-a^2u_{xx}=0,x \in R,t>0\\ u(x,0)=\varphi(x),x \in R \end{cases} {ut​−a2uxx​=0,x∈R,t>0u(x,0)=φ(x),x∈R​ 方法： Flourier变换法 step1：对方程和初始条件关于x进行Flourier变
【Python-CFD】02：线性对流方程
weixin_43863651的博客
02-27 2390
线性对流方程 差分与扩散 一维常系数对流方程是CFD中最简单最基本的方程。然而这个简单的方程能够告诉我们很多关于CFD方面的东西。模型方程如下式： 下面在Python中实现该求解过程。 代码如下所示： import numpy as np import matplotlib.pyplot as plt nx = 41 #x方向网格节点数量 dx = 2/(nx-1) #空间网格步长，x方向...
一阶电路误差分析_PDE有限差分方法(12)——对流方程数值格式的分析方法
weixin_39900180的博客
10-27 845
今天直接把张强书上的第6章大部分内容完成掉. 主要讲对流方程数值格式的相容性、稳定性以及数值振荡的分析方法。参考书：(1) J.W. Thomas - Numerical Partial Differential Equations_ Finite Difference Methods (1995, Springer)(2) 张强《偏微分方程的有限差分方法》科学出版社，2019年1月版。(3) K...
几种简单差分格式
lk040384的博客
05-05 1万+
以对流方程为例： 一阶迎风 （时间向前空间向前） （时间向前空间向后） 二阶迎风 三阶迎风 二阶中心 （无条件不稳定） 将上式的取成平均值，就得到了Lax-Friedrichs格式 在时间和空间上都用二阶中心差分离散，可以得到蛙跳（leap-frog）格式 将在进行Taylor级数展开 根据对流方程有 代入Taylor展开式得 用中...
一维常系数对流方程的学习——来自流沙公众号
代码
09-25 7975
一维对流扩散方程的学习——来自流沙公众号 认识： 1、时间步长越小越靠后移动；网格越小则波的形状越一致,波形失真在减小，引出Courant数：u*dt/dx<sigma； 2、对numpy中的ones（）、linspace（）、zeros（）有基础认识； 3、再次练习了matplotlib中的pyplot。 4、明白了一维方程编程的简单运算，基础 import numpy as np # n...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

Galerkin码农选手

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

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

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

打赏作者

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

抵扣说明：

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

余额充值