×

一种求解制导中心问题的后向半拉格朗日方法的有效轨迹跟踪算法。 (英语) Zbl 07506195号

总结:在本文中,我们开发了一种有效的数值算法,该算法可以跟踪使用反向半拉格朗日方法求解制导中心模型所需的轨迹。在数值计算方面,设计了两种适用于出发点的快速算法。一个是每个柯西问题的离散系统数值解的完全显式公式。该公式的特点是数值因子小于通常高斯消去中使用的乘法数的一半。另一种方法使用插值方法找到所需的出发点,对于严重的柯西问题,该方法的成本至少降低30%。最后,我们提出了一种改进解决方案的方法,以改进物理量的估计,例如质量守恒,这些物理量可能会在出发点计算的插值解中丢失。结果表明,该方法不仅节省了大量的计算时间,而且比传统方法更好地保留了质量和总动能等物理量。为了证明数值证据,我们使用该方法模拟了几个问题,如不可压缩Euler方程、Kelvin-Helmholtz不稳定性、Diocotron不稳定性和三维导向中心模型。

理学硕士:

70-XX岁 粒子和系统力学
76倍 流体力学

软件:

FFTW公司;艾根
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Bermúdez,A。;Nogueiras,M.R。;Vázquez,C.,具有高阶特征的对流-扩散-反应问题的数值分析/有限元。第一部分:时间离散化,SIAM J.Numer。分析。,44, 1829-1853 (2006) ·Zbl 1126.65080号
[2] Bermúdez,A。;Nogueiras,M.R。;Vázquez,C.,具有高阶特征的对流-扩散-反应问题的数值分析/有限元。第二部分:全离散格式和求积公式,SIAM J.Numer。分析。,44, 1854-1876 (2006) ·Zbl 1126.65081号
[3] Boukir,K。;Maday,Y。;Métivet,B.,不可压缩Navier-Stokes方程的高阶特征线方法,计算。方法应用。机械。工程,116211-218(1994)·Zbl 0824.76060号
[4] Boukir,K。;Maday,Y。;Métivet,B。;Razandrakoto,E.,《不可压缩Navier-Stokes方程的高阶特征/有限元方法》,国际期刊Numer。液体方法,251421-1454(1997)·Zbl 0904.76040号
[5] 蔡晓凤。;郭,W。;邱,J.M.,半拉格朗日间断Galerkin格式线性和非线性输运模拟的比较,arXiv预印本
[6] 蔡,X.F。;郭,W。;邱,J.M.,二维不可压Euler方程和无算子分裂的Guding中心Vlasov模型的高阶半拉格朗日间断Galerkin方法,J.Sci。计算。,79, 1111-1134 (2019) ·Zbl 1419.65055号
[7] 北卡罗来纳州克鲁塞尔斯。;拉图,G。;Sonnendrücker,E.,用于并行求解Vlasov-Poisson方程的贴片上的Hermite样条插值,国际数学杂志。计算。科学。,17, 335-349 (2007) ·Zbl 1159.65016号
[8] 法尔科内,M。;Ferretti,R.,一类高阶半拉格朗日平流格式的收敛性分析,SIAM J.Numer。分析。,35909-940(1998年)·Zbl 0914.65097号
[9] 费尔贝特,F。;Proufur,C.,反向半拉格朗日方法的高阶时间离散化,J.Compute。申请。数学。,303, 171-188 (2016) ·Zbl 1334.65149号
[10] Fourestey,G。;Piperno,S.,一种应用于风工程和桥梁剖面控制的二阶时间精度ALE-Lagrange-Galerkin方法,计算。方法应用。机械。工程,193,4117-4137(2004)·Zbl 1175.76101号
[11] 弗里戈,M。;Johnson,S.,FFTW3的设计与实现,Proc。IEEE,93,2,216-231(2005)
[12] Guennebaud,G。;Jacob,B.,Eigen v3(2010年)
[13] Kelvin,L.,《流体动力学溶液和观察》,Philos。Mag.,42,362-377(1871)
[14] Kim,P。;Kim,D。;Piao,X。;Bak,S.,用于求解Navier-Stokes方程的BSLM中Cauchy问题的完全显式格式,J.Compute。物理。,401,第109028条pp.(2020)·Zbl 1453.76034号
[15] Kim,P。;Piao,X。;Kim,S.D.,基于切比雪夫配置解决刚性问题的误差修正欧拉方法,SIAM J.Numer。分析。,49, 2211-2230 (2011) ·Zbl 1273.65095号
[16] Knauer,W.,等离子体和气体放电中的二极管不稳定性,J.Appl。物理。,37, 2, 602-611 (1996)
[17] Piao,X。;Kwon,J。;Yi,D。;Kim,S.D。;Kim,P.,指导中心问题的无迭代向后半拉格朗日格式,SIAM J.Numer。分析。,53619-643(2015)·Zbl 1318.65054号
[18] Piao,X。;Bu,S。;贝克·S。;Kim,P.,求解不可压缩Navier-Stokes方程的无迭代反向半拉格朗日格式,J.Compute。物理。,283, 189-204 (2015) ·Zbl 1351.76179号
[19] Piao,X。;Kim,P。;Kim,D.,向后半拉格朗日格式的一步\(L(\alpha)\)稳定时间积分及其在指导中心问题中的应用,J.Comput。物理。,366327-340(2018)·Zbl 1406.65098号
[20] Pollock,D.S.G.,《关于Kronecker积、张量积和矩阵微分学》,国际计算机杂志。数学。,90, 2462-2476 (2013) ·Zbl 1282.15022号
[21] Notsu,H。;Tabata,M.,不可压缩Navier-Stokes方程的二阶时间单步特征曲线有限元格式,J.Sci。计算。,38, 1-14 (2009) ·Zbl 1203.76086号
[22] 邱建明。;Christlieb,A.,Vlasov方程的保守高阶半拉格朗日WENO方法,J.Compute。物理。,229, 1130-1149 (2010) ·Zbl 1188.82069号
[23] 邱建明。;Shu,C.W.,不可压缩流平流的守恒高阶半拉格朗日有限差分WENO方法,J.Compute。物理。,230, 863-889 (2011) ·Zbl 1391.76489号
[24] Robert,A.,原始气象方程的稳定数值积分格式,Atmos-《海洋》,19,35-46(1981)
[25] 鲁伊·H。;Tabata,M.,对流扩散问题的二阶特征有限元格式,数值。数学。,92, 161-177 (2002) ·Zbl 1005.65114号
[26] Smolarkiewicz,P.K。;Pudykiewicz,J.A.,流体的一类半拉格朗日近似,J.Atmos。科学。,49, 22, 2082-2096 (1992)
[27] Sonnendrücker,E。;罗氏,J。;Bertrand,P。;Ghizzo,A.,Vlasov方程数值求解的半拉格朗日方法,J.Compute。物理。,149201-220(1999年)·Zbl 0934.76073号
[28] Staniforth,A。;Cóté,J.,《大气模型的半拉格朗日积分方案——综述》,周一。《天气评论》,1192206-2223(1991)
[29] 坦普顿,C。;Staniforth,A.,一个有效的两时间级半拉格朗日半隐式积分方案,Q.J.R.Meteorol。Soc.,1131025-1039(1987)
[30] 托达,K。;Ogata,Y。;Yabe,T.,浅水方程特征CIP的多维保守半拉格朗日方法,J.Comput。物理。,2284917-4944(2009年)·Zbl 1166.76045号
[31] Wiin-Nielsen,《轨道方法在数值预报中的应用》,Tellus XI,2180-196(1959)
[32] 秀,D。;Karniadakis,G.,Navier-Stokes方程的半拉格朗日高阶方法,J.Compute。物理。,172, 658-684 (2001) ·Zbl 1028.76026号
[33] 秀,D。;Sherwin,S.J。;Dong,S。;Karniadakis,G.E.,《不可压缩流动的半拉格朗日方法的强形式和辅助形式》,《科学杂志》。计算。,25, 323-346 (2005) ·Zbl 1203.76122号
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不声称其完整性或完全匹配。