Hermite WENO 重构格式

最新推荐文章于 2022-11-26 22:01:56 发布
最新推荐文章于 2022-11-26 22:01:56 发布
阅读量703 收藏 1
点赞数
原文链接：http://www.cnblogs.com/li12242/p/5821922.html
版权

Hermite WENO 单元重构

本文主要介绍采用 Hermite WENO 重构方法作为斜率限制器应用于二维或高维单元中。

1.简介1

ENO格式最早由 Harten 等2提出，ENO格式避免在高梯度地区进行插值，其重构过程通过多个不同位置模板，并且选取其中最光滑模板上的解进行重构，由此保证在重构过程中具有较高的精度。然而这种方法会导致收敛问题并且在连续区域降低精度，随后Liu等3提出了 weighted ENO 格式。WENO 格式采用不同的权重系数组合各个模板，而非只采用一个判断光滑程度的系数。Qiu 和 Shu 首先将 WENO 格式作为一种斜率限制器应用于DG方法中4，随后他们又根据 Hermite 多项式构造了一系列的 WENO 格式，称为 HWENO 格式并也将其作为限制器应用于 DG 方法中5

最低0.47元/天 解锁文章
R芮R
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
三次Hermite插值
littlely_ll的博客
05-03 4万+
设f(x)f(x)在节点a≤x0,x1,⋯,xn≤ba\le x_0, x_1,\cdots,x_n\le b处的函数值为f0,f1,...,fnf_0,f_1,...,f_n，设P(x)为f(x)P(x)为f(x)在区间[a,b][a,b]上的具有一阶导数的插值函数 （1）若要求P(x)P(x)在[a,b][a,b]上具有一阶导数（一阶光滑度） P(xi)=f(xi)=fiP′(xi)=f′(
weno格式源代码
10-06
weno代码
weno格式构造流程
06-04
高精度TVD格式的够早的一般流程，给出了一个简单的算例，欢迎大家参考借鉴学习
阅读笔记 > 重构的格式
孤独的梦1012
10-14 232
重构的格式每个重构手法一般都有如下五个部分：首先是名称（name) 建造一个重构词汇表，名称是很重要的。名称之后是一个简短概要（summary） 简单介绍此重构手法的适用情景，以及它所做的事情。这部分可以帮助你更快找到你所需要的重构手法。动机（motivation） 为你介绍「为什么需要这个重构」和「什么情况下不该使用这个重构」。作法（mechanics） 简明扼要地一步一步介绍如何进行此一重构。范
WENO格式自动推导
李文韬的博客
11-26 1536
WENO格式是CFD中的一种高精度的数值格式。如果函数光滑，使用r个模板可以在空间上达到2r−1。如果出现间断，那么WENO格式退化为ENO格式。理论上WENO可达任意阶精度，但是推导过程比较繁琐。本文使用Mathematica软件完成WENO格式的自动推导，可以实现任意阶精度。
一维sod激波管问题weno格式
02-17
一维sod激波管问题weno格式
五阶weno格式的m文件
06-07
五阶精度weno格式matlab代码 时间推进三阶rk格式 求解激波稀疏波问题
分段Hermite插值推导
weixin_44378835的博客
11-30 2117
分段Hermite插值 分段线性插值多项式S(x)S(x)S(x)在插值区间[a,b][a,b][a,b]上只能保证连续性，而不光滑。要想得到在插值区间上光滑的分段线性插值多项式，可采用分段埃尔米特（Hermite）插值，这里我们考虑在整个[a,b][a,b][a,b]上用分段三次埃尔米特插值多项式来逼近f(x)f(x)f(x)。一般的将带有导数的插值多项式称为Hermite插值多项式。 如果已知函数y=f(x)y = f(x)y=f(x)在节点a=x0<x1<…<xn=ba = x_0&
Hermite 插值多项式——JAVA实现
qq_43799743的博客
04-02 1739
Hermite 插值多项式计算近似值——JAVA实现 import java.util.Scanner; public class Hermite { //x x //y y //m m // x0 插入的值 // y0 近似值 public static double[] aaa(double x[],double y[], double m[], double x0[]) { ...
双曲守恒律的ENO格式和WENO格式.zip
12-14
双曲守恒律的ENO格式和WENO格式的MATLAB代码，求解的是Burgers方程的初值问题，分别使用有限体积法4阶ENO格式，有限体积法3阶和5阶WENO方法求解。时间方向用三阶TVD Runge-Kutta方法。研究格式在解光滑情况下和有间断的情况下的数值精度并作图。说明文档参考博客https://blog.csdn.net/lusongno1/article/details/79070049。
快速Hermite径向基函数曲面重构
02-09
在这篇论文中,我们提出了一个使用局部支撑径向基函数对三维散乱点进行Hermite插值或逼近的快速曲面重构方法.通过构造给定数据点集的一个层次结构,采用逐层精化的方式实现了全局曲面重构的效果,解决了因使用局部支撑...
Hermite_Hermite曲线_
09-29
hermite曲线绘制，hermite曲线的变化是受四个基函数控制的，相当于四个基函数的加权后累加。
MFC实现Hermite曲线
12-01
*MFC实现Hermite曲线，可以运行,有说明文档。
Hermite函数
03-09
定义Hermite插值函数的matlab代码，Hermite插值，matlab。
18.Hermite插值
Sweeping_Robot的博客
06-13 209
hermite插值matlab代码,hermite插值matlab
weixin_36358109的博客
03-18 854
数值计算方法与MATLAB 包能胜 汕头大学工学院 第四讲 数据插值与MATLAB 多项式与插值 ? 多项式插值的主要目的是用一个多项式 拟合离散点上的函数值,使得可以用该 ......多项式、插值与数据拟合 ? 多项式MATLAB命令 ? 插值– –––– Lagrange插值 Hermite插值 Runge现象和分段插值 分段插值 样条插值的MATLAB表示 ? 数据拟合–多项式......M...
三次hermite插值多项式例题_hermite插值例题
weixin_36403128的博客
01-11 6288
证毕。对于齐次问题,可知节 x 点 例题 讨论一个节点处的 Hermite 插值问题? 解 假定在一个节点 x 处满足 k+1个插值条件: 0 p (j) ( x ) ? c (0 ......关键词:龙格现象 分段差值 三次 Hermite 进行插值 1、实验目的 1) 通过对分段三次 Hermite 插值算法程序的编写,提高自己编写程序的 能力 2) 体会分段三次 ......(x x1)( ...
分段三次Hermite插值
04-09
分段三次Hermite插值
HWENO fortran 代码
最新发布
10-14
这是一个简单的HWENO（High-Order Weighted Essentially Non-Oscillatory）Fortran代码示例： ``` subroutine HWENO(u, uHWENO, N, dx) implicit none integer, intent(in) :: N real, intent(in) :: dx real, dimension(0:N-1), intent(in) :: u real, dimension(0:N-1), intent(out) :: uHWENO real, dimension(5) :: alpha real, dimension(5) :: w real, dimension(5) :: uL real, dimension(5) :: uR real :: eps = 1e-6 integer :: i do i = 2, N-3 alpha = [1/10.0, 6/10.0, 3/10.0, 0.0, 0.0] w = alpha / (eps + (u(i-2:i+2)-u(i))**2) w = w / sum(w) uL = [u(i-2), u(i-1), u(i), u(i+1), u(i+2)] uR = [u(i-1), u(i), u(i+1), u(i+2), u(i+3)] uHWENO(i) = w(1)*(-uL(1)+5*uL(2)+2*uL(3))/6.0 & + w(2)*(2*uL(2)+5*uL(3)-uL(4))/6.0 & + w(3)*(2*uL(3)+5*uL(2)-uL(1))/6.0 & + w(4)*(11*uL(3)-7*uL(2)+2*uL(1))/6.0 & + w(5)*(11*uL(3)-13*uL(2)+5*uL(1))/6.0 end do uHWENO(0) = u(0) uHWENO(1) = u(1) uHWENO(N-2) = u(N-2) uHWENO(N-1) = u(N-1) end subroutine HWENO ``` 其中，`u`是输入的原始数据，`uHWENO`是输出的经过HWENO处理后的数据，`N`是数据点数，`dx`是网格间距。该代码使用了五阶HWENO格式来计算每个内部网格点的值，边界点直接复制原始数据的值。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值