《应用化学》31卷(2014年)投稿简则

<正>《应用化学》是经国家科委批准向国内外公开发行的学报类期刊。《应用化学》由中国科学院主管,中国化学会和中国科学院长春应用化学研究所主办,科学出版社出版。着重报道化学及交叉学科有应用前景的创新性基础科学研究和创造性科研技术成果,介绍该领域中的新发现、新理论、新方法、新技术、新产品及相关科技信息,为推动应用化学学科的     

广义非线性Schr(o)dinger方程的一个新的守恒差分格式

对广义非线性Schr(o)dinger方程提出了一种新的差分格式.揭示了该差分格式满足两个守恒律,并证明该格式的收敛性和稳定性.数值实验结果表明,新的差分格式优于Crank-Nicolson格式以及Zhang Fei等人提出的格式.     

二维变系数反应扩散方程的紧交替方向差分格式

1 引言 在研究热传导过程、气体扩散现象和电磁场的传播等问题时,常常遇到抛物型偏微分方程。用有限差分方法研究这类问题的数值解法目前已经有了许多工作。对于二维、三维抛物方程的数值求解比较理想的方法是交替方向法。     

一种基于块雅可比迭代的高阶FR格式隐式方法

最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.      

基于子结构解耦的连接特性识别新方法

连接部件动态特性的准确辨识对预测装配式机械结构的动力学行为有重要意义. 针对传统基于子结构解耦的连接结构动力学特性识别方法难以直接利用可测量数据进行辨识及辨识结果受噪声影响显著等问题, 本文提出了一种新方法. 首先, 提取子结构解耦基本方程在测试自由度上的分量, 并经矩阵变换得到显含连接动刚度矩阵的形式, 而后由真实连接动刚度矩阵分解为已知的初始矩阵与待求的增量矩阵, 推导了具有收敛性质的增量型方程以增强界面自由度较多时辨识的数值稳定性, 并采用多项式拟合动刚度将其转化为了拟合系数的频域估计方程, 按给定准则选取合适的频率点联立方程后, 得到了只需装配体测试自由度上的频响函数来辨识连接特性的迭代公式. 最后, 以若干算例说明了算法的具体流程. 对10自由度弹簧?质量块系统进行了数值仿真, 验证了所提方法的正确性及抗噪性; 对包含一处胶接连接的T形梁结构和包含两处螺栓连接的L形梁结构进行了试验, 所辨识连接结构与残余结构重组的装配体有限元模型计算的频响函数与测量值在较宽频带内吻合较好, 表明了该方法能有效识别实际装配体结构中的连接特性.      

改进的五阶WENO-Z+格式

WENO-ZWENO-Z$+\!$格式的性能提升依赖于新增项的作用,该项的作用是在WENO-Z格式的基础上进一步增大欠光滑子模板上的权重. 系数$\lambda$被设置用来调控该项的作用, 以避免负耗散. 本文指出了WENO-Z$+\!$格式的缺陷,其所采用$\lambda $的取值方式既不能充分发挥格式的潜力, 也未完全消除负耗散;提出$\lambda $的值应随当地流场数据变化,方能充分发挥新增项在降低数值耗散、提高分辨率上的潜力. 基于此,本文重新设计了$\lambda $的计算公式,该公式能自适应地调控新增项的作用: 只在欠光滑子模板上的权重容易过度增大的地方削弱该项的作用,以避免负耗散; 在其他地方则充分发挥新增项的作用,最大限度增大欠光滑子模板上的权重, 提高格式的分辨率.将使用该系数公式的新格式命名为WENO-Z++, 并对其数值性能进行了系统的研究.理论分析表明, 新格式在间断处具有基本无振荡(essentially non-oscillatory,ENO)特性和更低的数值耗散. 对近似色散关系(approximate dispersion relation,ADR)的研究表明,新格式有效地避免了因过度增大欠光滑子模板上的权重而带来的负耗散,其频谱特性也得到了显著提升.本文还推导了使新格式在极值点处也能保持最优阶的精度的参数设置.一系列求解Euler方程的数值试验表明,新格式的激波捕捉能力和对复杂流场结构的分辨率都显著优于原WENO-Z$+\!$格式.}     

一种提高极值点处收敛精度的三阶WENO-Z格式

为了提高三阶WENO-Z格式在极值点处的计算精度,通过理论推导给出三阶WENO格式满足收敛精度的充分条件。采用泰勒级数展开的方式,推导给出所构造格式非线性权重的计算公式,并综合权衡计算精度和计算稳定性确定所构造格式的参数。通过两个典型的精度测试,验证了改进格式在光滑流场极值点区域逼近三阶精度。进一步选用激波与熵波相互作用和Richtmyer-Meshkov不稳定性等经典算例,证实了本文提出的改进格式WENO-PZ3相较其他格式（WENO-JS3和WENO-Z3）不仅具有较高的精度,而且降低了格式的耗散,提高了对流场结构的分辨率。     

一种求解带底部浅水波方程组的平衡保正的MUSCL-Hancock中心迎风格式

本文改进了基于MUSCL方法求解一维浅水波方程组的中心迎风格式.利用MUSCL-Hancock方法在每个时间步的中间时刻求解Riemann问题,并分别在tn时刻和tn+1/2时刻对自由面的斜率进行了修正,得到了MUSCL-Hancock中心迎风格式.在合理的CFL条件下,我们证明了该格式的平衡性与保正性.数值结果表明,...     

N-S方程基于投影法的特征线算子分裂有限元求解

提出了一种求解非定常不可压缩纳维-斯托克斯方程（N-S方程）的新型有限元法:基于投影法的特征线算子分裂有限元法.在每一个时间层上将N-S方程分裂成扩散项、对流项、压力修正项.对流项采用多步显式格式,且在每一个对流子时间步内采用更加精确的显式特征线-伽辽金法进行时间离散,空间离散采用标准伽辽金法.应用此算法对平面泊肃叶流、方腔流和圆柱绕流进行数值模拟,所得结果与基准解符合良好.尤其对于Re=10000的方腔流,给出了方腔中分离涡发展和运动的计算结果,并发现在该雷诺数下存在周期解,表明该算法能较好地模拟流体流动中的小尺度物理量以及流场中分离涡的运动.      

基于拓扑描述函数的连续体结构拓扑优化方法

提出了一种利用拓扑描述函数(TDF)作为拓扑设计变量求解连续体结构拓扑优化问题 的新方法. 优化问题的目标函数是结构的整体柔顺性，约束条件为对于可利用材料的体积限 制. 这种方法不仅可以消除拓扑优化中经常出现的棋盘格式等数值不稳定现象，而且能够有 效地抑制传统算法处理此类优化问题时所引发的边界扩散效应. 与其它的基于水平集描述函 数的拓扑优化方法相比，所提出的算法不仅无需求解控制水平集函数演化的双曲守恒方 程，而且合理地考虑了目标函数的拓扑导数信息，因而使得算法的计算效率有了显著的提高.