《水动力噪声理论》一书中Lighthill声类比方程与FW-H方程的联系与区别
时间: 2023-06-19 13:08:52 浏览: 101
Lighthill声类比方程和Ffowcs Williams-Hawkings (FW-H)方程都是用来描述流体中声波传播的方程,但它们的物理假设和适用范围略有不同。
Lighthill声类比方程是由Lighthill于1952年提出的,它基于Navier-Stokes方程、连续性方程和能量方程,将流体力学和声学联系在了一起。该方程假设声波的传播速度远大于流体中的平均速度,所以可以将流体中的扰动分为快速变化的声波和慢速变化的湍流扰动。在此基础上,Lighthill声类比方程将声波扰动和湍流扰动联系在了一起,从而描述了声波随时间和位置的演变。
FW-H方程是由Ffowcs Williams和Hawkings于1969年提出的,它基于Euler方程和Kirchhoff积分定理,假设声源是小振幅、高频率的点源,并且远离声源时流体中的声波可以看作是平面波。在此基础上,FW-H方程描述了声波在流体中的传播方式,能够计算声源产生的声功率和声辐射方向。
因此,Lighthill声类比方程和FW-H方程的适用范围略有不同。Lighthill声类比方程适用于低频声波和湍流扰动相互作用的情况,而FW-H方程适用于高频声波和点源附近的声波传播情况。不同的物理假设和适用范围决定了它们的计算方法和结果的精度。
相关问题
用matlab语言生成一个带有周期性边界条件的,初始时刻车流密度为0.25,迭代公式为z(j,t+2)=z(j,t+1)-r^2*V*(p*(z(j+1,t)-z(j,t))+(1-p)*(z(j+2,t)-z(j+1,t)))+lambda2*r^2*V*(z(j+1,t+1)-2*(z(j,t+1)+z(j-1,t+1)));的时空演化图代码
好的,下面是一个 MATLAB 代码,可以生成一个带有周期性边界条件的时空演化图,初始时刻车流密度为 0.25,迭代公式为:
$$
z_{j,t+2} = z_{j,t+1} - r^2 V \left[ p (z_{j+1,t} - z_{j,t}) + (1-p) (z_{j+2,t} - z_{j+1,t}) \right] + \lambda_2 r^2 V (z_{j+1,t+1} - 2(z_{j,t+1}+z_{j-1,t+1}))
$$
其中 $z_{j,t}$ 表示时刻 $t$ 时第 $j$ 个车道上的车辆密度,$r$ 表示空间步长,$V$ 表示车辆的平均速度,$p$ 表示车辆前进时会看到前一个车辆的概率,$\lambda_2$ 表示车辆的加速度,时间步长为 $1$。
```matlab
% 设置参数
L = 100; % 道路长度
r = 1; % 空间步长
V = 1; % 车辆平均速度
p = 0.5; % 前车占有率
lambda2 = 1; % 加速度
p0 = 0.25; % 初始时刻车流密度
time = 100; % 模拟时长
dt = 1; % 时间步长
% 初始化车流密度
z = zeros(L, time+2);
for i = 1:L
if rand() < p0
z(i, 1) = 1;
end
end
% 开始模拟
for t = 1:time
for j = 1:L
% 计算迭代公式
z(j, t+2) = z(j, t+1) - r^2 * V * (p * (z(mod(j,L)+1, t) - z(j,t)) + (1-p) * (z(mod(j+1,L)+1, t) - z(mod(j,L)+1, t))) + lambda2 * r^2 * V * (z(mod(j,L)+1, t+1) - 2 * (z(j,t+1) + z(mod(j-1,L)+1, t+1)));
end
% 绘图
plot(1:L, repmat(t+2, L, 1), 'w', 'LineWidth', 3);
hold on;
plot(1:L, repmat(t+2, L, 1), 'k--', 'LineWidth', 1);
plot(1:L, repmat(t+1, L, 1), 'w', 'LineWidth', 3);
hold on;
plot(1:L, repmat(t+1, L, 1), 'k--', 'LineWidth', 1);
plot(1:L, repmat(t, L, 1), 'w', 'LineWidth', 3);
hold on;
plot(1:L, repmat(t, L, 1), 'k--', 'LineWidth', 1);
plot(1:L, z(:, t), 'bo', 'MarkerSize', 5, 'MarkerFaceColor', 'b');
axis([0 L 0 time+3]);
xlabel('Position');
ylabel('Time');
drawnow;
end
```
这个代码使用了 Lighthill-Whitham-Richards (LWR) 模型来模拟车辆的运动,在 LWR 模型的基础上加入了加速度项。在每个时间步长中,程序首先计算车辆间距,并根据间距和迭代公式来更新车辆密度。然后,程序绘制时空演化图。
需要注意的是,由于这个模型是比较复杂的,因此运行程序可能需要一定的时间。同时,这个模型也比较抽象,可能无法完全反映真实的交通流动情况,仅供参考。
matlab交通流量模型
在 Matlab 中可以使用以下交通流量模型:
1. LWR 模型:Lighthill-Whitham-Richards 模型是一个宏观交通流模型,其基本假设是交通流是连续的。
2. CTM 模型:Cell Transmission Model 是一个微观交通流模型,它将路段分成若干个小单元,计算每个单元内的车辆数量和速度。
3. PDM 模型:Partial Differential Equation Model 是一个基于偏微分方程的交通流模型,它使用偏微分方程来描述车流量和车速的变化。
4. GHR 模型:Greenberg-Hastings-Rose 模型是一个基于离散动力学的交通流模型,它模拟了车流量和车速的非线性变化。
这些模型在 Matlab 中都有相应的实现方法和工具箱,你可以根据自己的需求选择合适的模型进行交通流量模拟。