// Pipe flow cross-section animation
(function() {
const el = document.getElementById('pipeAnimCanvas');
const ctx = el.getContext('2d');
const NPART = 80;
const particles = [];
let animT = 0;
function getParams() {
const D = (parseFloat(document.getElementById('diam').value) || 50) / 1000;
const mu_slider = parseFloat(document.getElementById('mu').value) || 1.0;
const rho_slider = parseFloat(document.getElementById('rho').value) || 998;
// Get velocity from stat card or slider
const useQ = document.getElementById('useQ').checked;
let U;
if (useQ) {
const qLog = parseFloat(document.getElementById('flow').value) || -2;
const Q = Math.pow(10, qLog);
const A = Math.PI * D * D / 4;
U = Q / A;
} else {
U = parseFloat(document.getElementById('velocity').value) || 1.0;
}
const mu = mu_slider * 1e-3;
const Re = rho_slider * U * D / mu;
return { D, U, Re };
}
function initParticles(W, H, pipeR) {
particles.length = 0;
const cx = W / 2;
const cy = H / 2;
for (let i = 0; i < NPART; i++) {
const angle = Math.random() * Math.PI * 2;
const r = Math.sqrt(Math.random()) * pipeR * 0.98;
particles.push({
x: Math.random() * W * 0.6 + W * 0.05,
y: cy + Math.cos(angle) * r,
r: r,
angle: angle
});
}
}
function resize() {
const dpr = window.devicePixelRatio || 1;
const w = el.parentElement.clientWidth - 36;
const H = Math.round(w * 0.42);
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);
const pipeR = Math.min(H * 0.4, w * 0.18);
initParticles(w, H, pipeR);
}
return { W: w, H };
}
function uProfile(r, R, Re) {
// Laminar: parabolic u/U_max = 1 - (r/R)^2
// Turbulent: 1/7 power u/U_max = (1 - r/R)^(1/7)
const rFrac = r / R;
if (rFrac >= 1) return 0;
if (Re < 3000) return Math.max(0, 1 - rFrac * rFrac);
else return Math.pow(Math.max(0, 1 - rFrac), 1/7);
}
function frame() {
const { W, H } = resize();
const { D, U, Re } = getParams();
const pipeR = Math.min(H * 0.4, W * 0.18);
const cx = W / 2;
const cy = H / 2;
const profileX = cx + pipeR + 40; // profile drawn to right of pipe
ctx.clearRect(0, 0, W, H);
// Dark background
ctx.fillStyle = '#001F3F';
ctx.fillRect(0, 0, W, H);
// Pipe walls
ctx.strokeStyle = '#007BFF';
ctx.lineWidth = 3;
ctx.beginPath();
ctx.moveTo(20, cy - pipeR);
ctx.lineTo(cx - 10, cy - pipeR);
ctx.stroke();
ctx.beginPath();
ctx.moveTo(20, cy + pipeR);
ctx.lineTo(cx - 10, cy + pipeR);
ctx.stroke();
// Pipe opening (ellipse)
ctx.strokeStyle = '#007BFF';
ctx.lineWidth = 2;
ctx.beginPath();
ctx.ellipse(cx - 10, cy, 10, pipeR, 0, 0, Math.PI * 2);
ctx.stroke();
// Fill pipe interior with gradient
const pipeGrad = ctx.createLinearGradient(20, cy - pipeR, 20, cy + pipeR);
pipeGrad.addColorStop(0, 'rgba(0,123,255,0.08)');
pipeGrad.addColorStop(0.5, 'rgba(0,123,255,0.18)');
pipeGrad.addColorStop(1, 'rgba(0,123,255,0.08)');
ctx.fillStyle = pipeGrad;
ctx.fillRect(20, cy - pipeR, cx - 30, pipeR * 2);
// 速度 profile arrows at exit
const nArrows = 12;
const uMaxDisplay = pipeR * 0.95;
ctx.fillStyle = '#003875';
ctx.fillRect(profileX - 15, cy - pipeR - 4, uMaxDisplay + 30, pipeR * 2 + 8);
ctx.strokeStyle = 'rgba(0,180,216,0.2)';
ctx.lineWidth = 1;
ctx.beginPath();
ctx.moveTo(profileX - 15, cy - pipeR);
ctx.lineTo(profileX - 15, cy + pipeR);
ctx.stroke();
for (let i = 0; i <= nArrows; i++) {
const rNorm = (i / nArrows) * 2 - 1; // -1 to 1
const r = Math.abs(rNorm) * pipeR;
const uFrac = uProfile(r, pipeR, Re);
const arrowLen = uFrac * uMaxDisplay;
const yPos = cy + rNorm * pipeR * 0.95;
const hue = Math.round(200 + uFrac * 40);
const lum = Math.round(40 + uFrac * 40);
ctx.strokeStyle = `hsl(${hue},90%,${lum}%)`;
ctx.lineWidth = 1.5;
ctx.beginPath();
ctx.moveTo(profileX - 15, yPos);
ctx.lineTo(profileX - 15 + arrowLen, yPos);
ctx.stroke();
if (arrowLen > 6) {
ctx.fillStyle = `hsl(${hue},90%,${lum}%)`;
ctx.beginPath();
ctx.moveTo(profileX - 15 + arrowLen, yPos);
ctx.lineTo(profileX - 15 + arrowLen - 5, yPos - 3);
ctx.lineTo(profileX - 15 + arrowLen - 5, yPos + 3);
ctx.closePath();
ctx.fill();
}
}
// Profile curve
ctx.beginPath();
ctx.strokeStyle = '#00B4D8';
ctx.lineWidth = 2;
for (let i = 0; i <= 40; i++) {
const rNorm = (i / 40) * 2 - 1;
const r = Math.abs(rNorm) * pipeR;
const uFrac = uProfile(r, pipeR, Re);
const arrowLen = uFrac * uMaxDisplay;
const yPos = cy + rNorm * pipeR * 0.95;
const x = profileX - 15 + arrowLen;
if (i === 0) ctx.moveTo(x, yPos);
else ctx.lineTo(x, yPos);
}
ctx.stroke();
// Labels
ctx.fillStyle = 'rgba(255,255,255,0.7)';
ctx.font = '11px Roboto Mono, monospace';
ctx.fillText('U_max', profileX - 15 + uMaxDisplay - 2, cy - pipeR - 8);
ctx.fillText(Re < 3000 ? '層流 (抛物线)' : '乱流 (1/7乗則)', profileX - 10, cy + pipeR + 16);
ctx.fillStyle = Re < 3000 ? '#2ecc71' : '#e17055';
ctx.font = 'bold 11px Noto Sans JP';
ctx.fillText(`Re = ${Re.toFixed(0)}`, 25, 18);
// Animate particles inside pipe
animT += 1;
for (let p of particles) {
const r = p.r;
const uFrac = uProfile(r, pipeR, Re);
const spd = uFrac * 2.5;
p.x += spd;
const alpha = 0.4 + uFrac * 0.5;
const g = Math.round(100 + uFrac * 120);
ctx.beginPath();
ctx.arc(p.x, p.y, 2, 0, Math.PI * 2);
ctx.fillStyle = `rgba(50,${g},255,${alpha})`;
ctx.fill();
if (p.x > cx - 12) {
p.x = 22 + Math.random() * 10;
const angle2 = Math.random() * Math.PI * 2;
p.r = Math.sqrt(Math.random()) * pipeR * 0.97;
p.y = cy + Math.cos(angle2) * p.r;
}
}
requestAnimationFrame(frame);
}
frame();
})();
什么是管道流动压力损失
🙋
水在管子里流得好好的，为什么还会“损失”压力呢？是什么东西在阻碍水流啊？
🎓
简单来说，压力损失就是流体为了克服管道内壁的摩擦和内部涡流而消耗的能量。你可以想象一下，水流过一条布满沙砾的河床，速度肯定会变慢，对吧？在实际工程中，比如给高楼供水，如果没算准这个损失，顶楼的水压就会不足。你试着在模拟器里把“粗糙度”滑块从光滑的“玻璃管”拖到粗糙的“混凝土管”，看看压力损失ΔP会怎么变化，马上就能感受到管壁粗糙的影响了。
🙋
诶，真的吗？我试了一下，粗糙度变大，压力损失果然飙升！那除了粗糙度，流量大小也会有影响吗？
🎓
当然有，而且影响很大！这里的关键是流态。流量小的时候，水流是分层平滑流动的“层流”，阻力主要来自流体内部的粘性。流量大到一定程度，就会变成混乱的“湍流”，阻力会急剧增加。这个转折点就是“雷诺数”。你可以在模拟器里把“流量Q”从很小慢慢调大，观察左下角“摩擦系数f”的变化曲线。在某个点之后，f值会突然变得平缓，但压力损失却因为流速增加而变得更大，这就是进入湍流区了。
🙋
原来如此！那工程师是怎么算出这个复杂的摩擦系数f的呢？总不能每次都做实验吧？
🎓
问得好！这就是这个模拟器的核心了。对于湍流，我们用一个著名的“Colebrook-White方程”来隐式求解f，它把粗糙度和雷诺数的影响都考虑进去了。这个方程没有直接的公式解，所以工具在后台帮你进行迭代计算。你改变参数后，看到的Moody图上的红点位置和对应的f值，就是方程的解。工程现场常见的是直接查Moody图，但CAE软件和我们的工具都是实时计算的。你试试换不同的流体（比如从水换成空气），看看在相同流量下，摩擦系数和压力损失有什么不同，非常直观！
物理模型与关键公式
计算压力损失的核心是Darcy-Weisbach方程，它表明压降与摩擦系数、管长、流速的平方成正比，与管径成反比。
$$\Delta P = f \cdot \frac{L}{D}\cdot \frac{\rho U^2}{2}$$
其中，$\Delta P$是压力损失（Pa），$f$是Darcy摩擦系数（无量纲），$L$是管长（m），$D$是管径（m），$\rho$是流体密度（kg/m³），$U$是平均流速（m/s）。
摩擦系数$f$取决于流态（雷诺数$Re$）和管道相对粗糙度$\epsilon/D$。对于湍流，由Colebrook-White方程确定：
$$\frac{1}{\sqrt{f}}= -2 \log_{10}\left( \frac{\epsilon/D}{3.7}+ \frac{2.51}{Re \sqrt{f}}\right)$$
其中，$Re = \frac{\rho U D}{\mu}$为雷诺数，$\mu$为动力粘度（Pa·s），$\epsilon$为管壁绝对粗糙度（m）。这是一个隐式方程，需要迭代求解。层流时公式简化为$f=64/Re$。
现实世界中的应用
城市供水与排水系统：设计市政管网时，必须精确计算从水厂到千家万户沿途的压力损失，以确保远端用户有足够的水压。同时，排水管道需要保证足够的坡度（由压降决定）来使污水顺利流动，防止堵塞。
石油与天然气输送：长达数千公里的输油/气管线，泵站或压缩机站的间距和功率完全由管道摩擦压力损失决定。计算不准会导致输送成本激增或流量不达标。
化工与过程工业：在化工厂的反应器、换热器、蒸馏塔之间，流体通过复杂的管道网络输送。准确的压力损失计算对于泵的选型、系统平衡和能耗优化至关重要。
建筑暖通空调（HVAC）：中央空调系统需要将冷/热水或空气输送到建筑的各个房间。风道和水管系统的设计必须计算沿程阻力（摩擦损失）和局部阻力（弯头、阀门等），以选择合适的风机和水泵。
常见误解与注意事项
模型假设：本模拟器所用数学模型基于线性、均质、各向同性等简化假设。在将计算结果直接用于设计决策之前，务必确认实际系统是否满足这些假设。
单位与量纲：许多计算错误源于单位换算错误或数量级判断失误。请时刻注意各参数输入框旁标注的单位。
结果验证：始终将模拟器输出结果与物理直觉或手算结果进行核对。若结果出乎意料，请检查输入参数或采用独立方法进行验证。
进阶学习指引
深化理论：在本工具的简化模型基础上，进一步研究非线性效应、三维行为和时间依赖现象。阅读专业教材和学术论文，掌握严格的数学推导，是提升工程解题能力的关键。
数值方法：系统学习有限元法（FEM）、有限差分法（FDM）和有限体积法（FVM），理解商业CAE求解器的内部运行机制，这将显著提升您设置有效仿真的能力。
实验验证：理论和仿真结果必须通过实验数据加以验证。养成将计算结果与测量值进行对比的习惯，这正是V&V（验证与确认）的精髓所在。
CAE工具：准备好后，可进一步探索Ansys、Abaqus、OpenFOAM、COMSOL等业界主流工具。通过本模拟器培养的物理直觉，将帮助您更有效地配置和使用这些工具。