条件设置
主要公式
Blasius: $\delta = 5x/\sqrt{Re_x}$
湍流(1/7次): $\delta = 0.37x/Re_x^{0.2}$
$C_{f,lam}= 0.664/\sqrt{Re_x}$
$C_{f,turb}= 0.0592/Re_x^{0.2}$
// Boundary layer flow animation
(function() {
const el = document.getElementById('blAnimCanvas');
const ctx = el.getContext('2d');
const NPART = 60;
const particles = [];
let animT = 0;
function getParams() {
const Uinf = parseFloat(document.getElementById('uinf').value) || 5;
const L = parseFloat(document.getElementById('length').value) || 1.0;
const fs = document.getElementById('fluidSel').value;
const nuMap = { water: 1.004e-6, air: 1.516e-5, oil: 4.6e-5, custom: 15e-6 };
const nu = nuMap[fs] || 1.516e-5;
return { Uinf, L, nu };
}
function deltaLam(x, Uinf, nu) {
return x > 0 ? 5 * x / Math.sqrt(Uinf * x / nu) : 0;
}
function initParticles(W, H) {
particles.length = 0;
for (let i = 0; i < NPART; i++) {
particles.push({
x: Math.random() * W,
y: Math.random() * H * 0.95,
speed: 0
});
}
}
function resize() {
const dpr = window.devicePixelRatio || 1;
const w = el.parentElement.clientWidth - 36;
const H = Math.round(w * 0.38);
if (Math.abs(el.width - w * dpr) > 2 || el.height !== H * dpr) {
el.width = w * dpr;
el.height = H * dpr;
el.style.height = H + 'px';
ctx.setTransform(dpr, 0, 0, dpr, 0, 0);
initParticles(w, H);
}
return { W: w, H };
}
function frame() {
const { W, H } = resize();
const { Uinf, L, nu } = getParams();
const padL = 30, padR = 20;
const drawW = W - padL - padR;
const wallY = H - 22;
const showTurb = document.getElementById('showTurb').checked;
ctx.clearRect(0, 0, W, H);
// Background gradient (fluid domain)
const bgGrad = ctx.createLinearGradient(0, 0, 0, wallY);
bgGrad.addColorStop(0, '#001F3F');
bgGrad.addColorStop(1, '#002d5a');
ctx.fillStyle = bgGrad;
ctx.fillRect(padL, 0, drawW, wallY);
// Wall
ctx.fillStyle = '#003875';
ctx.fillRect(padL, wallY, drawW, H - wallY);
ctx.fillStyle = '#007BFF';
ctx.fillRect(padL, wallY, drawW, 2);
// Wall label
ctx.fillStyle = 'rgba(255,255,255,0.5)';
ctx.font = '10px Roboto Mono, monospace';
ctx.fillText('平板', padL + 4, wallY + 14);
// Boundary layer profile line
const N = 80;
ctx.beginPath();
ctx.setLineDash([5, 3]);
ctx.strokeStyle = '#00B4D8';
ctx.lineWidth = 1.5;
for (let i = 0; i <= N; i++) {
const xFrac = i / N;
const xPhys = xFrac * L;
const px = padL + xFrac * drawW;
let dPx;
if (showTurb && xPhys > 0 && Uinf * xPhys / nu > 5e5) {
const delta = 0.37 * xPhys / Math.pow(Uinf * xPhys / nu, 0.2);
dPx = Math.min(delta / L * drawW * 2.5, wallY * 0.85);
} else {
const delta = deltaLam(xPhys, Uinf, nu);
dPx = Math.min(delta / L * drawW * 2.5, wallY * 0.85);
}
const py = wallY - dPx;
if (i === 0) ctx.moveTo(px, py);
else ctx.lineTo(px, py);
}
ctx.stroke();
ctx.setLineDash([]);
// Label for boundary layer
ctx.fillStyle = '#00B4D8';
ctx.font = '11px Noto Sans JP, sans-serif';
ctx.fillText('δ(x) 境界層', padL + 8, 18);
// Free stream arrow
ctx.strokeStyle = 'rgba(255,255,255,0.6)';
ctx.lineWidth = 1;
ctx.beginPath();
ctx.moveTo(padL + 8, H * 0.15);
ctx.lineTo(padL + 55, H * 0.15);
ctx.stroke();
ctx.beginPath();
ctx.moveTo(padL + 55, H * 0.15);
ctx.lineTo(padL + 48, H * 0.15 - 4);
ctx.lineTo(padL + 48, H * 0.15 + 4);
ctx.closePath();
ctx.fillStyle = 'rgba(255,255,255,0.6)';
ctx.fill();
ctx.fillStyle = 'rgba(255,255,255,0.6)';
ctx.font = '10px Roboto Mono, monospace';
ctx.fillText('U∞', padL + 58, H * 0.15 + 4);
// Animate particles
animT += 1;
for (let p of particles) {
const xFrac = (p.x - padL) / drawW;
const xPhys = Math.max(0.001, xFrac * L);
let delta;
if (showTurb && Uinf * xPhys / nu > 5e5) {
delta = 0.37 * xPhys / Math.pow(Uinf * xPhys / nu, 0.2);
} else {
delta = deltaLam(xPhys, Uinf, nu);
}
const deltaScaled = Math.min(delta / L * drawW * 2.5, wallY * 0.85);
const distFromWall = wallY - p.y;
// 速度 profile: Blasius-like, 1 in free stream, ~0 near wall
let uFrac;
if (distFromWall <= 0) {
uFrac = 0;
} else if (deltaScaled < 1 || distFromWall > deltaScaled * 1.5) {
uFrac = 1.0;
} else {
const eta = distFromWall / deltaScaled;
// Approximate Blasius: u/U ≈ tanh(2.5*eta)
uFrac = Math.tanh(2.5 * eta);
}
const spd = uFrac * Uinf * 0.3;
p.x += spd;
// Fade color by layer: blue near wall, white in free stream
const alpha = 0.55 + uFrac * 0.35;
const r = Math.round(uFrac * 180);
const g = Math.round(100 + uFrac * 100);
const b = 255;
ctx.beginPath();
ctx.arc(p.x, p.y, 2, 0, Math.PI * 2);
ctx.fillStyle = `rgba(${r},${g},${b},${alpha})`;
ctx.fill();
// Reset particle when it exits right or is in the wall
if (p.x > padL + drawW || p.y >= wallY) {
p.x = padL + Math.random() * 20;
p.y = Math.random() * (wallY - 4);
}
}
requestAnimationFrame(frame);
}
frame();
})();
什么是平板边界层
🎓
简单来说，就是流体（比如空气或水）流过物体表面时，紧贴着物体表面那一层“被拖慢”的流体。比如你把手伸出快速行驶的汽车窗外，感觉风很大，但紧贴你手表面的空气其实几乎是静止的，这个从静止到高速的过渡层就是边界层。在这个模拟器里，你试着拖动“主流速度”的滑块，就能看到边界层厚度如何随着流速变化。
🙋
诶，真的吗？那为什么有时候叫“层流”，有时候又叫“湍流”呢？
🎓
这取决于流体微团的运动是否规则。层流时，流体像一层层光滑的薄片平行滑动；湍流时，流体微团会剧烈地横向掺混。在实际工程中，比如飞机机翼，前段通常是层流以减少阻力，后段可能变成湍流。你可以在模拟器里打开“显示转捩点”开关，它会标记出从层流变为湍流的位置（通常是雷诺数 $Re_x = 5 \times 10^5$ 的地方），然后改变平板长度L，看看这个点怎么移动。
🙋
原来是这样！那这个厚度和摩擦力对工程师有什么用呢？
🎓
太有用了！边界层越厚，对物体形状的“掩盖”越明显，会影响后续的流动。而摩擦力直接决定了阻力大小。比如在设计船舶外壳时，工程师需要精确计算摩擦阻力来预估燃油消耗。你可以在模拟器里把流体从“空气”换成“油”，它的运动粘度ν更大，然后观察图表中摩擦系数 $C_f$ 曲线的变化，就能直观感受到粘度对阻力的巨大影响。
物理模型与关键公式
边界层厚度的计算依赖于流动状态（层流或湍流），核心是局部雷诺数 $Re_x = \frac{U_\infty x}{\nu}$，其中 $U_\infty$ 是来流速度，$x$ 是距平板前缘的距离，$\nu$ 是流体的运动粘度。
$$ \text{层流 (Blasius解)}: \quad \delta(x) = \frac{5x}{\sqrt{Re_x}}$$
$\delta(x)$：距离前缘 $x$ 处的边界层厚度。层流边界层厚度与 $\sqrt{x}$ 成正比，增长较慢。
当流动转变为湍流后，边界层厚度增长会显著加快，通常采用基于实验的1/7次方速度剖面律进行估算。
$$ \text{湍流 (1/7次方律)}: \quad \delta(x) = \frac{0.37x}{Re_x^{0.2}}$$
$\delta(x)$：湍流边界层厚度。湍流边界层厚度与 $x^{0.8}$ 成正比，增长远快于层流。
工程师更关心的是表面摩擦带来的阻力，这由局部摩擦系数 $C_f$ 和平均摩擦系数 $C_{f,avg}$ 描述。
$$ \text{局部摩擦系数}: \quad C_f = \begin{cases}\dfrac{0.664}{\sqrt{Re_x}}, & \text{层流}\\[10pt]
\dfrac{0.0592}{Re_x^{0.2}}, & \text{湍流}\end{cases}$$
$C_f$：局部摩擦系数，用于计算单位面积上的摩擦剪应力 $\tau_w = C_f \cdot \frac{1}{2}\rho U_\infty^2$。
现实世界中的应用
航空航天：飞机机翼和机身表面的边界层控制至关重要。通过计算边界层增长和转捩位置，可以优化机翼外形，设计翼梢小翼等装置，甚至采用层流翼型来延迟转捩，从而显著降低飞行阻力，节省燃油。
船舶与潜艇设计：船体在水下部分受到巨大的摩擦阻力。工程师使用此类计算来预估不同航速下的阻力，从而优化船体线型、评估不同防污涂层的效果（涂层会影响表面粗糙度，进而影响边界层发展），是船舶动力性能预报的核心环节。
能源与动力工程：在燃气轮机叶片、风力涡轮机叶片以及各种换热器管道的内流道中，边界层的发展直接影响气动效率、热交换效率和流动损失。准确预测边界层是进行高效叶型设计和流道优化的基础。
汽车工业：在汽车外流场分析中，车身表面（尤其是引擎盖、车顶、行李箱盖）的边界层发展会影响气动升力、阻力以及后部涡流的形成。这是进行CFD风洞仿真，降低风噪和油耗的关键输入之一。
常见误解与注意事项
模型假设：本模拟器所用数学模型基于线性、均质、各向同性等简化假设。在将计算结果直接用于设计决策之前，务必确认实际系统是否满足这些假设。
单位与量纲：许多计算错误源于单位换算错误或数量级判断失误。请时刻注意各参数输入框旁标注的单位。
结果验证：始终将模拟器输出结果与物理直觉或手算结果进行核对。若结果出乎意料，请检查输入参数或采用独立方法进行验证。
进阶学习指引
深化理论：在本工具的简化模型基础上，进一步研究非线性效应、三维行为和时间依赖现象。阅读专业教材和学术论文，掌握严格的数学推导，是提升工程解题能力的关键。
数值方法：系统学习有限元法（FEM）、有限差分法（FDM）和有限体积法（FVM），理解商业CAE求解器的内部运行机制，这将显著提升您设置有效仿真的能力。
实验验证：理论和仿真结果必须通过实验数据加以验证。养成将计算结果与测量值进行对比的习惯，这正是V&V（验证与确认）的精髓所在。
CAE工具：准备好后，可进一步探索Ansys、Abaqus、OpenFOAM、COMSOL等业界主流工具。通过本模拟器培养的物理直觉，将帮助您更有效地配置和使用这些工具。