Bernoulli Equation

news/2025/3/17 16:34:49/文章来源:https://www.cnblogs.com/jxsme/p/18777110

伯努利方程

伯努利方程并非一个独立的定律,而是在不同条件下由 Navier-Stokes 动量方程(式1)和能量方程(式2)推导而来的。

\[\rho\left(\frac{\partial u_{j}}{\partial t} + u_{i}\frac{\partial u_{j}}{\partial x_{i}}\right)=-\frac{\partial p}{\partial x_{j}}+\rho g_{j}+\frac{\partial}{\partial x_{i}}\left[\mu\left(\frac{\partial u_{j}}{\partial x_{i}}+\frac{\partial u_{i}}{\partial x_{j}}\right)+\left(\mu_{v}-\frac{2}{3}\mu\right)\frac{\partial u_{m}}{\partial x_{m}}\delta_{ij}\right] \tag{1} \]

\[\rho\frac{De}{Dt}=-p\frac{\partial u_{m}}{\partial x_{m}} + 2\mu\left(S_{ij}-\frac{1}{3}\frac{\partial u_{m}}{\partial x_{m}}\delta_{ij}\right)^{2}+\mu_{v}\left(\frac{\partial u_{m}}{\partial x_{m}}\right)^{2}+\frac{\partial}{\partial x_{i}}\left(k\frac{\partial T}{\partial x_{i}}\right) \tag{2} \]

考虑无粘流动时\(\mu = \mu_v = 0\),重力是唯一的体积力,因此式1就简化为欧拉方程

\[\frac{\partial u_{j}}{\partial t}+u_{i}\frac{\partial u_{j}}{\partial x_{i}}=-\frac{1}{\rho}\frac{\partial p}{\partial x_{j}}-\frac{\partial \varPhi}{\partial x_{j}} \tag{3} \]

式中:\(\varPhi = gz\) 表示体积势函数;\(g\)为重力加速度;z轴为垂直方向,如果流动是正压流动,那么\(\rho = \rho(p)\),并且:

\[\frac{1}{\rho}\frac{\partial p}{\partial x_j}=\frac{\partial }{\partial x_j}\int_{p_0}^{p}\frac{dp'}{\rho(p')} \tag{4} \]

式中,\(dp/\rho\)是全微分,\(p_0\)是参考压力,\(p'\)是积分变量。在上式中,积分置于端点有关,而与积分路径无关。恒密度、等温以及等熵流动都是正压流动(barotropic)。除此之外式3中对流加速度可以改写成速度-涡量的向量积和单位质量动能梯度的形式

\[u_{i}\frac{\partial u_{j}}{\partial x_{i}}=-(\mathbf{u} \times \boldsymbol{\omega})_{j}+\frac{\partial}{\partial x_{j}}\left(\frac{1}{2}u_{i}^{2}\right) \tag{5} \]

将式(5)和式(4)带入到式(3)中:

\[\frac{\partial u_{j}}{\partial t}+\frac{\partial }{\partial x_{j}}\left[\frac{1}{2}u_{i}^{2}+\int_{p_0}^{p}\frac{dp'}{\rho(p')}+gz\right]=(\mathbf{u} \times \boldsymbol{\omega})_{j} \tag{6} \]

所有的梯度项都被合并成了伯努利函数B,即中括号中的内容。式(6)可以用于推导无粘正压流动伯努利方程的演化。

首先考虑定常流动(\(\partial u_j/\partial t = 0\),这样式(6)简化成

\[\nabla B=\mathbf{u}\times\boldsymbol{\omega} \tag{7} \]

左端项是正交于常数曲面B的向量,右端项是同时垂直于\(\boldsymbol{u}\)\(\boldsymbol{\omega}\)的向量。如图所示,
image

这是一个由流线和涡线定义的曲面,在该曲面上,伯努利函数是一个常量,注意流线和涡线可以为任一方位。

由此可知B为常数的曲面必定包含流线和涡线,因此,无粘、定常、正压流动满足:

\[\frac{1}{2}u_{i}^{2}+\int_{p_0}^{p}\frac{dp'}{\rho(p')}+gz = 沿流线和涡线为常数值 \tag{8} \]

这是伯努利几种可能方程的第一种形式,另外,如果流动无旋(\omega = 0) 的,那么式(7)就有

\[\frac{1}{2}u_{i}^{2}+\int_{p_0}^{p}\frac{dp'}{\rho(p')}+gz=\text{constant everywhere} \tag{9} \]

可以证明,存在包含流线和涡线的曲面的充分条件是流动为正压流动,顺便说一下,这种类型的曲面称为 Lamb surfaces ,是为了纪念杰出的英国应用数学家和流体力学家霍勒斯・兰姆。

在一般的非正压流动中,可以在流场中的任意两点之间绘制由流线和涡线段组成的路径。那么式(8)是有效的,条件是在所选的特定路径上计算积分。如前文所述,(式(8)要求流动是定常的、无粘性的,并且只有重力(或其他保守)体积力作用于其上。

我们在此仅指出一个要点:在非旋转参考系中,如果粘性效应可以忽略不计,正压无旋流动将保持无旋。考虑绕固体物体的流动,比如翼型(如图所示)。
image

在靠近物体表面的薄粘性层(称为边界层)之外的所有点处,流动都是无旋的。这是因为在粘性层之外的流线上的粒子 P 从某点 S 出发,在 S 点处流动是均匀的,因此是无旋的。所以,在这个例子中,伯努利方程式(9)在粘性层之外的所有地方都成立。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/900435.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

社区演讲-基于.NET 技术栈的研发过程管理和智能化探索

Hi all 2024年11月23日,我作为【项目管理实践探索者大会】专题讲师,给社区分享了主题《基于.NET 技术栈的研发过程管理和智能化探索》 现场参与人数:100+ 演讲议题介绍: 体系建设:如何设计一个合适的研发过程管理体系 平台提效:通过.NET技术构建任务驱动+流程驱动,释放研…

halcon 深度学习教程(二)halcon基于深度学习的OCR检测

原文作者:aircraft 原文链接:https://www.cnblogs.com/DOMLX/p/18777081深度学习教程目录如下,还在继续更新完善中 深度学习系列教程目录 本篇非常简单,但是非常好用的OCR识别,传统的OCR识别方法一大堆各种各样的操作,麻烦的要死。halcon 22版本开始支持深度学习OCR,并且…

[ SpringMVC ] SpringMVC如何通过是否有@RestController注解来判断返回ModelAndView还是Json

引言之前在面试的遇到面试官问我SpringMVC的执行流程,我那时候回答的是SpringMVC的DispatcherServlet的dodispatch方法找到ControllerMethod之后将返回值通过convert成Json返回响应体,事后想了一下回答的其实并不正确,因为SpringMVC之前学习的时候有使用ModelAndView返回视图,我…

从HR+AI到AI+HR,企业人力资源AI进程已过半

一、人力资源管理智能化应势而上,核心价值受企业管理层肯定 过往各项研究表明,AI 已经被广泛应用于企业经营的各个环节中。根据易路于2023 年发布的《AI 在企业人力资源中的应用白皮书1.0》(以下简称《白皮书1.0》),我们可以明确感受到:AI 已不同程度应用于招聘管理、员工…

multi-object tracking in the dark

创新点:构建LMOT数据集:开发了一种双摄像头系统,同步采集明暗视频帧,形成高度对齐的低光视频对,并提供高质量的多目标跟踪标注。该数据集包含大量城市户外场景视频,涵盖多种动态物体,为模型训练和评估提供了丰富的数据支持。 提出LTrack方法:引入自适应低通下采样模块(…

LGP11831_1 [UPTS 2025] 追忆 学习笔记

LGP11831_1 [UPTS 2025] 追忆 学习笔记 Luogu Link 前言 又幻想了……唉! 幻想自己场切这道题,最后标准分上升至 \(\text{598pts}\),翻掉了 \(\text{yyz}\),不至于一点脸不要。 本题解基本借鉴这篇题解。 题意简述 给定一个 \(n\) 点 \(m\) 边的简单有向图 \(G\),有 \(m\)…

sql 在两个数据表中,A表存在字段以逗号分隔存储B表的多id对象,进行关联查询

A 表:B表:关联查询 需求为,查询出A表的数据列表,需要将A表关联B表的数据id,概要通过B表的 name 进行输出显示 SELECT A.id,A.name,A.creator,A.created_at,GROUP_CONCAT(B.name SEPARATOR , ) AS B_names -- 将 c_name 合并为逗号分隔的字符串 FROM xf_service_type A LE…

算力市场何以拥有巨大潜力

算力市场未来确实具有巨大的潜力,这一判断基于多个方面的因素: 一、算力成为经济增长的主要驱动力 随着全球数字化转型的迅猛推进,算力已成为推动经济增长的关键引擎。各国纷纷加大在算力基础设施方面的投入,以期为经济发展注入新的活力。例如,欧盟委员会批准了一项名为“…

算法备案拟公示内容编写指南

除了自评估报告,算法备案复审中的拟公示内容也是难度颇大的一份材料,导致很多开发者的算法备案申请被驳回。今天我就提供一份简易模板供大家学习参考(请结合实际情况撰写,不要照抄,不要买模板,否则会判定真实性存疑或高度雷同,影响备案)。(各类文件套模板都会判定该真…

CH585 RF_Basic例程讲解含单向和双向发送

CH585_RF基础通讯例程见下图路径:1、RF初始化参数配置/******************************************************************************** @fn RFRole_Init** @brief RF应用层初始化** @param None.** @return None.*/ void RFRole_Init(void) {rfTaskID = TMOS_…

uniapp整合SQLite(Android)

一、勾选SQLite数据库选项 (1)HBuilder工具打开项目 (2)项目/manifest.json =>App模块配置 => 勾选SQLite(数据库)二、封装sqlite.ts 在项目根目录下创建sqlite/sqlite.ts// 数据库名称 const dbName = scan/*** 数据库地址* @type {String} 推荐以下划线为开头 _d…

No.68 Vue---vue3新特性

一、vue3新特性 1.1 六大亮点二、组合API(setup)2.1 ref或者reactive 1、创建项目 vue create vue-demo5 2、进入文件,启动服务。 3、 2.2 methods中定义的方法写在setup() 2.3setup()中使用props和context 在2.x中,组件的方法中可以通过this获取到当前组件的实例,并执…