CN113312797A - 一种熔体超声空化强度计算方法及系统 - Google Patents

一种熔体超声空化强度计算方法及系统 Download PDF

Info

Publication number
CN113312797A
CN113312797A CN202110710107.5A CN202110710107A CN113312797A CN 113312797 A CN113312797 A CN 113312797A CN 202110710107 A CN202110710107 A CN 202110710107A CN 113312797 A CN113312797 A CN 113312797A
Authority
CN
China
Prior art keywords
noise
cavitation
cavitation noise
steady
sound pressure
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.)
Granted
Application number
CN202110710107.5A
Other languages
English (en)
Other versions
CN113312797B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110710107.5A priority Critical patent/CN113312797B/zh
Publication of CN113312797A publication Critical patent/CN113312797A/zh
Application granted granted Critical
Publication of CN113312797B publication Critical patent/CN113312797B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10Noise analysis or noise optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明公开了一种熔体超声空化强度计算方法及系统。该方法包括:对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。本发明普适性强、可准确描述超声空化强度的表征。

Description

一种熔体超声空化强度计算方法及系统
技术领域
本发明涉先进材料制备及加工技术领域,特别是涉及一种熔体超声空化强度计算方法及系统。
背景技术
在金属及合金凝固、化学反应、表面清洁过程中施加振动是改善其组织结构、调控化学反应过程、增强清洗效果最有效的方法之一。常规机械振动因存在振动频率低、能量密度受限、振动过程不易控制等问题,而逐渐被功率超声所取代。功率超声是物体超高机械振动能量的传播形式,具有频率高、功率密度大、束射性和方向性良好等显著物理特性,与液体、金属或合金熔体相互作用时存在压力场、热效应、声空化和声流等一系列非线性超声效应,显著影响晶体的形核和长大、晶粒的细化和粗化、溶质的扩散、固相粒子的弥散、熔体的流动、空化泡的运动等,表现出晶粒细化、组织均匀化与净化等宏观效果,可大大提高金属凝固、化学反应与表面清洁的性能与效果。
在上述过程中,超声空化强度直接决定了功率超声施加所获得的性能提升程度,是描述超声空化场的关键参数,也是开展金属凝固、化学反应调控、表面清洁中功率超声作用机理这一基础科学问题的重要参数。因此,建立描述超声空化场强弱的方法,即超声空化强度表征,对于超声空化场作用机理研究工作的开展至关重要。
相应的,目前超声空化强度表征常用形式及其局限性为:(1)以超声导入装置消耗的所消耗的电功率为表征方法,该方式与超声导入装置的具体结构、具体电路、超声换能器、变幅杆、施加方式等密切相关,即便是同一批装置,也会由于电子元器件、换能器之间的差异导致测量结果大相径庭,普适性差;(2)以空蚀法、化学反应法、声致发光法为代表的间接测量法,这类方法一般对被测环境要求比较严格,如化学反应法、空蚀法均要求被测环境为常温状态,声致发光法要求被测环境为透明状态等;(3)以水听器测量得到的超声空噪声信号的时域总能量为表征方法,该方式普适性较好,但空化场辐射声信号的时域总能量中同时包含了换能器噪声、稳态空化噪声、瞬态空化噪声三部分噪声的分能量。其中,换能器噪声与超声换能器状态、溶液声阻抗、电声转换效率有关,与空化效应无关,不应计入超声空化场强度,同时,稳态空化和瞬态空化所造成的物理效应完全不同,将二者相加无法准确表征空化强度;(3)以超声空化噪声信号中的某些分量的频域能量为表征方法,该方法通过测量空化噪声中基频、基频谐波与超谐波等线状谱分量的能量来实现空化强度的表征,但空化泡振动与溃灭时除了产生超谐波信号之外,通常还伴有强烈的宽带噪声,在频域上表现为可达MHz级的宽带连续谱信号,该方法忽略了这一分量的能量,其表征方法并不完善。
鉴于现有超声空化强度表征方法存在普适性差、无法准确反应空化场大小的现状,研发普适性强、可准确描述超声空化强度的表征方法成为当务之急。
发明内容
基于此,本发明的目的是提供一种熔体超声空化强度计算方法及系统。
为实现上述目的,本发明提供了如下方案:
一种熔体超声空化强度计算方法,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;
基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
进一步地,基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声频域能量;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声频域能量;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量以及所述稳态空化噪声频域能量,计算瞬态空化噪声频域能量。
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
频域能量计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;
能量比值计算模块,用于基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
一种熔体超声空化强度计算方法,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
进一步地,基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声声压级;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声声压级;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级。
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声压级差值计算模块,用于基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
一种熔体超声空化强度计算方法,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强;
基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
进一步地,基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声声压级;
基于所述驱动噪声声压级以及参考声压计算驱动噪声声强;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声声压级;
基于所述稳态空化噪声声压级以及参考声压计算稳态空化噪声声强;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级;
基于所述瞬态空化噪声声压级以及参考声压计算瞬态空化噪声声强。
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声强计算模块,用于基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强;
声强比值计算模块,用于基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明将超声空化场噪声按照空化泡行为分为驱动噪声、稳态空化噪声、瞬态空化噪声三部分,对应不同的空化效应作用机制,具有明确的物理意义。超声空化强度表征以能量占比这一相对物理量作为表征参数,其比值与超声导入装置之间的差异无关,表征方法具有普适性。
超声空化强度表征以声强占比作为表征参数,其比值为一相对量,与传感器灵敏度无关,表征方法不仅具有普适性,而且无需校准传感器。
超声空化强度表征以声压级之差作为表征参数,其差值消除了超声导入装置环境噪声影响,且使用了标准物理量表示空化强度,超声导入装置之间的差异转换为可量化的物理量,其表征方法同样具有普适性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一熔体超声空化强度计算方法的流程图;
图2为本发明实施例二熔体超声空化强度计算方法的流程图;
图3为本发明实施例三熔体超声空化强度计算方法的流程图;
图4为典型空化噪声频谱图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一:
如图1所示,一种熔体超声空化强度计算方法包括以下步骤:
步骤101:对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声。
步骤102:基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率。
步骤103:基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量。具体的:
步骤1031:确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
步骤1032:基于所述第一最高能量值计算驱动噪声频域能量;
步骤1033:确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
步骤1034:基于所述第二最高能量值计算稳态空化噪声频域能量;
步骤1035:基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量以及所述稳态空化噪声频域能量,计算瞬态空化噪声频域能量。
步骤104:基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
详细过程如下:
功率超声施加于金属熔体时所产生的声信号称之为空化噪声,本发明将按照熔体内空化泡行为将空化噪声分为三类:驱动噪声、稳态空化噪声、瞬态空化噪声。经模数转换之后,对空化噪声做傅里叶变换,得到典型频谱图(即空化噪声的频域表现形式),如图4所示,f0为功率超声的主频率。其中,驱动噪声与稳态空化噪声在频域上分别表现为中心频率为nf0和(n+0.5)f0的线谱(n为整数,其取值从1到15),瞬态空化噪声在频域上表现为除前述线谱之外的宽带连续谱。
驱动噪声频域能量ED、稳态空化噪声频域能量ES、瞬态空化噪声频域能量ET的计算方法为:
Figure BDA0003133281870000081
Figure BDA0003133281870000082
Figure BDA0003133281870000083
式中S(f)表示空化噪声s(t)的离散傅里叶变换的模值,即空化噪声的总能量;,f0表示施加于金属熔体中的功率超声主频率,是与超声换能器和功率放大器(二者功率超声施加装置)有关的已知固定值;nf0表示f0的整数倍,即f0的谐波成分,二者来自于驱动噪声,与空化场无关,其频谱形状呈线状,称为线谱;(n+0.5)f0为稳态空化噪声对应的频率,来自于空化场中空化泡参数振动与非线性脉动,与稳态空化有关;S(f)的剩余部分则来自于空化泡的混沌运动与溃灭,与瞬态空化有关。
在金属凝固过程中施加功率超声可显著提升金属性能,其原因在于空化场中大量空化泡的混沌运动和溃灭可分别产生除气与枝晶断裂的作用,其作用效果与空化泡的瞬态空化行为有关。针对上述特点,本发明提出一种用于描述金属熔体中空化强度的方法,该方法以瞬态空化噪声频域能量ET为基础,以ET占总能量比值以及ET与ED、ES的比值为表征方法,其计算公式为:
Figure BDA0003133281870000091
Figure BDA0003133281870000092
Figure BDA0003133281870000093
式中,QE表示瞬态空化噪声能量与空化噪声总能量的比值;QES表示瞬态空化噪声能量与稳态空化噪声能量的比值;QED分别表示瞬态空化噪声能量与驱动噪声能量的比值;三者从不同角度描述了瞬态空化噪声的强弱,即作为空化强度表征。
由ED、ES、ET的计算公式可知,nf0和(n+0.5)f0的取值决定了其各部分能量计算的准确性。然而由于超声换能器与功率放大器的为高功率器件,每次启动时,已知的固定主频率f0均会产生微小的随机频偏Δf,即主频率变为f0+Δf。对于f0的谐波nf0和(n+0.5)f0来说,由于倍数n的存在,会导致Δf被放大n倍,使得ED、ES、ET的计算出现误差。因此,本发明提出了一种新的计算方法用于消除误差,即在nf0和(n+0.5)f0的邻域内取能量最高值作为其能量值,算法如下:
Figure BDA0003133281870000101
Figure BDA0003133281870000102
Figure BDA0003133281870000103
其中,Max表示取域内最大值,Δf的取值与具体仪器有关,一般取不大于0.2f0的值即可满足要求。则相应的空化表征方法修正为
Figure BDA0003133281870000104
Figure BDA0003133281870000105
计算公式为:
Figure BDA0003133281870000106
Figure BDA0003133281870000107
Figure BDA0003133281870000108
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
频域能量计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;
能量比值计算模块,用于基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
本发明实施例一所述超声空化强度表征方法,将超声空化场噪声按照空化泡行为分为驱动噪声、稳态空化噪声、瞬态空化噪声三部分,对应不同的空化效应作用机制,具有明确的物理意义。将空化噪声能量按其对应的不同噪声类型分为驱动噪声频域能量ED、稳态空化噪声频域能量ES、瞬态空化噪声频域能量ET,并依据三者在频域上的差异给出了其计算公式。
与稳态空化相比,瞬态空化会产生局部高温高压与微射流,与熔体凝固后的组织均匀细化密切相关,即瞬态空化越强其空化效应愈明显。故以瞬态空化噪声能量与空化噪声总能量、稳态空化噪声能量、驱动噪声能量的比值QE、QES、QED作为超声空化强度表征,可更为准确的刻画熔体中空化效应的强弱。超声空化强度表征以能量占比这一相对物理量作为表征参数,其比值与超声导入装置之间的差异无关,表征方法具有普适性。
实施例二:
如图2所示,一种熔体超声空化强度计算方法包括以下步骤:
步骤201:对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声。
步骤202:基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
步骤203:基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级。具体的:
步骤2031:确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
步骤2032:基于所述第一最高能量值计算驱动噪声声压级;
步骤2033:确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
步骤2034:基于所述第二最高能量值计算稳态空化噪声声压级;
步骤2035:基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级。
步骤204:基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
详细过程如下:
功率超声施加于金属熔体时所产生的声信号称之为空化噪声,本发明将按照熔体内空化泡行为将空化噪声分为三类:驱动噪声、稳态空化噪声、瞬态空化噪声。经模数转换之后,对空化噪声做傅里叶变换,得到典型频谱图,如图4所示,f0为功率超声的主频率。其中,驱动噪声与稳态空化噪声在频域上分别表现为中心频率为nf0和(n+0.5)f0的线谱,瞬态空化噪声表现为除前述线谱之外的宽带连续谱。
驱动噪声声压级SPLD、稳态空化噪声声压级SPLS、瞬态空化噪声声压级SPLT的计算方法为:
Figure BDA0003133281870000121
Figure BDA0003133281870000122
Figure BDA0003133281870000123
式中Lm为传感器灵敏度,Lr为传感器调理放大电路的增益,二者单位均为dB;S(f)表示空化噪声s(t)的离散傅里叶变换,f0表示施加于金属熔体中的功率超声主频率,是与超声换能器和功率放大器有关的已知固定值;nf0表示f0的整数倍,即f0的谐波成分,二者来自于驱动噪声,与空化场无关,其频谱形状呈线状,称为线谱;(n+0.5)f0为稳态空化噪声对应的频率,来自于空化场中空化泡参数振动与非线性脉动,与稳态空化有关;S(f)的剩余部分则来自于空化泡的混沌运动与溃灭,与瞬态空化有关。
在金属凝固过程中施加功率超声可显著提升金属性能,其原因在于空化场中大量空化泡的混沌运动和溃灭可分别产生除气与枝晶断裂的作用,其作用效果与空化泡的瞬态空化行为有关。针对上述特点,本发明提出一种用于描述金属熔体中空化强度的方法,该方法以瞬态空化噪声声压级SPLT为基础,以SPLT本身、SPLT与总声压级SPLTotal之差QSPL、以及SPLT与SPLD之差QSPLD、以及SPLT与SPLS之差QSPLS为表征方法,其计算公式为:
Figure BDA0003133281870000131
Figure BDA0003133281870000132
QSPLD=SPLT-SPLD
QSPLS=SPLT-SPLS
由SPLD、SPLS、SPLT的计算公式可知,nf0和(n+0.5)f0的取值决定了其各部分声压级计算的准确性。然而由于超声换能器与功率放大器的为高功率器件,每次启动时,已知的固定主频率f0均会产生微小的随机频偏Δf,即主频率变为f0+Δf。对于f0的谐波nf0和(n+0.5)f0来说,由于倍数n的存在,会导致Δf被放大n倍,使得SPLD、SPLS、SPLT的计算出现误差。故,本发明提出了一种新的计算方法用于消除误差,即在nf0和(n+0.5)f0的邻域内取能量最高值作为其能量值,算法如下:
Figure BDA0003133281870000133
其中,
Figure BDA0003133281870000134
Figure BDA0003133281870000135
其中,
Figure BDA0003133281870000136
Figure BDA0003133281870000137
其中,Max表示取域内最大值,Δf的取值与具体仪器有关,一般取不大于0.2f0的值即可满足要求。则相应的空化表征方法修正为
Figure BDA0003133281870000138
Figure BDA0003133281870000139
计算公式为:
Figure BDA0003133281870000141
Figure BDA0003133281870000142
Figure BDA0003133281870000143
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声压级差值计算模块,用于基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
声压级之差作为表征参数,其差值消除了超声导入装置环境噪声影响,且使用了标准物理量声压级来表示空化强度,超声导入装置之间的差异转换为可量化的物理量,其表征方法同样具有普适性。
实施例三:
如图3所示,一种熔体超声空化强度计算方法包括以下步骤:
步骤301:对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声。
步骤302:基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率。
步骤303:基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级。
步骤304:基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强。具体包括:
步骤3041:确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
步骤3042:基于所述第一最高能量值计算驱动噪声声压级;
步骤3043:基于所述驱动噪声声压级以及参考声压计算驱动噪声声强;
步骤3044:确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
步骤3045:基于所述第二最高能量值计算稳态空化噪声声压级;
步骤3046:基于所述稳态空化噪声声压级以及参考声压计算稳态空化噪声声强;
步骤3047:基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级;
步骤3048:基于所述瞬态空化噪声声压级以及参考声压计算瞬态空化噪声声强。
步骤305:基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
详细过程如下:
功率超声施加于金属熔体时所产生的声信号称之为空化噪声,本发明将按照熔体内空化泡行为将空化噪声分为三类:驱动噪声、稳态空化噪声、瞬态空化噪声。经模数转换之后,对空化噪声做傅里叶变换,得到典型频谱图,如图4所示,f0为功率超声的主频率。其中,驱动噪声与稳态空化噪声在频域上分别表现为中心频率为nf0和(n+0.5)f0的线谱,瞬态空化噪声表现为除前述线谱之外的宽带连续谱。
驱动噪声声强ID、稳态空化噪声声强IS、瞬态空化噪声声强IT的计算方法为:
Figure BDA0003133281870000161
Figure BDA0003133281870000162
Figure BDA0003133281870000163
式中Lm为传感器灵敏度,Lr为传感器调理放大电路的增益,二者单位均为dB;Pref为参考声压,此处取为1uPa。SPLD表示驱动噪声声压级,SPLS表示稳态空化噪声声压级,SPLT表示瞬态空化噪声声压级。S(f)表示空化噪声s(t)的离散傅里叶变换,f0表示施加于金属熔体中的功率超声主频率,是与超声换能器和功率放大器有关的已知固定值;nf0表示f0的整数倍,即f0的谐波成分,二者来自于驱动噪声,与空化场无关,其频谱形状呈线状,称为线谱;(n+0.5)f0为稳态空化噪声对应的频率,来自于空化场中空化泡参数振动与非线性脉动,与稳态空化有关;S(f)的剩余部分则来自于空化泡的混沌运动与溃灭,与瞬态空化有关。
在金属凝固过程中施加功率超声可显著提升金属性能,其原因在于空化场中大量空化泡的混沌运动和溃灭可分别产生除气与枝晶断裂的作用,其作用效果与空化泡的瞬态空化行为有关。针对上述特点,本发明提出一种用于描述金属熔体中空化强度的方法,该方法以瞬态空化噪声声强IT为基础,以声强值IT本身、IT与总声强ITotal的比值QI、IT与ID、IS的比值QID、QIS为表征方法,其计算公式分别为:
Figure BDA0003133281870000164
Figure BDA0003133281870000165
其中
Figure BDA0003133281870000166
Figure BDA0003133281870000171
Figure BDA0003133281870000172
由ID、IS、IT的计算公式可知,nf0和(n+0.5)f0的取值决定了其各部分能量计算的准确性。然而由于超声换能器与功率放大器的为高功率器件,每次启动时,已知的固定主频率f0均会产生微小的随机频偏Δf,即主频率变为f0+Δf。对于f0的谐波nf0和(n+0.5)f0来说,由于倍数n的存在,会导致Δf被放大n倍,使得ID、IS、IT的计算出现误差。故,本发明提出了一种新的计算方法用于消除误差,即在,nf0和(n+0.5)f0的邻域内取能量最高值作为其能量值,算法如下:
Figure BDA0003133281870000173
其中,
Figure BDA0003133281870000174
Figure BDA0003133281870000175
其中,
Figure BDA0003133281870000176
Figure BDA0003133281870000177
Figure BDA0003133281870000178
其中,
Figure BDA0003133281870000179
其中,Max表示取域内最大值,Δf的取值与具体仪器有关,一般取不大于0.2f0的值即可满足要求。则相应的空化表征方法修正为
Figure BDA00031332818700001710
Figure BDA00031332818700001711
计算公式为:
Figure BDA0003133281870000181
其中
Figure BDA0003133281870000182
Figure BDA0003133281870000183
Figure BDA0003133281870000184
本发明还提供了一种熔体超声空化强度计算系统,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声强计算模块,用于基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强;
声强比值计算模块,用于基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
本发明所述超声空化强度表征方法,将超声空化场噪声按照空化泡行为分为驱动噪声、稳态空化噪声、瞬态空化噪声三部分,对应不同的空化效应作用机制,具有明确的物理意义。将空化噪声总声强按其对应的不同噪声类型分为驱动噪声声强ID、稳态空化噪声声强IS、瞬态空化噪声声强IT,并依据三者在频域上的差异,利用帕萨瓦尔定理给出了计算公式。与稳态空化相比,瞬态空化会产生局部高温高压与微射流,与熔体凝固后的组织均匀细化密切相关,即瞬态空化越强其空化效应愈明显。故以瞬态空化噪声声强IT、IT与总声强的比值QI、IT与ID的比值QID、IT与IS的比值QIS为表征方法,可更为准确的刻画熔体中空化效应的强弱。超声空化强度表征以声强绝对值为表征参数,为一意义明确的绝对物理量,具有普适性。超声空化强度表征以声强占比作为表征参数,其比值为一相对量,与传感器灵敏度无关,表征方法不仅具有普适性,而且无需校准传感器。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种熔体超声空化强度计算方法,其特征在于,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;
基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
2.根据权利要求1所述的熔体超声空化强度计算方法,其特征在于,基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声频域能量;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声频域能量;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量以及所述稳态空化噪声频域能量,计算瞬态空化噪声频域能量。
3.一种熔体超声空化强度计算系统,其特征在于,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
频域能量计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声频域能量、稳态空化噪声频域能量以及瞬态空化噪声频域能量;
能量比值计算模块,用于基于傅里叶变换后的空化噪声总能量、所述驱动噪声频域能量、所述稳态空化噪声频域能量以及所述瞬态空化噪声频域能量,计算用于表征空化强度的能量比值。
4.一种熔体超声空化强度计算方法,其特征在于,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
5.根据权利要求4所述的熔体超声空化强度计算方法,其特征在于,基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声声压级;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声声压级;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级。
6.一种熔体超声空化强度计算系统,其特征在于,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声压级差值计算模块,用于基于总声压级、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级,计算用于表征空化强度的声压级差值。
7.一种熔体超声空化强度计算方法,其特征在于,包括:
对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强;
基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
8.根据权利要求7所述的熔体超声空化强度计算方法,其特征在于,基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强,具体包括:
确定所述驱动噪声中心频率的邻域内的最高能量值,为第一最高能量值;
基于所述第一最高能量值计算驱动噪声声压级;
基于所述驱动噪声声压级以及参考声压计算驱动噪声声强;
确定所述稳态空化噪声中心频率的邻域内的最高能量值,为第二最高能量值;
基于所述第二最高能量值计算稳态空化噪声声压级;
基于所述稳态空化噪声声压级以及参考声压计算稳态空化噪声声强;
基于所述傅里叶变换后的空化噪声总能量、所述驱动噪声声压级以及所述稳态空化噪声声压级,计算瞬态空化噪声声压级;
基于所述瞬态空化噪声声压级以及参考声压计算瞬态空化噪声声强。
9.一种熔体超声空化强度计算系统,其特征在于,包括:
频谱图确定模块,用于对空化噪声做傅里叶变换,得到频谱图;所述空化噪声为功率超声施加于金属熔体时产生的声信号;所述空化噪声包括驱动噪声、稳态空化噪声和瞬态空化噪声;
中心频率确定模块,用于基于所述频谱图确定所述驱动噪声以及稳态空化噪声在频域上的中心频率;
声压级计算模块,用于基于傅里叶变换后的空化噪声总能量、驱动噪声中心频率以及稳态空化噪声中心频率分别计算驱动噪声声压级、稳态空化噪声声压级以及瞬态空化噪声声压级;
声强计算模块,用于基于参考声压、所述驱动噪声声压级、所述稳态空化噪声声压级以及所述瞬态空化噪声声压级分别计算驱动噪声声强、稳态空化噪声声强以及瞬态空化噪声声强;
声强比值计算模块,用于基于总声强、所述驱动噪声声强、所述稳态空化噪声声强以及所述瞬态空化噪声声强,计算用于表征空化强度的声强比值。
CN202110710107.5A 2021-06-25 2021-06-25 一种熔体超声空化强度计算方法及系统 Active CN113312797B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110710107.5A CN113312797B (zh) 2021-06-25 2021-06-25 一种熔体超声空化强度计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110710107.5A CN113312797B (zh) 2021-06-25 2021-06-25 一种熔体超声空化强度计算方法及系统

Publications (2)

Publication Number Publication Date
CN113312797A true CN113312797A (zh) 2021-08-27
CN113312797B CN113312797B (zh) 2022-11-25

Family

ID=77380211

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110710107.5A Active CN113312797B (zh) 2021-06-25 2021-06-25 一种熔体超声空化强度计算方法及系统

Country Status (1)

Country Link
CN (1) CN113312797B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060262993A1 (en) * 2005-05-18 2006-11-23 Lite-On Technology Corporation Method and a device for reducing image noises
CN101530320A (zh) * 2009-03-31 2009-09-16 西安交通大学 聚焦超声空化及其微汽泡实时提取装置与检测方法
CN102590338A (zh) * 2012-01-19 2012-07-18 清华大学 基于共振峰的超声空化状态识别方法
CN104157295A (zh) * 2014-08-22 2014-11-19 中国科学院上海高等研究院 用于检测及抑制瞬态噪声的方法
JP6307684B1 (ja) * 2017-04-19 2018-04-11 本多電子株式会社 高音圧音場の音圧分析装置及び方法、超音波洗浄機、超音波処理機
CN110031083A (zh) * 2018-12-31 2019-07-19 瑞声科技(新加坡)有限公司 一种噪音总声压级测量方法、系统及计算机可读存储介质
CN111415408A (zh) * 2020-04-14 2020-07-14 西安交通大学 一种超声空化的微秒级多尺度时空成像及特征图谱计算方法与系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060262993A1 (en) * 2005-05-18 2006-11-23 Lite-On Technology Corporation Method and a device for reducing image noises
CN101530320A (zh) * 2009-03-31 2009-09-16 西安交通大学 聚焦超声空化及其微汽泡实时提取装置与检测方法
CN102590338A (zh) * 2012-01-19 2012-07-18 清华大学 基于共振峰的超声空化状态识别方法
CN104157295A (zh) * 2014-08-22 2014-11-19 中国科学院上海高等研究院 用于检测及抑制瞬态噪声的方法
JP6307684B1 (ja) * 2017-04-19 2018-04-11 本多電子株式会社 高音圧音場の音圧分析装置及び方法、超音波洗浄機、超音波処理機
CN110031083A (zh) * 2018-12-31 2019-07-19 瑞声科技(新加坡)有限公司 一种噪音总声压级测量方法、系统及计算机可读存储介质
CN111415408A (zh) * 2020-04-14 2020-07-14 西安交通大学 一种超声空化的微秒级多尺度时空成像及特征图谱计算方法与系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LUHOVSKYI, A. F等: "The Impact of Static Pressure on the Intensity of Ultrasonic Cavitation in Aqueous Media", 《JOURNAL OF WATER CHEMISTRY AND TECHNOLOGY》 *
梁召峰等: "空化噪声谱的分离", 《声学技术》 *

Also Published As

Publication number Publication date
CN113312797B (zh) 2022-11-25

Similar Documents

Publication Publication Date Title
JP5070993B2 (ja) 音処理装置、位相差補正方法及びコンピュータプログラム
US7269098B2 (en) Ultrasonic transmitter, ultrasonic transceiver and sounding apparatus
JP6335168B2 (ja) インテリジェント電子ホーンおよびその実施方法
US8848931B2 (en) Method and device for testing and calibrating electronic semiconductor components which convert sound into electrical signals
CN113312797B (zh) 一种熔体超声空化强度计算方法及系统
TW201526669A (zh) 揚聲器的驅動訊號的產生裝置及方法
JP2009502466A (ja) 空気に基づく液体の工業的消泡用マクロ音波発生器
JP2003262623A (ja) ガスの音響速度測定装置及び方法
EP3712571A1 (en) Signal processing circuit, related chip, flowmeter, and method
CN101296533A (zh) 扬声器振膜等效质量及等效力顺的测量方法
TWI811768B (zh) 具散熱和前饋式主動噪音控制功能之電子系統
JP3456742B2 (ja) 変圧器の騒音レベル予測方法
CN108304617A (zh) 浮筏结构宽频线谱振动噪声快速预报方法
JPH0723492A (ja) ディジタルマイクロホン
Uhlendorf et al. Ultrasonic absorption and velocity measurements at mW peak power level in the range 50-500 MHz with variable pathlength cell
JPS6119297A (ja) 低域補正回路
CN111998934A (zh) 一种声源声功率测试方法
KR20230056762A (ko) 사운드 필드 감지를 기반으로 한 저전력 음향학의 근접 경보장치 및 경보방법
JP4191812B2 (ja) 超音波診断装置
JPS59187132A (ja) 防振制御装置
Kuroyama Viscosity dependence of harmonics and ultra-harmonics of acoustic cavitation noise
Yokoyama et al. 3P4-5 Transient Characteristics of Acoustic Cavitation Noise after Starting Ultrasound Irradiation
JP2714153B2 (ja) 利得制御装置
CN118102190A (zh) 一种非线性位移系统的建模方法、芯片和电子设备
JP2005011858A (ja) プラズマを用いた半導体製造におけるμ波パワー設定方法並びに当該設定方法を適用した半導体装置の製造装置及び上記設定方法を用いた半導体装置の製造方法

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