CN111611649A - 一种列车客室病原体传播机理数值仿真方法 - Google Patents
一种列车客室病原体传播机理数值仿真方法 Download PDFInfo
- Publication number
- CN111611649A CN111611649A CN202010392390.7A CN202010392390A CN111611649A CN 111611649 A CN111611649 A CN 111611649A CN 202010392390 A CN202010392390 A CN 202010392390A CN 111611649 A CN111611649 A CN 111611649A
- Authority
- CN
- China
- Prior art keywords
- train
- flow field
- model
- discrete
- transmission mechanism
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/80—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for detecting, monitoring or modelling epidemics or pandemics, e.g. flu
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Public Health (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Automation & Control Theory (AREA)
- Mathematical Analysis (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Pathology (AREA)
- Aviation & Aerospace Engineering (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种列车客室病原体传播机理数值仿真方法,包括:获得列车三维模型;将列车三维模型导入网格离散软件,对列车三维模型中的列车客室内流场空间和列车客室内饰进行网格划分,获得列车内流场离散模型;将列车内流场离散模型导入CFD仿真软件并进行边界条件设置,获得列车内流场数学计算模型;针对列车内流场数学计算模型,设定内流场送风模式,在列车客室内流场的气体连续相达到稳定之后,利用离散相模型从液体入口处加入液体离散相,在内流场的不同位置注入不同粒径的水滴,开启蒸发模型,实时追踪水滴粒子的轨迹及水滴粒子的粒径变化。本发明通过数值仿真方法可得到不同条件下的客室内带病原体微团传播机理,耗时短,成本低。
Description
技术领域
本发明属于试验方法领域,特别涉及一种列车客室病原体传播机理数值仿真方法。
背景技术
当前,病原体微团在封闭列车客室内的传播机理还没有完全掌握,全部通过试验去完成传播机理研究代价太大,因此进一步通过数值仿真的方法来研究病原体的传播机理是十分有必要的。
列车客室内的环境相对是一个比较特殊的环境,因此病原体在列车客室内的传播机理与飞机、邮轮等交通工具有所不同。另外对于不同型号的列车,其病原体的传播机理也存在着一定的差异,因此在数值仿真病原体的传播机理时,应该系统地研究不同型号和不同工况下运行的列车。
目前,现有技术中没有关于高速列车客室内病原体传播机理的仿真方法,一般是通过试验方法获得列车客室内病原体传播机理,耗时长,成本高。
因此,为了快速经济的对列车客室内病原体微团传播机理进行研究,亟待进行相关的仿真工作。
发明内容
本发明的目的在于,针对现有技术中使用试验方法研究列车客室内病原体传播机理耗时长、成本高的不足,提供一种列车客室病原体传播机理数值仿真方法,耗时短,成本低。
为解决上述技术问题,本发明所采用的技术方案是:
一种列车客室病原体传播机理数值仿真方法,其特点是包括以下步骤:
步骤一,对列车进行三维建模,获得列车三维模型;
步骤二,将列车三维模型导入网格离散软件,并利用网格离散软件对列车三维模型中的列车客室内流场空间和列车客室内饰进行网格划分,获得列车内流场离散模型;
步骤三,将从网格离散软件中导出的列车内流场离散模型导入 CFD仿真软件,并在CFD仿真软件内对列车内流场离散模型进行边界条件设置,获得列车内流场数学计算模型;
步骤四,针对列车内流场数学计算模型,设定内流场送风模式,在列车客室内流场的气体连续相达到稳定之后,利用离散相模型从液体入口处加入液体离散相,在内流场的不同位置注入不同粒径的水滴,开启蒸发模型,实时追踪水滴粒子的轨迹及水滴粒子的粒径变化。
进一步地,步骤四中,在内流场的不同位置注入不同粒径的水滴后,对列车内流场数学计算模型中的扩散项采用二阶精度的中心差分格式,对列车内流场数学计算模型中流项离散的差分格式采用具有二阶精度的QUICK格式。
作为一种优选方式,步骤四中,应用雷诺时均N-S方程和RNG k-ε湍流模型计算不同水滴粒子粒径大小、不同产生位置的水滴在列车客室内的运动规律。
作为一种优选方式,步骤三中,边界条件设置规则为:内流场进口采用速度进口条件,内流场出口设为压力出口,内流场近壁面处采用标准壁面函数,内流场固壁面上的各节点采用绝热无滑移壁面边界条件。
作为一种优选方式,步骤四中的送风模式为新风模式、回风风量变化模式、内循环模式、无回风模式、冬季送风模式或夏季送风模式。
作为一种优选方式,所述步骤4中,在列车客室内流场的气体连续相达到稳定之后,首先,通过UDF程序调用DPM飞沫离散相模型,在内流场的不同位置注入不同粒径的水滴模拟病原体微团,并形成高速列车内流场飞沫运动模型;然后,通过UDF程序控制蒸发相关参数并开启蒸发模型,进入飞沫蒸发模式,并且通过UDF程序跟踪记录飞沫蒸发特性和飞沫蒸发规律。
借由上述方法,本发明将飞沫蒸发模型与高速列车客室内流场数值仿真计算结合,高度还原液态飞沫在客室复杂环境中,在蒸发相关参数(如气流、温度、气压等)作用下的运动、蒸发特性与蒸发规律。
进一步地,所述步骤4中还包括,飞沫蒸发形成飞沫核后,利用 UDF程序连续追踪飞沫核运动,计算飞沫和飞沫核的运动轨迹。
作为一种优选方式,所述蒸发相关参数包括气流、温度、气压中的一种或多种。
与现有技术相比,本发明通过数值仿真方法可得到不同条件下的客室内带病原体微团传播机理,从而可避免大量相应的试验研究,对早日理清列车客室内带病原体微团传播机理,提高列车客室内环境的公共卫生安全等级研究有很大的时间和经济效益。
附图说明
图1为水滴粒子的轨迹及水滴粒子的粒径变化追踪原理图。
具体实施方式
以型号为CR400BF的中国标准动车组为例,对本发明所述的列车客室病原体传播机理数值仿真方法做进一步说明如下:
本发明所述的列车客室病原体传播机理数值仿真方法包括以下步骤:
步骤一,利用三维CAD软件依据真实车型全尺寸数据对列车进行三维建模,整车长L=25.0m,车宽W=3.4m,车高H=4m。从三维CAD软件导出列车三维模型。
步骤二,将列车三维模型导入网格离散软件,并利用网格离散软件,采用非结构(密集)/笛卡尔网格方法对列车三维模型中的列车客室内流场空间和列车客室内饰(为了真实模拟车内情况,对车内的座椅等设施也要进行精细的网格划分)进行非结构网格离散划分。通过非结构网格划分,获得接近于真实情况的列车内流场离散模型。离散后的非结构网格总数为1.2×108。
步骤三,将从网格离散软件中导出的列车内流场离散模型导入 CFD仿真软件ANSYS Fluent,并在CFD仿真软件内对列车内流场离散模型进行边界条件设置,获得列车内流场数学计算模型。边界条件设置规则为:内流场进口采用速度进口条件,内流场出口设为压力出口,内流场近壁面处采用标准壁面函数,内流场固壁面上的各节点采用绝热无滑移壁面边界条件。从而将列车内流场问题简化为不可压缩、非定常流动来处理。
步骤四,针对列车内流场数学计算模型,设定内流场送风模式,在列车客室内流场的气体连续相达到稳定之后,利用离散相模型从液体入口处加入液体离散相,在内流场的不同位置注入不同粒径的水滴 (1微米-400微米),模拟带病原体微团,通过不同粒径的水滴模拟带病原体微团,然后开启蒸发模型,实时追踪水滴粒子的轨迹及水滴粒子的粒径变化。对于离散相,由于气液体积比达到1000∶1,可通过增加射流压强提高液体横向射流贯穿深度。
离散相模型可以模拟不同大小的水滴(并且水滴之间没有相互关联,独自运动)在空间运动的规律,从而模拟仿真得出病原体的传播机理。
步骤四中,在内流场的不同位置注入不同粒径的水滴后,对列车内流场数学计算模型中的扩散项采用二阶精度的中心差分格式,对列车内流场数学计算模型中流项离散的差分格式采用具有二阶精度的 QUICK格式。
步骤四中,应用雷诺时均N-S方程和RNG k-ε湍流模型计算不同水滴粒子粒径大小、不同产生位置的水滴在列车客室内的运动规律。在CFD软件中,基于雷诺时均N-S方程和RNG k-ε湍流模型建立客室内流场的三维数学计算模型,时间迭代方法采用非稳态算法,且基于非正交曲线坐标下求解可压缩SIMPLE算法。
步骤四中的送风模式为新风模式、回风风量变化模式、内循环模式、无回风模式、冬季送风模式或夏季送风模式。
利用离散相模型从液体入口处加入液体离散相,在内流场的不同位置注入不同粒径的水滴;在CFD软件中设置各种送风模式均属于现有技术,在此不作赘述。
对于特殊车型,需要仿真计算冬季送风模式和夏季送风模式。
基于本发明所述的数值仿真方法,可获得不同送风模式下的液滴运动规律,得到不同条件下的列车客室内带病原体微团的传播机理。
如图1所示,所述步骤4中,在列车客室内流场的气体连续相达到稳定之后,首先,通过UDF程序调用DPM飞沫离散相模型,在内流场的不同位置注入不同粒径的水滴模拟病原体微团,并形成高速列车内流场飞沫运动模型;然后,通过UDF程序控制蒸发相关参数 (气流、温度、气压)并开启蒸发模型,进入飞沫蒸发模式,并且通过UDF程序跟踪记录飞沫蒸发特性和飞沫蒸发规律。飞沫蒸发形成飞沫核后,利用UDF程序连续追踪飞沫核运动,计算飞沫和飞沫核的运动轨迹。
本发明将飞沫蒸发模型与高速列车客室内流场数值仿真计算结合,高度还原液态飞沫在客室复杂环境中,在气流、温度、气压共同作用下的运动、蒸发特性与蒸发规律。同时,针对飞沫蒸发后形成飞沫核的现象,具备时空连续追踪计算飞沫与飞沫核运动轨迹的功能。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是局限性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护范围之内。
Claims (8)
1.一种列车客室病原体传播机理数值仿真方法,其特征在于,包括以下步骤:
步骤一,对列车进行三维建模,获得列车三维模型;
步骤二,将列车三维模型导入网格离散软件,并利用网格离散软件对列车三维模型中的列车客室内流场空间和列车客室内饰进行网格划分,获得列车内流场离散模型;
步骤三,将从网格离散软件中导出的列车内流场离散模型导入CFD仿真软件,并在CFD仿真软件内对列车内流场离散模型进行边界条件设置,获得列车内流场数学计算模型;
步骤四,针对列车内流场数学计算模型,设定内流场送风模式,在列车客室内流场的气体连续相达到稳定之后,利用离散相模型从液体入口处加入液体离散相,在内流场的不同位置注入不同粒径的水滴,开启蒸发模型,实时追踪水滴粒子的轨迹及水滴粒子的粒径变化。
2.如权利要求1所述的列车客室病原体传播机理数值仿真方法,其特征在于,步骤四中,在内流场的不同位置注入不同粒径的水滴后,对列车内流场数学计算模型中的扩散项采用二阶精度的中心差分格式,对列车内流场数学计算模型中流项离散的差分格式采用具有二阶精度的QUICK格式。
3.如权利要求1所述的列车客室病原体传播机理数值仿真方法,其特征在于,步骤四中,应用雷诺时均N-S方程和RNG k-ε湍流模型计算不同水滴粒子粒径大小、不同产生位置的水滴在列车客室内的运动规律。
4.如权利要求1所述的列车客室病原体传播机理数值仿真方法,其特征在于,步骤三中,边界条件设置规则为:内流场进口采用速度进口条件,内流场出口设为压力出口,内流场近壁面处采用标准壁面函数,内流场固壁面上的各节点采用绝热无滑移壁面边界条件。
5.如权利要求1所述的列车客室病原体传播机理数值仿真方法,其特征在于,步骤四中的送风模式为新风模式、回风风量变化模式、内循环模式、无回风模式、冬季送风模式或夏季送风模式。
6.如权利要求1所述的列车客室病原体传播机理数值仿真方法,其特征在于,所述步骤4中,在列车客室内流场的气体连续相达到稳定之后,首先,通过UDF程序调用DPM飞沫离散相模型,在内流场的不同位置注入不同粒径的水滴模拟病原体微团,并形成高速列车内流场飞沫运动模型;然后,通过UDF程序控制蒸发相关参数并开启蒸发模型,进入飞沫蒸发模式,并且通过UDF程序跟踪记录飞沫蒸发特性和飞沫蒸发规律。
7.如权利要求6所述的列车客室病原体传播机理数值仿真方法,其特征在于,所述步骤4中还包括,飞沫蒸发形成飞沫核后,利用UDF程序连续追踪飞沫核运动,计算飞沫和飞沫核的运动轨迹。
8.如权利要求6所述的列车客室病原体传播机理数值仿真方法,其特征在于,所述蒸发相关参数包括气流、温度、气压中的一种或多种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010392390.7A CN111611649A (zh) | 2020-05-11 | 2020-05-11 | 一种列车客室病原体传播机理数值仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010392390.7A CN111611649A (zh) | 2020-05-11 | 2020-05-11 | 一种列车客室病原体传播机理数值仿真方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111611649A true CN111611649A (zh) | 2020-09-01 |
Family
ID=72200102
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010392390.7A Pending CN111611649A (zh) | 2020-05-11 | 2020-05-11 | 一种列车客室病原体传播机理数值仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111611649A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765905A (zh) * | 2021-01-08 | 2021-05-07 | 西南交通大学 | 一种用于检测通风管道病毒消杀质量的方法 |
CN113191095A (zh) * | 2021-04-02 | 2021-07-30 | 南京航空航天大学 | 一种发动机吞水风险的评估方法及装置 |
CN113361215A (zh) * | 2021-06-21 | 2021-09-07 | 西南交通大学 | 一种亚音速真空管道列车周围气流参数计算方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080570A1 (en) * | 2003-09-15 | 2005-04-14 | Acosta Edward P. | Predicting probabilities of achieving a desired minimum trough level for an anti-infective agent |
CN101650221A (zh) * | 2009-09-14 | 2010-02-17 | 中南大学 | 铁路列车车内噪声测量分析方法及其系统 |
CN106573118A (zh) * | 2014-06-30 | 2017-04-19 | Syqe医药有限公司 | 用于活性剂的肺部递送的方法、装置及系统 |
US20170190764A1 (en) * | 2011-09-12 | 2017-07-06 | New York University | Peptides mimicking hiv-1 viral epitopes in the v2 loop for the gp120 surface envelope glycoprotein |
CN107220482A (zh) * | 2017-05-09 | 2017-09-29 | 清华大学 | 呼吸道传染病风险评估系统及评估方法 |
CN109470496A (zh) * | 2018-10-11 | 2019-03-15 | 中南大学 | 列车车体瞬态剧烈振动致振动舒适性的评估方法及系统 |
CN109753716A (zh) * | 2018-12-28 | 2019-05-14 | 山东大学 | 基于流场仿真的核/火电汽轮机组流体激励数值计算方法及系统 |
-
2020
- 2020-05-11 CN CN202010392390.7A patent/CN111611649A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080570A1 (en) * | 2003-09-15 | 2005-04-14 | Acosta Edward P. | Predicting probabilities of achieving a desired minimum trough level for an anti-infective agent |
CN101650221A (zh) * | 2009-09-14 | 2010-02-17 | 中南大学 | 铁路列车车内噪声测量分析方法及其系统 |
US20170190764A1 (en) * | 2011-09-12 | 2017-07-06 | New York University | Peptides mimicking hiv-1 viral epitopes in the v2 loop for the gp120 surface envelope glycoprotein |
CN106573118A (zh) * | 2014-06-30 | 2017-04-19 | Syqe医药有限公司 | 用于活性剂的肺部递送的方法、装置及系统 |
CN107220482A (zh) * | 2017-05-09 | 2017-09-29 | 清华大学 | 呼吸道传染病风险评估系统及评估方法 |
CN109470496A (zh) * | 2018-10-11 | 2019-03-15 | 中南大学 | 列车车体瞬态剧烈振动致振动舒适性的评估方法及系统 |
CN109753716A (zh) * | 2018-12-28 | 2019-05-14 | 山东大学 | 基于流场仿真的核/火电汽轮机组流体激励数值计算方法及系统 |
Non-Patent Citations (5)
Title |
---|
XUE-YI YOU 等: "The Study of Pathogen Particles Spread around Buildings by Considering the Pathogen Survival Time in Air", 《2009 2ND INTERNATIONAL CONFERENCE ON BIOMEDICAL ENGINEERING AND INFORMATICS》 * |
ZHOU HAI-PING: "Delay included dynamic virus spreading process simulating", 《2013 25TH CHINESE CONTROL AND DECISION CONFERENCE (CCDC)》 * |
余勇: "基于DPM的客车室内生物源性污染物扩散及控制仿真研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
姬鹏 等: "高速列车过隧道压力波动数值模拟中的计算网格及湍流模型研究", 《JOURNAL OF CENTRAL SOUTH UNIVERSITY》 * |
陈小明: "空调室内生物气溶胶扩散分布的数值模拟", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765905A (zh) * | 2021-01-08 | 2021-05-07 | 西南交通大学 | 一种用于检测通风管道病毒消杀质量的方法 |
CN112765905B (zh) * | 2021-01-08 | 2022-09-02 | 西南交通大学 | 一种用于检测通风管道病毒消杀质量的方法 |
CN113191095A (zh) * | 2021-04-02 | 2021-07-30 | 南京航空航天大学 | 一种发动机吞水风险的评估方法及装置 |
CN113361215A (zh) * | 2021-06-21 | 2021-09-07 | 西南交通大学 | 一种亚音速真空管道列车周围气流参数计算方法 |
CN113361215B (zh) * | 2021-06-21 | 2022-02-18 | 西南交通大学 | 一种亚音速真空管道列车周围气流参数计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111611649A (zh) | 一种列车客室病原体传播机理数值仿真方法 | |
CN104268943B (zh) | 一种基于欧拉‑拉格朗日耦合方法的流体仿真方法 | |
CN107284442B (zh) | 一种用于自动驾驶车辆的弯道行驶纵向控制方法 | |
CN112215350B (zh) | 一种基于强化学习的智能体控制方法及装置 | |
CN109709956A (zh) | 一种自动驾驶车辆速度控制多目标优化的跟驰算法 | |
CN114168796B (zh) | 一种建立飞行器高空气动力数据库的方法 | |
CN108983605B (zh) | 一种基于深度强化学习进行流体导向的刚体控制的方法 | |
CN107729613A (zh) | 飞行模拟机数据包架构的构建方法 | |
CN113191095A (zh) | 一种发动机吞水风险的评估方法及装置 | |
Uddin et al. | Artificial intelligence-based Monte-Carlo numerical simulation of aerodynamics of tire grooves using computational fluid dynamics | |
CN103049931A (zh) | 基于粒子系统的云团快速三维仿真算法 | |
CN104504170A (zh) | 车辆动画仿真方法及系统 | |
Pakkam | High Downforce Aerodynamics for Motorsports. | |
CN202869772U (zh) | 气流模拟实验装置 | |
CN104133969B (zh) | 一种真空羽流实验系统液氦深冷低温泵吸附捕获面积计算方法 | |
CN117078876A (zh) | 一种三维流场模拟方法、系统及电子设备 | |
CN116644666A (zh) | 基于策略梯度优化算法的虚拟装配路径规划引导方法 | |
CN105183965A (zh) | 用于预测雾化过程的大涡模拟方法 | |
CN108804791A (zh) | 一种适用于埋入式进气道布局的飞行器参数化方法 | |
CN104951592A (zh) | 一种基于cfd技术优化车体外形的分析方法 | |
CN106934119A (zh) | 一种利用轮地接触压力评价行车振动噪声的方法 | |
De Vanna et al. | Hypersonic intake design using a CFD data-driven multi-objective optimisation strategy | |
JP2021179992A (ja) | コンピュータ支援設計定義のジオメトリにおけるメッシュ空洞空間特定および自動シード検出 | |
CN106919050B (zh) | 多旋翼无人机高速开环动作自适应学习方法 | |
Kobayashi et al. | CFD application in automotive industry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200901 |
|
RJ01 | Rejection of invention patent application after publication |