一种非线性声场的基频和高次谐波频域有限差分计算方法技术

本发明专利技术提供一种非线性声场的基频和高次谐波频域有限差分计算方法，包括：得到频域上的非线性波动方程；计算出频域非线性波动方程的通解；基频采用线性近似表达式，求解得到二次谐波与基频的对应关系；利用得到的二次谐波与基频的对应关系，替代掉基频表达式中卷积项的二次谐波，对基频表达式进行修正；采用黎曼和近似，对二次谐波通解积分项写成便于数值求解的形式；基于声波波长确定空间步长，将基频总解代入到黎曼和，计算出黎曼和；将黎曼和代入到二次谐波通解中，求解出二次谐波在空间中的声场。本发明专利技术考虑了介质中基频和二次谐波声速及声衰减的不平等，精确地描述了各阶谐波在频散、衰减介质中的非线性传播，且计算量较小。

【技术实现步骤摘要】
一种非线性声场的基频和高次谐波频域有限差分计算方法
本专利技术涉及一种非线性声场的基频和高次谐波频域有限差分计算方法，利用傅立叶变换实现从时域向频域的转变，并能充分利用各阶谐波不等声速及声衰减的物理条件，真实的模拟具有强衰减、强频散介质中的非线性声传播。
技术介绍
自20世纪80年代开始，非线性效应的影响越来越受到人们的重视，尤其近二十多年，非线性声学已经渗透到各种场合，应用也越来越广，诸如物理学领域、生物医学领域、水声学领域等。由于实际的传播介质大都不是均匀介质，如水声领域中的含气泡水介质，医学中使用得造影剂等，均具有强衰减、强频散的特性，即声波的相速度、衰减系数随声波的频率变化而变化，这就使得一列声波在这种介质传播时，由于非线性效应产生的谐波相速度、衰减系数是不相等的。而典型的非线性波动方程，如Westervelt、KZK等，均是将不同频率声波的相速度和衰减系数等看做相等的加以处理，这与实际问题存在较大偏差。如何对不等声速和声衰减条件下的非线性声场进行模拟和计算成为亟需解决的问题。目前，主流的数值计算方法主要有以下两种：时域有限差分、频域有本文档来自技高网...

【技术保护点】
1.一种非线性声场的基频和高次谐波频域有限差分计算方法，其特征在于：步骤如下：/n(a)以空间z方向为主要声传播方向，对Westervelt方程进行时间维度以及空间x、y维度三维傅立叶变换，得到频域上的非线性波动方程；/n(b)利用一维格林函数，得到频域非线性波动方程的通解；/n(c)考虑弱非线性近似，基频采用线性近似表达式，求解得到二次谐波与基频的对应关系；/n(d)利用得到的二次谐波与基频的对应关系，替代掉基频表达式中卷积项的二次谐波，对基频表达式进行修正；/n(e)采用黎曼和近似，对二次谐波通解积分项写成便于数值求解的形式；/n(f)基于声波波长确定空间步长，将基频总解代入到黎曼和，计算...

【技术特征摘要】
1.一种非线性声场的基频和高次谐波频域有限差分计算方法，其特征在于：步骤如下：
(a)以空间z方向为主要声传播方向，对Westervelt方程进行时间维度以及空间x、y维度三维傅立叶变换，得到频域上的非线性波动方程；
(b)利用一维格林函数，得到频域非线性波动方程的通解；
(c)考虑弱非线性近似，基频采用线性近似表达式，求解得到二次谐波与基频的对应关系；
(d)利用得到的二次谐波与基频的对应关系，替代掉基频表达式中卷积项的二次谐波，对基频表达式进行修正；
(e)采用黎曼和近似，对二次谐波通解积分项写成便于数值求解的形式；
(f)基于声波波长确定空间步长，将基频总解代入到黎曼和，计算出黎曼和；
(g)将黎曼和代入到二次谐波通解中，求解出二次谐波在空间中的声场。


2.根据权利要求1所述的一种非线性声场的基频和高次谐波频域有限差分计算方法，其特征在于：(a)中的频域非线性波动方程为：



式中：w是角频率，ρ是介质的密度，β是介质的非线性系数，p(r,t)是时域声压，P(kx,ky,z,w)是时域声压经三维傅立叶变换后的频域复声压，卷积项与kx,ky,w相关，kx,...

【专利技术属性】
技术研发人员：时洁，邓安定，时胜国，张昊阳，刘宇林，田鸿兴，李松，李洪道，孙闻伯，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙;23