族谱网 头条 人物百科

开尔文船波

2020-10-16
出处:族谱网
作者:阿族小谱
浏览:659
转发:0
评论:0
多鞍点函数积分IntegrandofKelvinWakeIntegralKelvinShipWakeIntegrandcontourMapleplot当船只以速度V驶过深水湖面,波形的幅度在相对于船只为静止的极坐标(ρρ-->,ϕϕ-->{displaystylerho,phi}中在船只的速度

多鞍点函数积分

开尔文船波

Integrand of Kelvin Wake Integral

开尔文船波

Kelvin Ship Wake Integrand contour Maple plot

当船只以速度V驶过深水湖面,波形的幅度在相对于船只为静止的极坐标(ρ ρ -->,ϕ ϕ -->{\displaystyle \rho ,\phi }中在船只的速度矢量方向,ϕ ϕ -->=0{\displaystyle \phi =0}),由下列公式表示

K(ϕ ϕ -->,ρ ρ -->)=∫ ∫ -->− − -->π π -->/2π π -->/2cos(ρ ρ -->cos(θ θ -->+ϕ ϕ -->)cos2θ θ -->dθ θ -->{\displaystyle K(\phi ,\rho )=\int _{-\pi /2}^{\pi /2}cos(\rho {\frac {cos(\theta +\phi )}{cos^{2}\theta }}d\theta }

其中ρ ρ -->=gr/V2{\displaystyle \rho =gr/V^{2}}

1ρ ρ -->=V2gr{\displaystyle {\frac {1}{\rho }}={\frac {V^{2}}{gr}福禄数福禄数的平方Fr2{\displaystyle Fr^{2}}

g{\displaystyle g}为重力常数l{\displaystyle l}为船的长度。

上列K函数是下列多鞍点积分的正数部分:

K(ϕ ϕ -->,ρ ρ -->)=Re(∫ ∫ -->− − -->∞ ∞ -->∞ ∞ -->exp⁡ ⁡ -->(i∗ ∗ -->ρ ρ -->∗ ∗ -->f(θ θ -->,ρ ρ -->)dθ θ -->){\displaystyle K(\phi ,\rho )=Re(\int _{-\infty }^{\infty }\exp(i*\rho *f(\theta ,\rho )d\theta )} 其中,多鞍点积分的核函数为

f(θ θ -->,ϕ ϕ -->)=− − -->cos(θ θ -->+ϕ ϕ -->)cos2θ θ -->{\displaystyle f(\theta ,\phi )=-{\frac {cos(\theta +\phi )}{cos^{2}\theta }}}

此核函数是一个多鞍点函数,振荡剧烈如图

求其极点,

df(θ θ -->,ϕ ϕ -->)dθ θ -->=sin(θ θ -->+ϕ ϕ -->)cos(θ θ -->)2− − -->2∗ ∗ -->cos(θ θ -->+ϕ ϕ -->)∗ ∗ -->sin(θ θ -->)cos(θ θ -->)3=0{\displaystyle {\frac {df(\theta ,\phi )}{d\theta }}={\frac {sin(\theta +\phi )}{cos(\theta )^{2}}}-{\frac {2*cos(\theta +\phi )*sin(\theta )}{cos(\theta )^{3}}}=0}

解之,得

θ θ -->1=arctan((1/4)∗ ∗ -->(1+(1− − -->8∗ ∗ -->tan(ϕ ϕ -->)2))tan(ϕ ϕ -->))=− − -->arctan((1/4)∗ ∗ -->(− − -->1+(1− − -->8∗ ∗ -->tan(ϕ ϕ -->)2))tan(ϕ ϕ -->)){\displaystyle \theta _{1}=arctan({\frac {(1/4)*(1+{\sqrt {(1-8*tan(\phi )^{2}))}}}{tan(\phi )}})=-arctan({\frac {(1/4)*(-1+{\sqrt {(}}1-8*tan(\phi )^{2}))}{tan(\phi )}})}

由此

ϕ ϕ -->1=19.47{\displaystyle \phi _{1}=19.47}度,

ϕ ϕ -->2=− − -->19.47{\displaystyle \phi _{2}=-19.47}度

这就是凯尔文船波的V型波包线的夹角,最早由凯尔文男爵发现,而且角度与船速无关.至于波纹本身则与船速矢量的夹角为

θ θ -->=π π -->− − -->19.47=35.3{\displaystyle \theta =\pi -19.47=35.3}°

开尔文驻相法

开尔文船波

Kelvin Wake (Maple density plot)

开尔文船波

开尔文船波波形

开尔文船波积分K(ϕ ϕ -->,ρ ρ -->){\displaystyle K(\phi ,\rho )}必须通过数值积分计算。开尔文男爵根据被积分函数在积分区间内剧烈震荡的特点,提出了驻相法(Method of Stationary Phase)。

原理:当被积分函数剧烈震荡时,除了在极点外,震荡的被积分函数正负相抵消,因此可以将此被积分函数在极点的值作为整个积分的近似,驻相法乃是拉普拉斯方法的推广。

被积分函数 f(θ θ -->,ϕ ϕ -->)=− − -->cos(θ θ -->+ϕ ϕ -->)cos2θ θ -->{\displaystyle f(\theta ,\phi )=-{\frac {cos(\theta +\phi )}{cos^{2}\theta }}} 的两个极点是:

θ θ -->p=arctan((1/4)∗ ∗ -->(1+(1− − -->8∗ ∗ -->tan(ϕ ϕ -->)2))tan(ϕ ϕ -->)){\displaystyle \theta _{p}=arctan({\frac {(1/4)*(1+{\sqrt {(1-8*tan(\phi )^{2}))}}}{tan(\phi )}})}

θ θ -->m=− − -->arctan((1/4)∗ ∗ -->(− − -->1+(1− − -->8∗ ∗ -->tan(ϕ ϕ -->)2))tan(ϕ ϕ -->)){\displaystyle \theta _{m}=-arctan({\frac {(1/4)*(-1+{\sqrt {(}}1-8*tan(\phi )^{2}))}{tan(\phi )}})}

fm=f(θ θ -->m,ϕ ϕ -->)=sin((1/2)∗ ∗ -->ϕ ϕ -->− − -->(1/2)∗ ∗ -->arcsin(3∗ ∗ -->sin(ϕ ϕ -->)))sin((1/2)∗ ∗ -->ϕ ϕ -->+(1/2)∗ ∗ -->arcsin(3∗ ∗ -->sin(ϕ ϕ -->))){\displaystyle f_{m}=f(\theta _{m},\phi )={\frac {sin((1/2)*\phi -(1/2)*arcsin(3*sin(\phi )))}{sin((1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}}}

fp=f(θ θ -->p,ϕ ϕ -->)=cos((1/2)∗ ∗ -->ϕ ϕ -->+(1/2)∗ ∗ -->arcsin(3∗ ∗ -->sin(ϕ ϕ -->)))cos(− − -->(1/2)∗ ∗ -->ϕ ϕ -->+(1/2)∗ ∗ -->arcsin(3∗ ∗ -->sin(ϕ ϕ -->))){\displaystyle f_{p}=f(\theta _{p},\phi )={\frac {cos((1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}{cos(-(1/2)*\phi +(1/2)*arcsin(3*sin(\phi )))}}}

fbar:=1/2∗ ∗ -->(fp+fm){\displaystyle fbar:=1/2*(f_{p}+f_{m})}

D2F=d2F(θ θ -->,ϕ ϕ -->)dθ θ -->2{\displaystyle D2F={\frac {d^{2}F(\theta ,\phi )}{d\theta ^{2}}}}

D2Fp=D2F(θ θ -->p,ϕ ϕ -->){\displaystyle D2F_{p}=D2F(\theta _{p},\phi )}

D2Fm=D2F(θ θ -->m,ϕ ϕ -->){\displaystyle D2F_{m}=D2F(\theta _{m},\phi )}

Δ Δ -->:=(3/4∗ ∗ -->(fm− − -->fp))(2/3){\displaystyle \Delta :=(3/4*(f_{m}-f_{p}))^{(}2/3)}

u=Δ Δ -->1/22∗ ∗ -->(1D2Fp+1− − -->D2Fm){\displaystyle u={\sqrt {\frac {\Delta ^{1/2}}{2}}}*({\frac {1}{\sqrt {D2F_{p}}}}+{\frac {1}{\sqrt {-D2F_{m}}}})}

v=2Δ Δ -->1/2∗ ∗ -->(1D2Fp− − -->1− − -->D2Fm){\displaystyle v={\sqrt {\frac {2}{\Delta ^{1/2}}}}*({\frac {1}{\sqrt {D2F_{p}}}}-{\frac {1}{\sqrt {-D2F_{m}}}})}

K(ϕ ϕ -->,ρ ρ -->)≈ ≈ -->2∗ ∗ -->π π -->∗ ∗ -->(u∗ ∗ -->cos(ρ ρ -->∗ ∗ -->fbar)∗ ∗ -->AiryAi(− − -->ρ ρ -->(2/3)∗ ∗ -->Δ Δ -->)/ρ ρ -->(1/3)+v∗ ∗ -->sin(ρ ρ -->∗ ∗ -->fbar)∗ ∗ -->AiryAi(1,− − -->ρ ρ -->(2/3)∗ ∗ -->Δ Δ -->)/ρ ρ -->(2/3)){\displaystyle K(\phi ,\rho )\approx 2*\pi *(u*cos(\rho *fbar)*AiryAi(-\rho ^{(}2/3)*\Delta )/\rho ^{(}1/3)+v*sin(\rho *fbar)*AiryAi(1,-\rho ^{(}2/3)*\Delta )/\rho ^{(}2/3))}

开尔文船波的波峰,由下列两个参数方程式描述

x:=X∗ ∗ -->sin(β β -->)∗ ∗ -->(1− − -->(1/2)∗ ∗ -->sin(β β -->)2){\displaystyle x:=X*sin(\beta )*(1-(1/2)*sin(\beta )^{2})}

y:=X∗ ∗ -->sin(β β -->)2∗ ∗ -->cos(β β -->)/(2∗ ∗ -->M){\displaystyle y:=X*sin(\beta )^{2}*cos(\beta )/(2*M)}

参考文献

Frank J. Oliver, NIST Handbook of Mathematical Functions, 2010, Cambridge University Press

Jame Lighthill Waves in Fluids, Cambridge University Press 1979


免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。

——— 没有了 ———
编辑:阿族小谱
发表评论
写好了,提交
{{item.label}}
{{commentTotal}}条评论
{{item.userName}}
发布时间:{{item.time}}
{{item.content}}
回复
举报
点击加载更多
打赏作者
“感谢您的打赏,我会更努力的创作”
— 请选择您要打赏的金额 —
{{item.label}}
{{item.label}}
打赏成功!
“感谢您的打赏,我会更努力的创作”
返回

更多文章

更多精彩文章
打赏
私信

推荐阅读

· 开尔文
历史开尔文勋爵,该单位的同名人物重新定义的提议2005年,国际计量委员会着手于重新定义开尔文,让它有一个更严格的定义。目前(2010)的定义并不能满足于在低于20K和高于1300K处的使用。委员会提出,应重新定义开尔文,以使得波尔兹曼常数是一个定值1.3806505×10J/K。委员会原本希望,这个项目会在2011年的国际度量衡大会上通过,但到了2011年,这个提议又被推后到2014年,成为一个更大的项目中的一部分。从科学的观点来看,这个新定义会使开尔文和其它国际标准单位连接起来,并且可以不再依赖于某种特定的物质。从实际的观点来看,这个新定义不会造成任何区别,水在一个大气压下的凝固点仍然是0°C(32°F,273.15K)。实际使用色温开尔文常用于测量光源的色温。色温基于黑体发出的光的颜色取决于辐射体的温度的原理。温度低于约4000K的黑体为浅红色,温度高于约7500K为浅蓝色。色温在投影...
· 阿波罗飞船
指令/服务舱阿波罗飞船的主要部分是能够承载三人的指令/服务舱,用于在地球轨道上、地球与月球之间、以及在月球轨道上飞行,并能返回地球。该部分是由北美航空(即之后的洛克威尔)所负责建造。指令舱(CM)阿波罗指令舱及其在土星5号火箭上所处的位置指令舱是阿波罗飞船的主要控制中心以及三名宇航员的生活住处。其中包含加压的主船员舱、宇航员的卧椅、控制仪表板、光学电子导航系统(英语:ApolloPGNCS)、通讯系统、环境控制系统、电池、防热盾、反推力系统、前端对接舱口、侧舱门、五个窗口、以及降落伞回收系统。指令舱是整体阿波罗飞船及土星运载火箭中唯一完好返回地球的部分。服务舱(SM)阿波罗服务舱未加压的服务舱包含一个主要的服务推进引擎以及进出月球轨道所需的推进器、一个能进行姿态控制及平移能力的反推力系统、含有氢氧反应物的燃料槽、发散余热至太空中的散热器、以及一个高增益天线。燃料槽除了含有供人呼吸的氧气外,...
· 开尔文环流定理
其中Γ{\displaystyle\Gamma}为材料围线C(t){\displaystyleC(t)}的环流。用更简单的话来说,这条定理所指的是,若观察闭合围线并注意它一段时间(注意所有流体元的运动)的话,则始终两者间的环流相等。本定理在有黏性应力、非保守彻体力(例如科里奥利力)或非正压的压力-密度关系的情况下并不成立。数学证明参见:欧拉方程(流体动力学)材料围线C(t){\displaystyleC(t)}的环流Γ{\displaystyle\Gamma}的定义为:其中u为速度矢量,ds为沿着闭合围线的单元。彻体力保守的非黏性流体的主宰方程为其中D/Dt为实质导数,ρ为流体密度,p为密度,以及Φ为彻体力的势。上式为带彻体力的欧拉方程。正压性条件意味着密度是压力的函数,且为其唯一自变量,即ρ=ρ(p){\displaystyle\rho=\rho(p)}。取环流的实质导数,得:把主宰方程...
· 开尔文-亥姆霍兹机制
开尔文-亥姆霍兹收缩产生的能量在理论上曾推论,来自于收缩释放出的重力位能是太阳的能量来源。计算在这种历程中太阳能释放出多少的能量(假设密度是均匀的),他是依个接近理想的同心圆球壳,重力位能是对所有球壳,从中心到最外层半径,积分的结果。由牛顿力学得知重力位能的形式为:U=−−-->Gm1m2r{\displaystyleU=-{\frac{Gm_{1}m_{2}}{r}}}此处G是万有引力常数,两个质量分别是每一层半径为r厚度为dr的球壳所拥有的质量,从0到所有球壳半径的一次积分。这个陈述(转换)的结果是:U=−−-->G∫∫-->0Rm(r)4ππ-->r2ρρ-->rdr{\displaystyleU=-G\int_{0}^{R}{\frac{m(r)4\pir^{2}\rho}{r}}\,dr}此处R是球体最外层的半径,m(r)是在半径为r之处以内的总质量。将m(r)以体积和密度来表示...
· 波谲云诡,神来之笔——《草船借箭》赏析
“草船借箭”(《三国演义》第46回)是紧接“蒋干中计”之后的又一个精彩的斗智故事。历史上并无诸葛亮用计“借箭”的史实。与这个故事略有瓜葛的记载见于《三国志·吴书·吴主传》注引《魏略》,说建安十八年(213年)孙权与曹操相持于濡须,孙权乘大船去观察曹军营寨,曹操下令乱箭射之;船的一面受了许多箭,偏重将覆,孙权沉着应付,命令将船掉头,让另一面受箭,等“箭均船平,乃还”。这只是被动的“受箭”,而不是主动的“借箭”。在元代的《三国志平话》中,周瑜挂帅出兵后,与曹操在江上打话,曹军放箭,周瑜让船接满箭支而回。但这也只是随机应变的“接箭”,同样不是有计划的“借箭”。由此可见,“草船借箭”完全是《三国演义》的一段杰出创造。作者对事件的主角、时间、地点、原因、过程都进行了根本性的改造,把它纳入“斗智”的范畴,从而写出了这一脍炙人口的篇章。《演义》第45回末尾,在写毕蒋干中计,导致曹操误杀蔡瑁、张允之后...

关于我们

关注族谱网 微信公众号,每日及时查看相关推荐,订阅互动等。

APP下载

下载族谱APP 微信公众号,每日及时查看
扫一扫添加客服微信