[go: up one dir, main page]

CN106290065B - A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor - Google Patents

A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor Download PDF

Info

Publication number
CN106290065B
CN106290065B CN201610609356.4A CN201610609356A CN106290065B CN 106290065 B CN106290065 B CN 106290065B CN 201610609356 A CN201610609356 A CN 201610609356A CN 106290065 B CN106290065 B CN 106290065B
Authority
CN
China
Prior art keywords
shear rate
viscosity
reactor
critical point
relationship curve
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
Application number
CN201610609356.4A
Other languages
Chinese (zh)
Other versions
CN106290065A (en
Inventor
曹秀芹
尹伟齐
赵振东
袁海光
丁浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Civil Engineering and Architecture
Original Assignee
Beijing University of Civil Engineering and Architecture
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Civil Engineering and Architecture filed Critical Beijing University of Civil Engineering and Architecture
Priority to CN201610609356.4A priority Critical patent/CN106290065B/en
Publication of CN106290065A publication Critical patent/CN106290065A/en
Application granted granted Critical
Publication of CN106290065B publication Critical patent/CN106290065B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N11/00Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
    • G01N2011/0026Investigating specific flow properties of non-Newtonian fluids

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Mixers Of The Rotary Stirring Type (AREA)

Abstract

The present invention relates to the method and devices that a kind of determining non-newtonian fluid stirs dead zone in reactor.The described method includes: obtaining the relation curve of the shear rate of the viscosity and reactor agitating device of material in reactor;The shear rate effect critical point of the relation curve is obtained according to preset algorithm;The shear rate of reactor agitating device is stirring dead zone when being less than the shear rate effect critical point.Described device is realized based on method as discussed above.The present invention can determine different stirring dead zones according to different non-newtonian fluids, so as to provide technical support for reactor design, parameter optimization and saving energy consumption.

Description

一种确定非牛顿流体在反应器内搅拌死区的方法及装置A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor

技术领域technical field

本发明涉及一种确定非牛顿流体在反应器内搅拌死区的方法及装置。The invention relates to a method and a device for determining a non-Newtonian fluid stirring dead zone in a reactor.

背景技术Background technique

反应器内的物料从流体性质而言分为牛顿流体和非牛顿流体两大类,牛顿流体流体特征比较清晰,但非牛顿流体的流动规律比较复杂,对于非牛顿流体,黏度对其流动行为有较大影响,黏度越大,越易黏着反应器壁上,越不易流动,进而在反应器内形成死区和短流。很难直接精确全面的预测物料在反应器内的流动行为和运动特性,造成难以改善和优化反应器内物料的搅拌混合效果。反应器经常会出现物料混合不均而出现明显分层、上浮和干化结壳的现象。The materials in the reactor are divided into two categories: Newtonian fluid and non-Newtonian fluid in terms of fluid properties. The fluid characteristics of Newtonian fluid are relatively clear, but the flow law of non-Newtonian fluid is more complicated. For non-Newtonian fluid, viscosity has an influence on its flow behavior. The greater the impact, the greater the viscosity, the easier it is to stick to the wall of the reactor, and the less likely it is to flow, thereby forming a dead zone and short flow in the reactor. It is difficult to directly, accurately and comprehensively predict the flow behavior and movement characteristics of materials in the reactor, making it difficult to improve and optimize the mixing effect of materials in the reactor. The reactor often has the phenomenon of uneven material mixing, resulting in obvious stratification, floating and dry crusting.

在搅拌操作时,死区是反映物料是评价反应器搅拌混合效果的重要指标,其能够反映物料在反应器内的混合均匀程度以及物料在反应器内的传质传热效果。国内外不同的研究学者对反应器内死区有不同的认识和判断。在直斜错位桨搅拌槽内流场的探究中,有学者将死区定义为速度小于等于0.02m/s的区域;在厌氧消化反应器的数值模拟中将速度小于0.001m/s的区域定义为死区;在高黏体系中最大叶片式搅拌桨直径的CFD优化中将速度小于0.01v(叶段速率)的区域定义为搅拌死区。可见学者对死区的概念还没有统一的科学定义,仅仅从流体或者搅拌桨叶速度的角度来考虑问题,并没有结合反应器内物料自身流动性质来考虑问题,因此往往具有局限性。During the stirring operation, the dead zone reflects the material and is an important indicator for evaluating the stirring and mixing effect of the reactor. It can reflect the uniformity of mixing of the material in the reactor and the mass and heat transfer effect of the material in the reactor. Different researchers at home and abroad have different understandings and judgments on the dead zone in the reactor. In the exploration of the flow field in the stirred tank with straight inclined misaligned impellers, some scholars define the dead zone as the area where the velocity is less than or equal to 0.02m/s; in the numerical simulation of the anaerobic digestion reactor, the area where the velocity is less than 0.001m/s It is defined as the dead zone; in the CFD optimization of the maximum diameter of the blade-type stirring impeller in the high-viscosity system, the area where the velocity is less than 0.01v (blade velocity) is defined as the stirring dead zone. It can be seen that scholars do not have a unified scientific definition of the concept of dead zone. They only consider the problem from the perspective of fluid or stirring blade speed, and do not consider the problem in conjunction with the flow properties of the material in the reactor itself, so it often has limitations.

发明内容Contents of the invention

针对现有技术中的缺陷,本发明提供一种确定非牛顿流体在反应器内搅拌死区的方法及装置,用于解决现有技术中未考虑非牛顿流体在特殊情况而导致的搅拌不均匀的问题,可以根据不同非牛顿流体确定不同的搅拌死区,从而可以为反应器设计、参数优化以及节省能耗提供技术支撑。Aiming at the defects in the prior art, the present invention provides a method and device for determining the non-Newtonian fluid stirring dead zone in the reactor, which is used to solve the non-uniform stirring caused by non-Newtonian fluid under special circumstances in the prior art According to different non-Newtonian fluids, different stirring dead zones can be determined, which can provide technical support for reactor design, parameter optimization and energy saving.

第一方面,本发明提供了一种确定非牛顿流体在反应器内搅拌死区的方法,所述方法包括:In a first aspect, the present invention provides a method for determining a non-Newtonian fluid stirring dead zone in a reactor, the method comprising:

获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;Obtain the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device;

根据预设算法获取所述关系曲线的剪切速率效应临界点;Obtaining the shear rate effect critical point of the relationship curve according to a preset algorithm;

反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区。When the shear rate of the stirring device of the reactor is lower than the critical point of the shear rate effect, it is a stirring dead zone.

可选地,所述预设算法为分割法或者黏度倍数法。Optionally, the preset algorithm is a division method or a viscosity multiple method.

可选地,所述分割法包括:Optionally, the segmentation method includes:

获取物料黏度μ和剪切速率γ的取值范围;Obtain the value range of material viscosity μ and shear rate γ;

将所述取值范围分割成n等份,n=1,2,……;Divide the value range into n equal parts, n=1, 2, ...;

当剪切速率γ满足:时,剪切速率γi为剪切速率效应临界点。When the shear rate γ satisfies: When , the shear rate γ i is the critical point of shear rate effect.

可选地,所述黏度倍数法包括:Optionally, the viscosity multiple method includes:

获取物料黏度的最大黏度值μmax、极限黏度值μmin和参考黏度μ1;所述参考黏度μ1=m×μmin;m为非零整数;Obtain the maximum viscosity value μ max , limit viscosity value μ min and reference viscosity μ 1 of the material viscosity; the reference viscosity μ 1 = m×μ min ; m is a non-zero integer;

当关系曲线首次落入直线y=μ1和直线y=μmin形成的范围内时,对应的剪切速率为剪切速率效应临界点。When the relationship curve falls into the range formed by the straight line y= μ1 and the straight line y=μmin for the first time, the corresponding shear rate is the critical point of the shear rate effect.

可选地,获取所述关系曲线的步骤包括:Optionally, the step of obtaining the relationship curve includes:

采用HAAKE Viscotester-550型旋转黏度计对预设含固率的物料进行流变测试得到物料黏度与剪切速率的实验曲线;The HAAKE Viscotester-550 rotary viscometer is used to conduct rheological tests on materials with preset solid content to obtain the experimental curves of material viscosity and shear rate;

采用Sisko模型方程对所述实验曲线进行拟合,得出物料-剪切应力方程即关系曲线。Sisko model equations were used to fit the experimental curves to obtain the material-shear stress equation, namely the relationship curve.

第二方面,本发明实施例还提供了一种确定非牛顿流体在反应器内搅拌死区的装置,所述装置包括:In the second aspect, the embodiment of the present invention also provides a device for determining the non-Newtonian fluid stirring dead zone in the reactor, the device comprising:

关系曲线获取模块,用于获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;The relationship curve acquisition module is used to acquire the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device;

临界点获取模块,用于根据预设算法获取所述关系曲线的剪切速率效应临界点;A critical point acquisition module, configured to obtain the shear rate effect critical point of the relationship curve according to a preset algorithm;

搅拌死区获取模块,用于反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区。The stirring dead zone acquisition module is used for the stirring dead zone when the shear rate of the stirring device of the reactor is lower than the critical point of the shear rate effect.

可选地,所述临界点获取模块采用分割法或者黏度倍数法获取剪切速率效应临界点。Optionally, the critical point acquisition module acquires the shear rate effect critical point by a segmentation method or a viscosity multiple method.

可选地,所述临界点获取模块用于执行以下步骤:Optionally, the critical point acquisition module is configured to perform the following steps:

获取物料黏度μ和剪切速率γ的取值范围;Obtain the value range of material viscosity μ and shear rate γ;

将所述取值范围分割成n等份,n=1,2,……;Divide the value range into n equal parts, n=1, 2, ...;

当剪切速率γ满足:时,剪切速率γi为剪切速率效应临界点。When the shear rate γ satisfies: When , the shear rate γ i is the critical point of shear rate effect.

可选地,所述临界点获取模块用于执行以下步骤:Optionally, the critical point acquisition module is configured to perform the following steps:

获取物料黏度的最大黏度值μmax、极限黏度值μmin和参考黏度μ1;所述参考黏度μ1=m×μmin;m为非零整数;Obtain the maximum viscosity value μ max , limit viscosity value μ min and reference viscosity μ 1 of the material viscosity; the reference viscosity μ 1 = m×μ min ; m is a non-zero integer;

当关系曲线首次落入直线y=μ1和直线y=μmin形成的范围内时,对应的剪切速率为剪切速率效应临界点。When the relationship curve falls into the range formed by the straight line y= μ1 and the straight line y=μmin for the first time, the corresponding shear rate is the critical point of the shear rate effect.

可选地,所述关系曲线获取模块通过以下步骤获取所述关系曲线:Optionally, the relationship curve acquisition module acquires the relationship curve through the following steps:

采用HAAKE Viscotester-550型旋转黏度计对预设含固率的物料进行流变测试得到物料黏度与剪切速率的实验曲线;The HAAKE Viscotester-550 rotary viscometer is used to conduct rheological tests on materials with preset solid content to obtain the experimental curves of material viscosity and shear rate;

采用Sisko模型方程对所述实验曲线进行拟合,得出物料-剪切应力方程即关系曲线。Sisko model equations were used to fit the experimental curves to obtain the material-shear stress equation, namely the relationship curve.

由上述技术方案可知,本发明通过获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线,然后采用预设算法可到该关系曲线的剪切速率效应临界点,将剪切速度小于上述剪切速率效应临界点的区别作为搅拌死区。可见,本发明可以根据不同的物料确定其不同死区,从而可以为反应器设计时提供支撑,使反应器工作在效率最佳的状态,节省能耗。It can be seen from the above technical scheme that the present invention obtains the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device, and then adopts a preset algorithm to reach the critical point of the shear rate effect of the relationship curve, and the shear rate The difference between the velocity less than the critical point of the above-mentioned shear rate effect is regarded as the dead zone of stirring. It can be seen that the present invention can determine different dead zones according to different materials, so as to provide support for the design of the reactor, make the reactor work in the state of optimal efficiency, and save energy consumption.

附图说明Description of drawings

通过参考附图会更加清楚的理解本发明的特征和优点,附图是示意性的而不应理解为对本发明进行任何限制,在附图中:The features and advantages of the present invention will be more clearly understood by referring to the accompanying drawings, which are schematic and should not be construed as limiting the invention in any way. In the accompanying drawings:

图1是本发明实施例提供的一种确定非牛顿流体在反应器内搅拌死区的方法流程示意图;Fig. 1 is a schematic flow chart of a method for determining a non-Newtonian fluid stirring dead zone in a reactor provided by an embodiment of the present invention;

图2是分割方法原理示意图;Figure 2 is a schematic diagram of the principle of the segmentation method;

图3是极限黏度倍率法原理示意图;Fig. 3 is a schematic diagram of the principle of the limiting viscosity multiplier method;

图4是本发明一实施例中图2与图3所示的死区定义示意图;Fig. 4 is a schematic diagram of the dead zone definition shown in Fig. 2 and Fig. 3 in an embodiment of the present invention;

图5是本发明另一实施例中图2与图3所示的死区定义示意图;Fig. 5 is a schematic diagram of the dead zone definition shown in Fig. 2 and Fig. 3 in another embodiment of the present invention;

图6是本发明实施例提供的一种确定非牛顿流体在反应器内搅拌死区的装置框图。Fig. 6 is a block diagram of a device for determining a non-Newtonian fluid stirring dead zone in a reactor provided by an embodiment of the present invention.

具体实施方式Detailed ways

为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments It is a part of embodiments of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without creative efforts fall within the protection scope of the present invention.

第一方面,本发明实施例提供了一种确定非牛顿流体在反应器内搅拌死区的方法,如图1所示,所述方法包括:In the first aspect, an embodiment of the present invention provides a method for determining a non-Newtonian fluid stirring dead zone in a reactor, as shown in FIG. 1 , the method includes:

S1、获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;S1. Obtain the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device;

S2、根据预设算法获取所述关系曲线的剪切速率效应临界点;S2. Obtain the shear rate effect critical point of the relationship curve according to a preset algorithm;

S3、反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区。S3. When the shear rate of the stirring device of the reactor is lower than the critical point of the shear rate effect, it is a stirring dead zone.

需要说明的是,本发明实施例中反应器内的物料都是具有一定含固体率的非牛顿流体。其中,含固率指某一混合物中固体成分含量的百分比。It should be noted that the materials in the reactor in the embodiment of the present invention are all non-Newtonian fluids with a certain solid content. Wherein, the solid content refers to the percentage of solid content in a certain mixture.

本发明实施例中预设算法为分割法或者黏度倍数法。The preset algorithm in the embodiment of the present invention is the segmentation method or the viscosity multiple method.

分割法:如图2所示,实验测得反应器内物料黏度随剪切速率的数据,观察知其分布是离散的点。由实验数据可知,物料黏度为μ∈[μminmax],剪切速率为γ∈[γminmax]。结合离散数据相关知识,对黏度μ∈[μminmax]进行分割成n等份,即μ12,...μn(n=1,2,……),与之对应的剪切速率为γ12,...γn。当剪切速率γ满足:Segmentation method: As shown in Figure 2, the data of the viscosity of the material in the reactor with the shear rate is measured experimentally, and the distribution is a discrete point. According to the experimental data, the material viscosity is μ∈[μ minmax ], and the shear rate is γ∈[γ minmax ]. Combining the relevant knowledge of discrete data, the viscosity μ∈[μ minmax ] is divided into n equal parts, namely μ 12 ,...μ n (n=1,2,…), corresponding to The shear rate is γ 1 , γ 2 ,...γ n . When the shear rate γ satisfies:

则此时的γ就是剪切速率效应临界点γ0Then γ at this time is the critical point of shear rate effect γ 0 .

黏度总数法:如图3所示,μmin为物料黏度的极限黏度值,μmax为物料黏度的最大黏度值,μ1=m×μmin为参考黏度,m为非零整数。在关系曲线图中形成直线y=μmin和直线y=μ1。当关系曲线首次落入两直线形成的范围内时,此时的剪切速率即是剪切速率效应临界点γ0Viscosity sum method: as shown in Figure 3, μ min is the limit viscosity value of the material viscosity, μ max is the maximum viscosity value of the material viscosity, μ 1 =m×μ min is the reference viscosity, and m is a non-zero integer. A straight line y=μ min and a straight line y=μ 1 are formed in the relation diagram. When the relationship curve falls into the range formed by the two straight lines for the first time, the shear rate at this time is the critical point γ 0 of the shear rate effect.

下面就以上述两种分法对本发明实施例提供的确定非牛顿流体在反应器内搅拌死区的方法进行详细描述。The method for determining the non-Newtonian fluid stirring dead zone in the reactor provided by the embodiment of the present invention will be described in detail below by using the above two methods.

实施例一Embodiment one

本实施例中非牛顿流体为新鲜猪粪,采用HAAKE Viscotester-550型旋转黏度计对含固率为24.92%下的新鲜猪粪进行流变测试,得到黏度与剪切速率的实验曲线。根据上述实验曲线,采用Sisko模型方程对黏度-剪切速率实验曲线进行拟合,得出猪粪黏度-剪切应力方程为:In this example, the non-Newtonian fluid is fresh pig manure. A HAAKE Viscotester-550 rotational viscometer was used to conduct a rheological test on the fresh pig manure with a solid content of 24.92%, to obtain an experimental curve of viscosity and shear rate. According to the above experimental curve, the Sisko model equation was used to fit the viscosity-shear rate experimental curve, and the pig manure viscosity-shear stress equation was obtained as:

μ=0.0046+9.21·γ-0.6949μ=0.0046+9.21·γ −0.6949 .

上述新鲜猪粪的黏度μ∈[0.10,2.70],即猪粪的极限黏度为0.10Pa·s。The viscosity of the above-mentioned fresh pig manure is μ∈[0.10,2.70], that is, the limit viscosity of pig manure is 0.10Pa·s.

如图4中a图所示,本发明实施例中对黏度范围进行等分,n=10,第十个点处剪切速率γ10=160.90s-1,其坐标为(160.90,0.2744),第九个点处剪切速率γ9=59.35s-1,其坐标为(59.35,0.5441),第八个点处剪切速率γ8=33.11s-1,其坐标为(33.11,0.8139)。剪切速率按照(γii-1)/γi≥0.5的法则,将第十、第九和第八点代入计算,经计算验证其满足As shown in Figure a in Figure 4, in the embodiment of the present invention, the viscosity range is equally divided, n=10, the shear rate γ 10 =160.90s -1 at the tenth point, and its coordinates are (160.90, 0.2744), The shear rate at the ninth point is γ 9 =59.35s -1 , and its coordinates are (59.35, 0.5441), and the shear rate at the eighth point is γ 8 =33.11s -1 , and its coordinates are (33.11, 0.8139). According to the law of (γ ii-1 )/γ i ≥ 0.5, the shear rate is calculated by substituting the tenth, ninth and eighth points into the calculation, and it is verified by calculation that it satisfies

从而得出剪切速率效应临界点γ0为160.90s-1Thus, the critical point of shear rate effect γ 0 is 160.90s -1 .

需要说明的是,本发明实施例中上述数据0.50是为了说明实验曲线或者关系曲线中离散的剪切速率与相邻的剪切速率的变化情况而设定的。当然本领域技术人员可以根据具有使用场景进行设置,本发明不作限定。It should be noted that the above data of 0.50 in the embodiment of the present invention is set to illustrate the variation of discrete shear rates and adjacent shear rates in the experimental curve or relationship curve. Of course, those skilled in the art can make settings according to the usage scenario, which is not limited in the present invention.

如图4中b图所示,采用黏度倍率法计算剪切速度效应临界点γ0。在采用同样条件新鲜猪粪的情况下,HAAKE Viscotester-550型旋转黏度计的剪切速率范围0~600s-1。流变实验亦采用含固率为24.92%的新鲜猪粪。m值取为2.5,直线y=2.5×μmin与关系曲线μ=0.0046+9.21·γ-0.6949首次相交于点(165.76,0.27),即判定此交点即为剪切速率效应临界点γ0为165.76s-1As shown in figure b in Figure 4, the critical point γ 0 of the shear rate effect is calculated by the viscosity multiplier method. Under the same condition of fresh pig manure, the shear rate of HAAKE Viscotester-550 rotational viscometer ranges from 0 to 600s -1 . The rheological experiment also used fresh pig manure with a solid content of 24.92%. The value of m is taken as 2.5, the straight line y=2.5×μmin and the relationship curve μ= 0.0046 +9.21·γ -0.6949 intersect for the first time at the point ( 165.76 , 0.27), that is to say, this intersection point is judged to be the critical point of shear rate effect γ0 is 165.76s -1 .

实施例二Embodiment two

本实施例中非牛顿流体为生污泥,采用HAAKE Viscotester-550型旋转黏度计进行含固率为6.01%的生污泥的流变测试。In this example, the non-Newtonian fluid is raw sludge, and a HAAKE Viscotester-550 rotational viscometer is used to conduct a rheological test of the raw sludge with a solid content of 6.01%.

如图5中a图所示,生污泥黏度μ∈[0.10,4.09],即生污泥的极限黏度为0.10Pa·s。对黏度μ∈[0.10,4.09]范围,进行了n=10等分,第十个点处剪切速率γ10=177.34s-1,其坐标为(177.34,0.2778);第九个点处剪切速率γ9=80.38s-1,其坐标为(80.38,0.9030);第八个点处剪切速率γ8=50.60.11s-1,其坐标为(50.60,1.3025)。剪切速率按照(γii-1)/γi≥0.5的法则,将第十、第九和第八点代入计算,经计算验证其满足:As shown in figure a of Figure 5, the raw sludge viscosity μ∈[0.10,4.09], that is, the limiting viscosity of raw sludge is 0.10Pa·s. For the range of viscosity μ∈[0.10,4.09], n=10 equal divisions are carried out, the shear rate at the tenth point is γ 10 =177.34s -1 , and its coordinate is (177.34, 0.2778); Shear rate γ 9 =80.38s -1 , and its coordinates are (80.38, 0.9030); at the eighth point, shear rate γ 8 =50.60.11s -1 , and its coordinates are (50.60, 1.3025). According to the law of (γ ii-1 )/γ i ≥ 0.5, the shear rate is calculated by substituting the tenth, ninth and eighth points into the calculation, and it is verified by calculation that it satisfies:

从而得出剪切速率效应临界点γ0为177.34s-1Thus, the critical point of shear rate effect γ 0 is 177.34s -1 .

如图5中a图所示,采用极限黏度m倍率法则来判定剪切速率效应临界点γ0。HAAKEViscotester-550型旋转黏度计的剪切速率范围0~400s-1。流变实验亦采用含固率为6.01%的生污泥新鲜猪粪。m值取为5,直线y=5×μmin与猪粪黏度-剪切速率曲线μ=0.10+37.28·γ-0.8760首次相交于点(169.31,0.50),即判定此交点即为剪切速率效应临界点γ0为169.31s-1As shown in figure a of Fig. 5, the critical point γ 0 of the shear rate effect is determined by the limiting viscosity m multiplying rule. The shear rate range of HAAKEViscotester-550 rotational viscometer is 0~400s -1 . The rheological experiment also used raw sludge and fresh pig manure with a solid content of 6.01%. The value of m is taken as 5, and the straight line y=5× μmin and the pig manure viscosity-shear rate curve μ=0.10+37.28·γ-0.8760 intersect for the first time at the point (169.31, 0.50), that is, it is determined that the intersection point is the shear rate The effect critical point γ 0 is 169.31s -1 .

第二方面,本发明实施例还提供了一种确定非牛顿流体在反应器内搅拌死区的装置,如图6所示,所述装置包括:In the second aspect, the embodiment of the present invention also provides a device for determining the non-Newtonian fluid stirring dead zone in the reactor, as shown in Figure 6, the device includes:

关系曲线获取模块M1,用于获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;The relationship curve acquisition module M1 is used to acquire the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device;

临界点获取模块M2,用于根据预设算法获取所述关系曲线的剪切速率效应临界点;The critical point obtaining module M2 is used to obtain the shear rate effect critical point of the relationship curve according to a preset algorithm;

搅拌死区获取模块M3,用于反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区。Stirring dead zone acquisition module M3 is used for stirring dead zone when the shear rate of the reactor stirring device is lower than the critical point of the shear rate effect.

可选地,所述临界点获取模块M2采用分割法或者黏度倍数法获取剪切速率效应临界点。Optionally, the critical point acquisition module M2 acquires the shear rate effect critical point by using a division method or a viscosity multiple method.

可选地,所述临界点获取模块M2用于执行以下步骤:Optionally, the critical point acquisition module M2 is configured to perform the following steps:

获取物料黏度μ和剪切速率γ的取值范围;Obtain the value range of material viscosity μ and shear rate γ;

将所述取值范围分割成n等份,n为自然数;Divide the range of values into n equal parts, where n is a natural number;

当剪切速率γ满足:时,剪切速率γi为剪切速率效应临界点。When the shear rate γ satisfies: When , the shear rate γ i is the critical point of shear rate effect.

可选地,所述临界点获取模块M2用于执行以下步骤:Optionally, the critical point acquisition module M2 is configured to perform the following steps:

获取物料黏度的最大黏度值μmax、极限黏度值μmin和参考黏度μ1;所述参考黏度μ1=m×μmin;m为非零整数;Obtain the maximum viscosity value μ max , limit viscosity value μ min and reference viscosity μ 1 of the material viscosity; the reference viscosity μ 1 = m×μ min ; m is a non-zero integer;

当关系曲线首次落入直线y=μ1和直线y=μmin形成的范围内时,对应的剪切速率为剪切速率效应临界点。When the relationship curve falls into the range formed by the straight line y= μ1 and the straight line y=μmin for the first time, the corresponding shear rate is the critical point of the shear rate effect.

可选地,所述关系曲线获取模块M1通过以下步骤获取所述关系曲线:Optionally, the relationship curve acquisition module M1 acquires the relationship curve through the following steps:

采用HAAKE Viscotester-550型旋转黏度计对预设含固率的物料进行流变测试得到物料黏度与剪切速率的实验曲线;The HAAKE Viscotester-550 rotary viscometer is used to conduct rheological tests on materials with preset solid content to obtain the experimental curves of material viscosity and shear rate;

采用Sisko模型方程对所述实验曲线进行拟合,得出物料-剪切应力方程即关系曲线。Sisko model equations were used to fit the experimental curves to obtain the material-shear stress equation, namely the relationship curve.

由上可以看出,本发明实施例提供的装置基于上文所述的方法实现,因而可以解决同样的技术问题,并取得相同的技术效果,在此不再一一赘述。It can be seen from the above that the device provided by the embodiment of the present invention is implemented based on the method described above, so it can solve the same technical problem and achieve the same technical effect, which will not be repeated here.

在本发明中,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。术语“多个”指两个或两个以上,除非另有明确的限定。In the present invention, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance. The term "plurality" means two or more, unless otherwise clearly defined.

虽然结合附图描述了本发明的实施方式,但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型,这样的修改和变型均落入由所附权利要求所限定的范围之内。Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present invention. within the bounds of the requirements.

Claims (4)

1.一种确定非牛顿流体在反应器内搅拌死区的方法,其特征在于,所述方法包括:1. a method for determining non-Newtonian fluid stirring dead zone in reactor, it is characterized in that, described method comprises: 获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;Obtain the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device; 根据预设算法获取所述关系曲线的剪切速率效应临界点;Obtaining the shear rate effect critical point of the relationship curve according to a preset algorithm; 反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区;When the shear rate of the reactor stirring device is less than the critical point of the shear rate effect, it is the dead zone of stirring; 所述预设算法为分割法或者黏度倍数法;The preset algorithm is a segmentation method or a viscosity multiple method; 其中,所述分割法包括:Wherein, the segmentation method includes: 获取物料黏度μ和剪切速率γ的取值范围;Obtain the value range of material viscosity μ and shear rate γ; 将所述取值范围分割成n等份,n=1,2,……;Divide the value range into n equal parts, n=1, 2, ...; 当剪切速率γ满足:时,剪切速率γi为剪切速率效应临界点;When the shear rate γ satisfies: , the shear rate γ i is the critical point of the shear rate effect; 所述黏度倍数法包括:Described viscosity multiple method comprises: 获取物料黏度的最大黏度值μmax、极限黏度值μmin和参考黏度μ1;所述参考黏度μ1=m×μmin;m为非零整数;Obtain the maximum viscosity value μ max , limit viscosity value μ min and reference viscosity μ 1 of the material viscosity; the reference viscosity μ 1 = m×μ min ; m is a non-zero integer; 当关系曲线首次落入直线y=μ1和直线y=μmin形成的范围内时,对应的剪切速率为剪切速率效应临界点。When the relationship curve falls into the range formed by the straight line y= μ1 and the straight line y=μmin for the first time, the corresponding shear rate is the critical point of the shear rate effect. 2.根据权利要求1所述的方法,其特征在于,获取所述关系曲线的步骤包括:2. The method according to claim 1, wherein the step of obtaining the relationship curve comprises: 采用HAAKE Viscotester-550型旋转黏度计对预设含固率的物料进行流变测试得到物料黏度与剪切速率的实验曲线;The HAAKE Viscotester-550 rotary viscometer is used to conduct rheological tests on materials with preset solid content to obtain the experimental curves of material viscosity and shear rate; 采用Sisko模型方程对所述实验曲线进行拟合,得出物料-剪切应力方程即关系曲线。Sisko model equations were used to fit the experimental curves to obtain the material-shear stress equation, namely the relationship curve. 3.一种确定非牛顿流体在反应器内搅拌死区的装置,其特征在于,所述装置包括:3. A device for determining non-Newtonian fluid stirring dead zone in reactor, is characterized in that, described device comprises: 关系曲线获取模块,用于获取反应器内物料的黏度与反应器搅拌装置的剪切速率的关系曲线;The relationship curve acquisition module is used to acquire the relationship curve between the viscosity of the material in the reactor and the shear rate of the reactor stirring device; 临界点获取模块,用于根据预设算法获取所述关系曲线的剪切速率效应临界点;A critical point acquisition module, configured to obtain the shear rate effect critical point of the relationship curve according to a preset algorithm; 搅拌死区获取模块,用于反应器搅拌装置的剪切速率小于所述剪切速率效应临界点时即为搅拌死区;The stirring dead zone acquisition module is used for the stirring dead zone when the shear rate of the reactor stirring device is less than the critical point of the shear rate effect; 所述临界点获取模块采用分割法或者黏度倍数法获取剪切速率效应临界点;The critical point acquisition module obtains the shear rate effect critical point by a segmentation method or a viscosity multiple method; 其中,所述临界点获取模块采用分割法时用于执行以下步骤:Wherein, when the critical point acquisition module adopts the segmentation method, it is used to perform the following steps: 获取物料黏度μ和剪切速率γ的取值范围;Obtain the value range of material viscosity μ and shear rate γ; 将所述取值范围分割成n等份,n=1,2,……;Divide the value range into n equal parts, n=1, 2, ...; 当剪切速率γ满足:时,剪切速率γi为剪切速率效应临界点;When the shear rate γ satisfies: , the shear rate γ i is the critical point of the shear rate effect; 或者所述临界点获取模块采用黏度倍数法时用于执行以下步骤:Or when the critical point acquisition module adopts the viscosity multiple method, it is used to perform the following steps: 获取物料黏度的最大黏度值μmax、极限黏度值μmin和参考黏度μ1;所述参考黏度μ1=m×μmin;m为非零整数;Obtain the maximum viscosity value μ max , limit viscosity value μ min and reference viscosity μ 1 of the material viscosity; the reference viscosity μ 1 = m×μ min ; m is a non-zero integer; 当关系曲线首次落入直线y=μ1和直线y=μmin形成的范围内时,对应的剪切速率为剪切速率效应临界点。When the relationship curve falls into the range formed by the straight line y= μ1 and the straight line y=μmin for the first time, the corresponding shear rate is the critical point of the shear rate effect. 4.根据权利要求3所述的装置,其特征在于,所述关系曲线获取模块通过以下步骤获取所述关系曲线:4. The device according to claim 3, wherein the relationship curve obtaining module obtains the relationship curve through the following steps: 采用HAAKE Viscotester-550型旋转黏度计对预设含固率的物料进行流变测试得到物料黏度与剪切速率的实验曲线;The HAAKE Viscotester-550 rotary viscometer is used to conduct rheological tests on materials with preset solid content to obtain the experimental curves of material viscosity and shear rate; 采用Sisko模型方程对所述实验曲线进行拟合,得出物料-剪切应力方程即关系曲线。Sisko model equations were used to fit the experimental curves to obtain the material-shear stress equation, namely the relationship curve.
CN201610609356.4A 2016-07-28 2016-07-28 A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor Active CN106290065B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610609356.4A CN106290065B (en) 2016-07-28 2016-07-28 A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610609356.4A CN106290065B (en) 2016-07-28 2016-07-28 A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor

Publications (2)

Publication Number Publication Date
CN106290065A CN106290065A (en) 2017-01-04
CN106290065B true CN106290065B (en) 2019-08-20

Family

ID=57662940

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610609356.4A Active CN106290065B (en) 2016-07-28 2016-07-28 A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor

Country Status (1)

Country Link
CN (1) CN106290065B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115479867B (en) * 2022-08-02 2024-05-28 广东工业大学 Method for rapidly distinguishing non-Newtonian fluid

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101072991A (en) * 2004-12-10 2007-11-14 罗狄亚化学公司 Method and device for determining rheological properties of fluids and corresponding identification method
CN101782494A (en) * 2010-03-16 2010-07-21 同济大学 Loxiong viscosity behavior measuring method of high-viscosity modified asphalt and application thereof

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070056358A1 (en) * 2005-09-12 2007-03-15 Liu James Z Micro-rheometer

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101072991A (en) * 2004-12-10 2007-11-14 罗狄亚化学公司 Method and device for determining rheological properties of fluids and corresponding identification method
CN101782494A (en) * 2010-03-16 2010-07-21 同济大学 Loxiong viscosity behavior measuring method of high-viscosity modified asphalt and application thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
猪粪流变特性与表观粘度模型研究;石惠娴等;《农业机械学报》;20140228;第45卷(第2期);188-193

Also Published As

Publication number Publication date
CN106290065A (en) 2017-01-04

Similar Documents

Publication Publication Date Title
Kazemzadeh et al. Intensification of mixing of shear-thinning fluids possessing yield stress with the coaxial mixers composed of two different central impellers and an anchor
Bakker et al. Numerical modelling of non-Newtonian slurry in a mechanical flotation cell
Jegatheeswaran et al. Process intensification in a chaotic SMX static mixer to achieve an energy-efficient mixing operation of non-newtonian fluids
Ameur Mixing of complex fluids with flat and pitched bladed impellers: Effect of blade attack angle and shear-thinning behaviour
CN204177704U (en) A kind of High Temperature High Pressure drag reducer evaluating device
Liu et al. Adjusting the rheological properties of corn-straw slurry to reduce the agitation power consumption in anaerobic digestion
Fernandes del Pozo et al. Hydrodynamic analysis of an axial impeller in a non‐Newtonian fluid through particle image velocimetry
Ameur 3D hydrodynamics involving multiple eccentric impellers in unbaffled cylindrical tank
CN112100944B (en) Visualization method and application of anaerobic digestion flow field under multi-scale conditions based on CFD simulation and PIV measurement
CN101650292B (en) Method for measuring mixing effect of fluid
CN105043905A (en) Method for determining long-term strength parameter of rock based on steady-flow variable-rate tangent line
Li et al. Rheological analysis of Newtonian and non‐Newtonian fluids using Marsh funnel: Experimental study and computational fluid dynamics modeling
CN106290065B (en) A method and device for determining a non-Newtonian fluid stirring dead zone in a reactor
Luan et al. Effect of the 6PBT stirrer eccentricity and off-bottom clearance on mixing of pseudoplastic fluid in a stirred tank
Annas et al. A Scale‐up Strategy for the Fluid Flow in Biogas Plants
CN111524035A (en) Hydrodynamic regulation and control threshold determination method for improvement of river network water environment in plain city
Liu et al. Effect of bubble morphology and behavior on power consumption in non-Newtonian fluids’ aeration process
Chen et al. CFD simulation of two-phase flows in anaerobic digester
Liang et al. Numerical analysis in the effects of blade’s arrangement on the torque load characteristics of the three-blade planetary mixer
Núñez‐Ramírez et al. Mixing analysis for a fermentation broth of the fungus Beauveria bassiana under different hydrodynamic conditions in a bioreactor
Conti et al. Numeric Simulation‐Based Analysis of the Mixing Process in Anaerobic Digesters of Biogas Plants
Wang et al. Power consumption of inner–outer helical ribbon impellers in viscous Newtonian and non-Newtonian fluids
CN104462816A (en) Method for determining separation size of hydrocyclone
Ghorbanian et al. Mixing analysis of PCS slurries in a horizontal scraped surface bioreactor
Youcefi et al. On the influence of liquid elasticity on mixing in a vessel agitated by a two-bladed impeller

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant