节点不连续伽辽金方法在求解线性和非线性平流方程中的一维实现(Matlab代码实现)

💥💥💞💞欢迎来到本博客❤️❤️💥💥

****🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。

⛳️**座右铭:**行百里者,半于九十。

📋📋📋++本文目录如下:++🎁🎁🎁

目录

[💥1 概述](#💥1 概述)

[📚2 运行结果](#📚2 运行结果)

[🎉3 参考文献](#🎉3 参考文献)

[🌈4 Matlab代码实现](#🌈4 Matlab代码实现)


💥1 概述

本文提出了节点不连续伽辽金方法在求解线性和非线性平流方程中的一维实现,无需任何滤波器或限幅器。本公式旨在向CFD从业者介绍该方法,因此它意味着可读性而不是非常有效的实施。

该方法的基本思想是将平流方程离散化为一维节点不连续伽辽金方程,然后使用迭代方法求解该方程。

在一维情况下,平流方程可以写为:

∂u/∂t + c ∂u/∂x = 0

其中,u是要求解的变量,t是时间,x是空间坐标,c是平流速度。

为了将平流方程离散化为节点不连续伽辽金方程,我们首先将时间和空间离散化为有限的网格点。然后,我们使用中心差分方法来近似平流方程中的导数项。具体来说,我们使用以下公式来近似时间导数和空间导数:

∂u/∂t ≈ (u_i^n+1 - u_i^n) / Δt

∂u/∂x ≈ (u_i+1^n - u_i-1^n) / (2Δx)

其中,u_i^n表示在第n个时间步长和第i个空间节点处的变量值,Δt是时间步长,Δx是空间步长。

将这些近似代入平流方程,我们得到节点不连续伽辽金方程:

(u_i^n+1 - u_i^n) / Δt + c (u_i+1^n - u_i-1^n) / (2Δx) = 0

为了求解该方程,我们使用迭代方法,如迭代法或牛顿法。迭代方法的基本思想是将方程转化为一个求根问题,并通过迭代逼近根的解。

总之,节点不连续伽辽金方法是一种求解线性和非线性平流方程的一维实现方法,它不需要滤波器或限幅器。然而,由于其可读性和实施的简单性,它可能不是非常高效的方法。

📚 2 运行结果

部分代码:

% Build 1d mesh

xgrid = mesh1d([0 2*pi],nE,'LGL',K);

dx = xgrid.elementSize; J = xgrid.Jacobian;

x = xgrid.nodeCoordinates; w = xgrid.weights';

xc = xgrid.elementCenter;

% Load DG tools

tool = DGtools(xgrid.solutionPoints);

V = tool.Vandermonde2;

invM = tool.nodalInvMassMatrix;

Dr = tool.nodalCoefDiffMatrix;

% Build Lift Operator

Emat = zeros(K+1,2); % array of element's shape function

Emat(1,1)=1; Emat(K+1,2)=1;

Lift = V*(V'*Emat);

% IC

u0 = IC(x,2);

% Set plot range

plotrange = [xgrid.range(1),xgrid.range(2),...

min(min(min(0.9*u0)),min(min(1.1*u0))),1.1*max(max(u0))];

%% Solver Loop

% Set initial time & load IC

t=0; u=u0; it=0;

% Using a 3rd Order 3-stage SSPRK time integration

while t < tEnd

uo = u;

% update time

dt = cfl*dx/max(max(abs(dflux(u)))); t = t+dt;

% iteration counter

it = it+1;

% 1st stage

dF = residual(u,flux,dflux,Lift,Dr);

u = uo-dt*dF/J;

% 2nd Stage

dF = residual(u,flux,dflux,Lift,Dr);

u = 0.75*uo+0.25*(u-dt*dF/J);

% 3rd stage

dF = residual(u,flux,dflux,Lift,Dr);

u = (uo+2*(u-dt*dF/J))/3;

% build cell averages

u_bar = w*u/2;

% Plot u

subplot(1,2,1); plot(x,u,x,u0,'-+'); axis(plotrange); grid off;

subplot(1,2,2); plot(xc,u_bar,'ro'); axis(plotrange); grid off;

%if rem(it,10) == 0

drawnow;

%end

end

%% Final Plot for IC 2

subplot(1,2,1); plot(x,u,x,u0,'-+'); axis(plotrange);

title('Nodal DG','interpreter','latex','FontSize',18);

xlabel('\\it{x}','interpreter','latex','FontSize',14);

ylabel({'\\it{u(x)}'},'interpreter','latex','FontSize',14);

subplot(1,2,2); plot(x,u0,'k-',xc,u_bar,'ro'); axis(plotrange);

title('Cell Averages','interpreter','latex','FontSize',18);

xlabel('\\it{x}','interpreter','latex','FontSize',14);

ylabel({'\\it{u(x)}'},'interpreter','latex','FontSize',14);

🎉3 参考文献

文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。

1\]陈长胜,纪立人,陈嘉滨,等.JFNK方法在求解全隐式一维非线性平流方程中的应用\[J\].大气科学, 2007, 31(5):10.DOI:10.3878/j.issn.1006-9895.2007.05.19. \[2\]赵晓旭.勒让德-伽辽金法求解两类非线性积分微分方程\[D\].哈尔滨师范大学\[2023-08-11\]. \[3\]王俊.时域不连续伽辽金法在计算电磁学中的应用\[D\].西安电子科技大学,2015.DOI:10.7666/d.D01068543. ## [🌈](https://mp.weixin.qq.com/mp/appmsgalbum?__biz=Mzk0MDMzNzYwOA==&action=getalbum&album_id=2591810113208958977#wechat_redirect "🌈")****4 Matlab代码实现****

相关推荐
艾莉丝努力练剑9 分钟前
【LeetCode&数据结构】单链表的应用——反转链表问题、链表的中间节点问题详解
c语言·开发语言·数据结构·学习·算法·leetcode·链表
_殊途2 小时前
《Java HashMap底层原理全解析(源码+性能+面试)》
java·数据结构·算法
橡晟4 小时前
深度学习入门:让神经网络变得“深不可测“⚡(二)
人工智能·python·深度学习·机器学习·计算机视觉
墨尘游子4 小时前
神经网络的层与块
人工智能·python·深度学习·机器学习
Leah01054 小时前
机器学习、深度学习、神经网络之间的关系
深度学习·神经网络·机器学习·ai
珊瑚里的鱼5 小时前
LeetCode 692题解 | 前K个高频单词
开发语言·c++·算法·leetcode·职场和发展·学习方法
xingshanchang5 小时前
Matlab的命令行窗口内容的记录-利用diary记录日志/保存命令窗口输出
开发语言·matlab
秋说6 小时前
【PTA数据结构 | C语言版】顺序队列的3个操作
c语言·数据结构·算法
lifallen6 小时前
Kafka 时间轮深度解析:如何O(1)处理定时任务
java·数据结构·分布式·后端·算法·kafka
liupenglove7 小时前
自动驾驶数据仓库:时间片合并算法。
大数据·数据仓库·算法·elasticsearch·自动驾驶