CN1138967C - 控制连续介质的边界或壁层的方法和装置 - Google Patents

控制连续介质的边界或壁层的方法和装置 Download PDF

Info

Publication number
CN1138967C
CN1138967C CNB961969989A CN96196998A CN1138967C CN 1138967 C CN1138967 C CN 1138967C CN B961969989 A CNB961969989 A CN B961969989A CN 96196998 A CN96196998 A CN 96196998A CN 1138967 C CN1138967 C CN 1138967C
Authority
CN
China
Prior art keywords
continuous medium
fluid
wall
active force
effect
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.)
Expired - Fee Related
Application number
CNB961969989A
Other languages
English (en)
Other versions
CN1196790A (zh
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of CN1196790A publication Critical patent/CN1196790A/zh
Application granted granted Critical
Publication of CN1138967C publication Critical patent/CN1138967C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F28—HEAT EXCHANGE IN GENERAL
    • F28F—DETAILS OF HEAT-EXCHANGE AND HEAT-TRANSFER APPARATUS, OF GENERAL APPLICATION
    • F28F13/00—Arrangements for modifying heat-transfer, e.g. increasing, decreasing
    • F28F13/02—Arrangements for modifying heat-transfer, e.g. increasing, decreasing by influencing fluid boundary
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F15—FLUID-PRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
    • F15D—FLUID DYNAMICS, i.e. METHODS OR MEANS FOR INFLUENCING THE FLOW OF GASES OR LIQUIDS
    • F15D1/00—Influencing flow of fluids
    • F15D1/10—Influencing flow of fluids around bodies of solid material
    • F15D1/12—Influencing flow of fluids around bodies of solid material by influencing the boundary layer
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64C—AEROPLANES; HELICOPTERS
    • B64C2230/00—Boundary layer controls
    • B64C2230/12—Boundary layer controls by using electromagnetic tiles, fluid ionizers, static charges or plasma
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05D—INDEXING SCHEME FOR ASPECTS RELATING TO NON-POSITIVE-DISPLACEMENT MACHINES OR ENGINES, GAS-TURBINES OR JET-PROPULSION PLANTS
    • F05D2270/00—Control
    • F05D2270/01—Purpose of the control system
    • F05D2270/17—Purpose of the control system to control boundary layer
    • F05D2270/172—Purpose of the control system to control boundary layer by a plasma generator, e.g. control of ignition
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T50/00—Aeronautics or air transport
    • Y02T50/10—Drag reduction

Abstract

本发明涉及空气流体动力学以及热和质量的传送，特别是涉及一种控制在表面(1)附近连续介质边界层或壁层的方法和装置，连续介质由气体、液体和/或其混合物组成，通过对流体施加作用改变连续介质粒子的速度，以改变流体结构，湍流程度，冲量的传送，热和/或混合物的传送。

Description

控制连续介质的边界或壁层的方法和装置
发明领域
本发明涉及空气流体动力学，以及热和质量的传送，特别是涉及控制连续介质气体和/或液体的边界或壁层的问题。连续介质随受热的或等温的自由的流线型表面运动，还随压力表面和/或能量传送通道的表面而运动。所说的表面是由该连续介质加热的。本发明还特别涉及按权利要求1开头给出的方法，以及按权利要求8的开头给出的装置。
众所周知，当流体通过一个表面时，边界或壁层是空气流体动力学阻力的一个主要来源，参见Schlichting G.的“边界层原理”(“Theory of BoundaryLayer”，莫斯科NAUKA出版社(1974年)，第712页)，同样它也是在从加热表面到流体核心或从受热流体到流体流动的表面的热传递中确定热阻的因素。
1904年，L.Prandtl曾提出一种控制边界层的方法，即通过吸走移动的使流动的壁层减速的介质粒子(参见Prandtl L.，Titiens O.“流体和空气动力学”“Hydro-and Aerodynamic”IIV.1，2.M.，ONTI，1935年)。这种方法的本质是众所周知的，用在层流和湍流的情况有很实用的效果。
通过所谓的二维人工表面粗糙方法强化热和物质传送而得到的主要结果是可能在实际上任何几何形状的通道中使热传送增加1.5～3倍。在二维凸起具有三角形或矩形、或者是它们构成围绕流线表面的绕线形状时，为了达到这种热传送的强化，要求粗糙单元(roughness elements)的绝对高度是传递通道等效直径的0.02～0.03倍。对于一个圆柱形流线表面，二维横槽的适当深度是传送通道等效流体直径的0.03～0.05倍。通过二维人工粗糙表面而获得热和物质传送增强的方法，对于热载体流动的瞬态模型和续发湍流模型是有效的。但是，使用这种方法要大幅增加流体阻力，超过光滑通道的4～10倍。流体阻力的大幅增加是限制能量传送装置中为增强热传送而使用二维人工表面粗糙方法的主要原因。
专利GB 2 263 525给出一个流体控制装置，包含由磁流体动力学执行器列阵组成的承受流体通过的体表面。执行器在与导致预先确定的边界层状态的表面相平行的方向上施加一个力在导电流体上。WO-A-94/00342也是有关在流体流动时控制边界层的装置和方法，该流体具有预先确定的电导率，并相对于一表面运动。该装置包括：改变电导的装置，产生磁场的装置和产生电流的装置。
“单相流动的热传送管”(EP A-85 107 484.9)的作者提出了热交换管的权利要求，该热交换管的内表面至少具有一排沿螺旋曲线形的凸起(凹槽)，而这些凸起的横截面的轮廓线在凸起的任何高度(包括底部)均是光滑的曲线。此时，朝着凸起最高点方向横截面面积单调地减少，而凸起的高度是0.45～0.6mm。选择螺旋曲线，使得沿圆周变化的间距(圆周间距，方位角间距)是3.5～5m，沿轴向变化的间距是5～15mm。然而，专利的作者并没有指出凸起的尺寸与间距之间的关系，这些凸起的尺寸和间距可以说明凸起的位置，管的直径和热载体流动情况等特征。除此以外，作者也没有说明截面曲线的半径，在该截面上，管表面光滑部分与凸起表面相共轭。
同时，从Kahnin E.K.，Draitser G.A.和Yarko S.A.等人在“通道中热传送的强化”(“Intensification of Heat Transfer in Channels”，莫斯科MACHINOSTROYENIE出版社，1981出版)一书第206页可知，这些曲线半径决定了流体阻力值，从而也确定了热和物质的有效性。
在引作参考的PCT/RU 92 00 106中，通过详细的方程式定义了包含有很特殊形状波纹的流线形表面，该详细的方程式至少花费了申请文件的一页半篇幅。在这个方程式里，为确定一个静态表面形状，至少涉及8个不同的变量。特别是对有关流动的介质与速度的关系方面仍然没有给出什么提示。
本发明目的是提供一种控制连续介质边界或壁层的方法和装置。介质包括了气体、液体和/或其混合物。该方法很容易应用于改变所说边界层或壁层中流体的结构，湍流程度，动量传递以及热量和/或混合物的传送。该装置的优点是减少流体阻力，增加热和动量的传送。因此，可以在大量技术应用中使用。
本发明更进一步的目的是在不需要任何复杂而明确确定的表面形状的情况下，给出有关形成二阶类龙卷风涡流的一种非常简单和基本的方法。
上述目的可以通过以下方法和装置来达到：
一种控制在表面(1)附近连续介质边界或壁层的方法，连续介质由气体、液体和/或其混合物组成，通过对流体施加作用改变连续介质粒子的速度导致类龙卷风涡流的形成，以改变流体结构，湍流程度、冲量的传送、热和/或混合物的传送，
连续介质流受到力场的作用，该力场至少延表面(1)法线方向yn距离范围之内的壁区里，
连续介质粒子的速度矢量的转向是由于所述力场按空间和/或时间上重复作用引起的，
所说的作用力引起了所说的转向，在α角范围内向表面(1)和离开表面(1)交替转向；在β角范围内，相对于近壁流体的连续介质粒子的速率矢量方向，向左和向右交替转向，
其特征在于，所述yn是边界层厚度δ的0.005～0.3倍，或压力通道的等效流体直径，或近壁流体的特征流体尺寸；
所说α角为0.02～0.5弧度；
所说β角为0.02～0.3弧度；
所说作用的强度或所说力的强度是使得所说粒子流轨迹的最小曲率半径Rmin为沿流线形壁的垂线到粒子弯曲轨迹的平均距离(S)的2～30倍，
下面的任何一个或全部特征a)和/或b)都是有效的：
a)沿壁流流动方向，所说作用力的空间重复范围λ‖＝(3～30)yn，在垂直于壁流流动方向，所说作用力的空间重复范围λ⊥为(1～10)yn，
b)时间重复周期是距离yn的(3～30)倍除以边界层或壁层上的平均速度 v。
其中，所说类龙卷风涡流的特征是一个具有非零螺旋性v·rot(v)≠0的二阶类龙卷风涡流。
至少一部分作用在壁流上的作用力是利用空间和/或时间上交变的磁场实现的，所述磁场集中在壁区域里。
至少一部分作用在壁层上的作用力是通过磁场和电场共同实现的，磁场和电场按时间和/或空间并替作用，所述磁场和电场集中在壁区域里。
至少一部分作用在流体上的作用力是利用所说表面(1)的形状按空间和/或时间改变并在流动的连续介质里产生数值和方向均改变的压力梯度实现的。
且至少一部分作用在近壁流上的作用力是在所说表面(1)的多个位置上同时注入和吸出一定量的连续介质实现的，该注入和吸出是按空间和/或时间变化的。
至少一部分作用在近壁流上的作用力是利用以下两种方法组合实现的：一种方法是在该表面(1)的多个位置上同时注入和吸出一定量的连续介质，另一种方法是利用该表面(1)的形状按空间和/或时间变化，使流动的连续介质里产生数值和方向均改变的压力梯度。
另外，本发明提供一种控制包含气体、液体和/或其混合物的连续介质边界层或壁层的装置，该装置利用上述的方法，通过对流体施加作用力和改变连续介质粒子的速度矢量以改变所说层里流体的结构，湍流水平、冲量的传送、热和/或混合物的传送，该装置包含了表(1)，对连续介质(3)的流动施加作用力的源(2)，所说作用加力场组成，
连续介质流至少在其近壁区域内受到所说力场的作用，作用范围是沿表面(1)的法线方向距离yn的范围内，
其特征在于：
所说范围yn是边界层厚度δ的0.005～0.3倍，或压力通道等效流体直径，或壁流的特征流体尺寸，
所说作用力按空间和/或时间上重复作用，引起连续介质粒子速度矢量的转向，所说作用引起所说的转向，在角α为0.02～0.5弧范围内，向表面(1)和离开表面(1)之间交替转向，在角β为0.02～0.3弧底范围内，相对于连续介质近壁流的连续介质粒子速度矢量的方向上向左和向右交替转向，
所说作用的强度是使所说粒子流轨迹的最小曲率半径Rmin为沿流线型壁的垂线到粒子弯曲轨迹的平均距离S的2～30倍，
以下特征a)和b)不论是其中一个还是两者均是有效的：
a)沿壁流方向时，所说作用力的空间重复范围λ‖＝(3～30)yn，在垂直于壁流体方向，所说作用力的空间重复范围λ⊥为(1～10)yn，
b)时间重复周期为距离yn的3～30倍除以边界层或壁层处的平均速度 v，引起类龙卷风涡流的形成。
其中，所说的作用源(2)以确定区域形式形成，它产生了磁和/或电场，和/或这种作用源是表面形变，和/或多人流体注入和吸出的位置，
所说作用源(2)分别是线性和/或区域单元，
在表面形变的情况，至少有一个形变是边界层厚度δ的0.005～0.3倍，或压力通道等效流体直径，或壁流的特征流体尺寸。
且所说的源(2)是由凹面和/或凸面形成的在表面上的作用源，
凹面的深度或凸面的高度是边界层厚度δ的0.005～0.3倍，或压力通道的等效流体直径，或壁流的特征流体尺寸，凹面或凸面主要部分的最小曲率半径Rmin是(0.1～1.0)δ。
另外，所说的源是以空洞的方式使作用力作用在流体上，所说空洞上是沿空洞的周边固定的薄膜，空洞直径是边界层厚度δ的0.01～0.6倍，或是压力通道的等效流体直径，或是壁流特征流体尺寸，薄膜中心部分的位移范围是(0.005～0.3)δ。
所说作用源设计成空洞的形式，空洞的直径是边界层厚度δ的0.005～0.1倍，或是压力通道的等效流体直径，或是壁流的特征流体尺寸，这些空洞按照2～10个空洞直径的间距排列，注入或吸出所说空洞的连续介质速度矢量是限定在与所说表面成0.1～1.0弧度的角度范围。
值得注间的是发明者找到了一个形成二阶类龙卷风涡流很基本的解。它是流体介质速度的函数，通过许多不同的物理力和作用，可以得到这个解。而各从属权利要求主要是对较佳实施例一种选择的详细说明。该基本解是针对流体结构本身，特别针对具有非零螺旋性v·rot(v)≠0的二阶类龙卷风涡流的形成。因此在基本技术教诲的广泛意义上，本发明不受任何独特的装置或较佳实施例的限制。
下面进而更详细地给出在不同形体的影响下如何获得其流体结构的完整的方法。
附图的简单描述
本发明连同几个较佳实施例可以从以下的结合附图的详细说明中得到了很好的了解。
图1：流体穿过简单的凹槽形成涡流，特别是
图1-1：流体层状区内，对称涡流结构的照片和图；
图1-2：流体湍流区内，反对称涡流结构照片和图；
图1-3：图1-2所描述的涡流结构图。
图2：流体穿过凹型格栅的流体结构的目视图。
图3：由平均速度向量和第一表面法线所确定平面的横截面，更一般的发明实施例也表明了在沿壁流方向(λ‖)上的作用具有空间重复性。
图4：表面的正视图，源自图3。
图5～图8：在表面附近产生电场和磁场的导体图。
图9和图10：由包含有对构成薄膜的可变形表面单元的作用源的表面的正视图。
图11：包含有通过可注入和吸走部分流动介质的空洞的作用源头的表面，图为穿过这一表面的截面图。
图12：流体阻力以cd-系数度量，图中以流动介质的速度为函数，将本发明的表面3的cd-系数同已有技术的人工表面4相比较。
图13：流体阻力以cd-系数度量，图中以表面相对于流体介质的倾斜角(攻角)为函数，将图12中的两表面的cd-系数作一比较。
较佳实施例的详细说明
在以下讨论中，黑体字描述了矢量。有下划线的数值是各个物理量的平均数值。为了清楚和更好地理解以下发明概念，给出以下的定义。重要表述的定义
表面(1)：在本发明的实施例中“表面(1)”可以是任何封闭或敞开容器的内表面或外表面。特别指，任何四周是流动气体，液体和/或其混合物的管道或容器的内表面，或是飞机、船舶或汽车的外表面。
粒子：在本发明的实施例中“粒子”可以是连续介质流的元(小的)体积(单元，部分)或在气体，液体和/或其混合物中的固体粒子流。
近壁流：连续介质的近壁流是指在流线形壁(表面)的直接近邻(最近的)连续介质流动。
壁区：是表面(1)附近的空间区域。该区域内，连续介质流动受表面(1)存在的影响。
边界层厚度δ：是沿流线形壁表面法线方向上的一段距离。在该距离，流速达到外部潜在流速的0.99(参照莫斯科NAUKA出版社1974年出版的Schlichting G.所著的“边界层理论”(“Theory of Boundary Layer”)第712页，在此引作参考)。
压力通道的等效流体直径：是该通道四倍横截面面积与横截面周长之比。详见莫斯科NAUKA出版社1974年出版的Schlichting G.“边界层理论”(“Theoryof Boundary Layer”)第712页。
近壁流的特征流体尺寸：是对应于流速较大变化壁区的距离。
平均距离S：沿着流线形壁的法线从壁到粒子曲线轨迹的平均距离S是沿流线形壁的法线从该壁到沿曲线轨迹运动的粒子之间的最大和最小距离之和的一半。
基本物理概述：形成本发明的重要方面是，研制一种控制流体阻力、热和质量传送过程以及在流体边界层或壁层上的混合物的蒸发、沉积等的方法和装置的问题。流体成份可以是气体、液体或它们的两相和多成份混合物，其运动状态可以是理想的无压力，或者处于压力通道状态。后者是因开始形成大规模涡流结构以及因各种动力(power)作用而导致大规模涡流结构进一步发展而引起的动力作用包括：压力、由于表面形变(静态或非静态)而产生的流体通过表面流动、在流体壁区集中的交变电磁场、连续介质的注入以及吸走。
通过对连续介质流施加作用来解决确定的任务。该作用是利用沿流线形表面法线上距离yn之内的流体壁区域里的力场。该距离yn是边界层厚度δ的0.005～0.3或压力通道的等效流体直径或者壁流的特征流体尺寸。该作用力使连续介质粒子的速率矢量在空间和/或时间上交替改变。速度矢量的改变的方式是以角度α＝0.02～0.5弧度向流线形表面和从该表面离开，以及相对于连续介质壁流运动方向在角度β＝0.02～0.3弧度范围之间左右转向。在这种情形，作用的强度是使得在离开壁指定距离范围内力场作用下的连续介质粒子轨迹的最小曲率半径Rmin是沿流线形壁的垂线到粒子曲线轨迹的平均距离S的2到30倍。而沿壁流动方向作用力的空间重复范围为λ‖＝(3～30)yn。在穿过壁流的方向，  作用力的空间重复范围为λ⊥＝(1～10)yn。时间重复周期T为距离yn3到30倍除以边界层或壁层处的平均速率v。这为二阶类龙卷风涡流的形成提供了依据。该涡流产生了边界层和壁层的特定结构，确定了强度等级、动量、热和混合物等传送的等级。
确定的问题还可以通过在流动区内及时对壁流相关的重复作用来解决。沿壁流方向，流动区域大小为L‖＝(15～150)yn；在穿过壁层方向，  流动区域大小为L⊥＝(5～50)yn。众所周知，连续介质的边界层或壁层的状态可以通过这些层的速度场v，通过雷诺胁强 v′ i v′ j；混合物湍流 v′ i c′和热湍流 v i ′T来表征。这里，v是一段时间间隔内的平均速度矢量；vi′是速度分量I的脉动；c′，T′是混合物和热的脉动。自然，下面各脉动值的乘积表示在时间上的平均。 v′i， v′i v j ′， v j ′ T′( v j ′c′)的动态值是用合适的方程式来描述的。该方程式确定了这些数值平均速度场正切线上的变化。详见列宁格勒ENERGOATOMIZDAT出版社1987年出版的Belov I.A.，KudryavtsevN.A.，所著“热传递和管壁的阻力”(“Heat Trans ferand Resistanceof Tube Banks”)第223页，其内容列作参考。在这些方程式中，主要的几项，也就是用于预计 v i ′v j ′和 v i ′T′值的变化方向的项，被称为雷诺胁强和热流产生项，用这些值和平均速度相对于坐标的导数的乘积来表示。按这种方法，平均速度矢量转向，会引起雷诺胁强项和热流产生项的变化。连续介质粒子轨迹的曲率半径越小，该项就越大，它们对于 v i′ v j ′和 v i ′T′值变化的累加作用是与平均速度矢量的转向角成比例。
作为规律，当平均速度矢量转向流线形表面，从流动的连续介质到流体经过的表面之间的冲量和热传送减少；当上述矢量从流体所经过的表面转过来时，冲量和热传送增加。速度矢量相对于壁流体方向，向左或向右转动，涉及到越过上述流体并与流体流过表面的法线正交的冲量传送。
速度矢量的转动对平均速度的偏移产生作用，即对于平均速度绝对值在垂直于平均速度矢量方向上的导数产生作用。改变雷诺胁强也会使速率分量相对于座标的导数变化。在平均速率按三维改变的条件下，随着流体管伸长，这些因素导致各类涡流结构形成，其中包括类龙卷风涡流，这些涡流结构在它们转向时又影响冲量，热和混合物的传送。
在从壁算起的yn距离内，力场会对于连续介质流起作用，同时还会引起连续介质粒子速度矢量转向。距离yn对应于在壁湍流机理中起着重要作用的相干大规模结构的形成和转变区域。在连续介质的湍流中，这个距离通常处在边界层厚度δ的0.005～0.3倍内，或压力通道的等效流体半径内，或壁流特征流动尺寸大小范围内。
通过及时地对相关壁流重复作用，可以控制连续介质边界或壁层。对壁流作用的轨迹可延伸至接近20yn的流体动力学长度的范围。在这种情况下，依靠最初作用力的种类和强度，在轨迹区域内，可以观察到冲量、热、或混合物的湍流所具有的单调行为。对壁流重复作用的相关性，允许增加对流体单一作用所希望的影响，也允许减少不希望的影响。
对流体的作用
对流体的作用可以通过在空间和/或时间上集中在壁区yn＝(0.005～0.3)δ的距离范围之内的交变的磁场或者电场和磁场联合实施。
对流体的作用可以通过改变流线形表面的形状来实现。如同圆周自身固定的可变形的薄膜，依靠空间和/或在时间上交替变化，导致数值和方向改变的压力梯度下实现变形。
通过连续介质按空间和/或时间交替变化地在流体所经过表面的各个断面上吹进和吸出，这也会对流体产生作用。
为了更好地理解这项发明的实质，请参见附图。
图3表示壁流区域图。图中示出 v的速率分布；连续介质粒子轨迹A；从流体流经的壁算起的距离yn的范围；壁流的边界层厚度δ或壁流体的特征流体尺寸；在力场F的作用下，连续介质粒子速度矢量朝流线形表面B或从B离开的转角α；从流线形壁到粒子弯曲轨迹A之间平均距离S；连续介质粒子轨迹A的最小曲率半径值Rmin，沿壁流体方向λ‖的作用的空间重复范围。
图4示出了在流线形表面B的凸起上的连续介质粒子的两条轨迹A1和A2。并表明了在力场F的作用下，连续介质粒子速度矢量 v相对于壁流方向向左和向右转角β和在穿过壁流方向λ⊥上的作用的空间重复范围。
图5、图6、图7和图8示出了导体a、b、c、d、e、f、g、h的示意图。这些导体是磁和/或电场的源。按顺序，导体a、b和d中的电流值分别是I1、I2和I3。对于技艺熟练的人都知道，为了产生电或磁力，要在导体a、b、c、d、e、f、g、h上施加能量(激励)。进一步将其连接到控制系统上，为了清楚起见省去了计算机控制系统。
图9和图10显示了流线形表面B上的区域C。该区域的形状按空间和/或时间交替变化。表面C可以是一个周边保持密封的弹性膜，从而被表面1另一面上传送流体的压力所激励。该流体压力可以由技艺熟练的人所熟知的方法对整个或部分薄膜进行控制。
图11示出流线形表面B。B面上带有将连续介质注入的空洞D，  还有将这些介质吸出的空洞E。
这里必须提到，以下方法可以实现控制连续介质的边界或壁层。如图3和图4所示，在壁区内，沿从流线形表面B的法线方向距离yn范围内，连续介质流体受到力场F的作用。这个范围是边界层厚度δ的0.005～0.3倍，或是压力通道的等效流体直径，或是壁流体的特征流体尺寸。借助于这种作用，连续介质粒子速度矢量按空间和/或时间随即交替转向，转向的角度α是在0.02～0.5弧度的范围穿过流线形表面，或者从该表面穿出，以及相对于连续介质壁流体的v或 v的方向，在β角0.02～0.3弧度范围，向左或向右转向。在这种情况下，作用强度是使从壁算起指定距离范围内施加力F场的作用之下的连续介质粒子轨迹A的最小曲率半径Rmin为沿从流线形表面B的法线从表面B到粒子弯曲轨迹A之间平均距离的2～30倍。而该作用的空间重复性在沿着流体方向为λ‖＝(3～30)yn，穿过流体方向为λ⊥＝(1～10)yn。时间重复周期T是3～30yn距离被边界层或壁层处平均速率v除。这样就为二阶类龙卷风涡流的形成作了准备，而这种涡流形成了边界层或壁层结构。这个结构又决定了湍流的等级、冲量、热和混合物的传送。它们所产生的类龙卷风三维结构，其特征是非零螺旋性v·rot(v)≠0。
按现代知识参见Levich.E.，shilman L.，Tur A.V.等人所著“流体动力学湍流中粘滞性的起因”(“The Origin of Coherence in HydrodynamicalTurbulence”)physici A 176，Tg，9I，第241～296页，或Belyan A.V.，Moiseev S.S.，Chkhetiani O.G.等人所著“在螺旋性湍流中湍流动的粘滞性”(科学院报告(俄罗斯)1994年334卷44～46页“About TurbulentViscosity in Helicity Turbulence”)，非零螺旋性涡流体区域诸如类龙卷风结构，导致了湍流范围能量传送的异常效应，导致了负值湍流粘度，导致了在热传送方向上雷诺模拟的紊乱。所以，类龙卷见涡流体控制和形成了边界或壁的流体的结构并产生螺旋性湍流。
通过图5～图11所示的装置可以作用于流体。最终，电流通过导体a、b、c、d、e、f、g(图5～图8)给它们加上电位(特别是导体h)。在这种情况下，穿过导体的电流产生了磁场，由于在连续介质的壁层中感应的电流流动，磁场作用于壁流体，其中包括由于电位差而形成的作用。根据上述方法，源动力F引起了连续介质粒子速度矢量的转动，F＝σ[E×B+B×(B×V)]。电流导体有各种构形，特别有可能是线性或区域性的。在这种情形下，导体的尺寸至少有一个是(0.005～0.3)δ。
凹状和/或凸状(c)，不论哪一种作用于流体(图9，图10)。此时，凹面(凸面)的深度(高度)为(0.005～0.3)δ，凹面(凸面)主要部分的最小曲率半径Rmin为(0.1～1.0)δ。
利用固定在任何形状空洞周边的可移动薄膜(图7和图8示出的薄膜C示意图)实现对流体施加作用。在上述情况，空洞直径(0.01～0.6)δ。薄膜中央部分的位移是(0.005～0.3)δ。
将连续介质注入空洞内和吸出空洞外，实现对流体施加作用。在这种情况下，空洞直径为(0.005～0.1)δ，这些空洞以2～10倍空洞直径的间距排列。注入或吸出的连续介质的速度矢量与流线形表面构成一个夹角，该角度为0.1～1.0弧度。
提供凹面、凸面、薄膜，注入和吸出流体的孔洞等的装置，操作很简单。当流体通过凹面、凸面、薄膜、空洞时，形成了压力梯度场，根据上述方法，梯度又引起了连续介质粒子速率矢量在流体壁区域转向。
为了实践本发明的方法和实现本发明装置，可以使用上述不同的物理作用。此外，还有一项供检测和优化条件的试验和误差技术，该技术在本文的权利要求中详述。
对于给定的表面，可以从测量流体的速率矢量，这对技术娴熟的人是共知的；还可通过大家都知道的技术将该矢量演示，如激光风速计、测量热交换、流体中粒子位置目视记录等等。
首先，在表面1中或表面1附近，可从确定出作用源2的位置。通过测试流体结构可以提供速率矢量v方向改变的信息，也就是平均速度矢量v的方向改变的信息，可以提供与表面平行的延伸平面上和法线方向上的流体速度角分布信息，详见图1-1～图1-2。
在所得结果的基础上，可以修正所施加场的强度，即感应力的强度或形变的深度。从本质上看，这是在表面1法线方向上的分布的问题。
在测得的流体速度矢量的横向分布的基础上，第一源2作用的横向大小可以与权利要求的定义相适应。
在第一源2优化之后，可以在如权利要求中所确定的距离上确定出第二和更多源的位置。这样，对于几乎任何形状的任何表面以及对于几乎所有的专门的相关装置，都可以逐步地获得合适源的网格。这种网格的实施例参见图2。另外，源2的分布是在权利要求中专门确定范围之内随机的或按统计规律分布。显然，如船体、飞机翼或车体等，在不同地方有不同的局部速度。因此，源2的最终位置也将相应变化。
利用封闭容器，透明的初样模型可以在容器自身之内获得上述结果。
然而，关键在于发明的范围，不仅是使用位于上面分别描述的各个距离上的源2的场，而且使用源2的很致密的场。例如磁线圈、电位盘、薄膜表面形变、和/或流体注入吸出的空洞。当然这些并非同时采用。如果将个别错位源2的激励忽略不计，分别控制对不同速度的局部作用是可能的。
进一步的发明实施例是基于使用的计算机控制系统。如在船体的情况下，控制适当区域的冷却水的注入和喷出，该船体也可以有合适的电板或表面形变。
在飞机外表面，可以使用一个进气和出气口作为作用源。这些作用源与涡轮机的几个不同等级的压力相连接。进、出口可单独使用，也可以同合适的表面形变组合使用。
从图12和图13可以看到一个特殊发明实施例所产生的结果。图12示出了一个发明表面的cd-系数(流体阻力的重要度量)与流体介质速度的函数关系。图中还将它与以前的人工表面相比较，两种表面均放入风洞，按照技艺熟练的专业人员所共知方法测量cd-系数。本发明的表面包含了球形的弯凹，即相对于平面，凹面有最大深度0.5mm，半径约为4.5mm。弯凹边界之间的纵向间距，即流体方向的间距约为2.0mm。几排弯凹以一个排在另一个下面而构成交错列阵，排间距离大致为9.5mm。当介质是空气时，流体介质的速度按马赫数给出。当0.75马赫时，可以看到，先前的飞机表面曲线4的流体阻力比曲线3的阻力高出15％，曲线3也就是本发明表面的流体阻力。
为了示出倾角对发明效果的影响，将发明表面和先前的人工表面相对流体介质的速度倾斜。正如所预计的那样，当倾角为0度时，可以得到最低阻力，也就是最低cd-系数。此时与图12中示出的当介质流动基本上与表面相平行的情况一样。很清楚，在很宽倾角范围内，本发明对减少流体阻力的作用是有效的，倾角亦称攻角。

Claims (12)

1.一种控制在表面(1)附近连续介质边界或壁层的方法，连续介质由气体、液体和/或其混合物组成，通过对流体施加作用改变连续介质粒子的速度导致类龙卷风涡流的形成，以改变流体结构，湍流程度、冲量的传送、热和/或混合物的传送，
连续介质流受到力场的作用，该力场至少延表面(1)法线方向yn距离范围之内的壁区里，
连续介质粒子的速度矢量的转向是由于所述力场按空间和/或时间上重复作用引起的，
所说的作用力引起了所说的转向，在α角范围内向表面(1)和离开表面(1)交替转向；在β角范围内，相对于近壁流体的连续介质粒子的速率矢量方向，向左和向右交替转向，
其特征在于，所述yn是边界层厚度δ的0.005～0.3倍，或压力通道的等效流体直径，或近壁流体的特征流体尺寸；
所说α角为0.02～0.5弧度；
所说β角为0.02～0.3弧度；
所说作用的强度或所说力的强度是使得所说粒子流轨迹的最小曲率半径Rmin为沿流线形壁的垂线到粒子弯曲轨迹的平均距离(S)的2～30倍，
下面的任何一个或全部特征a)和/或b)都是有效的；
a)沿壁流流动方向，所说作用力的空间重复范围λ‖＝(3～30)yn，在垂直于壁流流动方向，所说作用力的空间重复范围λ⊥为(1～10)yn，
b)时间重复周期是距离yn的(3～30)倍除以边界层或壁层上的平均速度 v。
2.如权利要求1所述的方法，所说类龙卷风涡流的特征是一个具有非零螺旋性v·rot(v)≠0的二阶类龙卷风涡流。
3.如权利要求1所述的方法，其特征是至少一部分作用在壁流上的作用力是利用空间和/或时间上交变的磁场实现的，所述磁场集中在壁区域里。
4.如权利要求1所述的方法，其特征是至少一部分作用在壁层上的作用力是通过磁场和电场共同实现的，磁场和电场按时间和/或空间交替作用，所述磁场和电场集中在壁区域里。
5.如权利要求1所述的方法，其特征是至少一部分作用在流体上的作用力是利用所说表面(1)的形状按空间和/或时间改变并在流动的连续介质里产生数值和方向均改变的压力梯度实现的。
6.如权利要求1所述的方法，其特征是至少一部分作用在近壁流上的作用力是在所说表面(1)的多个位置上同时注入和吸出一定量的连续介质实现的，该注入和吸出是按空间和/或时间变化的。
7.如权利要求1所述的方法，其特征是至少一部分作用在近壁流上的作用力是利用以下两种方法组合实现的：一种方法是在该表面(1)的多个位置上同时注入和吸出一定量的连续介质，另一种方法是利用该表面(1)的形状按空间和/或时间变化，使流动的连续介质里产生数值和方向均改变的压力梯度。
8.一种控制包含气体、液体和/或其混合物的连续介质边界层或壁层的装置，该装置利用权利要求1-7所述的方法，通过对流体施加作用力和改变连续介质粒子的速度矢量以改变所说层时流体的结构，湍流水平、冲量的传送、热和/或混合物的传送，该装置包含了表面(1)，对连续介质(3)的流动施加作用力的源(2)，所说作用力由力场组成，
连续介质流至少在其近壁区域内受到所说力场的作用，作用范围是沿表面(1)的法线方向距离yn的范围内，
其特征在于：
所说范围yn是边界层厚度δ的0.005～0.3倍，或压力通道的等效流体直径，或壁流的特征流体尺寸，
所说作用力按空间和/或时间上重复作用，引起连续介质粒子速度矢量的转向，所说作用引起所说的转向，在角α为0.02～0.5弧度范围内，向表面(1)和离开表面(1)之间交替转向，在角β为0.02～0.3弧度范围内，相对于连续介质近壁流的连续介质粒子速度矢量的方向上向左和向右交替转向，
所说作用的强度是使所说粒子流轨迹的最小曲率半径Rmin为沿流线型壁的垂线到粒子弯曲轨迹的平均距离S的2～30倍，
以下特征a)和b)不论是其中一个还是两者均是有效的：
a)沿壁流方向时，所说作用力的空间重复范围λ‖＝(3～30)yn，在垂直于壁流体方向，所说作用力的空间重复范围λ⊥为(1～10)yn，
b)时间重复周期为距离yn的3～30倍除以边界层或壁层处的平均速度v，引起类龙卷风涡流的形成。
9.如权利要求8所述的装置，其特征是所说的作用源(2)以确定区域形式形成，它产生了磁和/或电场，和/或这种作用源是表面形变，和/或多个流体注入和吸出的位置，
所说作用源(2)分别是线性和/或区域单元，
在表面形变的情况，至少有一个形变是边界层厚度δ的0.005～0.3倍，或压力通道等效流体直径，或壁流的特征流体尺寸。
10.如权利要求8或9所述的装置，其特征是所说的源(2)是由凹面和/或凸面形成的在表面上的作用源，
凹面的深度或凸面的高度是边界层厚度δ的0.005～0.3倍，或压力通道的等效流体直径，或壁流的特征流体尺寸，凹面或凸面主要部分的最小曲率半径Rmin是(0.1～1.0)δ。
11.如权利要求10所述的装置，其特征是所说的源是以空洞的方式使作用力作用在流体上，所说空洞上是沿空洞的周边固定的薄膜，空洞直径是边界层厚度δ的0.01～0.6倍，或是压力通道的等效流体直径，或是壁流特征流体尺寸，薄膜中心部分的位移范围是(0.005～0.3)δ。
12.如权利要求8或9所述的装置，其特征是所说作用源设计成空洞的形式，空洞的直径是边界层厚度δ的0.005～0.1倍，或是压力通道的等效流体直径，或是壁流的特征流体尺寸，这些空洞按照2～10个空洞直径的间距排列，注入或吸出所说空洞的连续介质速度矢量是限定在与所说表面成0.1～1.0弧度的角度范围。
CNB961969989A 1995-07-19 1996-07-19 控制连续介质的边界或壁层的方法和装置 Expired - Fee Related CN1138967C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP95111337.2 1995-07-19
EP95111337 1995-07-19

Publications (2)

Publication Number Publication Date
CN1196790A CN1196790A (zh) 1998-10-21
CN1138967C true CN1138967C (zh) 2004-02-18

Family

ID=8219446

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB961969989A Expired - Fee Related CN1138967C (zh) 1995-07-19 1996-07-19 控制连续介质的边界或壁层的方法和装置

Country Status (7)

Country Link
US (1) US6119987A (zh)
EP (1) EP0839309B1 (zh)
CN (1) CN1138967C (zh)
AU (1) AU6701596A (zh)
DE (1) DE69620185T2 (zh)
ES (1) ES2172673T3 (zh)
WO (1) WO1997004280A1 (zh)

Families Citing this family (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE9703217L (sv) * 1997-09-08 1998-12-21 Vattenfall Ab Diffusor
AU2002323881A1 (en) * 2001-07-05 2003-01-21 Vida, Nikolaus Tornado-type wind turbine
US6644921B2 (en) 2001-11-08 2003-11-11 General Electric Company Cooling passages and methods of fabrication
UA89470C2 (uk) * 2002-04-18 2010-02-10 Эйрбас Дойчланд Гмбх Пристрій керування ламінарною течією і літальний апарат, який використовує такий пристрій
US20060099122A1 (en) * 2002-11-25 2006-05-11 Vida Nikolaus M Method and device for generating mixtures of fluids in a boundary layer
US20040107718A1 (en) * 2002-12-06 2004-06-10 Michael Bowman Method, system and apparatus for cooling high power density devices
ATE413535T1 (de) 2003-03-19 2008-11-15 Vida Nikolaus Dr Dreidimensionale oberflächenstruktur für reduzierten strömungswiderstand und verbesserten wärmeübergang
EP1606512A1 (en) 2003-03-21 2005-12-21 Vida, Nikolaus, Dr. Tornado-type wind or water turbine
DE10347022A1 (de) * 2003-10-07 2005-05-04 Nikolaus Vida Oberfläche mit reduzierter Partikelablagerung und reduzierter Eisbildung
US7510149B2 (en) * 2004-08-02 2009-03-31 Lockheed Martin Corporation System and method to control flowfield vortices with micro-jet arrays
US20060102801A1 (en) * 2004-11-01 2006-05-18 The Boeing Company High-lift distributed active flow control system and method
US20060099073A1 (en) * 2004-11-05 2006-05-11 Toufik Djeridane Aspherical dimples for heat transfer surfaces and method
US20090090423A1 (en) * 2005-03-04 2009-04-09 Gennady Iraklievich Kiknadze Method of forming a current that generates Tornado Like Jets (TLJ) embedded into the flow, and the surface for its implementation
US7635107B2 (en) * 2005-08-09 2009-12-22 The Boeing Company System for aerodynamic flows and associated method
US8033510B2 (en) * 2005-08-09 2011-10-11 The Boeing Company Lift augmentation system and associated method
US7874525B2 (en) * 2006-05-04 2011-01-25 Lockheed Martin Corporation Method and system for fully fixed vehicle control surfaces
EP2103818B1 (en) * 2006-08-31 2013-09-18 Gennady Iraklievich Kiknadze Friction reducing surface and a mass and heat transfer enhancing surface
GB0803730D0 (en) * 2008-02-29 2008-04-09 Airbus Uk Ltd Shock bump array
US10352171B2 (en) 2008-11-01 2019-07-16 Alexander J. Shelman-Cohen Reduced drag system for windmills, fans, propellers, airfoils, and hydrofoils
US20100219296A1 (en) * 2008-11-01 2010-09-02 Alexander J. Shelman-Cohen Reduced drag system for windmills, fans, propellers, airfoils, and hydrofoils
WO2010005337A1 (ru) 2008-12-29 2010-01-14 Kiknadze Gennady Iraklievich Преобразователь и способ преобразования энергии, смерчевой нагнетатель и турбина
DE102011013572A1 (de) * 2011-03-10 2012-09-13 Norma Germany Gmbh Durchflussoptimierte Fluidleitung
GB2497139B (en) * 2011-12-02 2015-11-11 Vkr Holding As Phase change material pack
US9845902B2 (en) 2012-05-13 2017-12-19 InnerGeo LLC Conduit for improved fluid flow and heat transfer
WO2015178988A2 (en) 2014-03-03 2015-11-26 The Florida State University Research Foundation, Inc. Swirling jet actuator for control of separated and mixing flows
JP6060236B1 (ja) * 2015-09-30 2017-01-11 富士重工業株式会社 インストルメントパネル用空気流通装置
JP6560975B2 (ja) * 2015-12-17 2019-08-14 三菱重工機械システム株式会社 自流型境界層制御装置及びこれを用いる風洞試験装置
US10683076B2 (en) 2017-10-31 2020-06-16 Coflow Jet, LLC Fluid systems that include a co-flow jet
US10718362B2 (en) * 2017-11-09 2020-07-21 The Florida State University Research Foundation, Inc. Systems and methods for actively controlling a vortex in a fluid
US11261890B2 (en) * 2017-11-29 2022-03-01 Khaled Abdullah Alhussan High speed rotating bodies with transverse jets as a function of angle of attack, reynolds number, and velocity of the jet exit
US11149766B2 (en) 2018-08-24 2021-10-19 Quest Engines, LLC Controlled turbulence system
GB2600584B (en) 2019-07-23 2024-03-06 Coflow Jet Llc Fluid systems and methods that address flow separation
DE102021004807A1 (de) 2020-10-07 2022-04-07 Mathias Herrmann Antriebskonzept zur Kombination von konventionellen Raketentriebwerken und luftatmenden Triebwerken (Heber-Konzept)
CN113090946B (zh) * 2021-04-06 2022-04-05 西南石油大学 一种促进掺氢天然气管道内天然气与氢气混合的蚌式管道结构
DE102021004784A1 (de) 2021-09-22 2023-03-23 Mathias Herrmann Adaptives Triebwerkskonzept (z.B. Raketen, Überschallraketenbrennkammern, luftatmende Triebwerke wie z.B. Düsentriebwerke, Unterschall- Staustahltriebwerken, Ramjets Scramjets, Dualmode, Pulsejets, Detonationsstriebwerke, Raketen-Staustrahl-Antriebe, Laval-/Düsen) und Strukturkonzept (z.B. für Flugzeuge und Raketen) für flexible und gezielt veränderliche Geometrien
CN114048661B (zh) * 2021-11-29 2022-09-16 中南大学 一种层状介质下直流点源电位和电场处理方法和装置
CN116395131A (zh) * 2023-03-09 2023-07-07 中国科学院力学研究所 一种超高速飞行器稀薄流动减阻方法及系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR1346271A (fr) * 1962-06-26 1963-12-20 Perfectionnement des échangeurs de chaleur
CH533289A (de) * 1970-11-10 1973-01-31 R Dr Lage James Verfahren zur Verbesserung der Wärmeübertragung und Einrichtung zur Durchführung des Verfahrens
DE3228939C1 (de) * 1982-08-03 1983-11-24 Messerschmitt-Bölkow-Blohm GmbH, 8000 München Verfahren und Einrichtung zur Beeinflussung der Grenzschicht von umstroemten Koerpern
US4932612A (en) * 1986-02-25 1990-06-12 Blackwelder Ron F Method and apparatus for reducing turbulent skin friction
US5386955A (en) * 1986-05-22 1995-02-07 Rolls-Royce Plc Control of fluid flow
EP0588788A1 (en) * 1990-02-08 1994-03-30 RUSSELL, Francis, Michael Fluid flow apparatus
GB9118646D0 (en) * 1991-08-30 1992-02-19 Ruston Gas Turbines Ltd Flow control arrangement
US5362179A (en) * 1991-11-19 1994-11-08 Lawrence Sirovich Method of and apparatus for controlling turbulence in a wall-bounded fluid flow field
US5320309A (en) * 1992-06-26 1994-06-14 British Technology Group Usa, Inc. Electromagnetic device and method for boundary layer control
US5346745A (en) * 1993-06-01 1994-09-13 The United States Of America As Represented By The Secretary Of The Navy Elastic micro-fabricated surface layer for reducing turbulence and drag on an object while it moves through a fluid medium

Also Published As

Publication number Publication date
CN1196790A (zh) 1998-10-21
WO1997004280A1 (en) 1997-02-06
US6119987A (en) 2000-09-19
DE69620185D1 (de) 2002-05-02
AU6701596A (en) 1997-02-18
EP0839309B1 (en) 2002-03-27
DE69620185T2 (de) 2002-11-07
EP0839309A1 (en) 1998-05-06
ES2172673T3 (es) 2002-10-01

Similar Documents

Publication Publication Date Title
CN1138967C (zh) 控制连续介质的边界或壁层的方法和装置
Tsukahara et al. Flow regimes in a plane Couette flow with system rotation
Ozgoren et al. Comparison of flow structures in the downstream region of a cylinder and sphere
Ozgoren et al. Investigation of flow characteristics around a sphere placed in a boundary layer over a flat plate
Lee et al. Aerodynamics of a translating comb-like plate inspired by a fairyfly wing
CN107972850A (zh) 一种高速绕流边界层的被动式降热减阻装置及方法
Koushan et al. Experimental and numerical study of pre-swirl stators PSS
Canpolat et al. Observation of the vortical flow over a yawed delta wing
Young et al. Numerical modeling of supercavitating propeller flows
Bakhuis et al. Controlling secondary flow in Taylor–Couette turbulence through spanwise-varying roughness
Ozturk et al. PIV measurements of flow past a confined cylinder
Ozen et al. Control of vortical structures on a flapping wing via a sinusoidal leading-edge
JP4824190B2 (ja) 乱流摩擦抵抗低減表面
Lin et al. Shape optimization and hydrodynamic simulation of a Magnus anti-rolling device based on fully parametric modeling
EP2103818A1 (en) Friction reducing surface and a mass and heat transfer enhancing surface
Ozalp et al. The effect of a high-speed rotational oscillating cylinder on a heated cylinder
RU2425260C2 (ru) Поверхность тела для уменьшения трения и поверхность тела для интенсификации теплообмена
Yayla Flow Characteristic of Staggered Multiple Slotted Tubes in the Passage of a Fin Tube Heat Exchanger.
Zharkova et al. Free stream turbulence effect on the flow structure over the finite span straight wing
WO2006098649A1 (fr) Procede de formation d'un courant de formation de jets tourbillonnants integres a un flux et surface conçue pour sa mise en oeuvre
Esfahani et al. COMPARATIVE INVESTIGATIONS IN THE EFFECT OF ANGLE OF ATTACK PROFILE ON HYDRODYNAMIC PERFORMANCE OF BIO-INSPIRED FOIL.
CN201661557U (zh) 减小运动阻力的钣金
García-Magariño et al. Micro/Bubble drag reduction focused on new applications
Go et al. Towed underwater PIV measurement of propeller wake in self-propelled condition
Li et al. Time-averaged flow field behind a transversely spinning sphere: An experimental study

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20040218

Termination date: 20100719