船舶水动力数值计算中的有限体积法
cfd有限体积法
cfd有限体积法CFD有限体积法CFD（Computational Fluid Dynamics）是指利用计算机模拟流体运动的科学技术。
而有限体积法（FVM，Finite Volume Method）是CFD中的一种数值方法，它将流域分割成许多小的控制体积，然后通过对每个控制体积内的物理量进行离散化，将偏微分方程转化为代数方程组，从而求解出流场的各个物理量。
1. FVM基本原理1.1 控制体积FVM方法将流域分割成许多小的控制体积，每个控制体积都是一个封闭区域。
在这个区域内，可以计算出各种物理量（如密度、速度、压力等），并且这些物理量在整个区域内都是均匀的。
1.2 通量通量是指单位时间内通过单位面积所传递的某种物理量。
在FVM中，通量是一个重要的概念。
通过对每个控制体积进行质量守恒和动量守恒方程进行离散化，可以得到通量在各个边界上的表达式。
1.3 离散化离散化是将偏微分方程转化为代数方程组的过程。
在FVM中，通过对控制体积内的物理量进行离散化，可以得到每个控制体积内的物理量与相邻控制体积内的物理量之间的关系式。
1.4 数值求解离散化后，可以得到代数方程组。
通过数值方法（如迭代法、高斯消元法等），可以求解出这个方程组，并得到流场各个物理量的数值解。
2. FVM优点2.1 适用性广FVM方法适用于各种复杂流动问题，如湍流、多相流、非牛顿流等。
2.2 精度高FVM方法是一种高精度的数值方法，能够准确地计算出流场各个物理量的分布情况。
2.3 稳定性好FVM方法具有良好的稳定性和收敛性，在计算过程中不会出现发散等问题。
3. FVM应用领域3.1 航空航天工业在航空航天工业中，FVM方法被广泛应用于飞行器气动力学、燃烧室燃烧过程模拟、液体火箭发动机喷注等领域。
3.2 汽车工业在汽车工业中，FVM方法被用于模拟气动力学、燃烧过程、发动机燃料喷射等问题。
3.3 能源领域在能源领域中，FVM方法被用于模拟火电厂锅炉内的流动和传热过程、风力发电机叶片的气动特性等问题。
船舶水动力学特性研究
船舶水动力学特性研究船舶作为人们日常生活和经济活动中的不可或缺的交通工具，其水动力学特性的研究对于提升航行安全性和经济效益至关重要。
因此，本文将从船舶水动力学研究的概念、基本原理、影响因素等方面进行探讨，希望能够对船舶水动力学研究有一定的了解。
一、船舶水动力学研究概念船舶水动力学研究是指对船舶在水中运动时的艏波、尾流、水阻、稳定性、操纵性等各种现象进行研究和分析的学科。
其研究对象是各种类型和尺寸的船舶，包括商船、军舰、游艇等。
二、船舶水动力学研究基本原理1. 流态分析原理流态分析原理是船舶水动力学研究的核心。
该原理是将船舶在水中运动时所遇到的现象分别按其流态性质分类，以便于研究和分析船舶的水动力学特性。
根据流态分析原理，船舶的水动力学特性主要包括以下三个方面：(1) 湍流：指水流具有乱流、波浪等不规则的流动状态。
船舶在湍流中行驶时，会产生较大的阻力和动荡现象，对航行安全和经济效益均有不利影响。
(2) 节流：指水流具有规则但收缩的流动状态。
船舶在节流中行驶时，会感受到明显的水阻，对船速和经济性能产生较大影响。
(3) 局部分离：指水流对船身表面产生局部脱离的流动状态。
局部分离会对船舶的水动力学性能产生很大影响，如增大阻力和降低操纵性能等。
2. 流场数值模拟原理流场数值模拟原理是指通过计算机对船舶在水中运动时产生的复杂流场进行模拟和分析，以便于更加精确地研究船舶的水动力学特性。
该原理主要采用有限体积法、有限元法、边界元法等方法进行计算。
流场数值模拟可以预测船舶在水中的运动轨迹、速度、阻力等物理量，为航行安全和经济效益的提升提供了重要依据。
三、船舶水动力学特性影响因素船舶水动力学特性的研究不仅需要掌握其基本原理，还需要考虑各种影响因素。
以下是船舶水动力学特性研究中的一些重要影响因素。
1. 船体形状船体形状是影响船舶水动力学性能的最重要因素之一。
船体的前后部形状、槽宽、吃水深度以及外型等参数对船舶的水阻、波浪、稳定性等性能均有影响。
有限体积法介绍
有限体积法介绍有限体积法1 有限体积法基本原理上⼀章讲到的有限差分法将数值⽹格的节点上定义为计算节点，并在⽹格节点上对微分形式的流体基本⽅程进⾏离散，⽤⽹格节点上的物理量的代数⽅程作为原PDE 的近似。
在本章所要学习的有限体积法则采⽤了不同的离散形式。
⾸先，有限体积法离散的是积分形式的流体⼒学基本⽅程：d q ds ds SSΩΩ+??Γ=?φφρφn n v(1)计算域⽤数值⽹格划分成若⼲⼩控制体。
和有限差分法不同的是，有限体积法的⽹格定义了控制体的边界，⽽不是计算节点。
有限体积法的计算节点定义在⼩控制体内部。
⼀般有限体积法的计算节点有两种定义⽅法，⼀种是将⽹格节点定义在控制体的中⼼，另⼀种⽅法中，相邻两个控制体的计算节点到公共边界的距离相等。
第⼀种⽅法的优点在于⽤计算节点的值作为控制体上物理量的平均值具有⼆阶的精度；第⼆种⽅法的好处是在控制体边界上的中⼼差分格式具有较⾼的精度。
积分形式的守恒⽅程在⼩控制体和计算域上都是成⽴的。
为了获得每⼀个控制体上的代数⽅程，⾯积分和体积分需要⽤求⾯积公式来近似。
2 ⾯积分的近似采⽤结构化⽹格，在⼆维情况下，每⼀个控制体有4个⾯，⼆维情况，每⼀个控制体有6个表⾯。
计算节点⽤⼤写字母表⽰，控制体边界和节点⽤⼩写字母表⽰。
为了保证守恒性，控制体不能重叠，每⼀个⾯都是相邻两个控制体的唯⼀公共边界。
控制体边界上的积分等于控制体个表⾯的积分的和：∑??=kkfds fdS(2)上式中，f 可以表⽰n u ρφ或nΓφ。
显然，为了获得边界上的积分，必须知道f 在边界上的详细分布情况，这是不可能实现的，由于只是计算节点上的函数值，因此必须采⽤近似的⽅法来计算积分。
整个近似过程分成两步第⼀步：⽤边界上⼏个点的近似积分公式第⼆步：边界点上的函数值⽤计算节点函数值的插值函数近似⾯积分可采⽤以下不同精度的积分公式：⼆阶精度积分：e e e e S e Sf S f fds F e≈==?(3)上式中e f 为边界中点出的函数值。
volume of fluid 有限体积法
volume of fluid 有限体积法
Volume of Fluid (VOF) method 和有限体积法（Finite Volume Method）都是计算流体力学中的数值方法，用于模拟和分析流体流动。
Volume of Fluid (VOF) method 是一种界面捕捉方法，利用流体体积函数处理界面破碎、融合以及大变形等问题。
这种方法通过求解体积分数的输运方程，实现多相流动界面形状及演化的计算。
体积分数的空间分布隐含着界面的位置和形状，通过求解体积分数的输运方程，可以计算多相流动界面形状及演化的计算。
有限体积法（Finite Volume Method）是一种常用的数值算法，着重从物理观点来构造离散方程。
每一个离散方程都是有限大小体积上某种物理量守恒的表示式，推导过程物理概念清晰，离散方程系数具有一定的物理意义，并可保证离散方程具有守恒特性。
这种方法将计算区域划分为一系列不重复的控制体积，每一个控制体积都有一个节点作代表，将待求的守恒型微分方程在任一控制体积及一定时间间隔内对空间与时间作积分。
总之，这两种方法都是计算流体力学中常用的数值方法，用于模拟和分析流体流动。
VOF方法更适合处理界面捕捉问题，而有限体积法更适合处理物理量守恒的问题。
有限容积法和有限体积法
有限容积法和有限体积法有限容积法和有限体积法是计算流体力学中常用的两种数值方法，它们在流体动力学的数值计算中占有非常重要的地位。
本文将从概念、原理、特点、应用等方面，对这两种方法进行详细介绍。
一、有限容积法1.概念有限容积法（Finite Volume Method，FVM）是一种离散化的数值方法，它将连续的物理量离散化为有限个体积元，在每个体积元内计算其平均值，进而求解整个流体系统的物理量。
FVM方法的核心是质量守恒原理，即物质的进出必须平衡，这种保证了物理量在每个体积元内的守恒关系，从而保证了数值计算的准确性。
2.原理FVM方法的数值计算是基于网格的，它将流体动力学问题离散化为一个由有限体积元组成的系统，将原问题转化为流量守恒方程的求解，即$$\frac{\Delta m}{\Delta t}=\Sigma_{faces}\rho uA$$其中，$\Delta m$是在$\Delta t$时间内通过一个表面的质量变化量，$\rho$是介质的密度，$u$是速度，$A$是面积。
对于每个有限体积元，上式可以写为其中，$F_{ij}^p$和$F_{ij}^n$分别是流向有限体积元内部和外部的通量，$i,j$是有限体积元的编号。
3.特点（1）FVM方法基于质量守恒原理，具有非常强的数值稳定性和保真性；（2）FVM方法的计算结果具有局部守恒性，能够准确反映流场内部的物理现象；（3）FVM方法可以处理非结构化网格，适用范围广泛；（4）FVM方法求解的是面积分，所需的时间和空间存储相对较少。
4.应用（1）流体力学领域，如空气动力学、水力学、燃烧问题等；（2）材料科学领域，如薄膜生长、材料变形等。
有限体积法（Finite Element Method，FEM）是一种离散化的数值方法，它将求解的物理场离散化为有限个单元，然后在每个单元内进行近似计算。
相比于FVM方法，FEM方法更加精确，适用于需要高精度计算的问题。
带附体潜艇尾流场的数值模拟与验证
带附体潜艇尾流场的数值模拟与验证李艳;姚震球【摘要】对美国DARPA潜艇模型SUBOFF尾流场进行了数值模拟,采用湍流模型k-ε、k-ω与壁面函数相结合,直接求解RANS方程,并将计算结果与试验结果进行了比较分析,验证了数值方法的可靠性.在此基础上将SUBOFF模型改变为"X"型尾翼,对其尾流场进行了数值模拟分析.【期刊名称】《江苏科技大学学报（自然科学版）》【年(卷),期】2006(020)002【总页数】6页(P7-12)【关键词】计算流体力学;潜艇;尾流场【作者】李艳;姚震球【作者单位】江苏科技大学,船舶与海洋工程学院,江苏,镇江,212003;江苏科技大学,船舶与海洋工程学院,江苏,镇江,212003【正文语种】中文【中图分类】U661.30 引言含附体潜艇周围粘性流场特别是尾流场不仅是其水动力性能的基本反映,而且是艇体水动力噪声的主要来源,因此潜艇尾流场流动特性的研究直接关系到新型安静型潜艇的设计。
一直以来模型试验是研究船舶流场的重要方法。
随着计算机技术的迅猛发展,数值计算方法逐渐成为研究船舶流场的另一有效方法。
对于潜艇流场的数值计算,国内外水动力学界正在进行研究。
中国船舶科学研究中心对潜艇流场也进行了大量的试验和计算研究。
赵峰、周连弟[1](1996)采用k-ε湍流模型对潜艇含指挥台附体区域周围粘性流场进行了计算。
张楠[2](2005)采用k-ε、RNG k-ε、k-ω 3种湍流模型预报了美国DARPA潜艇模型SUBOFF与CSSCR潜艇模型SM-x 的阻力与尾流场,并将计算结果与试验结果对比验证。
Huang[3](1992)对SUBOFF 艇型进行了系统的风洞尾流场试验,采用BL模型、修正的BL模型与k-ε模型详细计算了桨盘面处的轴向无量纲速度。
Bull[4](1996)通过变化网格的形式和数量、湍流模型、流动求解器,对SUBOFF AFF-1(裸艇体)与AFF-8(全附体)的尾流场的数值模拟进行了系统的验证。
某型船舶水力学性能计算及优化研究
某型船舶水力学性能计算及优化研究近年来，随着海洋经济的发展和对船舶安全性、运输效率等方面的要求越来越高，对船舶水动力性能的研究也越来越重要。
其中，船舶水力学性能在船舶设计与运行中占有重要地位。
本文旨在探讨某型船舶水力学性能的计算与优化。
一、简介某型船舶是一种大规模的海洋运输器械，主要用于散货、油料等大宗物资运输。
该型船整体长度180米，宽度32米，速度16节，排水量36000吨。
由于船舶特性、运输特性和环保要求等不同因素的影响，某型船舶的水动力学效能需要进行研究和计算。
二、计算方法船舶水动力学性能的计算方法多种多样，甚至有些深奥难懂，但总的来说主要可以分为两大类，即实验研究方法和计算机模拟方法。
实验研究方法主要是通过模型试验和原型试验等手段，测定船舶在不同流速和舵角条件下的阻力、推力、速度、操纵性能等参数，从而获取船舶水动力学特性的基础数据。
但是这种方法成本较高且需要较长时间。
计算机模拟方法可以通过数值模拟手段来计算船舶水动力特性参数，包括速度、流量、粘性系数、推力、抗力等。
主要可以分为数值流体计算（CFD）和势流计算（Potential Flow）等两种方式。
从实用角度出发，本文采用了CFD数值流体计算来进行水动力性能的研究。
三、数值模拟分析采用CFD数值流体计算方法，可以准确计算船舶船体的流速、压力、阻力、推力等水动力学参数，进而分析不同条件下船舶性能的变化规律。
1. 建模与网格化首先需要对某型船舶进行建模，包括将船舶几何形状、细节参数等输入计算软件，并进行网格化处理。
建模过程中需要注意对船舶舰形和细节特征的合理处理，确保计算结果的准确性。
2. 模拟条件设置模拟条件的设置直接关系到数值计算结果的准确性和适用性。
需要确定的条件有：（1）不同流速条件下船舶的阻力和推力分布。
（2）不同舵角条件下船舶的操纵能力。
（3）船舶在波浪中的性能表现。
通过对不同条件的模拟分析，可以了解某型船舶的水动力学性能特点，找到性能瓶颈，最终实现性能的优化。
流体力学中的数值模拟与应用
流体力学中的数值模拟与应用流体力学是一门研究流体运动规律和流体力学性质的科学，包括流体静力学和流体动力学两个方面。
流体力学应用广泛，如气体动力学、水力学、气象学等领域。
本文将着重介绍流体力学中的数值模拟及其应用。
一、数值模拟数值模拟是一种通过计算机模拟流体等物理现象的方法，通常采用计算流体动力学（CFD）方法来实现。
CFD模拟涉及到流体力学、数值分析、计算机科学和编程等多个领域。
1. 基本原理运用数值方法，将物理问题转换为数值问题，对流场进行数值求解。
其中最重要的是Navier-Stokes方程组，它是描述不可压缩流体动力学的基本方程。
Navier-Stokes方程组分为连续性方程和动量方程两个部分。
连续性方程描述了质量守恒，动量方程描述了动量守恒。
对于不同的物理问题，还需要加入适当的边界条件和物理参数等。
2. 数值方法常用的CFD数值方法有有限体积法、有限元法、谱方法、边界元法等。
其中，有限体积法是最常用的方法，它的特点是离散区域简单，应用广泛。
3. 软件工具目前使用最多的CFD软件工具包括ANSYS Fluent、OpenFOAM、Star-CCM+、COMSOL Multiphysics等。
这些工具包括了多种物理模型和数值方法，能够应用到各种工程领域中。
二、应用领域数值模拟在工程应用中具有重要的作用，可以帮助工程师和科学家更好地理解物理现象，预测实验结果，并进行优化。
以下是几个CFD在流体力学中的应用领域：1. 汽车工程CFD模拟可以帮助汽车设计师预测汽车性能和气动外观。
如在汽车的外形设计上，可以通过CFD模拟预测空气阻力对车辆速度和燃油经济性的影响。
在引擎设计中，CFD模拟可以预测燃烧过程中的热传导、气流动力学和挥发物扩散等。
2. 船舶工程在船舶工程中，CFD模拟可以用来预测船体的水动力性能，包括船速、抗拒力、侧向力和推力等。
同时，CFD模拟还可以分析船舶的安全性和稳定性，优化船体结构和舵型设计。
港口航道的水动力模型研究
港口航道的水动力模型研究一、引言港口航道作为海洋与内陆之间的重要连接通道，其水动力特性对于港口的运营、船舶的航行安全以及周边环境的保护都具有至关重要的意义。
水动力模型作为研究港口航道水流、波浪等水动力现象的有效工具，能够为港口的规划、设计和管理提供科学依据。
二、水动力模型的基本原理水动力模型通常基于流体力学的基本方程，如连续性方程、动量方程和能量方程等。
这些方程描述了水流的运动规律和物理特性。
在港口航道的水动力模型中，还需要考虑边界条件，如海岸线、港口建筑物、船舶等对水流的影响。
同时，模型还需要对波浪、潮汐等因素进行合理的模拟。
三、常见的水动力模型类型（一）二维水动力模型二维水动力模型主要考虑水平方向上的水流运动，适用于研究大面积的水域，如海湾、河口等。
它能够较好地模拟水流的平均状态和宏观趋势，但对于垂直方向上的水流变化和局部复杂地形的模拟能力相对较弱。
（二）三维水动力模型三维水动力模型能够更全面地考虑水流在空间三个方向上的运动，对于港口航道中复杂的水流结构、漩涡和分层现象等具有更好的模拟能力。
然而，三维模型的计算量较大，对计算资源和数据要求较高。
（三）浅水方程模型浅水方程模型是一种简化的水动力模型，适用于水深相对较浅的港口航道。
它在保证一定精度的前提下，能够大大提高计算效率。
四、水动力模型的构建过程（一）数据收集构建水动力模型首先需要收集大量的基础数据，包括地形数据、水文数据、气象数据等。
地形数据如海岸线、水深等对于准确模拟水流的流动路径至关重要；水文数据如潮位、流速、流向等能够为模型提供初始条件和验证依据；气象数据如风场、气压等则会影响波浪的生成和传播。
（二）网格划分根据研究区域的大小和复杂程度，将其划分为一系列的网格单元。
网格的大小和形状会直接影响模型的精度和计算效率。
在港口航道等重点区域，通常需要采用较精细的网格，以捕捉局部的水动力特征。
（三）参数设置模型中涉及到众多的参数，如糙率系数、涡粘系数等，这些参数的取值需要根据实际情况进行合理的估计或通过现场观测和实验数据进行率定。
计算流体力学中的有限体积法
计算流体力学中的有限体积法
有限体积法（FVM）是一种数值计算方法，用于模拟流体力学问题。
它是通过把流场分成很多个离散化的小体积来描述流体的运动的。
有限体
积法的基本思想是在每一个小体积中应用质量、动量和能量守恒方程，然
后将它们组合成一个离散化的形式，以便于数值计算。
其中，质量守恒方
程描述了流体的连续性，动量守恒方程描述了流体的运动，能量守恒方程
则描述了流体的温度和压力等性质随时间的变化。
有限体积法的计算流程一般包括以下步骤：
1.网格划分：将流场划分成若干个小体积，每个小体积称为一个网格
单元。
2.定义控制体：在每个网格单元内，定义一个控制体。
控制体是一个
虚拟的小体积，它可以是任意形状，但通常为正交体。
3.求解守恒方程：对于每个控制体，应用守恒方程，得到一个自由度
方程组。
4.数值求解：利用数值方法求解自由度方程组，得到解。
5.更新场变量：根据求解得到的解，更新场变量（如速度、压力等）。
6.考虑边界条件：在每个边界上，根据物理条件定义边界条件，用于
修正解。
7.重复以上步骤：对于每个时间步长，重复以上步骤，直到计算结束。
需要注意的是，有限体积法是一种局域方法，只考虑每个网格单元内
部的守恒方程，没有直接考虑两个网格单元之间的相互作用。
因此，在计
算边界处或流场中存在复杂流动结构的区域时，需要采用一些特殊的技术（如插值方法、外推方法等）来处理。
基于CFD的潜艇阻力及流场数值计算
基于CFD的潜艇阻力及流场数值计算涂海文;孙江龙【摘要】The article numerically simulated the resistance and flow field of bare submarine body with different appendages and submarine body full-appended with RNG k-e model. And solving solutions with RANS equations. The girds were divided by ICEM CFD of CFD pretreatment software. Through numerical simulation, it obtained the submarine surface pressure distribution and some characteristics of the flow field near appendages. This is further to lay the foundation for optimizing the boat type and analyzing noise of submarine that induced by flow. And resistance contrast verified the reliability of numerical simulation to some degree.%运用雷诺平均N-S方程,使用CFD前处理软件ICEM CFD划分流场网格,采用RNG k-ε湍流模型,实现了对裸潜体、带指挥台围壳艇体、带十字尾翼艇体、全附体潜艇4种模型的阻力及粘性流场的数值模拟.通过数值模拟,得到了潜艇表面压力分布情况和附体附近流场的一些特性,为进一步优化潜艇的艇型和分析潜艇的流噪声打下了基础.而阻力的对比在一定程度上验证了数值模拟的可靠性.【期刊名称】《舰船科学技术》【年(卷),期】2012(034)003【总页数】7页(P19-25)【关键词】潜艇;CFD;RNG k-ε模型;摩擦阻力【作者】涂海文;孙江龙【作者单位】华中科技大学船舶与海洋工程学院,湖北武汉430074;华中科技大学船舶与海洋工程学院,湖北武汉430074【正文语种】中文【中图分类】U674.76;TB5潜艇周围的流场特别是尾流场的流动特性，不仅对潜艇的水动力性能产生直接的影响，而且引发的流动是潜艇水动力噪声的主要来源之一，对潜艇隐蔽性有重要影响。
有限体积法(水力学问题及应用)
目录
目录 ..................................................................................................................................................1 第 8 章 有限体积法.........................................................................................................................2
8.1.1 几个术语
在进行数值计算时，要把计算区域划分成一系列互不重叠的离散小区域，然后在该小区 域上离散控制方程求解待求物理量。在有限差分法中只涉及到网格节点的概念，而有限体积 法因为物理解释需要，形成了以下几个常用几何要素的相关名词。
2
N
N3
n
N2
W w
P
单元
e
E
N4
sห้องสมุดไป่ตู้
S
N1
图 8-1 单元中心有限体积法中典型的矩形网格 如图 8-1 所示，为单元中心格式（cell-centered scheme）有限体积法中典型的矩形网格关 系。其包含以下几个要素： 控制体积（control volume）：图中阴影部分，方程积分离散时的小体积单元（二维为面积 单元）。 单元（cell）：控制体积的中心，常用形心来表示，为待求物理量的几何位置。图中用空心 园来表示，如点 P、W、E、N、S 等。常用单元来代表整个控制体积。以下如若不作特殊说明，则 用“单元”来代表控制体积。 网格线（grid lines）：用来分割计算区域内各控制体积的交错曲线簇，如图中折线 N1N2、 N2N3、N3N4、N4N1 等 网格节点（nodes）：网格线之间的交点， 图中用黑圆点来表示，如 N1、N2、N3、N4等 单元界面（cell faces）：相邻两个控制体积间的公共面(二维则可以认为是公共边，这样看
基于有限体积法对洪水汇合流流场的数值模拟
散得：
∑ ( ) ( ) Uin+=1
U in
−
∆t Vi
3
T −1 θi, j
j =1
F Ui, j
Li, j + ∆tSi
(5)
上式中源项采用中心差分进行离散，基于局部 Riemann 问题自相似解波系结构分析，采用 TVD-WAF 格式给出数值通量为：
∑ ( ) FWAF=
1 2
( c4
FR
2. 数学模型
2.1. 控制方程
通过建立蓄洪区退洪、长江河道行洪和分洪区分洪的数学物理方程，采用二维浅水波方程组研究多 股洪水汇合流问题。浅水波方程组采用以水位和单宽流量为变量的守恒形式，即：
∂U + ∂F + ∂G = S ∂t ∂x ∂y
(1)
式中 t 代表时间，x 和 y 分别表示空间变量。U 表示物理守恒变量， F (U ) 、G (U ) 分别为 x 和 y 方向的
q 2y
+
1
gH 2
−
gbH
h 2
0
0
S
=
Sb
+ Sf
=
− gH
∂b ∂x
+
−
g ρ
τ
x
− gH
∂b ∂y
−
g ρ
τ
y
其中 qx = hu ， qy = hv 表示 x，y 方向的单宽流量，u，v 分别表示 x，y 两个方向上的流速， H= h + b 为
fvm 有限体积法
fvm 有限体积法有限体积法（Finite Volume Method，FVM）是一种常用的数值计算方法，用于求解流体力学和热传导等守恒方程。
它将计算区域划分为离散的控制体，通过在控制体上应用平衡方程，将守恒方程转化为代数方程组，进而得到数值解。
在有限体积法中，计算区域被划分为若干个控制体，每个控制体代表一个小区域，用来计算相应的物理量。
这些控制体之间通过边界面相连，形成一个网格。
在每个控制体内，平均物理量被定义为该控制体上物理量的积分平均值。
通过在控制体上应用平衡方程，可以得到守恒方程的离散形式。
有限体积法的基本思想是，根据质量守恒、动量守恒和能量守恒等守恒方程，将守恒方程在每个控制体上进行积分，然后通过对积分方程进行离散化，得到代数方程组。
这个代数方程组可以通过数值方法求解，得到每个控制体上的物理量的数值解。
在有限体积法中，流体流动被描述为流体在控制体内的质量、动量和能量的变化。
通过在控制体上应用守恒方程，我们可以得到质量守恒方程、动量守恒方程和能量守恒方程的离散形式。
这些方程是以每个控制体为基础的，通过将守恒方程在控制体上进行积分，可以得到每个控制体上物理量的离散形式。
有限体积法的求解过程包括以下几个步骤：首先，将计算区域划分为离散的控制体，并在每个控制体上定义平均物理量。
然后，通过在控制体上应用守恒方程，将守恒方程转化为代数方程组。
接下来，使用数值方法求解代数方程组，得到每个控制体上物理量的数值解。
最后，根据数值解，可以得到流体流动的各种性质，如速度、压力、温度等。
有限体积法在工程领域得到了广泛的应用。
例如，在流体力学领域，有限体积法可以用来模拟流体的流动，计算流动的速度、压力分布等。
在热传导领域，有限体积法可以用来模拟热量的传输，计算温度的分布等。
在材料科学领域，有限体积法可以用来模拟材料的变形和应力分布等。
有限体积法是一种常用的数值计算方法，适用于求解守恒方程。
它通过在控制体上应用平衡方程，将守恒方程转化为代数方程组，并通过数值方法求解代数方程组，得到物理量的数值解。
船舶艏侧推器推力减缩试验与数值计算研究
（ １ ） 、（ ２ ）和 （ ３ ）计 算得
到 相应 的 推力 系数 Ｋｒ 、转 矩系 数 和 效率 ， ７ 。
Ｙ
＝
＝
Ｑ
（ ２ ）
Ｋ
。
Ｊ
（ ３ ）
黼 ７ ｋ｛ 验 幺 圭单 加 嘉
Ｉ昕 景
表 １ 敝 水 试验 结 果
进速 系数 流速／ （ １ ｎ ／ ｓ ） 桨模转速／ （ ｓ ） Ｋｒ Ｍ １ ０ Ｋ Ｑ Ｍ 『 ７
１ ． ２ ． ２ 试 验结 果
艏 侧 推 器螺旋 桨 为 可调 距 螺旋 桨 ， 试验 中， 仪１ ， ５ ／ ， ＋ １ ９ 。 桨 叶 为 试验 对 象 ， 桨 模直 径 ０ ． ２ ４ ｍ， 转速 ９ ． ７ ｌ ｒ ／ ｓ ，
４
国
造
船
学 术论 文
很 姑 螺 旋 采 牛 田 女 装 削 传 感 器 侍 削 螺 旋 采 瞅 水 Ｈ 削 惟 刀 利 转 旭 ， 由式
０ ． ０ ９ ０ ４ ６ ０ １ ３ ５ ７
０． ２ Ｏ． ２ ５ ０． ３ ０． ３ ５
０ ． ２ １ Ｏ ． ３ ２
０． ４ ７ ０． ５ ８ ０． ７ ０ ０． ８ ２
９ ． ７ ２ ９ ． ７ ２
９ ． ７２ ９ ． ７２ ９． ７２ ９ ． ７２
的存在而产生推力减缩 ，致使推力达不到设计要求。因此 ，开展侧推器水动力性能及推力减缩研究具
有 重 要 的理论 意义 和 工程 应 用价 值 。
Ｌ ｅ ｈ ｎ ｔ ｌ ￣ 和 Ｍｏ ｂ ｅ ｒ ｇ ［ ２ １ 对两 个螺 旋 桨呈 前后 布置 的形 式进 行 了试验 研 究 ，结果表 明：螺旋 桨之 间 的干 扰 强度 同螺旋桨 的距 离成 负相 关 ，同时 ，螺旋 桨 的尾 流影 响 区域较大 ，在 下游 １ ６倍直径 处 布置第 二个 螺 旋桨 ，其 推力 减 缩幅度 仍 有 可能 达 ２ ５ ％。Ｎｉ ｅ ｎ ｈ ｕ ｉ ｓ ｌ ３ ］ 对不 同角度来 流和 不 同距 离 的布 置方 案进 行 了试
船舶排水体体积计算公式
船舶排水体体积计算公式船舶排水体体积是指船舶在水中排水时所占据的体积，它是船舶设计和建造中非常重要的一个参数。
船舶排水体体积的计算可以通过以下公式进行：V = L x B x T x C。
其中，V表示船舶排水体体积，L表示船长，B表示船宽，T表示船舶的吃水深度，C表示船舶的垂线系数。
这个公式是根据船舶的几何形状和水线形状来推导的，通过这个公式可以比较准确地计算出船舶的排水体积。
船舶排水体体积是船舶设计和建造中的重要参数，它直接影响着船舶的稳定性、载重能力和航行性能。
因此，对船舶排水体体积的准确计算是非常重要的。
首先，我们来看一下船长、船宽和吃水深度这三个参数的定义和影响。
1. 船长（L），船长是指船舶的整体长度，它是船舶设计中的一个重要参数。
船长的大小直接影响着船舶的载重能力和稳定性。
一般来说，船长越长，船舶的载重能力就越大，但也会增加船舶的阻力和航行成本。
2. 船宽（B），船宽是指船舶的最大宽度，它也是船舶设计中的一个重要参数。
船宽的大小直接影响着船舶的稳定性和航行性能。
一般来说，船宽越宽，船舶的稳定性就越好，但也会增加船舶的阻力和航行成本。
3. 吃水深度（T），吃水深度是指船舶在水中的下沉深度，它是船舶设计中的一个重要参数。
吃水深度的大小直接影响着船舶的载重能力和稳定性。
一般来说，吃水深度越深，船舶的载重能力就越大，但也会增加船舶的阻力和航行成本。
接下来，我们来看一下船舶的垂线系数（C）。
垂线系数是一个无量纲的参数，它是根据船舶的几何形状和水线形状来计算的。
垂线系数的大小直接影响着船舶的排水体积和稳定性。
一般来说，垂线系数越小，船舶的排水体积就越小，但也会影响船舶的稳定性和航行性能。
通过上面的公式和参数的定义，我们可以看到船舶排水体体积的计算是一个比较复杂的过程，需要考虑到船舶的几何形状、水线形状和稳定性等多个因素。
因此，在实际的船舶设计和建造中，需要借助计算机辅助设计软件来进行排水体积的计算，以确保计算结果的准确性和可靠性。
船舶航行性能的数值模拟与分析
船舶航行性能的数值模拟与分析船舶作为重要的水上交通工具，其航行性能的优劣直接关系到航行的安全、效率和经济性。
随着计算机技术的飞速发展，数值模拟已成为研究船舶航行性能的重要手段。
通过数值模拟，可以在船舶设计阶段就对其航行性能进行预测和优化，从而减少试验次数、缩短研发周期、降低成本。
数值模拟的基本原理是基于流体力学和船舶动力学的相关理论，将船舶和周围的水流视为一个连续的流体场，通过求解一系列的控制方程来获得船舶周围流场的压力、速度等参数，进而计算船舶的航行性能。
在船舶航行性能的研究中，常用的数值模拟方法包括有限元法、有限体积法和边界元法等。
船舶航行性能主要包括阻力性能、推进性能、操纵性能和耐波性能等方面。
阻力性能是船舶航行性能中的一个关键指标，它直接影响船舶的动力需求和燃油消耗。
通过数值模拟，可以对船舶在不同速度、吃水和姿态下的阻力进行计算，并分析船体形状、附体布置等因素对阻力的影响。
例如，优化船体的首部形状可以减少兴波阻力，而合理设计船底的粗糙度可以降低摩擦阻力。
推进性能是衡量船舶动力系统效率的重要指标。
数值模拟可以用于研究螺旋桨的水动力性能，包括推力、扭矩和效率等。
通过模拟螺旋桨在不同转速和进流条件下的工作情况，可以优化螺旋桨的叶片形状和布置，以提高推进效率。
同时，还可以考虑螺旋桨与船体之间的相互干扰，从而更准确地评估推进系统的性能。
操纵性能是船舶在航行中改变航向和速度的能力。
数值模拟可以模拟船舶在不同舵角和螺旋桨转速下的运动响应，计算船舶的回转半径、转向时间等操纵性指标。
对于多桨多舵的船舶，还可以研究不同桨舵组合对操纵性能的影响。
此外，还可以通过模拟船舶在风浪中的操纵情况，评估船舶在恶劣海况下的操纵安全性。
耐波性能是船舶在波浪中航行时的性能表现，包括船舶的运动响应、砰击和上浪等。
数值模拟可以模拟船舶在规则波和不规则波中的运动，分析船舶的纵摇、横摇和垂荡等运动特性。
通过优化船舶的型线和结构，可以减少船舶在波浪中的运动幅度，提高船舶的舒适性和安全性。
船舶浅水水动力导数的数值计算
第34卷第3期 2011年9月中 国 航 海N AV IG AT ION O F CH IN AVol.34No.3 S ep.2011收稿日期:2011-05-16作者简介:石爱国(1956-),男,吉林德惠人,教授,从事舰船操纵性研究。
E -mail:AGShi56@.文章编号:1000-4653(2011)03-0069-05船舶浅水水动力导数的数值计算石爱国, 闻 虎, 李 理, 刘 可, 刘 博(海军大连舰艇学院,辽宁大连116018)摘 要:船舶浅水水动力导数对研究浅水中船舶操纵性有重要的意义。
以/M ariner 0船模为研究对象,采用Realiz -ablek-E 湍流模型来封闭RA N S 方程,运用SIM PL E 算法,对两种水深的浅水定漂角、定舵角、纯艏摇试验进行了数值模拟。
实现了浅水水动力导数的求取,并将计算结果与模型试验结果进行对比,验证了方法的有效性。
关键词:船舶,舰船工程;操纵性;计算流体力学;浅水水动力导数;数值模拟;模型试验中图分类号:U 661.1 文献标志码:AComputation of Hydrodynamic Derivatives for Ships in Shallow WaterShi A ig uo , Wen H u , L i L i, L iu K e, L iu Bo (Dalian N av al Academ y,Dalian 116018,China)Abstract:H ydro dy namic deriv atives of ships in shallow w ater ar e essential for study ing ship maneuv erability in sha-l lo w w ater.T aking t he ship model "M a riner"as the study ing object,this pa per carr ies o ut a series of numerical sim -ulatio n under stat ic dr ift,st atic rudder and pure y aw in tw o differ ent w ater depths by so lving the Reynolds -aver age N -S equations w ith Realizable turbulence model using SIM PL E algo rithm.T he numerical results are compared w ith ex per iment al results and calculated results published in literature to pro ve t he effect iveness of the numer ical method.Key words:ship,nav al engineer ing ;maneuv erability ;CFD;hy dr odynamic der ivatives f or shallow water ;numer ical simulation;mo del ex periment操纵性是船舶的重要航海性能,与航行的安全性和营运的经济性密切相关。
船体绕流场及流噪声的CFD模拟方法
船体绕流场及流噪声的CFD模拟方法刘波; 武国启【期刊名称】《《舰船科学技术》》【年(卷),期】2019(041)010【总页数】7页(P63-69)【关键词】绕流场; 流致噪声; CFD技术; 水下噪声; 数值模拟【作者】刘波; 武国启【作者单位】大连测控技术研究所辽宁大连116013【正文语种】中文【中图分类】U661.10 引言为了确定船体绕流场及其绕流发声问题，船模试验是一个非常重要的实验手段。
但是与CFD 方法比较，船模试验在时间和试验成本等消耗上是巨大的，且不能对船体周围流场及声场做出局部的特征估计和预报。
在2000 年召开的第四次国际船舶水动力学数值计算学术研讨会上，参加会议的学者讨论认为：采用计算流体动力学研究船舶绕流场已经达到很髙的精度要求，船舶各个部分的绕流情况基本可以利用CFD 进行精确的模拟。
对于物体的绕流问题，Zebib 等[1] 采用数值计算的方法研究了圆柱体的绕流扰动问题，结果表明圆柱体在对称情况下绕流扰动是稳定的，非对称情况下绕流不稳定。
Y.Ahmed 等[2] 利用CFD 数值模拟了VLCC 船体周围的自由液面情况，为船舶阻力确定提供了很好的预报；之后又利用CFX 模拟了不同弗劳德数下复杂船体绕流场的变化情况，模拟结果与实验数据非常接近[3]。
物体绕流发声问题首先在航空和航海领域得到高度关注，航空方面它会造成航空乘员感官和心理上的不适，航海方面对海洋生物短期和长期造成严重负面影响，特别是对海洋哺乳动物的影响更大。
对于物体绕流发声的研究，Lighthill[4-5]提出的声类比理论奠定了流致发声的理论基础，随后各国学者对绕流发声问题开展广泛的研究。
Paula Kellett[ 6 ]等基于CFD 研究了LNG 船水动力性能预报问题，利用FW-H方程计算船舶水下辐射噪声，并与实船噪声测量数据进行了比较，验证了船舶水下噪声预报的可行性。
曾文德[7]、魏英三[8] 等针对SUBOFF潜艇开展CFD 数值计算研究，获得了潜艇表面辐射声源的强度分布，指出水流动转捩区是主要的发声部位。
船舶水动力学模拟的数值方法与应用
船舶水动力学模拟的数值方法与应用随着计算机技术不断发展，数值计算在船舶水动力学研究中得到了广泛应用。
船舶水动力学模拟的数值方法能够模拟船舶在不同工况下的流场、压力、速度等水动力参数，并结合船体结构进行分析，从而为船舶设计、性能评估和运行控制等提供科学依据。
本文将从数值方法和应用两个方面进行阐述。
数值方法目前，船舶水动力学模拟的数值方法可大致分为两类：基于潜水体积法的方法和基于边界元法的方法。
潜水体积法是通过将流场分割成小的节点，并针对每个节点进行流场的密度和速度计算。
这种方法最适合模拟大规模的流场，具有高精度和高稳定性的特点，但计算时间较长。
相对而言，边界元方法的计算速度更快，但是只考虑流场的边界，并未将整个流场分割成节点。
这种方法的精度和稳定性略低于潜水体积法。
应用船舶模拟的数值方法在船舶设计、性能评估和船体结构分析等方面有着广泛的应用。
首先，船舶设计方面。
在船舶设计中，可通过船舶水动力学数值模拟来预估船舶的性能表现和稳定性，从而为船体设计提供依据。
同时，也可通过模拟船舶在不同风浪情况下的性能表现，评估船体的适应性。
其次，性能评估方面。
在船舶建造完成後，通过模拟船舶在不同工作状态下的性能表现，评估船舶的性能，如航行速度、耗油量和稳定性等，可为船舶的管理提供科学依据。
最后，船体结构分析方面。
在船舶使用过程中，各种外在因素可能会对船体产生损伤，如海浪、风、撞击等。
通过船舶水动力学模拟数值方法，可以分析船体在外力作用下的应力分布情况和变形情况，进而实现对船体结构的修复和保养。
总结在数值计算方面，船舶水动力学数值模拟主要分为潜水体积法和边界元法。
在应用方面，船舶水动力学模拟数值方法广泛应用于船舶设计、性能评估和船体结构分析等方面，为船舶的设计、建造、保养和性能管理等提供了一种有力的工具。
随着计算机技术的发展，船舶水动力学模拟数值方法也会不断地得到改进和完善，并为未来的船舶建造和维护提供更为准确和高效的技术支持。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
武 汉 造 船 No. 4. 1998
船舶水动力数值计算中的 有限体积法
张谢东 吴秀恒
( 武汉交通科技大学)
摘 要 有限体积法是一种较新的流体力学数值计算方法 ,其对于数值求解船舶水动力学中常常遇到 的偏微分方程十分有效 。本文对有限体积法作了初步的研究与探讨 ,并采用此方法对圆柱绕流问题进行了实 例计算 。数值计算结果与解析作了比较 ,吻合良好 。 关键词 船舶 流体力学 有限体积法
A B CD 内运用子域法及 Green 公式 ,可以得到
4
© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
1998 年第 4 期
武 汉 造 船 No. 4. 1998
1 前言
在现代化船舶研究中 ,随着船舶性能的不断提 高 ,对船舶水动力的研究也提出了更高的要求 。除 了要对船舶的宏观量如力和力矩进行精确预测外 , 还需要对船舶周围的流场细节有精确的了解 。要 做到这些 , 除了可以进行实船或船模的试验研究 外 ,还可以采用 CFD 技术进行理论计算研究 。众 所周知 ,船舶试验研究费时费力 , 且其精度还往往 受到试验条件和环境的影响 。因此 ,经济快速的理 论计算研究越来越受到人们的重视 。 采用理论计算方法 ,准确地预报出船舶周围的 流场及船体水动力 , 能够用于指导船舶设计 , 这对 于提高效率 , 节省人力和物力 , 设计出高效 、 节能 、 降噪的优良船型 , 具有及其重要的意义 , 并能产生 巨大的经济效益 。因此 ,船舶流体力学的数值计算 方法研究十分必要 ,而且具有广阔的前景 。 当前 ,在船舶流体力学中 , 通常采用的数值计 算方法有 : 有限差分法 、 有限元法 、 边界元法和有限 体积法 。有限差分法 ( Finite Difference Met hod) 是 流体力学方程数值解法中应用最广泛的方法之一 , 它基于 Taylor 级数展开 , 在离散的网格点上把各 偏导数项转化为差商 ,从而将微分方程按照一定的 格式离散为差分方程 ,然后来求它的数值解 。影响 差分方程求解精度的关键因素是差分格式的选取 。 有限元法 ( Finite Element Met hod) 是将偏微分方程 的求解区域划分为有限个子域或单元 ,将有限个单
5 5y 5 5y 5 5y
Δx j +1 , j + ΔyBcΔy j +1 , j ) / S Bc , PB C = (Δx B C
2 2 Q CD = (Δx CD + Δy CD ) / S CD ,
j +1/ 2
Δx B C + Δx CD +
PCD = (Δx CDΔx k +1 , k + Δy CDΔy k +1 , k ) / S CD ,
j +1 , k +1
+
+
j- 1, k j +1 , k
Q DA分 + 0 . 25 ( PCD - PA B ) Q B C + 0 . 25 ( PA B - PCD )
+ +
元
图3 求解圆柱绕流的计算区域
0 . 25 ( PDA - PA B ) 0 . 25 ( PAB - PB C)
j- 1 , k- 1
j , k +1 / 2
j- 1 , k
-
j , k) + PDA (
-
D) = 0
其中
Q AB = Δx 2 y2 AB + Δ A B ) / S AB , PA B = (Δx ABΔx k - 1 , k + Δy A BΔy k - 1 , k ) / S A B ,
2 2 Q B C = Δx B C + Δy B C ) / S B C ,
ε
图2 曲线坐标的有限体积
以同样的方法处理 fN 5 /5x 式 ( 4) 可得
Q AB (
j , k- 1 C
j +1/ 2 , k
等项 , 然后代入
j +1 , k D
B)
j , k)
+ PAB (
j , k- 1
B
-
A) j , k)
+ QBC ( + PCD (
C)
j , k)
+
图1 二维的有限体积
PB C ( Q DA (
-
+ Q CD (
A
-
+ ( 8)
方程 ( 2) 可近似地表示为 DA 5 5 ( Δy + Δx ) = 0 ( 3) 6 5 x 5y AB 上式可以写为 5 Δy 法-, 5 Δx + 能产 5 x j , k - 1/ 2 A B 5 y j , k - 1/ 2 A B 5 Δ yl 分法 y 5 x j +1/ 2 , k B C 5 5x 5 5x
-
j , k)
+ Δy k - 1 , k (
B
-
A)
S AB
当网格点确定之后 ,式 ( 11) 中的 PA B 和 Q AB 等就可 ы х 以唯一地求出 , 可以称之为网格参数 。
5
© 1994-2010 China Academic Journal Electronic Publishing House. All rights reserved.
+
j , k- 1 n
Q A B + 0 . 25 ( PDA - PB C )
j +1 , k - 1
+ ( 13 )
/
( Q A B + Q B C + Q CD + Q DA )
其求解控制方程为 Laplace 方程 ( 1 ) , 在采用 FVM 计算时 ,将计算区域 X Y Y′ X′ 划分为 m × m ( ) ( ) 个网格单元 ,利用式 12 及式 13 按超松弛迭代法 足 ( SOR) 进行迭代计算 ,直到收敛 。 这里 ,取圆柱半径 R a = 0 . 1 , 计算区域半径 R b
l(
AB CD
∫ 5 5 dy + dx = 0 ∫ 5x 5y
AB CD
52 52 ) d xd y = 2 + 5x 5 y2
( 2)
5 5y
j , k - 1/ 2
=
-
Δx AB (
j , k- 1
rS
of j, k)
+ Δx k - 1 , k ( S AB
B
-
A)
( 7)
方程 ( 2) 是一种守恒形式的方程 。与有限差分 法总是将控制方程变为差分形式不同 ,有限体积法 通常是将控制方程表述为积分形式 。
= =
1
S A′ B′ C′ D′
5 5y
1
S A′ B′ C′ D′
j , k - 1/ 2
κ κ
5 d x dy = 5x 5 d xd y = 5y
1
S A′ B′ C′ D′
1
S A′ B′ C′ D′
∫ dx ∫
dy
( 5)
( 8) 代入式 此 , 体积九点表达式 : 0 . 25 ( PCD - PDA ) 0 . 25 ( PB C - PCD )
+
+
j- 1, k j,k
Δy C′ D′ +
Δy D′ , A′
A′ B′ C′ D′
∫ d x 同理 。
Q DA + 0 . 25 ( PCD - PAB ) ( Q AB + Q B C + Q CD + Q DA ) Q B C + 0 . 25 ( PA B - PCD )
-
+ +
如果网格不是十分弯曲 ,则有 Δy A ′ y C″ yA B , B′≈ - Δ D′≈ Δ ΔyB′ y D′ yk - 1 , k 和 C′≈ - Δ A ′≈ Δ Δy k - 1 , k - Δy A BΔx k - 1 , k ( 6) S AB = S A ′ B′ C′ D′ = Δx A B 于是 ,式 ( 5) 变为
元的结点作为离散点 。如果单元内的函数值能通 过这些离散点上未知函数值唯一地确定 ,则通过在 整个求解区域内使泛函的变分为零 ,便能引出一组 联立代数方程组 。求解此方程组 ,即可确定所有结 点上 的 函 数 值 。边 界 元 法 ( Boundary Element
Met hod) 只将计算区域的边界进行离散 , 然后从边
它 的p ,于是可以得到 Laplace l
今后 方程 ( 1 ) 的有限
j - 1 , k +1
+
所 示
j , k +1
式中
j,k
A′ B′ C′ D′
∫
dy ≈
A
j , k- 1
Δy A ′ B′ +
B
ΔyB′ C′ +
Q CD + 0 . 25 ( PB C - PDA )
j +1 , k +1
1998 年第 4 期
武 汉 造 船 No. 4. 1998 示)
方程 ( 11) 能够非常方便地采用逐次超松弛迭 代方法 ( SOR) 进行求解 ,它的解转变为以迭代值的 形式
n +1 j,k
2 2 离散 Q DA = (Δ x DA + Δy DA ) / S DA ,
Δy CD 度的 Δy DA -