CN110083866A - 一种建筑物内部水淹的模拟系统及方法 - Google Patents
一种建筑物内部水淹的模拟系统及方法 Download PDFInfo
- Publication number
- CN110083866A CN110083866A CN201910232763.1A CN201910232763A CN110083866A CN 110083866 A CN110083866 A CN 110083866A CN 201910232763 A CN201910232763 A CN 201910232763A CN 110083866 A CN110083866 A CN 110083866A
- Authority
- CN
- China
- Prior art keywords
- water logging
- water
- building
- interior
- test
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 191
- 238000000034 method Methods 0.000 title claims abstract description 79
- 238000004088 simulation Methods 0.000 title claims abstract description 46
- 238000013178 mathematical model Methods 0.000 claims abstract description 40
- 239000002245 particle Substances 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 23
- 239000004744 fabric Substances 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims description 60
- 239000012530 fluid Substances 0.000 claims description 13
- 230000001133 acceleration Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 8
- 238000012805 post-processing Methods 0.000 claims description 7
- 238000013139 quantization Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 239000000523 sample Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 239000011521 glass Substances 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000001035 drying Methods 0.000 claims description 3
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 claims description 3
- 238000010998 test method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 12
- 238000002474 experimental method Methods 0.000 abstract description 11
- 238000011160 research Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 4
- 230000007812 deficiency Effects 0.000 abstract description 2
- 239000007788 liquid Substances 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种建筑物内部水淹的模拟系统及方法,该方法首先根据研究工况开展典型情况物理模型试验;基于SPH方法进行水动力数学模型的建模,采用物理模型试验结果对其进行率定和验证;然后采用经过验证的数学模型开展水淹事件数值模拟计算;分析和处理数值计算结果,获得水淹特征以及水体运动过程三维演示动画。本发明采用无网格的光滑粒子动力学方法模拟水淹三维动态过程,该方法适于计算液面大变形问题,避免了网格方法在水体翻转、破碎等复杂情况下追踪自由液面时易于产生失真的不足。核岛厂房建筑物结构复杂,计算量大,该方法通过并行计算解决计算效率问题,实现了核岛厂房建筑物内部水淹的精确计算。
Description
技术领域
本发明属于核电站核岛厂房设计技术,具体为一种针对核岛厂房所进行的建筑物内部水淹模拟系统及方法。
背景技术
核岛内部水淹是系统和设备的一种共模故障,核电站设计必须考虑这种潜在危害。包括福岛事件中因外部水淹导致的内部水淹也引起了国家核安全局和行业专家的重视。以我国目前在役的M310主力堆型为例,其内部水淹防护设计虽从确定论的角度进行了分析,但是缺少水淹发生后水体在建筑物内的漫延路径以及水淹高度的分析。因此,需要对水淹模型的精确度、数据准确性等方面进行试验验证,严格按照厂房实际条件进行水淹数值模拟,准确反映水淹时间和高度,为核岛内部水淹的分析研究提供数据支持及保障。
SPH(Smoothed Particle Hydrodynamics)是光滑粒子流体动力学方法的缩写,是在近20多年来逐步发展起来的一种无网格方法。该方法的基本思想是将连续的流体(或固体)用相互作用的质点组来描述,各个物质点上承载各种物理量,包括质量、速度等,通过求解质点组的动力学方程和跟踪每个质点的运动轨道,求得整个系统的力学行为。由于SPH方法中质点之间不存在网格关系,因此它可避免极度大变形时网格扭曲而造成的精度破坏等问题,并且也能较为方便的处理不同介质的交界面。
发明内容
本发明的目的在于提供一种针对核岛厂房所进行的基于SPH算法的建筑物内部水淹模拟系统及方法,并采用物理试验结果对数学模型进行验证,从而精确得出内部水淹分析的精确数据。
本发明的技术方案如下:一种建筑物内部水淹的模拟系统,包括:
物理试验模型,用于根据模拟条件开展物理试验,作为数学模拟结果的验证依据;
水动力建模模块,用于根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;
验证模块,用于根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
计算模块,用于采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算;
后处理模块,用于根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,并制作三维动态演示。
进一步,如上所述的建筑物内部水淹的模拟系统,其中,所述物理试验模型是依据建筑物结构特点和内部水淹源信息,按重力相似准则进行设计,结合试验结构尺度及设备性能确定模型长度比尺,设置满足试验条件的流量、粗糙度进行淹水试验;可选用有机玻璃作为模型材料,物理试验设备包括水流量控制设备和超声波水位测量系统,所述水流量控制设备包括蓄水池、水泵、电磁流量计、压力机、调剂管路。
进一步,如上所述的建筑物内部水淹的模拟系统,其中,所述水动力建模模块在建模中依据的建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯;水淹源信息包括水淹源的流量大小及水流方向。
进一步,如上所述的建筑物内部水淹的模拟系统,其中,所述计算模块在配备专业GPU科学运算加速设备的计算机上采用并行技术完成水淹事件模拟计算。
进一步,如上所述的建筑物内部水淹的模拟系统,其中,所述后处理模块采用ParaView软件直接读取计算模块的计算结果,并制作成三维演示动画。
一种建筑物内部水淹的模拟方法,包括如下步骤:
(S1)根据模拟条件开展物理试验,作为数学模拟结果的验证依据;
(S2)根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;
(S3)根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
(S4)采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算;
(S5)根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,并制作三维动态演示。
进一步,如上所述的建筑物内部水淹的模拟方法,步骤(S1)中物理试验针对每组试验工况进行多次重复,试验步骤如下:
(S1-1)调试供水系统,达到泄漏点出水要求;
(S1-2)根据泄漏点位置,架设管路;
(S1-3)在水淹主要路径的房间内调整水位测量探头;
(S1-4)对水位测量系统进行零点标定;
(S1-5)架设摄像装置,确定拍摄内容;
(S1-6)开始试验,记录不同房间内水位变化、最大水淹高度、水淹路径;
(S1-7)排空模型内水体,清洁、烘干模型,并返回步骤(S1-1)重复试验。
进一步,如上所述的建筑物内部水淹的模拟方法,步骤(S2)中所述的建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯,水淹源信息包括水淹源的流量大小及水流方向,采用DualSPHysics软件进行水动力数学模型的建模。
进一步,如上所述的建筑物内部水淹的模拟方法,步骤(S2)中以厂房和相关基座、设备结构作为模型固壁边界,计算中只考虑厂房内流体粒子,粒子运动由连续性方程、动量方程和状态方程控制,在笛卡尔坐标系下,控制方程如下:
其中:t为时间;v为流体速度;ρ为流体密度;ρ0为参考密度;m为流体质量;P为压力;Π为粘性项;W为核函数;g为重力加速度;C0为参考密度下的声速;γ为常系数;下标i表示第i个粒子,下标j表示i粒子作用域范围内的第j个粒子;下标ij表示第i个粒子与其作用范围内第j个粒子之间的相互作用关系;下标0为常数角标。
进一步,如上所述的建筑物内部水淹的模拟方法,步骤(S3)中对物理模型试验结果和水动力数学模型计算结果进行整理,对比两者的水体漫延状态以及水位高度数据,并进行数学模型参数的调整,最终达到两者淹水效果基本一致的情况。
进一步,如上所述的建筑物内部水淹的模拟方法,步骤(S5)中通过提取DualSPHysics软件的计算结果文件,量化水淹过程中水位随时间的变化过程,得到房间水淹高度、水淹高度的时程曲线以及水体漫延路径。
本发明的有益效果如下:本发明通过物理模型试验对数学模型进行率定和验证,从而能够获得可靠的数学模型。同时,采用无网格的光滑粒子动力学方法模拟水淹三维动态过程,该方法适于计算液面大变形问题,避免了网格方法在水体翻转、破碎等复杂情况下追踪自由液面时易于产生失真的不足。由于核岛厂房建筑物结构复杂,计算量大,本发明通过并行计算技术解决了计算效率问题,实现了核岛厂房建筑物内部水淹的精确计算,为核岛内部水淹的分析研究提供数据支持和保障。
附图说明
图1为本发明建筑物内部水淹的模拟系统结构示意图;
图2为本发明建筑物内部水淹的模拟方法流程图;
图3为本发明水流量控制设备的基本组成示意图;
图4为本发明内部水淹模拟系统的物理模型和数学模型试验数据对比图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提供的是一整套核岛厂房内部水淹情况量化分析的系统及方法。该核岛厂房内部水淹的模拟方法是采用数值模型进行全三维动态模拟的方式,利用SPH方法的优势,从而可以精确得出内部水淹分析中的关键数据——水淹高度和时程曲线。
如图1所示,基于SPH算法的核岛厂房内部水淹模拟系统,包括物理试验模型、水动力建模模块、验证模块、计算模块、后处理模块,其中,
物理试验模型,用于根据模拟条件开展物理试验,作为数学模拟结果的验证依据;物理试验模型是依据建筑物结构特点和内部水淹源信息,按重力相似准则进行设计,结合试验结构尺度及设备性能确定模型长度比尺,设置满足试验条件的流量、粗糙度进行淹水试验;可选用有机玻璃作为模型材料,物理试验设备包括水流量控制设备和超声波水位测量系统,所述水流量控制设备包括蓄水池、水泵、电磁流量计、压力计、调剂管路;
水动力建模模块,用于根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;水动力建模模块在建模中依据的建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯等,水淹源信息包括水淹源的流量大小及水流方向;
验证模块,用于根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
计算模块,用于采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算,为提高计算效率,在配备专业GPU科学运算加速设备的计算机上采用并行技术完成水淹事件模拟计算;
后处理模块,用于根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,并制作三维动态演示,从而可以直观监测水淹发生过程中,建筑物内构筑物的淹没情况,对核岛内部水淹的分析研究提供了重要的依据。
对应上述系统的核岛厂房内部水淹模拟方法的流程如图2所示,包括如下步骤:
(S1)根据模拟条件开展物理试验,作为数学模拟结果的验证依据,用于验证和改进数学模型对于各种流量水淹工况的模拟精度和效果;
(S2)根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯等,水淹源信息包括水淹源的流量大小及水流方向;本实施例可采用基于SPH方法的DualSPHysics软件进行水动力数学模型的建模;
(S3)根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
(S4)采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算;
(S5)根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,得到房间水淹高度、水淹高度的时程曲线以及水体蔓延路径,并制作三维动态演示。
下面对各方法步骤进行具体说明。
1.物理模型试验:
为确定数学模型计算结果的可靠性,开展物理模型试验,帮助调整模型参数,包括粒子的尺寸和时间步长等数据。通过分析厂房结构特点(包括建筑物结构、门类型、门槛、孔洞、楼梯等)和内部水淹源流量和位置信息(包括流量大小、水流方向),按重力相似准则设计物理模型试验,等比例建立核岛厂房的模型,结合试验结构尺度及设备性能等因素确定模型长度比尺,制作模型和开展试验,采用高精度动态水位监测设备对水淹时间中不同位置处水位变化进行测量。
建立物理模型试验主要包含:厂房模型制作、水流量控制设备安装和调试、水位测量仪和数据采集软件安装、模拟不同工况下厂房内部水淹模拟试验、试验数据整理分析。
可以选用高透明、高硬度的有机玻璃作为模型制作材料,保证模型尺度与原型几何相似。开展试验需要用到的仪器和设备包括水流量控制设备和超声波水位测量系统,可以动态定点测量水淹过程中水位时程变化曲线。通过确定试验流量,根据流量要求,定制高精度恒定水流量控制设备,水流量控制设备的基本组成如图3所示,包括蓄水池1、水泵2、电磁流量计5、压力计3、调节阀4、橡胶软管6、扩径/缩径短管7,通过超声波水位探头采集的数据及摄像机记录的影像进行提取,分析试验数据的合理性。
针对每组试验工况重复三次,具体试验步骤如下:
①调试供水系统,达到泄漏点出水要求;
②根据泄漏点位置,架设管路;
③在水淹主要路径的房间内调整水位测量探头;
④水位测量系统零点标定;
⑤架设摄像装置,确定拍摄内容;
⑥开始试验,记录不同房间内水位变化、最大水淹高度、水淹路径;
⑦排空模型内水体,清洁、烘干模型,重复试验。
试验结束,对超声波水位探头采集的数据及摄像机记录的影像进行提取并分析试验数据的合理性,将水位随时间的变化历程绘制成曲线,作为数学模型验证的依据。
2.水动力数学模型的建模:
(1)数值模型建立:
根据核岛厂房的建筑物结构、门类型、门槛、孔洞、楼梯等相关信息以及水淹源的流量大小及水流方向,采用基于SPH方法的DualSPHysics软件进行水动力建模,DualSPHysics软件为一款基于光滑粒子流体动力学(SPH)模型的开源软件。
以厂房和相关基座、设备结构等作为模型固壁边界,计算中只考虑厂房内流体粒子,粒子运动由连续性方程、动量方程和状态方程控制,在笛卡尔坐标系下,控制方程描述如下:
其中:t为时间;v为流体速度;ρ为流体密度;ρ0为参考密度;m为流体质量;P为压力;Π为粘性项;W为核函数;g为重力加速度;C0为参考密度下的声速;γ为常系数;下标i表示第i个粒子,下标j表示i粒子作用域范围内的第j个粒子;下标ij表示第i个粒子与其作用范围内第j个粒子之间的相互作用关系;下标0为常数角标。
(2)参数设定:
选用Quintic内核形式,相应的光滑尺度选为0.866;选用Symplectic时间积分格式,为了提高计算效率采用变时间步长方法,相应最大克朗数设置为0.3;粘性处理方式采用常用的人工粘性方式,粘性系数设置为0.01;重力加速度设置为9.81m/s2;水体密度设置为1000kg/s3;状态方程中常系数γ设置为7;使用Delta-SPH公式增加扩散项减小密度波动,Delta-SPH系数设置为0.1;固壁边界采用交错方式布置粒子。
3.对水动力数学模型进行率定和验证
数学模型计算中粒子尺寸、时间步长及模型参数的选取直接关系到计算精度、效率和流态的真实性。时间步长和粒子尺寸设置过大,粒子数量较少,但是计算结果精度低,直接导致水位高度结果和实际情况不符;时间步长和粒子尺寸设置过小,使计算所需粒子数量异常庞大,计算负荷远超现有计算机的能力。按照物理模型试验结果对上述数学模型的参数进行调整,可以确保后期进行厂房水淹时间的系列数学模拟时,能够在合理的时间范围内和比较经济的计算耗费下得到满足精度要求的结果。对物理模型试验结果和数学模型计算结果进行整理,对比两者的水体漫延状态以及水位高度数据,并进行数学模型参数的调整,最终达到两者淹水效果基本一致的情况。
图4为物理模型试验与水动力数学模型的试验数据对比图。
4.根据计算要求开展水淹事件模拟计算:
模型建好后,在配备专业GPU科学运算加速设备(2880个核心)的计算机上完成。DualSPHysics软件求解器采用先进、成熟的SPH算法,支持中央处理器(CPU)或图形处理器(GPU)并行计算,特别采用基于GPU的并行技术,实现了上千万规模粒子的高效流场模拟,满足实际工程中对计算精度和效率的要求。
5.计算结果的后处理:
(1)数据处理
通过已有的技术手段二次开发数据提取程序,这对于本领域的技术人员来说是可以自行实现的,以所开发的数据提取程序读取DualSPHysics的计算结果文件,量化水淹过程的水位随时间的变化过程,得到房间水淹高度、水淹高度的时程曲线以及水体蔓延路径。
(2)三维动态演示
DualSPHysics的计算结果可以直接用ParaView进行读取并制作成三维演示动画。通过制作成的三维演示动画,可以直观监测水淹发生过程中,建筑物内构筑物的淹没情况,对核岛内部水淹的分析研究提供了重要的依据。
综上所述,本发明描述了一整套核岛厂房内部水淹情况量化分析的方法。开展物理模型试验,用于数值模型进行率定和验证,获得可靠的数学模型。采用无网格粒子法来进行带自由面水体流动过程的模拟,粒子代表介质系统,适合模拟自由面大变形问题,不需要复杂的算法追踪自由表面、移动边界等运动特征,可以有效避免欧拉方法中模拟翻转、破碎等复杂自由表面变化时重构带来的失真。与传统网格方法相比计算精度较好,对于大量粒子的高负荷计算可以通过并行计算来提高计算效率。
本发明严格按照厂房实际水淹信息作为输入条件,进行数值模拟,准确反映水淹路径和水位高度,为核岛内部水淹的分析研究提供数据支持及保障,填补了核岛厂房内部水淹精确计算的空白。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其同等技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (12)
1.一种建筑物内部水淹的模拟系统,包括:
物理试验模型,用于根据模拟条件开展物理试验,作为数学模拟结果的验证依据;
水动力建模模块,用于根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;
验证模块,用于根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
计算模块,用于采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算;
后处理模块,用于根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,并制作三维动态演示。
2.如权利要求1所述的建筑物内部水淹的模拟系统,其特征在于所述物理试验模型是依据建筑物结构特点和内部水淹源信息,按重力相似准则进行设计,结合试验结构尺度及设备性能确定模型长度比尺,设置满足试验条件的流量、粗糙度进行淹水试验。
3.如权利要求2所述的建筑物内部水淹的模拟系统,其特征在于:选用有机玻璃作为模型材料,物理试验设备包括水流量控制设备和超声波水位测量系统,所述水流量控制设备包括蓄水池、水泵、电磁流量计、压力计、调剂管路。
4.如权利要求1所述的建筑物内部水淹的模拟系统,其特征在于:所述水动力建模模块在建模中依据的建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯;水淹源信息包括水淹源的流量大小及水流方向。
5.如权利要求1所述的建筑物内部水淹的模拟系统,其特征在于:所述计算模块在配备专业GPU科学运算加速设备的计算机上采用并行技术完成水淹事件模拟计算。
6.如权利要求1所述的建筑物内部水淹的模拟系统,其特征在于:所述后处理模块采用Para View软件直接读取计算模块的计算结果,并制作成三维演示动画。
7.一种建筑物内部水淹的模拟方法,包括如下步骤:
(S1)根据模拟条件开展物理试验,作为数学模拟结果的验证依据;
(S2)根据建筑物相关信息以及水淹源信息,基于SPH方法进行水动力数学模型的建模;
(S3)根据物理试验模型的试验结果,对水动力数学模型进行率定和验证;
(S4)采用经过验证的水动力数学模型根据计算要求开展水淹事件模拟计算;
(S5)根据模拟计算结果量化给出水淹过程中水位随时间的变化过程,并制作三维动态演示。
8.如权利要求7所述的建筑物内部水淹的模拟方法,其特征在于:步骤(S1)中物理试验针对每组试验工况进行多次重复,试验步骤如下:
(S1-1)调试供水系统,达到泄漏点出水要求;
(S1-2)根据泄漏点位置,架设管路;
(S1-3)在水淹主要路径的房间内调整水位测量探头;
(S1-4)对水位测量系统进行零点标定;
(S1-5)架设摄像装置,确定拍摄内容;
(S1-6)开始试验,记录不同房间内水位变化、最大水淹高度、水淹路径;
(S1-7)排空模型内水体,清洁、烘干模型,并返回步骤(S1-1)重复试验。
9.如权利要求7所述的建筑物内部水淹的模拟方法,其特征在于:步骤(S2)中所述的建筑物相关信息包括建筑物结构、门类型、门槛、孔洞、楼梯,水淹源信息包括水淹源的流量大小及水流方向,采用DualSPHysics软件进行水动力数学模型的建模。
10.如权利要求7或9所述的建筑物内部水淹的模拟方法,其特征在于:步骤(S2)中以厂房和相关基座、设备结构作为模型固壁边界,计算中只考虑厂房内流体粒子,粒子运动由连续性方程、动量方程和状态方程控制,在笛卡尔坐标系下,控制方程如下:
其中:t为时间;v为流体速度;ρ为流体密度;ρ0为参考密度;m为流体质量;P为压力;Π为粘性项;W为核函数;g为重力加速度;C0为参考密度下的声速;γ为常系数;下标i表示第i个粒子,下标j表示i粒子作用域范围内的第j个粒子;下标ij表示第i个粒子与其作用范围内第j个粒子之间的相互作用关系;下标0为常数角标。
11.如权利要求7所述的建筑物内部水淹的模拟方法,其特征在于:步骤(S3)中对物理模型试验结果和水动力数学模型计算结果进行整理,对比两者的水体漫延状态以及水位高度数据,并进行数学模型参数的调整,最终达到两者淹水效果基本一致的情况。
12.如权利要求9所述的建筑物内部水淹的模拟方法,其特征在于:步骤(S5)中通过提取DualSPHysics软件的计算结果文件,量化水淹过程中水位随时间的变化过程,得到房间水淹高度、水淹高度的时程曲线以及水体蔓延路径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910232763.1A CN110083866A (zh) | 2019-03-26 | 2019-03-26 | 一种建筑物内部水淹的模拟系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910232763.1A CN110083866A (zh) | 2019-03-26 | 2019-03-26 | 一种建筑物内部水淹的模拟系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110083866A true CN110083866A (zh) | 2019-08-02 |
Family
ID=67413657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910232763.1A Pending CN110083866A (zh) | 2019-03-26 | 2019-03-26 | 一种建筑物内部水淹的模拟系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110083866A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111008418A (zh) * | 2019-11-18 | 2020-04-14 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN111929427A (zh) * | 2020-07-01 | 2020-11-13 | 中国核电工程有限公司 | 橡胶膨胀节性能试验验证方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070043544A1 (en) * | 2004-12-27 | 2007-02-22 | Oh-Young Song | Method for simulating stable but non-dissipative water |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN104991999A (zh) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | 一种基于二维sph的溃坝洪水演进模拟方法 |
CN108334719A (zh) * | 2018-03-26 | 2018-07-27 | 四川理工学院 | 一种基于sph法的土质边坡稳定性及滑坡运动过程分析方法 |
-
2019
- 2019-03-26 CN CN201910232763.1A patent/CN110083866A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070043544A1 (en) * | 2004-12-27 | 2007-02-22 | Oh-Young Song | Method for simulating stable but non-dissipative water |
CN104143027A (zh) * | 2014-08-01 | 2014-11-12 | 北京理工大学 | 一种基于sph算法的流体热运动仿真系统 |
CN104991999A (zh) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | 一种基于二维sph的溃坝洪水演进模拟方法 |
CN108334719A (zh) * | 2018-03-26 | 2018-07-27 | 四川理工学院 | 一种基于sph法的土质边坡稳定性及滑坡运动过程分析方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111008418A (zh) * | 2019-11-18 | 2020-04-14 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN111008418B (zh) * | 2019-11-18 | 2023-06-06 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN111929427A (zh) * | 2020-07-01 | 2020-11-13 | 中国核电工程有限公司 | 橡胶膨胀节性能试验验证方法 |
CN111929427B (zh) * | 2020-07-01 | 2023-04-07 | 中国核电工程有限公司 | 橡胶膨胀节性能试验验证方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Savage et al. | Physical and numerical modeling of large headwater ratios for a 15 labyrinth spillway | |
Roland et al. | On the development and verification of a 2-D coupled wave-current model on unstructured meshes | |
Wu et al. | Numerical modeling of dam-break flood through intricate city layouts including underground spaces using GPU-based SPH method | |
CN110083866A (zh) | 一种建筑物内部水淹的模拟系统及方法 | |
Lin et al. | Investigation of the Reynolds number independence of cavity flow in 2D street canyons by wind tunnel experiments and numerical simulations | |
CN112464584A (zh) | 自由表面流的水位和流量推求方法 | |
Zhang et al. | An improved SWE model for simulation of dam-break flows | |
Thordal et al. | Engineering approach for a CFD inflow condition using the precursor database method | |
Echeverribar et al. | Analysis of the performance of a hybrid CPU/GPU 1D2D coupled model for real flood cases | |
Fan et al. | Development of a benchmark for drag correlations of nonspherical particles based on settling experiments of super-ellipsoidal particles | |
CN112485106B (zh) | 一种控制土体状态参数的物理模型分层制备与试验方法 | |
Ikinciogullari et al. | Comparison of scour properties of classical and Trapezoidal Labyrinth Weirs | |
Shi et al. | 3-D numerical simulation of curved open channel confluence flow with partially non-submerged rigid vegetation | |
CN108256266B (zh) | 一种一维水动力模型和二维水动力模型耦合方法及系统 | |
Fadaei-Kermani et al. | Numerical detection of cavitation damage on dam spillway | |
Majdalani et al. | Experimental investigation of loop rating curve on a small 3D printed laboratory channel | |
Kuiry et al. | A high-resolution shallow water model using unstructured quadrilateral grids | |
Liu et al. | A second-order treatment to the wet–dry interface of shallow water | |
Li et al. | Numerical simulation of random wave overtopping of rubble mound breakwater with armor units | |
Bolgov et al. | Computer hydrodynamic model of the Lower Volga | |
Saffar et al. | Flow-3D numerical modeling of converged side Weir | |
Berger Jr | Free surface flow over curved surfaces | |
Moussa et al. | MPI-Based Simulation of the Shallow Water Model using the Finite Volume Characteristics Scheme | |
Elkholy et al. | Investigations on the trajectory of large sandbags in open channel flow | |
Wilsnack | An Application of Computational Fluid Dynamics to the Hydraulic Analysis of a Water Intake Tower. |
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 |