跳转到内容

生理学计算汇总:修订间差异

来自osm&bio
武汉正宗糖三角留言 | 贡献
无编辑摘要
Raspierry留言 | 贡献
第93行: 第93行:
\end{bmatrix}</math>
\end{bmatrix}</math>


因为函数U(x,y,z)在方向l上的导数<math>{\partial U\over\partial l}</math>=<math>\nabla\!U\cdot\hat{l}</math>,其中<math>}\hat{l</math>是l方向的单位向量。
且函数 <math>U(x,y,z)</math> 在方向 <math>\vec l</math> 上的导数 <math>{\partial U\over\partial l}=\nabla\!U\cdot\hat{l}</math>,其中 <math>\hat l</math> 是 <math>\vec l</math> 方向的单位向量。


因此梯度向量指向当地函数值增加最快的方向,'''生理学中所谓的“顺浓度梯度”实际上是逆着浓度梯度向量方向,而像动物行为学当中的“动物顺浓度梯度寻找食物”才是顺浓度梯度向量方向。'''
因此梯度向量指向当地函数值增加最快的方向,'''生理学中所谓的“顺浓度梯度”实际上是逆着浓度梯度向量方向,而像动物行为学当中的“动物顺浓度梯度寻找食物”才是顺浓度梯度向量方向。'''

2026年3月21日 (六) 16:07的版本

虽然25联赛一题未考,但整理一下总归不是坏事

(生理还是要用goldman)

希望osm.bio能支持mhchem包的解析,这么多的\ce{}我可不想换成\mathrm{},毕竟后者不会自动把数字作为下标

后人:目前实在是看不了,还是手动改了,后面插入公式再看怎么引入mhchem扩展包......

二编:发现确实引入了mhchem扩展包,可以使用化学式输入数学公式,但是并不能使用数学公式输入化学式...

eg.

RTln([ClAA][ClAA])F

警告

请各位用户量力而行,以事实准确性为第一,尤其是物理与化学内容不容闪失。

请勿编写自己一知半解的内容或使用AI生成!否则即为对osm.bio的污染。

约定

  • 必须严格区分全微分(正体 d)和偏微分(
  • 必须对首次出现的不常见符号、物理量的含义给出解释
  • 对于梯度,给出三维形式,除非是离子通道等一维背景下的推导
  • 必须紧密结合物理学、物理化学原理与生物学现象、应用
  • 必须区分标量与向量,向量除了梯度算子等少数之外都需要加箭头A或加帽A^

综述

植物生理学和动物生理学是研究生命活动调节与稳态的科学,因此不可避免地会使用物理学、物理化学等方面的工具。此页面聚焦于CHSBO与CNBO的考试范围,给出常见生理学计算题所需使用的数学公式

植物生理学

离子跨膜

杜南平衡

引入

一半透膜将容器分隔为等体积的甲、乙两室。甲室盛有1L与蛋白质结合的Ca2+溶液,Ca2+浓度为0.15mol/L;乙室盛有1L CaCl2溶液,浓度为0.06mol/L。达到杜南平衡时,甲室中Ca2+的浓度为:
A.0.15mol/L
B.0.06mol/L
C.0.024mol/L
D.0.176mol/L

此题需要在微粒数、水量和电中性约束下满足两个平衡:杜南平衡和水势平衡。

化学势与杜南平衡

示例 甲乙两室被半透膜隔离,允许小分子量物质通过,阻挡蛋白质等大分子物质。

现以 CaClA2 为例进行计算。

对于ClA μCl=μ+RTln([Cl])Fφ甲侧 μCl=μ+RTln([Cl])Fφ乙侧ClA 扩散平衡时,两侧电化学势相等 RTln([Cl][Cl])+F(φ乙侧φ甲侧)=0Δφφ甲侧φ乙侧 则上式变为 RTln([Cl][Cl])FΔφ=0 同理,对于 CaA2+ 应用相同条件,得 RTln([Ca2+][Ca2+])+2FΔφ=0

注意离子电荷数

利用两式联立消除 Δφln([Ca2+][Ca2+])=2ln([Cl][Cl]) 同时作为 e 的指数 [Ca2+][Ca2+]=[Cl]2[Cl]2 即得杜南平衡条件

一般形式

Const=[Az]1z[Az]1z

电中性条件

多数情况下,应用杜南平衡条件可得 Δφ0 。但必须注意的是,两侧的溶液一定近似处于电中性状态。这一方面是因为实际上对电场建立有贡献的电荷都处于紧挨着膜的约 1nm 厚的薄层里,膜两侧的静正负电荷互相吸引,从而维持尽可能低的偶极矩;另一方面是因为膜两侧产生的电场在外侧近似抵消,使得 得 Δφ 被严格限制在膜附近极小的范围内。 倘若膜两侧溶液远离膜的部分不成近似电中性,将会产生很大的偶极矩 p 。如前所述,这是不可能存在的。尽管存在更为精确的理论,(如古伊-查普曼模型,Gouy-Chapman Theory),对经典吉布斯-杜南关系用界面理论修正,一般的杜南平衡+电中性条件也足够精确求解生理条件下的平衡了。

物质运输

菲克扩散定律

菲克第一定律:扩散法则

在理想扩散条件下:

J:2mAt=Dc(3维)

J=Dcx(1维)

  • 其中 J扩散通量密度D 为扩散常数(随溶剂和物质等不同而不同),c 为浓度梯度向量

因为梯度的定义是

:[xyz]

且函数 U(x,y,z) 在方向 l 上的导数 Ul=Ul^,其中 l^l 方向的单位向量。

因此梯度向量指向当地函数值增加最快的方向,生理学中所谓的“顺浓度梯度”实际上是逆着浓度梯度向量方向,而像动物行为学当中的“动物顺浓度梯度寻找食物”才是顺浓度梯度向量方向。

菲克第二定律:物质守恒

对于空间中任意一个固定的控制体,考虑由于扩散引起的质量变化。

从扩散通量密度出发,控制体内部物质的量变化率是:

nt|V控制体=SJn^dA

  • 其中 n 是物质的量,n^ 是垂直控制体表面 S 向外的单位法向量,dA 是面积微元

利用散度定理(高斯定理):

SJn^dA=VJdV

  • 其中 dτ 是控制体 V 内部的体积微元

如考虑在单位时间内物质的量改变一值Ct,则垂直控制体内物质的量变化率:

VCtdV

于是:

VCtdVVJdV

因为控制体 V 是任意的,必须有:

CtJ

即为散度形式的连续性方程(物质的量守恒/质量守恒)

带入菲克第一定律:

J=DC

得:

Ct=D2C

便是最常用的菲克第二定律表示形式。

哈根-泊肃叶定律

哈根 - 泊肃叶定律(Hagen-Poiseuille Law),简称泊肃叶定律,是流体力学中描述不可压缩牛顿流体在长直圆管中作稳态层流运动时,流量与压力差、几何尺寸及流体粘度之间关系的物理定律。哈根 - 泊肃叶定律并非适用于所有流体流动,必须严格满足以下条件:

  1. 层流(Laminar Flow):流体必须分层流动,无湍流混合。通常要求雷诺数(Reynolds number, Re )小于临界值(对于圆管,通常 Re<2000∼2300 )。如果流速过快导致湍流,该定律失效(此时压力差与流量的平方成正比)。
  2. 牛顿流体(Newtonian Fluid):流体的粘度 η 必须是常数,不随剪切速率变化。水、空气、低分子量油类符合;血液(非牛顿流体,但在大血管中可近似)、油漆、聚合物熔体在严格意义上不完全符合。
  3. 不可压缩流体:密度恒定。液体通常符合,气体在压差较小时可近似符合。
  4. 稳态流动(Steady Flow):流量不随时间变化。
  5. 长直圆管:
    • 截面必须是完美的圆形。
    • 管子必须足够长,以忽略入口效应(Entrance effect,即流体刚进入管道时速度分布尚未形成抛物线的区域)。通常要求 Lr
  6. 无滑移边界条件:管壁处的流体速度为零
  7. 需要充分发展流,即(v)v=0
  8. 轴对称流动。(仅存在轴向速度 vz(r)

牛顿粘滞定律

τ=ηdudy

严格来说本定律是实验定律,无法使用理论方法推导。

  • τ (Tau):剪切应力(Shear Stress)。
    • 定义:流体层之间单位面积上的内摩擦力。
    • 单位:帕斯卡 ( Pa 或 N/m2 )。
  • du/dy:速度梯度(Velocity Gradient),也称为剪切速率(Shear Rate, γ˙)。
    • 定义:流体速度 u 随垂直于流动方向的距离 y 的变化率。
    • 物理意义:描述流体层之间相对滑动的快慢程度。
    • 单位: s−1 (每秒)。
  • η (或 μ ):动力粘度(Dynamic Viscosity)。
    • 定义:比例常数,表征流体抵抗剪切变形的能力(即“粘稠度”)。
    • 单位: Pa⋅s (帕·秒) 或 N⋅s/m2 。常用单位还有泊 ( P ) 或厘泊 ( cP ),1 Pa⋅s=10 P=1000 cP 。
    • 关键特性:对于牛顿流体, η 是一个常数,只与流体种类和温度有关,不随剪切速率 du/dy​ 的变化而改变

圆管模型

在圆管内,笛卡尔坐标系可变为柱坐标系。此时,在满足基本条件下(见哈根-泊肃叶定律)纳维-斯托克斯方程:

ρ(vt+(v)v)=μ2vp+f

  • ρ 是流体密度。
  • v 是流体速度矢量。
  • t 是时间。
  • p是压力。
  • μ是流体动力粘度。
  • f是体积力。
  • 是梯度算子。
  • 2是拉普拉斯算子

变为平衡方程:

pz+μ(2vzr2+1rvzr)=0

再次变形可得:

1rr(rvzr)=1μpz

第一次积分:

rvzr=12μpzr2+C1

变形后再进行第二次积分:

vz(r)=(12μpzr+C1r)dr=14μpzr2+C1lnr+C2

  • 轴对称条件:在管道中心 (r=0) 处,速度梯度必须有限,因此 C1必须为零(否则 lnr 项在 r=0 处发散)。
  • 无滑移边界条件:在管壁 (r=R) 处,流体速度为零,即 vz(r)=0。此时:

C2=14μpzR2

得到速度剖面:

vz(r)=14μpz(R2r2)

为了得到流量Q,将速度剖面在整个管道横截面A上积分:

Q=Avz(r)dA=0Rvz(r)(2πr)dr

dA是剖面微元,dr是轴向微元。积得:

Q=πR48μpz

如果管道两端的压力差为 Δp=P1P2,管道长度为 L,则 pzΔpL​。此时代入Q得:

Q=πR4Δp4μL

如果将流体粘度μ换成η,就是课本上的哈根-泊肃叶定律。

杨-拉普拉斯方程

1.系统能量

  • 表面能: 弯曲界面的表面能 Es=γA,其中 A是界面面积,γ是表面张力。
  • 压力体积功: 当界面形状发生变化时,体积会发生变化,从而产生压力体积功。对于界面内部流体,其压力为 P内部,外部压力为 P外部​。当界面上的面积微元dA向外微小移动 dn 时,体积变化 dV=dAdn,做的功为 ​(P内部P外部)dV=ΔPdAdn

2.虚功原理:考虑界面发生一个虚位移 δn,即界面法向移动了无限小的距离。系统的总能量变化 δE应该为零,即δEsδW=0其中 δW是压力差产生的虚功。

  • 界面面积 A随法向位移 δn 的变化率为 δA=AHδn,其中H是平均曲率1R1+1R2。因此,表面能的变化为 δEs=γδA=γ(1R1+1R2)Aδn
  • 压力功变化: 压力差产生的虚功为 δW=ΔPAδn

3.平衡情形:将上述变形带入总能量为零的情形下 从而:

γ(1R1+1R2)AδnΔPAδn=0

因为:

Aδn0

所以:

ΔP=γ(1R1+1R2)

动物生理学

生物膜电生理

能斯特方程:电化学势平衡

弦电导方程:线性电阻模型

根据平衡电位的物理意义,当膜电位 Em 等于某离子的平衡电位 EX 时,离子流动的通量就为零了。因此,膜电位与离子平衡电位的差值 EmEX 是离子流动的推动力。离子顺着电势差流动形成电流,所以很容易想到,该离子流动贡献的电流与推动力之间的关系是(记外向电流为正):

IX=GX(EmEX)

稳态时,膜电位不变了,各离子虽然仍有可能流动,但总电流是为零的。于是:

IX=EmGXGXEX

从而:

Em=GXEXGX

特别地,一般我们只考虑钾、钠和氯三种离子通道,所以得到弦电导方程:

Em=Na,K,ClGXEXNa,K,ClGX

GHK方程:扩散-电迁移的精确建模

能斯特-普朗克方程:离子转运

能斯特-普朗克方程(Nernst-Planck equation)是描述带电粒子在电场和浓度梯度共同作用下输运的核心方程。它将离子的总通量分解为扩散、电迁移和对流三个基本分量。在生理学的多个领域都有作用。

  • 能斯特-普朗克方程是理解细胞膜上离子跨膜运动的基础。例如,神经元、心肌细胞等可兴奋细胞膜上钠离子(Na⁺)、钾离子(K⁺)、钙离子(Ca²⁺)和氯离子(Cl⁻)的跨膜流动,对于产生和传播电信号至关重要 。该方程可以定量描述这些离子在浓度梯度(扩散)、电场(电迁移)和对流驱动下的通量 。例如,在微流控通道中,非等温流体的能斯特-普朗克方程通过增加一个与温度和内部电势梯度相关的迁移项进行了修改,以解释热电化学迁移效应 。
  • 神经元的动作电位涉及离子通道的开放和关闭,导致离子通透性快速变化,进而引起膜电位的剧烈波动 。能斯特-普朗克方程及其衍生模型,如与泊松方程耦合的泊松-能斯特-普朗克(PNP)模型,被用于模拟离子在离子通道内的动态运输 。PNP模型在描述离子通道行为方面比GHK模型更精确,尤其是在通道半径小于两倍德拜长度的情况下,因为它考虑了离子在蛋白壁上诱导的表面电荷效应 。通过量化电化学驱动力,能斯特-普朗克方程有助于理解离子通道开放时离子流动的方向和强度,这对于神经信号传递至关重要 。例如,新的动作电位方程为神经电生理学提供了生物物理学见解,可以补充经典的GHK方程 。
  • 近年来,随着神经科学技术的发展,对神经元在完整回路中高时空分辨率记录的需求增加 。电缆-能斯特-普朗克模型被开发用于预测整个树突分支上具有微米特异性的突触活动 。此外,有研究表明,突触输入引起的离子浓度梯度可以增加树突棘中的电压去极化,即使是中等突触电流也能显著改变离子浓度,导致树突和棘头之间出现化学势梯度,产生可测量的电流 。分数阶能斯特-普朗克方程也被引入,用于模拟树突棘的异常次扩散,预测具有更大棘密度的树突上的突触后电位能更快到达胞体并保持在更高水平 。
爱因斯坦关系
  1. 扩散通量:根据菲克第一定律,由于浓度梯度引起的扩散通量Jdiff表示为: Jdiff=DC(详见菲克第一定律)
  2. 电迁移通量:带电离子在电场作用下会发生定向移动,产生电迁移通量Jelec。离子的平均漂移速度vd​与其电迁移率μ和电场强度成正比: vd=μE 。电迁移通量是离子浓度C与漂移速度的乘积: Jelec=Cvd=CμE。电场 E 可以通过电势Φ的梯度来表示:E=Φ。因此,电迁移通量可以写为: Jelec=CμΦ
  3. 热力学平衡:在热力学平衡条件下,系统中的离子没有净通量,即扩散通量与电迁移通量相互抵消: Jdiff=DC=(CμΦ)=Jelec,化简得DC=CμΦ
  4. 在热力学平衡时,离子的电化学势μ~​在整个体系中是均匀的,即μ~=0。离子的电化学势可以表示为: μ~=μ0+RTlnC+zFΦ 。其中,μ0是标准化学势,R 是理想气体常数,T 是绝对温度,z 是离子的化合价,F 是法拉第常数。 对电化学势求梯度: μ~=RTCC+zFΦ,由于在平衡时 μ~=0,则有: C=cFΦRT
  5. 最后代入上述热力学平衡条件,即可导出爱因斯坦关系:μ=DzFRT
电化学建模
  1. 扩散通量 (Ji,diff) 扩散是由于离子浓度不均匀而引起的粒子从高浓度区域向低浓度区域的净移动。这部分通量由菲克第一定律描述:Ji,diff=Dici(3维) 其中:
    • Di 是离子i的扩散系数(单位:m²/s)。
    • Ci​ 是离子i的浓度(单位:mol/m³)。
    • Ci​ 是离子i的浓度梯度。 负号表示通量方向与浓度梯度方向相反,即离子从浓度高的地方向浓度低的地方移动。
  2. 电迁移通量(Ji,elec)电迁移是由于离子带电并在电场作用下定向移动引起的通量。在电场 E的作用下,带电粒子会受到电场力,并以一定的漂移速度移动。根据上文所述(详见爱因斯坦关系),Ji,elec=DiziFRTCiΦ其中:
    • F 是法拉第常数(单位:C/mol)。
    • R 是理想气体常数(单位:J/(mol·K))。
    • T 是绝对温度(单位:K)
  3. 对流通量(Ji,conv)对流是指离子随溶液本身的宏观流动而发生的整体移动,此时有Ji,conv=Civ其中v是流体平均速度。
  4. 将以上三个通量分量相加,立得能斯特-普朗克方程:Ji=DiCiziFRTDiCiΦ+Civ

呼吸循环生理

心脏 pV 图与搏出功

理想条件搏出功

与理想气体 pV 图不同,心脏搏出功是指在1个心动周期内心脏某腔对血液做的功,也即血液对外界做的负功。因此,对于环路积分:

Wblood→out=pdV

可得:

Wheart→blood=Wblood→out=pdV

必须注意此处对环路积分取了相反数。热学中,我们知道气体 pV 图的顺时针环路积分表示体系对外做功;而在搏出功计算时,我们计算的是体系(血液)从外界(心脏)取得了多少正功。因此逆时针环路积分表示体系对外做负功,也即心脏做正功。

误差分析

显然,我们利用 pdV 计算心脏搏出功,暗含了准静态假设:假设此腔内血液是处于准静态平衡的,即在腔收缩与舒张的每一瞬间,腔内各处的压强 p 与腔体积 V 均是一致、确定的,而射血过程无限缓慢。这意味着我们忽略了血液在心腔内由于加速度产生的压强梯度、血液运动产生的动压和粘滞耗散。

倘若我们能测量心腔-血液交界处每一位置的外压总压 p(s) ,准静态假设就毫无必要。使用 0t周期Ap(s,t)v(s,t)dA dt 其中:
p(s,t)心脏壁受压强关于位置 s 和时间 t 的函数
v是血液垂直当地界面的速度分量,v<0 表示向血液侧运

即可精确计算心脏搏出功。可惜这根本无法实现。

弹性阻力

顺应性

顺应性是肺容积变化 dV 与压强变化 dp 的比值。

C顺应性=dVdP

比顺应性

顺应性没有考虑到肺自然状态下的容积差异,使用比顺应性能更客观描述肺的弹性阻力。比顺应性要将顺应性除以静息状态吸气开始的肺容积(功能余气量

C比顺应性=dVV功能余气 dP

万金油型

请此段原作者及时补充内容或将其删除 --Raspierry留言) 2026年3月13日 (五) 09:45 (CST)

一、顺应性

二、黄金代换:做功量

三、压力形成的推导

欧姆定律:

本质:(Flow流量;R阻力)

(S面积;v速度)

综上,