CN113033117B - 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 - Google Patents
一种运动荷电液滴诱导电场强度及电场力计算方法和系统 Download PDFInfo
- Publication number
- CN113033117B CN113033117B CN202110255713.2A CN202110255713A CN113033117B CN 113033117 B CN113033117 B CN 113033117B CN 202110255713 A CN202110255713 A CN 202110255713A CN 113033117 B CN113033117 B CN 113033117B
- Authority
- CN
- China
- Prior art keywords
- electric field
- discrete phase
- discrete
- charged liquid
- field strength
- 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.)
- Active
Links
- 230000005684 electric field Effects 0.000 title claims abstract description 135
- 239000007788 liquid Substances 0.000 title claims abstract description 84
- 230000006698 induction Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 50
- 239000012530 fluid Substances 0.000 claims abstract description 39
- 239000002245 particle Substances 0.000 claims description 31
- 230000003993 interaction Effects 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 4
- 238000005096 rolling process Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 14
- 230000008569 process Effects 0.000 description 8
- BWKQDFQFEBGOQJ-NKBQLWGUSA-N [(2E,6E)-9-[(2R)-3,3-dimethyloxiran-2-yl]-3,7-dimethylnona-2,6-dienyl] 2-diazoacetate Chemical compound C\C(CC\C=C(/C)CC[C@H]1OC1(C)C)=C/COC(=O)C=[N+]=[N-] BWKQDFQFEBGOQJ-NKBQLWGUSA-N 0.000 description 6
- 238000000889 atomisation Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005507 spraying Methods 0.000 description 2
- 241000135164 Timea Species 0.000 description 1
- 241000607479 Yersinia pestis Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002105 nanoparticle Substances 0.000 description 1
- 239000011236 particulate material Substances 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000007921 spray Substances 0.000 description 1
- 239000002341 toxic gas Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种运动荷电液滴诱导电场强度及电场力计算方法和系统,包括以下步骤:建立运动荷电液滴诱导电场强度及电场力的几何模型,几何模型内部为计算所需计算域;在离散单元法软件中,按需在计算域内生成相应离散相;编写荷电液滴诱导电场强度及相应电场力计算代码,嵌入计算流体力学软件中;对几何模型进行网格划分,得到网格文件;在计算流体力学软件中,进行边界条件设置,加载计算代码,选择求解器,计算运动荷电液滴诱导电场强度及电场力。本发明在计算流体力学软件中嵌入电场强度及相应电场力计算代码,可以在求解流体相的同时求解离散相荷电液滴的诱导电场强度及相应电场力,极大拓展了计算流体力学软件和离散单元法软件的应用范围。
Description
技术领域
本发明属于EHDA(Electrohydrodynamics atomization)数值模拟计算技术领域,具体涉及基于计算流体力学软件和离散单元法软件的运动荷电液滴诱导电场强度及相应电场力计算方法和系统。
背景技术
液体转变为液滴这一过程,即喷雾或者雾化,已经被人类广泛利用。雾化的方法有多种,其中基于电场方法,即荷电喷雾,因其成本低、污染小、生物相容性高已经被广泛应用于各个领域。例如,病虫防治领域水基农药的雾化;污染治理领域可溶性有毒气体和颗粒物的处理;生产制造领域单分散纳米颗粒材料的制备等。
荷电喷雾的基本单元是流体中运动的荷电液滴。荷电液滴在运动过程中,关键点是运动荷电液滴所诱导的电场强度以及在该电场强度中其它荷电离散相,例如荷电颗粒物,所受的电场力。
随着计算机技术的发展,数值模拟计算已经被广泛用于研究各种物理过程。然而EHDA(Electrohydrodynamics atomization)涉及多相流问题,相间作用及多场耦合机制复杂,目前关于EHDA(Electrohydrodynamics atomization)没有成熟商用代码。部分特定代码可以计算电场强度,但缺乏流体相计算能力。计算流体力学软件可以求解流体相,但没有电势、电场强度等计算代码。有学者基于计算流体力学软件中的热求解器近似计算电势、电场强度等,但涉及诸多假设:流体必须是电介质,流体必须各向同性且不可压缩,不能存在电荷源等。为了同时获得荷电液滴运动过程中的各种流体信息以及电场信息,例如流场、速度场、电势分布、电场强度分布等,必须基于计算流体力学软件和离散单元法软件同时耦合相关计算代码进行数值模拟计算。目前没有基于计算流体力学软件和离散单元法软件数值模拟计算运动荷电液滴诱导电场强度及电场力的研究。
发明内容
针对上述技术问题,本发明提供一种运动荷电液滴诱导电场强度及电场力计算方法和系统,求解出运动荷电液滴诱导电场强度及电场力,从而更好地模拟EHDA过程,拓展计算流体力学软件和离散单元法软件的应用范围。
为解决上述问题,本发明采用以下技术手段实现:
一种运动荷电液滴诱导电场强度及电场力数值计算方法,包括以下步骤:
建立几何模型:建立运动荷电液滴诱导电场强度及电场力的几何模型,几何模型内部为计算所需计算域;
生成离散相:在离散单元法软件中,按需在计算域内生成相应离散相,包括离散相荷电液滴和离散相荷电颗粒;
编写运动荷电液滴诱导电场强度及电场力计算代码:所述计算域中,有若干个离散相荷电液滴和离散相荷电颗粒,先计算单独荷电液滴诱导电势,如果存在多个荷电液滴,将每个荷电液滴诱导电势相互叠加得总诱导电势,总诱导电势求导得总电场强度;在总电场强度中,计算离散相荷电液滴受到的电场力,如果存在离散相荷电颗粒,计算离散相荷电颗粒受到的电场力,运动荷电液滴诱导电场强度及电场力的计算方程通过编写相应计算代码,嵌入计算流体力学软件中;
划分网格:对几何模型划分网格,得到网格文件;
计算:在计算流体力学软件中,设置边界条件,加载计算代码,选择求解器,计算运动荷电液滴诱导电场强度及电场力。
上述方案中,在计算流体力学软件中,通过以下基本方程求解所述计算域中不可压缩气相的速度和压力:
其中,ug为气相速度;ρg为气相密度;p为静压;μg为气相动力粘度;g为重力加速度。
进一步的,在离散单元法软件中,通过以下基本方程求解所述计算域中离散相的受力:
其中,mp为离散相质量;up为离散相速度;Fp,n为离散相所受法相接触力;Fp,t为离散相所受切向接触力;Fp,f为气相-离散相相互作用力;Fp,g为离散相所受重力;Fp,e为离散相所受电场力。Ip为离散相转动惯量;ωp为离散相角动量;Mp,t为离散相所受切向力矩;Mp,r为离散相所受滚动摩擦力矩。
进一步的,气相-离散相相互作用力主要包括曳力Fd和浮力Fb;
曳力基本方程如下:
Fd=mpfD(ug-up) 式五
其中,fD为单位质量曳力系数;
其中,dp为离散相直径;ρp为离散相密度;CD为曳力系数;Re为雷诺数
浮力基本方程如下:
上述方案中,所述总诱导电势求导得总电场强度的具体方程为:
其中,φ为单个荷电液滴的诱导电势;KE为库仑常数;Q为离散相荷电液滴荷电量;r为计算域内各处与荷电液滴质点间距;φsup为总诱导电势;E为总电场强度。
进一步的,在总电场强度中,离散相荷电液滴受到相应电场力Fp,e,d,具体方程为:
FP,e,d=EQ 式十四
上述方案中,在总电场强度中,离散相荷电颗粒受到相应电场力Fp,e,p,具体方程为:
Fp,e,p=Eq 式十五
其中,q是离散相荷电颗粒荷电量。
一种实现所述的运动荷电液滴诱导电场强度及电场力数值计算方法的系统,所述系统包括计算流体力学软件和离散单元法软件,所述计算流体力学软件需要额外嵌入荷电液滴诱导电场强度及电场力的计算代码。
与现有技术相比,本发明的有益效果是:本发明基于计算流体力学软件和离散单元法软件,在计算流体力学软件中嵌入电场强度及电场力计算代码,可以在求解流体相的同时求解离散相荷电液滴的诱导电场强度及相应电场力,从而更好地模拟EHDA过程,极大拓展了计算流体力学软件和离散单元法软件的应用范围。相比其它数值模拟计算方法,本发明仅涉及少量假设,模拟结果与实际物理过程更吻合。
附图说明
图1是本发明流程图;
图2是本发明运动荷电液滴诱导电场强度计算流程图;
图3是本发明运动荷电液滴诱导电场力计算流程图;
图4是本发明运动荷电液滴诱导电场强度及电场力计算域示意图;
图5是本发明离散相轨迹图;
图6是本发明计算域内流场速度场图;
图7是本发明离散相荷电液滴诱导电势图;
图8是本发明离散相荷电液滴诱导电场强度图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
实施例1
如图1所示,一种运动荷电液滴诱导电场强度及电场力数值计算方法,包括以下步骤:
建立几何模型:建立运动荷电液滴诱导电场强度及电场力的几何模型,几何模型内部为计算所需计算域;
生成离散相:在离散单元法软件中,按需在计算域内生成相应离散相,包括离散相荷电液滴和离散相荷电颗粒;
如图2和3所示,编写运动荷电液滴诱导电场强度及电场力计算代码:所述计算域中,有若干个离散相荷电液滴和离散相荷电颗粒,先计算单独荷电液滴诱导电势,如果存在多个荷电液滴,将每个荷电液滴诱导电势相互叠加得总诱导电势,总诱导电势求导得总电场强度;在总电场强度中,计算离散相荷电液滴受到的电场力,如果存在离散相荷电颗粒,计算离散相荷电颗粒受到的电场力。运动荷电液滴诱导电场强度及电场力的计算方程通过编写相应计算代码,嵌入计算流体力学软件中;
划分网格:对几何模型划分网格,得到网格文件;
计算:在计算流体力学软件中,设置边界条件,加载计算代码,选择求解器,计算运动荷电液滴诱导电场强度及电场力。
对于所述计算域中不可压缩的粘性气相,计算流体力学软件应用的基本方程如下:
式一和式二分别是气相的连续性方程和动量守恒方程。其中,ug为气相速度;ρg为气相密度;p为静压;μg为气相动力粘度;g为重力加速度。
对所述计算域中的离散相,离散单元法软件应用的基本方程如下:
其中,mp为离散相质量;up为离散相速度;Fp,n为离散相所受法相接触力;Fp,t为离散相所受切向接触力;Fp,f为气相-离散相相互作用力;Fp,g为离散相所受重力;Fp,e为离散相所受电场力。Ip为离散相转动惯量;ωp为离散相角动量;Mp,t为离散相所受切向力矩;Mp,r为离散相所受滚动摩擦力矩。
进一步的,气相-离散相相互作用力主要包括曳力Fd和浮力Fb。曳力基本方程如下:
Fd=mpfD(ug-up) 式五
其中,fD为单位质量曳力系数;
其中,dp为离散相直径;ρp为离散相密度;CD为曳力系数;Re为雷诺数
浮力基本方程如下:
边界条件:
对于非定常运动问题,需要考虑初始条件的设置。初始条件就是t=t0时各变量的分布,其表示如下:
ug=ug(x,y,t0)=ug,0(x,y) 式九
up=up(x,y,t0)=up,0(x,y) 式十
其中,ug为气相速度,初始条件ug是t0时刻的速度,ug,0(x,y)代表初始条件ug的具体数值。up为离散相速度,初始条件up是t0时刻的速度,up,0(x,y)代表初始条件up的具体数值。
在图4所示的计算域中,有若干个离散相荷电液滴和离散相荷电颗粒。关于荷电液滴的诱导电场强度,首先考虑单个荷电液滴的诱导电势,每个荷电液滴的诱导电势相互叠加得总诱导电势,总诱导电势求导得总电场强度,具体方程为:
其中,φ为单个荷电液滴的诱导电势;KE为库仑常数;Q为单个荷电液滴荷电量;r为计算域各处和荷电液滴间距;φsup为总诱导电势;E为总电场强度。
在总电场强度中,离散相荷电液滴受到相应电场力Fp,e,d,具体方程为:
FP,e,d=EQ 式十四
其中,Q是离散相荷电液滴荷电量。
在总电场强度中,离散相荷电颗粒受到相应电场力Fp,e,p,具体方程为:
Fp,e,p=Eq 式十五
其中,q是离散相荷电颗粒荷电量。
荷电液滴诱导电场强度及相应电场力的计算通过C语言编写相应代码,嵌入计算流体力学软件中。
如下为一个具体操作:
1.构建几何模型并对几何模型进行网格划分,得到网格文件。根据本实施例,优选的,所述运动荷电液滴诱导电场强度及电场力的几何模型示意图如图4所示,所述几何模型长度为L0,宽度为L1,高度为H0。几何模型内部为计算所需计算域。三维笛卡尔坐标系xyz固定在几何模型底部中心,x轴和y轴位于几何模型底部平面,相互垂直,z轴垂直几何模型底部平面,竖直向上。
2.离散单元法软件
2.1设定液滴和颗粒物材料属性。具体为液滴密度1000kg/m3,液滴形状优选为球型,液滴直径优选为2×10-3m;颗粒物密度优选为2200kg/m3,颗粒物形状优选为球型,颗粒物直径优选为1×10-5m。
2.2设定液滴和颗粒物生成区域。优选的,液滴生成区域中心坐标x=0,y=0.09m,z=0,液滴生成区域尺寸是1×10-2m×1×10-2m×1×10-2m;优选的,颗粒物生成区域中心坐标x=0, y=0.05m,z=0,颗粒物生成区域尺寸是2×10-2m×2×10-2m×0.1m。
2.3设定液滴和颗粒物生成参数。液滴生成数目优选为10,液滴生成时间优选为1×10-12s,液滴生成位置优选为随机;颗粒物生成数目优选为20000,颗粒物生成时间优选为1×10-12s,颗粒物生成位置优选为随机。
2.5设定重力加速。优选的,为-y方向,大小为9.81m/s2。
2.6设定时间步长参数。优选为1×10-9s。
2.7设定网格尺寸。优选为2×10-5m。
2.8开启和计算流体力学软件耦合接口。
3.计算流体力学软件
3.1.选择网格文件长度尺寸单位。优选为mm。
3.2.选择求解器。优选为瞬态计算。
3.3.设定重力加速度。优选的,为-y方向,大小为9.81m/s2。
3.4.选择湍流模型。优选为k-epsilon(2eqn)湍流模型。
3.5.设定边界条件。优选为固壁边界条件。
3.6.和离散单元法软件耦合。
3.7.求解器初始化。
3.8.加载计算代码文件。
3.9.设定求解时间参数。优选的,时间步长为1×10-4,时间步数为1000。
从图5可以看出,离散相荷电颗粒在离散相荷电液滴群附近有较明显聚集。随着荷电液滴群下落,其扫掠路径中形成了明显的颗粒物群流。原因一是受限于荷电液滴荷电量,荷电液滴群外部部分颗粒仅被荷电液滴群通过库仑力吸引至该区域,无法到达荷电液滴表面。原因二是荷电液滴群下落路径中会引发气体流动,形成相对低压区,位于该区域的颗粒在较短时间内难以扩散。从图6可以看出,仅在荷电液滴群附近及其下落路径中才有气体流动。从图7可以看出,荷电液滴群内部电势值最大,这是因为电势为标量,在荷电液滴群内部,各个荷电液滴诱导电势相互叠加而导致。距离荷电液滴群越远电势越小。从图8可以看出,荷电液滴群内部电场强度并非最大,这是因为电场强度是矢量,在荷电液滴群内部,各个荷电液滴诱导电场强度相互叠加,实际上互相削弱导致。
可以看到,通过本发明可以数值模拟计算多个荷电液滴诱导电场强度及相应电场力,从而更好地模拟EHDA(Electrohydrodynamics atomization)过程,极大拓展了计算流体力学软件和离散单元法软件应用范围。
实施例2
一种实现实施例1所述的运动荷电液滴诱导电场强度及电场力数值计算方法的系统,因此具有实施例1的有益效果,此处不再赘述。所述系统包括计算流体力学软件和离散单元法软件,所述计算流体力学软件需额外嵌入荷电液滴诱导电场强度及电场力的计算代码。
应当理解,虽然本说明书是按照各个实施例描述的,但并非每个实施例仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施例的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技艺精神所作的等效实施例或变更均应包含在本发明的保护范围之内。
Claims (8)
1.一种运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,包括以下步骤:
建立几何模型:建立运动荷电液滴诱导电场强度及电场力的几何模型,几何模型内部为计算所需计算域;
生成离散相:在离散单元法软件中,在计算域内生成相应离散相,包括离散相荷电液滴和离散相荷电颗粒;
编写运动荷电液滴诱导电场强度及电场力计算代码:所述计算域中,有若干个离散相荷电液滴和离散相荷电颗粒,先计算单独荷电液滴诱导电势,如果存在多个荷电液滴,将每个荷电液滴诱导电势相互叠加得总诱导电势,总诱导电势求导得总电场强度;在总电场强度中,计算离散相荷电液滴受到的电场力,如果存在离散相荷电颗粒,计算离散相荷电颗粒受到的电场力,运动荷电液滴诱导电场强度及电场力的计算方程通过编写相应计算代码,嵌入计算流体力学软件中;
划分网格:对几何模型划分网格,得到网格文件;
计算:在计算流体力学软件中,设置边界条件,加载计算代码,选择求解器,计算运动荷电液滴诱导电场强度及电场力。
2.根据权利要求1所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,在计算流体力学软件中,通过以下基本方程求解所述计算域中不可压缩气相的速度和压力:
▽·ug=0式一
其中,ug为气相速度;ρg为气相密度;p为静压;μg为气相动力粘度;g为重力加速度。
3.根据权利要求2所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,在离散单元法软件中,通过以下基本方程求解所述计算域中离散相的受力:
其中,mp为离散相质量;up为离散相速度;Fp,n为离散相所受法相接触力;Fp,t为离散相所受切向接触力;Fp,f为气相-离散相相互作用力;Fp,g为离散相所受重力;Fp,e为离散相所受电场力;Ip为离散相转动惯量;ωp为离散相角动量;Mp,t为离散相所受切向力矩;Mp,r为离散相所受滚动摩擦力矩。
4.根据权利要求3所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,所述气相与离散相相互作用力包括曳力Fd和浮力Fb;
所述曳力的基本方程如下:
Fd=mpfD(ug-up) 式五
其中,fD为单位质量曳力系数;
其中,dp为离散相直径;ρp为离散相密度;CD为曳力系数;Re为雷诺数;
浮力基本方程如下:
5.根据权利要求1所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,所述总诱导电势求导得总电场强度的具体方程为:
E=▽φsup 式十三
其中,φ为单个荷电液滴的诱导电势;KE为库仑常数;Q为离散相荷电液滴荷电量;r为计算域内各处与荷电液滴间距;φsup为总诱导电势;E为总电场强度。
6.根据权利要求5所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,在总电场强度中,离散相荷电液滴受到相应电场力Fp,e,d,具体方程为:
FP,e,d=EQ 式十四。
7.根据权利要求5所述的运动荷电液滴诱导电场强度及电场力数值计算方法,其特征在于,在总电场强度中,离散相荷电颗粒受到相应电场力Fp,e,p,具体方程为:
Fp,e,p=Eq 式十五
其中,q是离散相荷电颗粒荷电量。
8.一种实现权利要求1-7任意一项所述的运动荷电液滴诱导电场强度及电场力数值计算方法的系统,其特征在于,所述系统包括计算流体力学软件和离散单元法软件,所述计算流体力学软件需要额外嵌入荷电液滴诱导电场强度及相应电场力的计算代码。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110255713.2A CN113033117B (zh) | 2021-03-09 | 2021-03-09 | 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110255713.2A CN113033117B (zh) | 2021-03-09 | 2021-03-09 | 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113033117A CN113033117A (zh) | 2021-06-25 |
CN113033117B true CN113033117B (zh) | 2024-03-19 |
Family
ID=76467234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110255713.2A Active CN113033117B (zh) | 2021-03-09 | 2021-03-09 | 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113033117B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113625067B (zh) * | 2021-08-12 | 2022-06-10 | 华北电力大学 | 一种离子流场中悬浮液滴荷电特性测量装置及方法 |
CN113970662B (zh) * | 2021-10-19 | 2023-05-16 | 中山大学 | 一种基于单个囚禁离子的电场力探测系统 |
CN115294839B (zh) * | 2021-12-22 | 2024-06-11 | 浙江安防职业技术学院 | 模拟异极性荷电粒子在磁电凝并器中的运动情况的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005021151A1 (en) * | 2003-08-27 | 2005-03-10 | President And Fellows Of Harvard College | Electronic control of fluidic species |
CN107644132A (zh) * | 2017-09-18 | 2018-01-30 | 大唐环境产业集团股份有限公司 | 一种用于模拟带电除尘器除尘效率的计算方法 |
CN110287600A (zh) * | 2019-06-26 | 2019-09-27 | 中国人民解放军陆军装甲兵学院 | 磁控等离子体在圆筒内的流动及压力分布研究方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8308926B2 (en) * | 2007-08-20 | 2012-11-13 | Purdue Research Foundation | Microfluidic pumping based on dielectrophoresis |
-
2021
- 2021-03-09 CN CN202110255713.2A patent/CN113033117B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005021151A1 (en) * | 2003-08-27 | 2005-03-10 | President And Fellows Of Harvard College | Electronic control of fluidic species |
CN107644132A (zh) * | 2017-09-18 | 2018-01-30 | 大唐环境产业集团股份有限公司 | 一种用于模拟带电除尘器除尘效率的计算方法 |
CN110287600A (zh) * | 2019-06-26 | 2019-09-27 | 中国人民解放军陆军装甲兵学院 | 磁控等离子体在圆筒内的流动及压力分布研究方法 |
Non-Patent Citations (2)
Title |
---|
Characteristics and health impacts of particulate matter pollution in China (2001e2011);Zhen Cheng等;Atmospheric Environment;第65卷;第186-194页 * |
荷电液滴捕集细颗粒物的数值模拟研究;左子文等;工程热物理学报;第36卷(第3期);第559-563页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113033117A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113033117B (zh) | 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 | |
Eaton | Two-way coupled turbulence simulations of gas-particle flows using point-particle tracking | |
Pendar et al. | Numerical modeling of electrostatic spray painting transfer processes in rotary bell cup for automotive painting | |
Mitra et al. | Collision behaviour of a smaller particle into a larger stationary droplet | |
CN114168796B (zh) | 一种建立飞行器高空气动力数据库的方法 | |
CN112613246A (zh) | 一种固体火箭发动机在飞行过载下的两相流仿真方法 | |
Ellwood et al. | A finite-element model for an electrostatic bell sprayer | |
Li et al. | Hydrodynamic analysis of the energy dissipation of droplets on vibrating superhydrophobic surfaces | |
Wassen et al. | Simulation of cluster formation in gas–solid flow induced by particle–particle collisions | |
Lin et al. | Effect of droplet deformation and internal circulation on drag coefficient | |
Kosinski et al. | Modelling of dust lifting using the Lagrangian approach | |
Koolpiruck et al. | Numerical modeling of inhaled charged aerosol deposition in human airways | |
Vanderroost et al. | Modelling coating quality in fluidised bed coating: Spray sub-model | |
Zheng et al. | Analytical study on turbulent agglomeration in turbulence agglomerator | |
Kulkarni et al. | Development of a parallel lagrangian particle tracking code for 3d multi-block curvilinear grids | |
Frank et al. | Numerical Simulation of Gas-Droplet Flow Around a Nozzle in Cylindrical Chamber Using a Lagrangian Model Based on a Multigrid Navier-Stokes Solver | |
Abdulkadhim | Solid-particles deposition through a turbulent impinging jet using lattice Boltzmann method | |
Wang et al. | Critical comparison of electrostatic effects on hydrodynamics and heat transfer in a bubbling fluidized bed with a central jet | |
Hariharan | Large scale simulation of particle coating using coupled CFD-DEM | |
Jahannama et al. | Electrostatic effects on agricultural air-atomized sprays and deposition. Part II: A computational study | |
Pier et al. | CUDA-enabled particle-based 3D fluid haptic simulation | |
Matsuyama et al. | Lattice Boltzmann method study of Rayleigh instability of a charged droplet | |
Han et al. | Drag and lift force on an ellipsoid particle close to a wall | |
Goniva et al. | Modelling three-phase flow in metallurgical processes | |
Zhao et al. | Particle-scale simulation of fluidized bed with immersed tubes |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |