微重环境中三维液体非线性晃动的数值模拟

将ALE（任意的拉格朗日-欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkin加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界件的弱积分形式。模拟了考虑表面张力情况下圆筒形贮腔中液体的非线性晃动,揭示了考虑表面张力效应时液体非线性晃动的重要特征。     

城市大气环境的大涡模拟研究进展

本文回顾城市大气污染扩散的研究实践与进展, 介绍城市大气环境流动的特点和用现代计算流体力学手段开展城市大气环境的大涡模拟研究方法, 包括: 数学模型, 控制方程、亚网格湍流模式、定解条件及其数值方法. 其中, 重点介绍中尺度到微尺度的耦合模型与适用于复杂城市下垫面(满足大涡模拟分辨率要求且计算量较小)的组合模型;并给出实际算例和结果分析, 包括准确度的统计估算和湍流特性等. 最后, 讨论城市大气环境数值模拟方法进一步改进的方向和城市大气环境流动与污染物扩散数值研究应重点关注的几个问题.     

水平应力对顶板跨落角的影响及应用

为得到水平应力对顶板垮落角的影响,进而确定高位抽采钻孔布孔参数,采用理论分析及相似材料模拟实验的方法进行了研究。首先依据基本顶断裂位置围岩结构力学模型,将顶板岩层假设为梁结构,得到单一岩梁在不同水平应力作用下的断裂线位置变化规律。然后将各岩梁的断裂线位置叠加组合,得到水平应力对垮落角发育的影响。为验证理论分析结果,以水平应力值为影响因素设计了三组不同水平应力值条件下的双向加载相似材料模拟实验,得到了不同水平应力值对应的垮落角。研究结果表明:顶板垮落角由各岩梁断裂线位置决定。地基系数较小时,断裂线与地基初始接触位置距离随水平应力的增加表现为先减小后增加,垮落角变化规律同样为先减小后增加,但减小不明显;当地基系数增加时,随水平应力增加,断裂线与地基初始接触位置的距离减小趋势逐渐消除,仅表现为增加趋势,垮落角变化同样表现为增加趋势。现场应用表明,采用双向加载模型得到的垮落角对钻孔布置参数进行优化,优化后钻场抽采浓度及高浓度抽采时间均得到提高。     

水下爆炸冲击波的数值模拟研究

应用商业有限元程序MSC.DYTRAN数值模拟了球形药包在无限水域中爆炸产生的冲击波。通过和经验公式计算结果的比较,证明采用合理的计算参数和有限元模型能够较好地模拟水下爆炸冲击波的传播过程。通过调整水的状态方程参数,达到了提升冲击波峰值应力的效果,从而有效地降低了单元数量和计算时间。     

不同围压作用下非均匀岩石水压致裂过程的数值模拟

从岩石细观非均匀性的特点出发，提出一个描述非均匀材料渗流和破裂相互作用的数值模型。在这个数值模型中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流，损伤耦合迭代来实现。算例表明，该模型能较好地模拟出岩石类材料在水力压裂作用下，微结构非均匀分布和不同围压比对破裂模式、失稳压力的影响，非均匀性导致试件的开裂压力、失稳压力明显不同，裂纹扩展路径不规则发展，模拟结果和实验结果较为一致。     

基于MATLAB的GPS软件接收机捕获与跟踪算法实现

研究了GPS软件接收机的捕获和跟踪算法,并基于Matlab软件平台和射频前端在PC上实现了GPS软件接收机样机。介绍了GPS软件接收机的结构和数据采集硬件,讨论了GPSC/A码的特性、产生原理以及捕获过程。针对传统的串行搜索算法慢的缺点以及高动态GPS软件接收机的特点,在该样机中实现了快速的基于循环卷积的并行捕获算法,并联合使用超前滞后环和对相位反转不敏感的科斯塔斯锁相环分别对码相位和多普勒频偏进行跟踪,解调得到导航电文。仿真和测试结果表明,使用GPS软件接收机进行信号处理的思想使用户在算法处理和软件升级等方面具有更大的灵活性,可应用于下一代任何全球导航卫星定位系统(GNSS)和空基增强系统(SBAS)接收机的设计。     

磁驱动飞片的一维磁流体动力学数值研究

磁驱动高速飞片技术是近年来发展的一种新型实验技术,在冲击波物理领域得到应用。该过程伴随着磁场扩散,并由此引起焦耳加热,使得飞片加载面的相状态发生变化,这决定了飞片厚度的范围。基于拉格朗日坐标系,利用磁流体动力学方程组、电阻率方程和状态方程数据库,对磁驱动铝飞片进行了一维磁流体动力学数值计算,获得了不同时刻铝飞片密度、温度的剖面分布,得到了磁场扩散速率随加载电流密度的变化关系。文章所选取的电导率方程只考虑到汽化点为止,对于等离子体形成的过程无法描述,如果要精确描述更高电流密度下的驱动过程,需考虑更为普适的电导率方程。磁场扩散速率随加载电流密度的变化存在转折点,在转折点前后可分别用两个线性关系表达式加以刻画。利用这些关系和冲击波物理相关知识,对磁压加载等熵驱动飞片实验样品厚度的选择进行了研究。     

入水后爆炸问题的数值模拟技术

应用各种数值计算方法及计算处理技术,编制程序实现了对弹体入水后爆炸问题全过程的数值模拟.其中弹体与水,爆轰产物气体与弹壳之间的相互作用通过流固耦合技术来描述;水面与空气、爆轰产物与水、空气之间的相互作用采用VOF方法(Volume-of-Fluid)来描述;采用了刚体-柔体转换、单元失效删除等计算技术以更高效、更好地模...     

小波变换在湍流数值研究中的应用

小波变换具有时空双局部性特点,恰好适应了湍流特性。本文主要阐述了小波在湍流数值计算中的两大研究进展:一个是利用连续小波(高斯小波)可使L ap lace算子降阶的特性来求解N-S方程;另一个是正交小波与有限元法相结合的方法——相干涡模拟。指出它们的优缺点及其存在的问题,并对小波在湍流计算中的应用前景作了展望。     

环形浅池内硅熔体中的热流体波

为了了解工业Czochralski炉内硅熔体表面轮型的基本特征,对环形浅池内硅熔体的热毛细-浮力对流进行了三维数值模拟,硅熔池内径为15 mm,外径为50 mm,深度为3 mm,熔池外壁被加热,内壁被冷却,底部固壁和顶部自由表面均绝热或者允许一个垂直方向的传热。模拟结果表明,当径向温差较小时,熔池内会产生稳定的单胞热毛细-浮力流动,随着温差的增大,流动将转变为三维振荡流动,在熔体自由表面会出现沿周向运动的轮型,小的垂直方向的热流密度(3W/cm2)对这种振荡流动没有大的影响。同时,讨论了流动和温度波动的特征,并确定了振荡流动的临界条件。