边界层模式中参数反演问题研究

来源 :国防科技大学 | 被引量 : 0次 | 上传用户:hsgnln
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
一直以来,边界层都是大气科学中重点研究的课题,其主要的特点是湍流运动显著。边界层的湍流效应主要是通过模式中的湍流黏性系数来体现,准确的湍流黏性系数对于准确描述边界层中的动力学和热力学结构特征至关重要。湍流黏性系数难以通过直接的观测得到,以往的研究中大多采用半经验半理论的方式确定此项参数。在本文中,我们把获取湍流黏性系数视为反问题,通过模式输出值的观测反演计算得到参数的准确值。我们选取弱非线性普朗特模式和广义Ekman动量近似模式(General Ekman momentum approximation model,简记为GEM)进行研究,分别通过变分反演和统计反演计算出模式的湍流黏性系数。
  首先,利用变分反演方法反演弱非线性普朗特模式中的湍流黏性系数。采用WKB方法计算出模式正问题的解,避免了无穷边界的两点边值问题的求解。通过观测和模式计算值之差的平方构建泛函。利用模拟观测数据进行数值试验,并且尽可能使得模拟观测数据消除模式依赖性问题。试验结果表明,变分法能够很好的反演出湍流黏性系数,而且能够克服一定的观测误差带来的影响,证明了方法的有效性和鲁棒性。
  然后,利用广义混沌多项式结合EnKF方法(Polynomial Chaos-ensemble Kalman filter,简记为PC-EnKF)反演GEM的低维的湍流黏性系数。EnKF中是利用蒙特卡罗法求得似然函数的,此方法收敛速度慢。在本文中,利用混沌多项式展开求解似然函数,通过选择合适的多项式基函数能够加快收敛速度,仅需较小的计算量即可得到似然值。数值试验的结果证明混沌多项式展开结合回归法求解系数能快速且有效地求解随机广义Ekman动量近似模式的风速分布,计算时间仅需蒙特卡罗法的1/50。在此基础上构建出了PC-EnKF反演GEM中参数的计算流程,并进行数值试验。试验结果表明PC-EnKF是一个非常有效的求解低维随机变量的湍流黏性系数后验概率分布的方法,而且在试验中发现在风速先验分布不确定性较大的区域内选择的观测数据在反演时会有更好的效果,这对于观测点的位置的选择具有指导意义。
  最后,解决了高维湍流黏性系数的不确定性量化以及反演问题。高维参数在进行混沌多项式展开时会面临“维数灾难”的问题,我们提出了一种新的求解风速分布的方法:基于投影重采样和切片逆回归的控制变量方法Projective Resampling inverse regression uncertainty quantification-based control variate method,简记为:PRIRUQ-CV),方法利用投影重采样结合切片逆回归方法求得线性降维空间,并使之与蒙特卡罗法相结合,从而避免“维数灾难”。将此方法应用到PC-EnKF中去,提出了新的反演方法PRIRCV-EnKF。利用数值试验验证方法的可行性,结果表明经过线性降维后的空间能够捕捉到原模式的主要信息,但是精度不足以求出风速的统计信息;与蒙特卡罗法结合后能够准确求出统计信息,且结果明显好于相同样本数的蒙特卡罗法。说明PRIRUQ-CV能够避免混沌多项式展开方法在高维时计算量太大的问题,且保持着优于蒙特卡罗法的特点。在此基础上我们验证了PRIRCV-EnKF方法反演湍流黏性系数的有效性,结果表明此方法能有效反演出广义Ekman动量近似模式的参数,在先验分布基础上大大减小参数不确定性。
其他文献
一个囚犯被关在空间非常狭小的牢狱里。有一天,牢房里飞进一只苍蝇,嗡嗡叫个不停,他手脚并用,东扑一下,西扑一下,努力捕捉。捉了很久,还是无法捉到它,最后他慨叹说:我的小囚房不小啊!居然连一只苍蝇都捉不到!  他悟出一个道理:心中有事世间小,心中无事一床宽。  派了3名核查员赶赴他的农场,重新进行核查。核查员花了近一天时间,认定戴维斯确实饲养了10214头奶牛,少发了他一头奶牛的财政补贴。两天后,农业
期刊
一只老钟对一只小钟说:你一年里要摆525600下啊。小钟吓坏了:哇,这怎么可能?我怎么能完成那么重的任务呢!这时候,另一只老钟笑着说:不用怕,你只需一秒钟摆一下,坚持下来就可以了。小钟高兴了:一秒钟摆一下好像并不难啊,试试看吧。果然,很轻松地就摆了一下。不知不觉一年过去了,小钟已经摆了525600下!事实上,当我们面对大困难的时候,只要根据实际,分期制定小目标,一一完成就行了。
期刊
想当年,我的朋友队伍庞大。渐渐发现,真能心心相印、赤胆忠心的实在不多。不少人聚在一起,差不多是图个热闹。朋友就像竹林一样,需要适时挖笋和砍伐,竹林才能长得旺盛,不能太密互吞养料。有些全然没有真诚实情的内涵,有的只是利用和攀援,这种别有用心的朋友,还不如大刀阔斧地砍掉。还有的人往往初会时热情如火,求人中好话说尽,事成后却不思量。互相利用的朋友一列车,还不如肝胆相照的挚友一人!交朋友是一个大浪淘沙的行
期刊
荷兰有个名叫戴维斯的大农场主,饲养了1万多头奶牛。按照政府规定,每饲养一头奶牛可获得70多欧元的财政补贴。  一天,戴维斯在查询补贴款额时发现少了73欧元,于是,驱车前往农业局询问。农业局立即指
期刊
女友让我手写6万字的情书。神回复:亲爱的,以后我们在一起的每一天,我都会给你一封写有“我爱你”三个字的情书,这样,等到从今天算起的第55年,你一定能读懂我那6万字的情深。
期刊
最近,著名作家陈忠实的逝世引来如潮纪念,人们用各种美好语言来悼念他的离去,或曰巨星陨落,大树倾覆,或曰“中国文坛的天空塌了一个角”,这位当代最优秀的小说家代表当得起各种礼赞。  面对这种热闹场面,我们不会忘记, 陈忠实写《白鹿原》,从构思到写成,差不多用了十年时间,坐了十年冷板凳。其中很长一段时间,他把自己关在陕北农村的陋室里,远离尘嚣,倾情投入。吃饭是冷锅冷灶,睡觉是冷衾冷炕,面对的冷笔冷墨,冷
期刊
启功先生以教书为生,但却以书法知名,有人评价他的字“外若飞仙,飘逸洒脱;内似硬汉,钢筋铁骨”。面对他的书作,总似有一股挥洒不尽的朝霞之气扑面而来,令人神清气爽,尘杂俱灭。问先生是如何练就这一手超凡脱俗的书法,他总会笑着讲起年少时的一桩尴尬往事,正是那犹如被人打耳光般的火辣辣的疼痛,改变了他的命运。  自幼时起,启功便迷上了绘画。那是受到了祖父的影响,他经常看到祖父拿着笔蘸上墨彩,在扇面上涂抹几笔,
期刊
著名表演艺术家李雪健真正走进观众的视线是在空政话剧团排演的多幕话剧《九一三事件》,他因出演林彪一角而名声大噪的,不仅外形酷似而且气韵夺人,并因此夺得戏剧最高奖“梅花奖”。  可他并非是这部戏林彪扮演者的最初理想人选,当时团里初定的是李保田,导演认为当时凭借李保田瘦小的体型和娴熟的表演技巧,应该是演林彪的不二人选。但当时李保田业务繁忙,一直没办法抽身露面,这可急坏了导演和化妆师。  李雪健当时也在这
期刊
军人军属的心灵笔记 复转老兵的军旅记忆 热心军迷的寻梦之旅 文学前卫的精神领地  《前卫文学》是济南军区政治部主办的军事文化类双月刊,拥有国际标准刊号和国内统一刊号,面向国内外公开发行,邮发代号24-260,单本定价10.00元,每份年价60.00元。本刊尚有少量往年余刊,您可以随时联系本刊发行部邮购订阅,办法:从邮局汇款至山东省济南市经十一路80号前卫文学杂志社,邮编:250002,联系电话:0
期刊
战斗中,一架敌机向阵地俯冲下来。上尉并没有立刻卧倒,因为他发现离他四五米远处有一个小战士还愣在那儿。他顾不上多想,一个鱼跃将小战士压在了身下。此时一声巨响,飞溅起来的泥土纷纷落在他们的身上。上尉拍拍身上的尘土,回头一看,刚才自己所处的那个位置被炸成了一个大坑。  其实,在帮助别人的同时也是在帮助自己!
期刊