介绍重要的建模概念并展示相关应用案例

应用案例 2018-12-29 10:59:14

  physics® 软件及其附加的“声学模块”提供了线性纳维-斯托克斯接口,支持详细模拟对流动和声学之间复杂的相互作用。现在,当系统的声学属性可以由湍流背景流场改变或决定时,您可以对系统进行稳健的仿真;例如汽车的排气系统。在本文中,我们将介绍重要的建模概念,并展示相关应用案例。

  稳态背景流场和声场之间复杂的相互作用可以使用“声学模块”中的线性纳维-斯托克斯物理场接口来模拟。此接口支持详细分析流体流动——可以同时是湍流和非等温流——是如何影响不同系统中的声场的。这包括当背景流场与声场发生相互作用时,以及当流场改变声场时发生的所性效应。线性纳维-斯托克斯接口不包含流致噪声源项。这些方程基本上求解的是一般形式 CFD 方程的全线性扰动——质量、动量和能量守恒。

  对于许多行业和应用领域,模拟与仿真背景流场对声场的具体影响具有重要意义。在汽车工业中,流经的流体会改变排气和进气系统的声学属性,例如,旁背景流场的大小会影响消声器的传输损耗。在航空航天应用中,衬垫和穿孔板在系统引入流动时的声学性能是一个研究重点。子系统的具体声学属性(吸收、和反射系数)可以影响整个系统的性能,喷气发动机便是如此。

  在消声器和衬垫的示例中,线性纳维-斯托克斯方程也可以捕获背景流场中的湍流所导致的声信号衰减。此外,这些模型中的背景流场通常属于非等温流动。

  汽车应用的示例。上图显示了基于亥姆霍兹共振器的流场示例的结果。前侧的彩色表面图显示了声压级。后侧流线图显示了背景流场。

  线性纳维-斯托克斯接口提供了一个与结构相互耦合的内置多物理场,因此我们能够在频域(或者线性化的时域中)中现成地设置流-固耦合(FSI)模型。在许多应用中,流动、声学和结构振动的相互作用都是重要的考虑因素。一个应用案例是科里奥利流量计的流量感测功能。总而言之,这些接口适用于分析结构在背景流场的流体载荷作用下的振动特性变化。

  频域中的流-固耦合示例:在基频下驱动的科里奥利流量计的运动情况。表面显示了结构变形(为了加强可视化效果,我们特意夸大了相位和振幅),管道敞开的切口段显示了管道的内表面声压。

  线性纳维-斯托克斯接口还可用于研究燃烧不稳定性和一般的管道内声学,也可以研究更多学术类应用,例如分析流动不稳定性的起始点,或者研究易于发生啸叫的区域。

  现在,此接口增加了伽辽金最小二乘(Galerkin least squares,简称 GLS)稳定方案,助力提升仿真的稳健性。新增的默认设置能够更好地处理由控制方程的对流项和反应项引入的数值与物理不稳定性。此外,当使用迭代求解器求解模型时,非常适合使用重新推导后的滑移边界条件。对于必须求解的大型工业问题,这一点至关重要。

  线性纳维-斯托克斯方程对描述可压缩、粘性与非等温流体的整套控制方程(线性纳维-斯托克斯方程)的线性化进行了表征。此方程相当于由压力、速度、温度和密度(p0、u0、T0 和 ρ0)定义的稳态背景流场的一阶扰动方程。由此可推导出描述压力、速度和温度(p、u、和 T)——即因变量——的微小扰动的控制方程。在扰动理论中,下标 1 有时表示变量为一阶扰动项。控制方程(下标 0 的量表示背景场)写作:

  其中 Φ = ∇u : τ0 + u0 : τ 是粘性耗散函数;M、F 和 Q 代表可能的源项;κ 是传热系数(国际单位:W/m/K);αp 是(等压)热膨胀系数(国际单位:1/K);βT 是等温压缩率(国际单位:1/Pa);p 是恒定压力下的比热容(单位质量的热容)(国际单位:J/kg/K)。

  在频域内,iω 乘子表示时间导数。应力张量和线性化状态方程(密度扰动)的本构方程由下列公式给出:

  其中 τ是粘性应力张量(斯托克斯表达式),μ 是动力粘度(国际单位:Pa s),μB 是体积粘度(国际单位:Pa s)。

  我们将傅立叶导热定律应用到能量方程中。您可以在 Acoustics Module User’s Guide 中查阅方程的详细推导过程。然后可以利用线性纳维-斯托克斯,瞬态 接口或线性纳维-斯托克斯,频域 接口在时域或频域中求解方程。

  由于在接口中求解的是通用方程,这些方程默认模拟声学(可压缩)波、涡旋波和熵波的。后两种类型的波依赖于背景流场的速度进行对流,而不以声速。声波在过程中可以(通过反应项)与流体相互作用,声能量可以传递给声学模式,并从声学模式传递给涡旋模式和熵模式。控制方程中的反应项是引起类似的流动声学耦合的原因。这是因为涡旋波和熵波对背景流场的解产生了非声学(类似于 CFD)扰动,因此在某种程度上,反应项模拟的是 CFD 和声学之间的线性相互作用。

  许多气动声学公式忽略了反应项,因为它们也是 开尔文-亥姆霍兹不稳定性 的产生过程的背后原因。这些过程很难进行数值处理。另一方面,如果忽略这些项,则不能对声音衰减和放大进行精确建模。反应项完全包含在线性纳维-斯托克斯接口中。

  不稳定性增长在 COMSOL Multiphysics 中有两种处理方式。我们可以通过选择频域而非时域公式来解决随时间增长的不稳定性。如果不能正确地对涡旋模式进行解析,可能出现空间不稳定性,这时可以利用伽辽金最小二乘法稳定方案有效地进行处理。

  根据利用线性纳维-斯托克斯方程模拟的不同应用,我们可能需要解析声学、粘性和热边界层。如果存在无滑移和等温边界条件,则需要在固体表面上创建上述边界条件,从而对振荡流进行解析。通常情况下,在大型模型(同边界层厚度相比)中,没有必要考虑边界层的损耗。在液体中,我们通常也可以忽略热边界层,但是气体中一定要添加。通过在壁边界条件下勾选滑移或绝热选项,可以忽略这两种效应。

  值得一提的是,我们还能在背景流场和声场之间创建另一个间接耦合。当声波在涉及湍流背景流场的区域中时会衰减。只要将 CFD RANS 模型的湍流粘度耦合到声学模型中,即可将衰减效应引入模型。比如说,在分析存在流动现象的消声器系统的传输损耗时,衰减效应是一个重要的考虑因素。

  求解属于计算气动声学(computational aeroacoustics,简称 CAA)领域的线性纳维-斯托克斯方程时,我们需要仔细考虑、理解与处理数值挑战。如上所述,控制方程在物理(开尔文-亥姆霍兹)和数值方面具有不稳定性。由于接口应用了稳定性,那么剩下的关键数值难题就是避免在包含背景场变量(p0、u0、T0 和 ρ0)的项中引入数值噪音。如果反应项的变量存在梯度,尤其需要注意这一点。

  如果 CFD 和声学模型采用不同的网格,且/或背景流场和声学问题采用不同的离散化阶次,此问题发生的可能性更大。请注意,我们之所以使用不同的网格或离散化阶次,主要是因为这两个问题需要求解不同的物理场和长度尺度。为了防止此类情况,我们需要谨慎地将背景流场的数据从 CFD 映射到声学模型。这是计算气动声学建模中一个易于理解和描述的步骤。另外,映射步骤可用于平滑 CFD 数据,可以是整体平滑,也可以是特定细节——比如流体动力学边界层——的局部平滑,如果细节对于声学模型不重要的话。

  在 COMSOL Multiphysics 中,网格之间的映射由另外的研究步骤来完成。Acoustics Module User’s Guide 和线性纳维-斯托克斯物理场接口的应用教学模型描述了此步骤的详细信息。

  解析声学边界层:根据所模拟的物理效应和模型大小,判断是否需要对声学边界层进行解析。如果不需要,则将壁上默认的无滑移和等温条件改设为滑移和绝热条件。背景流场的解析细节也会影响这一选择。例如,如果背景流场边界层要实现全解析度,通常需要在声学问题中设置匹配的无滑移条件。

  网格应该解析 CFD 和声场:重要的几何特征、边界层和大梯度区域应该使用CFD 和声学仿真中的网格来解析。具体而言,声学仿真的网格(如果与 CFD 网格不同的话)应该解析波长和声学边界层(建模请参考上一条)等声学特征,以及背景流场特征。

  映射:使用映射步骤将 CFD 数据映射到声学问题,尤其当使用不同的网格或离散化阶次时。必要时,可以平滑解,(若声学研究应用了滑移条件)也可以平滑边界层。根据需要,可将背景流场的无滑移条件添加到映射中。

  离散化阶次:默认情况下,线性纳维-斯托克斯接口对因变量全部采用线性离散化,对于大多数模型这是一个合理的选择。但是,如果应用了无滑移和等温条件,则最好将速度和温度变量切换为二阶离散化。这可以增加壁附近的空间解析度,但是也引入了更多待求解的度。

  (排气系统中的)亥姆霍兹共振器会使特定的窄频带衰减。系统中的流体流动会改变共振器的声学属性和子系统的传输损耗。亥姆霍兹共振器教学模型研究了主管道中引入流动时的传输损耗(共振器位于主管道的侧分支)。

  计算平均流动时,采用马赫数为 Ma = 0.05 和 Ma = 0.1 的 SST 湍流模型。然后使用线性纳维-斯托克斯,频域 接口求解声学问题。接下来,将声学模型与平均流速、压力和湍流粘度耦合。损耗的仿真预测与期刊论文的发布数据高度吻合(Ref. 1)。为了准确找到共振,并传输损耗大小的正确性,模型必须适当地平衡对流项和扩散项。平衡是在模型中实现的。

  系统在频率为 100 Hz,马赫数为 Ma = 0.1 时的内部压力分布。平面波从流体上游左侧入射。

  在涉及切向背景流场的声学衬垫教学模型中,声学衬垫由八个带微缝的共振器组成,背景切向流场的马赫数为 0.3。衬垫上方的声压级计算结果与研究论文(Ref. 2)发表的数据高度一致。该示例使用“CFD 模块”的 SST 湍流模型来计算流动,并使用线性纳维-斯托克斯,频域 接口计算声。然后对声学边界层进行解析,并将默认的线性离散化选项修改为二阶离散化,从而改进壁附近的空间解析度。

  曲线显示了在四个不同的驱动频率下,衬垫上方的表面声压级。曲线的彩色部分突出显示了衬垫的范围。仿真结果与参考研究论文的实验结果呈现出高度一致。

  系统在频率为 100 Hz,马赫数为 Ma = 0.1 时的内部压力分布。平面波从流体上游左侧入射。

  在涉及切向背景流场的声学衬垫教学模型中,声学衬垫由八个带微缝的共振器组成,背景切向流场的马赫数为 0.3。衬垫上方的声压级计算结果与研究论文(Ref. 2)发表的数据高度一致。该示例使用“CFD 模块”的 SST 湍流模型来计算流动,并使用线性纳维-斯托克斯,频域 接口计算声。然后对声学边界层进行解析,并将默认的线性离散化选项修改为二阶离散化,从而改进壁附近的空间解析度。

  曲线显示了在四个不同的驱动频率下,衬垫上方的表面声压级。曲线的彩色部分突出显示了衬垫的范围。仿真结果与参考研究论文的实验结果呈现出高度一致。

  声速波动在衬垫上方作为平面波进行,动画显示前四个衬垫。驱动频率为 1000 Hz。彩色图显示速度大小,箭头显示速度矢量。在衬垫表面的小孔附近,流体与声学的相互作用产生了涡流。

  科里奥利流量计——又称质量流量计或惯性流量计——可以测量流经此流量计的流体质量流率。该装置还可以计算流体密度以及基于密度的体积流率。科里奥利流量计教学模型演示了如何利用弯曲的几何对通用的科里奥利流量计进行建模。

  当流体流经弹性结构(例如弯管)时,它会与弹性结构的震荡运动相互作用。科里奥利效应导致管道上变形的两点之间产生相位差,可用于计算质量流率。

  为了对此进行建模,我们借助内置的多物理场耦合,将线性纳维-斯托克斯,频域 接口耦合到固体力学 接口。然后使用湍流,SST 接口来模拟背景平均流动。通过这种方法,我们可以在频域中有效地模拟流-固耦合。

  上游点和下游点之间的相差(下面动画中的红点)。图中曲线表示运行科里奥利流量计所需的校准结果。

  三个不同的质量流率对应的科里奥利流量计的运动。流量计在结构的固有频率 fd = 163.5Hz 下被驱动。为了增强可视化效果,绘图夸大了变形幅度和相位。随着流率增加,上游和下游的相差随之增大。

  文章出处:【微信号:COMSOL-China,微信号:COMSOL】欢迎添加关注!文章转载请注明出处。

  在执行声学分析之前,研究人员将发动机外壳(发动机的外部表面)封闭在被笛卡尔完美匹配层(perfect....

  本文建立了基于svpwm变频调速的双电机偏差耦合控制算法的仿真模型,并用Matlab7.1仿真软件进....

  图1所示是双管阻容耦合放大器。这一多级放大器由两个单级放大器组成,两级放大器之间通过电容耦合,所以称....

  作者:Robert Kollman,仪器 (TI) 在这篇《电源设计小贴士》中,我们继续《电源设计小贴士#32-第 1 部分》...

  电子工程师指从事各类电子设备和信息系统研究、教学、产品设计、科技开发、生产和管理等工作的高级工程技术人才。一般分为硬...

  据瑞声科技上半年财报显示:上半年实现营收84.24亿元,同比下跌2.5%;净利润17.78亿元,同比....

  本文使用通用电仿真工具ANSYS有限元软件,对上述这两种耦合机构进行仿真。设定耦合机构中发射线圈....

  首先谈谈我们是如何测量ADC性能的吧,ADC芯片是AD9218,模拟信号是单端信号,所以参与AD831...

  描述 介绍了直流耦合和交流耦合应用中的宽带单端至差动转换参考设计。该设计评估 LMH5401 和 LMH6401 级联的性能并提供...

  856880是用于频带2的通用上行链声表面(SAW)滤波器。该滤波器是专门设计用于基站应用的3x3....

  看到ADN2807 DATASHEET第18页上讲到AC耦合模式下,可以1000个CID,即连续0或连续1,想问一下这是什么机理呢...

  由整流和逆变电组成的UPS变换电是UPS中的关键组成部分。根据其所采用的功率电子器件和UPS不同....

  要提高开关频率,提高开关电源产品的质量,电磁兼容性问题必须重点考虑。本文是在分析了干扰产生机理以及经....

  当来自两个不同电的电流流经一个公共时就会产生共耦合。上的压降由两个电决定,来自两个电....

  客户在不了解两个器件处能的这一参数时,就会向技术支持人员提出疑问。如果采用交流耦合设计,并且....

  串扰是信号完整性中最基本的现象之一,在板上走线密度很高时串扰的影响尤其严重。我们知道,线性无缘系统满....

  一直以来,设计中的电磁干扰(EMI)问题十分令人头疼,尤其是在汽车领域。为了尽可能的减小电磁干扰,设....

  使用子面板(SubPanel)技术,我们可以很方便地做出模块化的程序,并且模块之间不会产生耦合;对于....

  本文通过在地板上放置两个天线单元结构,单元之间用金属中和线进行连接,并在中和线上加载集总电感元件,对....

  本文可以让你掌握音频声学的基础理论和电\磁\机械学中与喇叭有关的基本知识,了解扬声器测试的要求和 T....

  该应用报告讨论了驱动高速ADCs与德克萨斯仪器LMH6521 DVGA高IFAC耦合应用。 采样数据....

  现如今有期望可以通过光子-磁子的强耦合体系得以很好地实现,并有希望在室温下的毫米尺寸的器件上得到突破....

  由于地线就是信号的回流线,因此当两个电共用一段地线时,彼此也会相互影响。一个电的地电位会受到另一....

  多回耦合输电方式是电网建设的必然趋势,以局部耦合多回线为主。受线间局部耦合影响,现有故障测距方法存....

  波暗室又叫无反射室、吸波暗室简称暗室。微波暗室由电磁屏蔽室、滤波与隔离、接地装置、通风波导、室内配电....

  一、电容 所谓电容,就是容纳和电荷的电子元器件。电容的基本工作原理就是充电放电,通交流,隔直流。....

  耦合指信号由第一级向第二级传递的过程,一般不加注明时往往是指交流耦合。退耦是指对电源采取进一步的滤波....

  为了方便优化高速干式切削工艺,对高速干式滚齿机切削进行热力耦合和工艺参数影响分析。首先结合高速干式滚....

  序言 70 年代中期发展起来的离子膜法烧碱工艺较传统的其他隔膜法相比具有能耗低,较环保,制得的产品...

  液压滑阀在工作过程中常常因黏性加热而出现阀芯热卡紧现象,基于流固耦合共轭传热方法,运用COMSOL软....

  以750 t履带起重机液压防后倾系统为例,基于多体动力学和虚拟仿真技术建立机液耦合模型,分析了在突然....

  前言 业界对高度、高性、低成本、易携带的生物传感器的研发兴趣有增无减。这些要求对于医疗、食品....

  针对TV模型去噪后图像容易产生阶梯效应的现象,提出一种全变分耦合图像去噪模型。首先,根据去噪过程中图....

  什么是耦合电容?什么是去耦电? 耦合指信号由第一级向第二级传递的过程,一般不加注明时往往是指....

  耦合是指两个或两个以上的电元件或电网络等的输入与输出之间存在紧密配合与相互影响,并通过相互作用从....

  耦合是指两个或两个以上的电元件或电网络等的输入与输出之间存在紧密配合与相互影响,并通过相互作用从一....

  集成电感对多相电压调节模块( VRM)稳态和动态特性有重要影响,合理的耦合度可以提高VRM的输出动态....

  针对增强图像中的弱边缘、细节纹理和消除二阶偏微分方程在图像平滑部分的阶梯效应问题,提出一种各向异性四....

  高速列车制动装置基本上采用盘形制动,在制动过程中制动盘最高温度和受到的最大应力关乎到制动装置工作的稳....

  基于高铁列车动态负载的动态耦合基础问题,提出一种发射端为单矩形线圈,接收端为多方形线圈级联的非对称耦....

  上海微系统所首次提出了微纳光纤耦合的SNSPD器件结构。该结构将SNSPD器件置于微纳光纤的倏逝场内....

  针对风电一氢储能与煤化工多能耦合系统( WP-HESCCMFCS)的经济性评估问题,提出了一种多能耦....

  影响单片机系统可靠安全运行的主要因素主要来自系统内部和外部的各种电气干扰,并受系统结构设计、元器件选....

  针对全球定位系统( GPS)中的多径误差消除问题,提出了一种基于完好性和深耦合结构相结合的多径误....

  超/特高压直流输电工程用接地极技术的应用日趋普遍,单极大地运行方式下直流接地极承载的入地电流大幅....

  消除公共耦合的途径有两个:一个是减小公共地线部分的,这样公共地线上的电压也随之减小,从而控制....

  零点漂移是怎样形成的:运算放大器均是采用直接耦合的方式,我们知道直接耦合式放大电的各级的Q点是相互....

  设计了一种新的耦合方法自变量斥耦合,并将其作用于斜帐篷映射来得到一个新的混沌系统。实验显示,这一系统....

  随着通讯技术不断发展,全波EM精确设计逐渐应用于很多滤波器设计中,但EM仿真的缺陷则在于仿真时间长,....

  示波器的输入耦合方式的意思是输入信号的传输方式。耦合是指两个或两个以上的电元件或电网络等的输入与输....

  受控源又称为非源。一般来说,一条支的电压或电流受本支以外的其它因素控制时统称为受控源。受控源....

  综合运用缝隙耦合馈电技术、双线馈线技术和引入空气层等方式展宽了天线的频带,设计并仿真出一种工作在Ku....

  耦合器的作用是把输入光信号和泵浦光耦合进掺铒光纤中,通过掺铒光纤作用把泵浦光的能量转移到输入光信号中....

  自从20世纪80年代初期第一片数字信号处理器芯片(DSP)问世以来,DSP就以数字器件特有的稳定性、....