跳转到内容

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

来自osm&bio
Raspierry留言 | 贡献
无编辑摘要
Raspierry留言 | 贡献
第69行: 第69行:
在理想扩散条件下:
在理想扩散条件下:


<math display="block">\vec{J}\equiv \frac{\partial^2m}{\partial A \partial t} = -D\nabla C \quad  \text{(3维)}</math>
<math display="block">\vec{J}:\equiv \frac{\partial^2m}{\partial A \partial t} = -D\nabla C \quad  \text{(3维)}</math>


<math display="block">J=-D{\partial C \over \partial x} \quad  \text{(1维)}</math>
<math display="block">J=-D{\partial C \over \partial x} \quad  \text{(1维)}</math>
第75行: 第75行:
其中 <math>\vec J</math> 为'''扩散通量密度''',<math>\nabla x</math> 为浓度梯度向量。因为梯度的定义是
其中 <math>\vec J</math> 为'''扩散通量密度''',<math>\nabla x</math> 为浓度梯度向量。因为梯度的定义是


<math display="block">\nabla \equiv \begin{bmatrix}\frac{\partial}{\partial x}\\ \frac{\partial}{\partial y}\\ \frac{\partial}{\partial z}
<math display="block">\nabla :\equiv \begin{bmatrix}\frac{\partial}{\partial x}\\ \frac{\partial}{\partial y}\\ \frac{\partial}{\partial z}
\end{bmatrix}</math>
\end{bmatrix}</math>



2026年3月18日 (三) 21:20的版本

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

(生理还是要用goldman)

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

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

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

eg.

RTln([ClAA][ClAA])F

综述

植物生理学和动物生理学是研究生命活动调节与稳态的科学,因此不可避免地会使用物理学、物理化学等方面的工具。此条目聚焦于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扩散通量密度x 为浓度梯度向量。因为梯度的定义是

:[xyz]

因此梯度向量指向当地函数值增加最快的方向,“顺浓度梯度”实际上是逆着浓度梯度向量方向。

菲克第二定律:物质守恒

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

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

nt|V控制体=SJn^dA

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

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

SJn^dA=VJdτ

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

而从浓度 c 对时间 t 的偏导出发,亦可得到控制体内部物质的量变化率

ddtVCdτ=VCtdτ

于是:

VCtdτVJdτ

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

CtJ

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


再运用质量守恒律:

● 流入速率:J(x,t)A

● 流出速率:J(x+Δx,t)

● 积累速率:微元内物质总量M对时间t的导数:

积累速率=Mt=AΔxC(x,t)t

将以上各式带入守恒方程:

AΔxCt=[J(x,t)A][J(x+Δx,t)A]

化简得质量守恒条件下的连续性方程:

Ct=Jx

带入非稳态条件下的菲克第一定律:

J=DCx

得:

Ct=D2Cx2

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

哈根-泊肃叶定律

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

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

牛顿粘滞定律

τ=η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​ 的变化而改变

圆管模型

杨-拉普拉斯方程

动物生理学

生物膜电生理

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

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

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

呼吸循环生理

心脏 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速度)

综上,