吸积过程的数值模拟与外流
一、研究外流的意义
1.黑洞系统的外流很普遍,解释这些外流的起源本身就是一个很重要的课题。
2.外流是活动星系核反馈的重要媒介。
补充
(1)外流对黑洞吸积盘的影响:比如说带走黑洞吸积盘的质量和能量,影响其密度以及光谱。
(2)活动星系核:河外天体中的一类中央核区活动性很强的系外星系的核心
(3)活动星系核反馈:黑洞的产出、辐射和外流会影响星际介质的性质,而黑洞的吸积原料来自于星际介质,当星际介质发生改变时,又会影响黑洞的增长,它们是一个相互反馈的。
外流如此重要,就要搞清楚其产生的机制、性质等等。
二、黑洞热吸积流外流的研究历史以及现状
黑洞热吸积流模型
(1)吸积流的能量方程:
其中,表示气体密度,T表示气体温度,S表示气体的熵,t表示时间,表示粘滞加热率,表示辐射制冷率。
(2)热吸积流的特点:温度高、盘厚、角动量小(气体压可部分抵抗引力)、径向速度大、吸积率低、光盘薄。
(3)热吸积流是一个吸积率系统的吸积模型。因为热吸积流的密度很低,辐射无效,辐射制冷可以忽略。此时,只有粘滞加热,粘滞加热产生的热量全部储存在气体中,作为气体的内能,径移到黑洞里面去,因此吸积流或者吸积盘的温度就很高;温度高,吸积流有热压支撑,因此其吸积流很厚;由于有热压的支撑,吸积流就不需要开普勒转动,所以其径向下落速度就很大,其角动量比较低,且吸积率很低,径向速度很大,所以热吸积流光盘很薄。
1.1对流主导的吸积流模型(CDAF)
1.吸积流的能量方程:
熵的增加等于加热率减去制冷率。当加热率大于制冷率时,吸积流的熵随半径的减小(吸积的进行)而增加;但加热率小于制冷率时,吸积流的熵随着半径的减小(吸积的进行)而减小。
2.热吸积流辐射无效,因此吸积率远大于制冷率,吸积流的熵随着半径的减小(吸积的进行)而增加。
3.对于无转动的吸积流,熵随着半径减小而增大时,吸积流是对流不稳定的(吸积流中会产生对流运动)。
4.对于转动的吸积流,产生对流变得困难:
无转动时:气体压平衡引力
有转动时,气体压和离心力平衡引力,气体压较小,产生对流更加困难。
5.基于一维的热吸积流的解析解,Narayan等人发现,热吸积流是对流不稳定的,吸积流中会产生对流。Narayan等人构建了对流主导的吸积流模型(convection dominated accretion flow;CDAF)。
这两种吸积流最大的区别:
(1)CDAF是对流主导的
(2)CDAF的吸积率是半径的函数,吸积率随着吸积的进行逐渐降低:
1.2绝热内流-外流模型(ADIOS)
1.伯努利参数(动能+焓+引力能):
伯努利参数大于0,气体不受黑洞引力束缚,有可能变成外流。
热吸积流辐射无效,因此,其伯努利参数可能大于0。基于此,Blandford&Payne 1999构建了绝热的内流-外流模型(adiabatic inflow-outflow solution,ADIOS)。
气体在向黑洞下落的过程中,不断产生外流。由于外流的产生,吸积率随着吸积的进行逐渐降低:
此模型只是假设外流存在
1.3CDAF模型与ADIOS模型的比较
(1)CDAF无外流,ADIOS有外流
(2)两个模型中,吸积率都随半径的减少而降低
1.4粒子轨迹线方法
(1)在数值模拟中,撒入“虚拟粒子”,虚拟粒子完全跟着流体元运动,通过追踪虚拟粒子的轨迹,来找出真正的外流。
(2)找到真正的外流后,即可计算其各自的物理性质。
(3)外流的质量流随半径的增大而增大:
(4)外流主要来自于大半径处
外流的驱动机制
(1)克服引力
(2)对于热吸积流,辐射无效,辐射压可忽略。
热吸积流外流的驱动
(1)热吸积流外流被气体压强梯度力、磁压梯度力和离心力共同驱动。
a.气体接近维理温度,气体压力接近引力,因此,气体压强梯度力重要。
b.磁压接近气体压力,因此,磁压梯度力重要。
c.气体转动稍稍亚开普勒,因此,离心力重要。
(2)热吸积流喷流,被磁压梯度力驱动。
1.5热吸积流外流的质量流到底有多大
(1)Yuan,Bu&Wu(2012)通过将“内”,“外”流的性质进行比较,证明热吸积流存在系统性的真正外流,但无法确定,正速度的流体中,真外流和湍流成分的比例。 Yuan…Bu(2015),发现外流质量流随半径的增大而增大
但还是无法给出外流的质量流到底有多大。
(2)Bondi半径:黑洞引力和气体内能相等的地方
a.物质以光速进入黑洞,而气体声速远小于光速,因此,黑洞附近气体的引力势能(绝对值)主导气体内能。
b.随着半径的增大(距离黑洞越大越远),引力势能(绝对值)越来越小,黑洞引力越来越弱。而气体总是有一定内能,因此,随着半径的增大,黑洞的引力势能的绝对值在某个半径将会等于气体内能。这个半径成为Bondi半径。
c.Bondi半径以外,引力势能(绝对值)继续减小,而气体内能几乎保持不变。Bondi半径外气体内能大于其引力能。
(3)Bondi半径一般认为是黑洞吸积系统的外边界。
(4)前面数值模拟的外边界远小于Bondi半径,因此公式:
中的“r”到底可以取多大值仍然未可知,前面的工作无法具体给出外流的质量流的定量值。
(5)解决这个问题需要开展大尺度的数值模拟,模拟的外边界需要大于Bondi半径。
(6)对于星系中心的超大黑洞的吸积系统来说,Bondi半径外,星系的恒星引力超过黑洞引力,数值模拟需要将黑洞引力和恒星引力统一考虑
(7)观测发现星系中恒星的弥散速度随半径的变化很小,因此,可以认为恒星弥散速度不随半径变化。
(8)弥散速度代表了恒星的引力,弥散速度不随半径变化,也就意味着恒星的开普勒转动速度不随半径变化。开普勒速度不随半径变化,说明恒星的引力反比于半径平方(而黑洞引力反比于半径的平方)。这样恒星的引力势能如下:
其中为恒星的弥散速度。不同星系其中心恒星的弥散速度不同,一般为几百公里每秒。典型值为200km/s.
黑洞的引力比恒星引力陡。
恒星引力和黑洞引力相对的地方和Bondi半径位置基本相同。
因此,Bondi半径外,恒星引力主导;Bondi半径以内,黑洞引力主导。
吸积过程的数值模拟与外流相关推荐
- wps 模拟分析 规划求解_【图书推荐】金属塑性加工过程有限元数值模拟及软件应用...
科学出版社 梅瑞斌 著 <金属塑性加工过程有限元数值模拟及软件应用>共6章,第1章主要讲述金属塑性加工工艺及求解方法:第2章主要讲述塑性力学及有限元理论基础:第3章分别采用有限元商用软件A ...
- 焊接过程计算机模拟研究,焊接过程的数值模拟
<焊接过程的数值模拟>是2017年科学出版社出版的图书,作者是金成. 书 名 焊接过程的数值模拟 别 名 金成出版社 科学出版社 ISBN 9787030551030 焊接过程 ...
- 复合材料专场 | ABAQUS车载四型复合材料气瓶固化过程的数值模拟分析攻略
复合材料气瓶固化的热场本质上可以认为包含两个阶段,复合材料气瓶表面和周围通过与空气的热对流换热,复合材料与内部塑料芯模以及金属接头的传热.在第一个阶段整体温度较低,热量从表面向复合材料层内部流入,此时 ...
- comsol计算机模拟过程,借助数值模拟分析多孔结构
储层.水坝.及其他户外结构都必须达到坚固.可靠的要求.压力变化可能会导致这些结构中的多孔材料被损坏,进一步引起流体流动及结构的逐渐垮塌和下沉.借助 COMSOL Multiphysics的多物理场仿真 ...
- 计算机cfd教学,CFD数值模拟过程
<CFD数值模拟过程>由会员分享,可在线阅读,更多相关<CFD数值模拟过程(8页珍藏版)>请在人人文库网上搜索. 1.CFD数值模拟过程,CFD 数值模拟简介 CFD软件简介 ...
- 计算机模拟在铸造领域的应用,铸造过程计算机模拟研究应用现状与发展
<铸造过程计算机模拟研究应用现状与发展>由会员分享,可在线阅读,更多相关<铸造过程计算机模拟研究应用现状与发展(3页珍藏版)>请在人人文库网上搜索. 1.铝合金铸造数值模拟的研 ...
- 深入Linux内核网络堆栈
前一段时间看到这篇帖子,确实很经典,于是翻出了英文原版再读,顺便再翻译出来供大家学习,这篇文章的中文版也早都有了,不过出于完全理解的目的,我还是将它翻译了出来,加进了自己的代码,虽然在上一周的翻译过程 ...
- 地理科学与计算机技术的关系,地理信息系统与地理学的关系
地理信息系统与地理学的关系表现在: (1)地理学是我们理解世界的基础科学,GIS使得地理科学活生生地应用到现实世界中,包括科学中的基础部分. (2)地理学与GIS密不可分,两者形成完美结合.地理学的进 ...
- 高斯旋转热源与双椭球热源_【干货】论焊接,3D打印模拟的热源模型——焊缝、3D打印高度变换模拟...
原标题:[干货]论焊接,3D打印模拟的热源模型--焊缝.3D打印高度变换模拟 焊接.3D打印是难度相对较高的数值模拟过程,在模拟温度场及应力场过程中,至关重要的是如何进行热源函数的建立和加载,同时配以 ...
最新文章
- Windows查找JDK的路径
- 0x00000000指令引用的内存不能为written_JVM03——对象实例化,内存布局,访问定位...
- 微软自带的防反编译工具dotfuscator.exe的使用
- SQL Server启动的几种方法
- websphere一直安装部署_构建独立部署系统,从部署规范开始
- 使用N2N,实现异地快速组建内网
- c++11并发与多线程
- Javascript特效:五彩小球
- 三十岁那年,我的梦想是年薪十万
- AV1代码学习:解码端aomdec.c的main函数
- Windows server 2016 安装sql server
- 高通 Camx debug log控制
- android 类似qq空间微博微信九宫格图片
- 奥鹏计算机应用基础考前,华东2020年秋季学期《计算机应用基础》在线考试
- DES加密算法的C++实现
- 张丽俊最新演讲:要像竹子一样扎根,你终会一飞冲天
- matlab读取声音文件
- vol.132 《流浪地球》带火了中国科幻?(嘉宾:郑军)
- 数据库课程设计:图书借阅系统(Java+MySQL)
- 10个优秀的Spring Boot开源项目整理分享
热门文章
- 泡在网上的日子 ---*****----开源项目*******--新功能
- 银行的east系统是什么
- Vue.js学习笔记—input-number:实战:开发一个数字输入框组件
- html+css+js之20个练手小项目(一)
- Mac M1 (Apple Silicon) CH340 驱动
- 支持Linux系统的加密狗
- wf工作流java_WF Workflow 状态机工作流 开发
- How to use Clang Static Analyzer
- 华为机试 经验(血泪史啊)
- 06 “eden没有发生minor gc, 对象直接分配在了old gen“ 的调试