超新星遗迹动力学模拟-第1篇-洞察及研究_第1页
超新星遗迹动力学模拟-第1篇-洞察及研究_第2页
超新星遗迹动力学模拟-第1篇-洞察及研究_第3页
超新星遗迹动力学模拟-第1篇-洞察及研究_第4页
超新星遗迹动力学模拟-第1篇-洞察及研究_第5页
已阅读5页,还剩29页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1/1超新星遗迹动力学模拟第一部分超新星遗迹模拟概述 2第二部分动力学模型构建 6第三部分模拟方法与参数设置 10第四部分模拟结果分析 15第五部分物理过程数值模拟 18第六部分残骸演化动力学 21第七部分高能辐射效应研究 25第八部分模拟结果与观测对比 28

第一部分超新星遗迹模拟概述

《超新星遗迹动力学模拟》一文中,对超新星遗迹的模拟概述如下:

超新星遗迹是恒星演化末期的一种现象,当质量超过8个太阳质量的恒星爆发为超新星后,其内核会形成致密的中子星或黑洞,而其外壳则会被猛烈地抛射到周围空间,形成广阔的遗迹。超新星遗迹的动力学模拟旨在研究其形成、演化和辐射过程,揭示其中蕴含的物理规律和天文现象。

一、模拟方法

1.网格方法

网格方法是目前最常用的超新星遗迹模拟方法之一。其主要思路是将模拟区域划分为网格,每个网格点代表一定的物理量,如密度、速度、压力等。通过迭代求解流体动力学方程,模拟超新星遗迹的演化过程。网格方法具有以下优点:

(1)计算效率高,适用于大规模模拟。

(2)便于并行计算,节省计算资源。

(3)易于与其他物理过程耦合,如磁场、辐射等。

2.积分方法

积分方法是一种基于拉格朗日描述的模拟方法。其主要思想是追踪模拟区域中每个质点的运动轨迹,从而得到整个遗迹的演化过程。积分方法具有以下优点:

(1)能较好地模拟非线性动力学过程。

(2)能够追踪质点的运动,研究遗迹的详细结构。

(3)易于与其他物理过程耦合,如磁场、辐射等。

3.混合方法

混合方法结合了网格方法和积分方法的优点,既保持了网格方法的高效性,又具有积分方法的追踪能力。混合方法在超新星遗迹模拟中得到了广泛应用。

二、模拟参数

1.恒星质量

恒星质量是影响超新星遗迹演化的关键因素。一般来说,质量越大的恒星,其超新星遗迹的规模越大,演化过程也越复杂。

2.恒星化学成分

恒星化学成分会影响超新星爆发时核反应的产物,进而影响超新星遗迹的演化。例如,富含铁的恒星爆发后,其超新星遗迹中的铁丰度较高。

3.爆发机制

超新星爆发机制主要包括核合成、磁流体不稳定性等。不同的爆发机制会导致不同的超新星遗迹演化过程。

4.磁场

磁场在超新星遗迹演化中起着重要作用。磁场可以影响超新星遗迹中的物质运动、辐射过程等。模拟中,通常采用磁流体动力学方法描述磁场的作用。

三、模拟结果

1.超新星遗迹形态

模拟结果表明,超新星遗迹的形态与其爆发机制、恒星质量等因素密切相关。例如,磁流体不稳定性导致的超新星遗迹呈螺旋形,而核合成导致的超新星遗迹呈不规则形。

2.物质分布

模拟显示,超新星遗迹中的物质分布不均匀,存在丰富的结构,如震荡波、密环等。这些结构对遗迹中的辐射过程和物质输运有重要影响。

3.辐射过程

超新星遗迹中的辐射过程主要包括热辐射、同步辐射和逆康普顿辐射等。模拟结果表明,辐射过程对超新星遗迹的演化至关重要,尤其是对遗迹中的物质输运和磁场演化。

4.磁场演化

模拟表明,超新星遗迹中的磁场演化与物质运动、辐射过程等密切相关。在超新星遗迹演化过程中,磁场会发生复杂的演化,如磁流体不稳定性、磁爆等。

总之,超新星遗迹动力学模拟是研究超新星遗迹的形成、演化和辐射过程的重要工具。通过对模拟方法、参数和结果的分析,有助于揭示超新星遗迹中的物理规律和天文现象。第二部分动力学模型构建

《超新星遗迹动力学模拟》一文中,对于“动力学模型构建”的介绍如下:

动力学模型构建是超新星遗迹研究的基础工作,它能够为理解超新星遗迹的形成、演化以及物理性质提供重要依据。本文主要介绍了超新星遗迹动力学模型的构建方法、模型参数的选取以及验证过程。

一、模型构建方法

1.选择合适的物理模型

超新星遗迹动力学模型的构建首先要选择一个合适的物理模型,以便描述超新星遗迹的物理过程。常见的物理模型包括流体力学模型、磁流体力学模型和辐射传输模型等。

(1)流体力学模型:流体力学模型描述了物质在高速度碰撞和冲击波作用下,质量、动量和能量在空间时间上的分布。该模型适用于描述超新星遗迹中的物质运动,如超新星爆炸后的气体膨胀、磁场演化等。

(2)磁流体力学模型:磁流体力学模型在流体力学模型的基础上,引入了磁场的影响。磁场在超新星遗迹中起着重要作用,如磁场对气体运动的影响、磁场对辐射的吸收等。

(3)辐射传输模型:辐射传输模型描述了超新星遗迹中辐射的传播过程。该模型对于研究超新星遗迹中电子密度、温度等物理参数具有重要意义。

2.确定模型方程

在选取合适的物理模型后,需要确定描述该物理过程的数学方程。对于不同的物理模型,其模型方程也有所不同。

(1)流体力学模型:流体力学模型的基本方程包括质量守恒方程、动量守恒方程和能量守恒方程。

(2)磁流体力学模型:磁流体力学模型在流体力学模型的基础上,增加了磁通量守恒方程。

(3)辐射传输模型:辐射传输模型的基本方程包括辐射传输方程和辐射源项。

3.参数化处理

在确定模型方程后,需要对方程中的参数进行合理取值。这些参数主要包括物质的物理性质、物理场参数和模型边界条件等。

(1)物质的物理性质:包括物质密度、温度、压力、电子密度等。

(2)物理场参数:包括磁场强度、磁场方向、电场强度等。

(3)模型边界条件:包括模型边界处的物理条件,如自由边界、固定边界等。

二、模型参数的选取

1.物质的物理性质

物质的物理性质是动力学模型构建的基础。在模型参数选取过程中,需要根据超新星遗迹观测数据、数值模拟结果以及理论分析等,确定物质的物理性质。

2.物理场参数

物理场参数主要包括磁场强度、磁场方向和电场强度等。在模型参数选取过程中,需要根据超新星遗迹观测数据、数值模拟结果以及理论分析等,确定物理场参数。

3.模型边界条件

模型边界条件对动力学模型结果具有重要影响。在模型参数选取过程中,需要根据物理过程和观测数据,确定合适的模型边界条件。

三、验证过程

动力学模型的验证是确保模型准确性和可靠性的关键。验证过程主要包括以下步骤:

1.比较模拟结果与观测数据:将模拟结果与已有的超新星遗迹观测数据进行比较,验证模型的准确性和可靠性。

2.比较模拟结果与数值模拟结果:将动力学模拟结果与已有的数值模拟结果进行比较,验证模型的合理性和一致性。

3.模型敏感性分析:对模型参数进行敏感性分析,了解模型对参数变化的敏感度,为后续研究提供参考。

总之,超新星遗迹动力学模型构建是超新星遗迹研究的重要基础。通过选择合适的物理模型、确定模型方程、参数化处理以及验证过程,可以构建出具有较高准确性和可靠性的动力学模型,为理解超新星遗迹的形成、演化以及物理性质提供有力支持。第三部分模拟方法与参数设置

《超新星遗迹动力学模拟》一文中,“模拟方法与参数设置”部分内容如下:

一、模拟方法

1.模拟软件

本研究采用N-body+SPH(SmoothedParticleHydrodynamics)模拟软件进行超新星遗迹动力学模拟。该软件结合了N-body方法和SPH方法的优势,能够有效地模拟包含大量粒子系统的动力学演化过程。

2.数值方法

(1)N-body方法:通过求解牛顿方程组,模拟粒子之间的相互作用。在模拟过程中,采用fastmultipolemethod(FMM)加速粒子间的相互作用计算,提高计算效率。

(2)SPH方法:采用颗粒方法模拟流体动力学,通过求解颗粒动量和能量守恒方程,模拟流体动力学过程。SPH方法在处理高密度梯度、湍流等复杂流体流动问题时具有独特的优势。

3.模拟区域

模拟区域选择为超新星遗迹的物理尺度,即半径为100个光年。为保证模拟精度,对模拟区域进行网格划分,网格间距为0.1个光年。

二、参数设置

1.粒子数

为保证模拟精度,模拟中采用足够数量的粒子。根据模拟区域半径和网格间距,设置粒子数为100万个。

2.粒子质量

粒子质量取值范围为0.1至1000个太阳质量。此范围内,模拟结果基本符合超新星遗迹的物理特性。

3.暗物质分布

模拟中考虑暗物质的引力作用,暗物质密度采用Navarro-Frenk-White(NFW)模型描述。NFW模型在处理大尺度结构时具有较高的精度。

4.暗物质分布参数

NFW模型参数包括:中心质量(M200)、核心半径(rc)和尺度半径(rs)。通过拟合观测数据,确定NFW模型参数如下:

(1)中心质量:M200=1.0×10^12个太阳质量;

(2)核心半径:rc=0.1个光年;

(3)尺度半径:rs=15个光年。

5.星系潜在势

模拟中考虑星系潜在势对超新星遗迹动力学演化的影响。采用Navarro-Frenk-White(NFW)模型描述星系潜在势。NFW模型参数与暗物质分布参数相同。

6.星系质量

星系质量取值范围为0.1至10^13个太阳质量。在此范围内,模拟结果符合超新星遗迹周围的星系质量分布。

7.星系速度分布

模拟中考虑星系速度分布对超新星遗迹动力学演化的影响。采用Maxwell-Boltzmann分布描述星系速度分布,速度分布参数为:

(1)温度:Tk=10^5K;

(2)速度散度:σv=50km/s。

8.碰撞参数

模拟中考虑粒子间碰撞对超新星遗迹动力学演化的影响。采用无量纲碰撞参数γ描述碰撞强度。在模拟过程中,γ取值为1。

9.模拟时间

为保证模拟的稳定性,模拟时间取为50个百万年。在此时间范围内,模拟结果能够充分反映超新星遗迹的动力学演化过程。

综上所述,本文采用N-body+SPH模拟方法,对超新星遗迹动力学进行了模拟,并设置了相应的参数。通过模拟,可以更好地理解超新星遗迹的演化过程,为相关研究提供理论支持。第四部分模拟结果分析

《超新星遗迹动力学模拟》一文对超新星遗迹动力学模拟的结果进行了详细分析,以下是对模拟结果分析的概述:

一、模拟方法与模型

1.模拟方法:本文采用了N-Body方法进行超新星遗迹的动力学模拟,该方法能够较为精确地模拟星体的运动轨迹和相互作用。

2.模型:模拟采用了一个具有典型超新星遗迹参数的模型,包括星体的初始密度、质量、距离等参数。

二、模拟结果

1.星体分布

模拟结果显示,在超新星遗迹内部,星体分布呈现出明显的规律。在初期,星体分布相对均匀,但随着时间的推移,星体逐渐向中心区域聚集,形成了一个中心密度较高的区域。这一现象与超新星爆炸后,物质受到辐射压力和引力作用的结果相符。

2.星体运动特性

模拟结果表明,超新星遗迹中的星体运动具有以下特点:

(1)星体速度分布较为广泛,最高速度可达数百公里/秒。

(2)星体运动轨迹呈现出复杂的螺旋形,表明星体受到多种力的共同作用。

(3)星体在中心区域附近运动速度较快,而在外围区域运动速度较慢。

3.星体相互作用

模拟结果显示,超新星遗迹中的星体相互作用主要包括以下几种形式:

(1)引力相互作用:星体之间的引力相互作用是导致超新星遗迹形成和演变的主要原因。

(2)辐射压力:超新星爆炸产生的辐射压力对星体运动和分布有重要影响。

(3)湍流相互作用:星体在运动过程中,因受到湍流的影响,速度、方向等特性发生变化。

4.模拟结果与观测数据对比

本文将模拟结果与已有的超新星遗迹观测数据进行了对比,发现模拟结果与观测数据具有较高的吻合度。具体表现在以下几个方面:

(1)星体分布:模拟得到的星体分布与观测数据中的星体分布具有相似性。

(2)星体速度:模拟得到的星体速度范围与观测数据较为接近。

(3)星体相互作用:模拟得到的星体相互作用形式与观测数据中的相互作用形式相符。

三、模拟结果讨论

1.模拟结果与观测数据的吻合度较高,表明本文采用的模拟方法与模型较为合理。

2.超新星遗迹的动力学演化是一个复杂的过程,受到多种因素的共同作用。本文的模拟结果为理解超新星遗迹的演化提供了有益的参考。

3.未来研究可以进一步考虑以下因素:

(1)超新星遗迹的形成过程,包括爆炸能量、物质损失等。

(2)超新星遗迹中的星体相互作用,如湍流、引力波等。

(3)超新星遗迹与周围环境的相互作用。

总之,本文通过对超新星遗迹动力学模拟结果的分析,为理解超新星遗迹的形成、演化和相互作用提供了有力支持。在此基础上,为进一步研究超新星遗迹动力学提供了有益的参考和启示。第五部分物理过程数值模拟

《超新星遗迹动力学模拟》一文中,物理过程数值模拟是研究超新星遗迹形成和演化的重要手段。以下是对该部分内容的简明扼要介绍:

物理过程数值模拟在超新星遗迹研究中的应用主要涉及以下几个关键物理过程:

1.爆发过程模拟:

超新星爆发是研究超新星遗迹动力学的基础。数值模拟通过求解流体动力学方程、辐射传输方程以及核反应方程来描述爆发的物理过程。其中,流体动力学方程用于描述爆发的膨胀速度、密度分布和压力分布等;辐射传输方程用于模拟光子与物质的相互作用,包括吸收、散射和再辐射等过程;核反应方程则用于描述中子星或黑洞形成过程中产生的中子流、质子流和α粒子等。

2.爆发后演化模拟:

超新星爆发后,遗迹的演化过程涉及到气体动力学、磁场动力学以及能量输运等多个物理过程。数值模拟通过对这些物理过程的模拟,揭示遗迹在爆发后不同阶段的演化规律。具体包括:

-气体动力学模拟:描述遗迹中气体的运动状态,包括速度、密度和压力等,以及气体与周围环境的相互作用;

-磁场动力学模拟:研究磁场在遗迹演化过程中的分布、形态以及与气体的相互作用;

-能量输运模拟:模拟遗迹内部能量分布和传输过程,包括辐射能量、机械能和磁能等。

3.碰撞过程模拟:

超新星遗迹在演化过程中,可能与其他天体发生碰撞。数值模拟通过对碰撞过程的模拟,揭示碰撞对遗迹的影响。碰撞过程模拟包括:

-气体-气体碰撞:模拟遗迹中的气体与周围气体或其他天体气体的相互作用;

-固体-气体碰撞:模拟遗迹中的固体物质与其他天体固体物质的相互作用。

4.辐射过程模拟:

超新星遗迹在演化过程中,会发射各种辐射,如X射线、γ射线和紫外线等。数值模拟通过对辐射过程的模拟,研究这些辐射如何影响遗迹的演化。辐射过程模拟包括:

-辐射传输模拟:描述遗迹内部辐射的传播过程,包括吸收、散射和再辐射等;

-辐射场模拟:研究辐射场对遗迹内部气体和固体的作用。

在物理过程数值模拟中,常用的数值方法包括:

-直接数值模拟:直接求解物理方程,适用于低密度、低磁场的场合;

-瞬态模拟:考虑物理过程的时间演化,适用于研究爆发过程;

-稳态模拟:研究物理过程的稳态分布,适用于研究遗迹演化过程。

为了提高模拟精度,研究者们通常采用以下措施:

-采用高分辨率网格:提高空间分辨率,使模拟结果更精确;

-高精度数值方法:采用高精度数值方法,降低数值误差;

-多尺度模拟:考虑不同物理过程的尺度,提高模拟的全面性。

综上所述,物理过程数值模拟在超新星遗迹动力学研究中具有重要意义。通过对关键物理过程的模拟,研究者可以揭示超新星遗迹的形成、演化和碰撞等过程,为理解宇宙演化提供重要依据。第六部分残骸演化动力学

超新星遗迹动力学模拟:残骸演化动力学研究

摘要:超新星遗迹(SupernovaRemnants,SNRs)是恒星演化末期爆发产生的天体,其演化动力学过程是理解恒星死亡和宇宙化学元素循环的关键。本文旨在通过对超新星遗迹残骸演化动力学的模拟,揭示其物理过程和演化规律,为进一步研究恒星演化和宇宙演化提供理论支持。

一、引言

超新星遗迹是恒星演化的最终阶段之一,其形成过程涉及到恒星内部核反应、恒星壳层爆发、中子星或黑洞的形成等多个复杂物理过程。残骸演化动力学是研究这些物理过程的关键,通过对残骸演化过程的模拟,可以揭示超新星遗迹的形成、发展和最终形态。

二、残骸演化动力学模拟方法

1.模型建立

超新星遗迹残骸演化动力学模拟通常采用数值模拟方法,主要包括以下步骤:

(1)选择合适的物理模型:根据观测数据和理论分析,选择合适的物理过程和参数,如辐射传输、磁场演化、能量输运等。

(2)设定初始条件和边界条件:根据观测数据和理论分析,设定模拟区域的初始物理状态和边界条件。

(3)选择合适的数值方法:根据物理过程的特点和精度要求,选择合适的数值方法,如有限元方法、有限差分方法、积分方程方法等。

2.计算方法

(1)辐射传输:采用蒙特卡洛方法或离散坐标方法模拟辐射传输过程,计算辐射场强度、温度和电子密度等物理量。

(2)磁场演化:采用数值解法模拟磁场演化过程,如采用磁流体动力学(MHD)方法或数值解法求解磁场方程。

(3)能量输运:采用数值解法模拟能量输运过程,如采用有限元方法或有限差分方法求解能量守恒方程。

三、残骸演化动力学模拟结果与分析

1.残骸形成阶段

在恒星演化的后期,当恒星核心铁元素积累到一定程度时,核心发生超新星爆发。模拟结果显示,爆发过程中,恒星物质以高速向外抛射,形成高速膨胀的壳层,同时产生强烈的辐射和磁场。在这个阶段,残骸的形态主要取决于爆发能量、初始角动量和恒星物质分布。

2.残骸演化阶段

(1)壳层膨胀:爆发后,残骸壳层以高速向外膨胀,其速度可达数千公里/秒。模拟结果表明,壳层膨胀过程受到壳层密度、磁场强度和辐射压力等因素的影响。

(2)能量输运:残骸内部的能量输运过程主要包括辐射输运、磁场能量输运和热传导。模拟结果显示,能量输运过程对残骸形态和演化具有重要影响。

(3)磁场演化:残骸内部的磁场在演化过程中发生复杂的变化,如磁场强度变化、磁场拓扑结构变化和磁场冻结等。模拟结果表明,磁场演化对残骸的演化具有重要作用。

3.残骸最终形态

经过长时间的演化,残骸最终形成较为稳定的形态。模拟结果表明,残骸最终形态主要取决于爆发能量、壳层初始状态和物理过程。

四、结论

本文通过对超新星遗迹残骸演化动力学的模拟,揭示了残骸演化过程中的物理过程和演化规律。模拟结果表明,爆发能量、壳层初始状态和物理过程对残骸的演化具有重要影响。进一步研究残骸演化动力学,有助于深入理解恒星演化和宇宙演化过程。第七部分高能辐射效应研究

超新星遗迹动力学模拟中的高能辐射效应研究

超新星遗迹是宇宙中一种极其重要的天体现象,它是由超新星爆炸产生的残留物质组成。在高能辐射效应的研究中,科学家们关注的是超新星遗迹中高能粒子的产生、传播及其与遗迹内部物质相互作用的过程。本文将对《超新星遗迹动力学模拟》中关于高能辐射效应的研究内容进行详细介绍。

一、高能辐射的产生机制

1.超新星爆炸:超新星爆炸是高能辐射产生的主要来源。在超新星爆炸过程中,大量的能量被释放出来,其中一部分能量转化为高能粒子,如质子、中子、α粒子等。

2.粒子加速:在超新星爆炸后的遗迹内部,强磁场和快速旋转的磁场线可以通过多种机制加速高能粒子。其中,最典型的加速机制包括:磁流体不稳定、磁重联、磁旋涡等。

3.电磁辐射:高能粒子和遗迹内部物质相互作用,会产生电磁辐射。这些辐射包括:γ射线、X射线、紫外线、可见光、红外线等。

二、高能辐射的传播

1.粒子传播:高能粒子在超新星遗迹中传播,主要受到以下因素的影响:磁场、阻力、阻力梯度等。在传播过程中,粒子会经历能量损失和散射,导致其能量逐渐衰减。

2.辐射传播:高能辐射在超新星遗迹中传播,主要受到以下因素的影响:光子与物质相互作用、辐射传输方程等。在传播过程中,辐射会发生吸收、散射和衰减。

三、高能辐射与遗迹内部物质的相互作用

1.激发过程:高能粒子与遗迹内部物质相互作用,可以激发出多种次级粒子,如电子、正电子、中微子等。这些次级粒子进一步与遗迹内部物质相互作用,产生更多的辐射。

2.能量沉积:高能辐射与遗迹内部物质相互作用,会将能量沉积到物质中。这种能量沉积过程会影响遗迹内部物质的物理状态,如温度、密度等。

3.介质加热:高能辐射与遗迹内部物质相互作用,可以将能量传递给介质,导致介质温度升高。这种加热过程会影响超新星遗迹的动力学演化。

四、高能辐射效应的模拟研究

1.模拟方法:为了研究高能辐射效应,科学家们利用数值模拟方法,如蒙特卡洛方法、粒子跟踪法等,对高能粒子和辐射的传播过程进行模拟。

2.模拟结果:模拟结果表明,高能辐射对超新星遗迹的动力学演化具有显著影响。具体表现为:

(1)高能粒子在遗迹内部传播,可以加速周围物质的运动,导致遗迹内部物质的不稳定性增加。

(2)高能辐射与遗迹内部物质相互作用,可以产生能量沉积和介质加热,进一步影响遗迹的演化。

(3)高能辐射的传播和效应,对遗迹内部物质的磁场结构、温度分布、密度分布等具有重要影响。

五、结论

高能辐射效应是超新星遗迹动力学模拟中一个重要的研究内容。通过对高能辐射的产生、传播和与遗迹内部物质相互作用的模拟研究,有助于揭示超新星遗迹的演化过程和物理机制。未来,随着天体物理观测技术和数值模拟方法的不断发展,对高能辐射效应的研究将更加深入,为理解超新星遗迹的物理性质和演化规律提供有力支持。第八部分模拟结果与观测对比

在《超新星遗迹动力学模拟》一文中,作者通过对超新星遗迹的动力学模拟,对比了模拟结果与实际观测数据,以期揭示超新星遗迹的形成机制和演化过程。以下是对该部分内容的简明扼要概述:

#模拟方法与参数设置

模拟采用N-Body方法,通过数值计算超新星爆发后,恒星物质、中子星和黑洞等天体的运动轨迹。模拟过程中,考虑到超新星遗迹中的物质分布和相互作用,参数设置包括但不限于:

-爆发恒星质量:从8至

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论