Chem-To-Math的博客

原创 泊松-玻尔兹曼方程的提出

泊松-玻尔兹曼Poisson–Boltzmann方程的提出1.简要历史回顾1.1 基本原理1.2 方程推导1.简要历史回顾泊松-玻尔兹曼方程是用来计算电解质溶液中离子浓度和电荷密度分布的一个微分方程，其基本形式为∇2ϕ(r)=−4πϵ∑ici0ziqe−βziqϕ(r)(1)\nabla^2\phi(\textbf r)=-\frac{4\pi}{\epsilon}\sum_{i} c_i^0z_iqe^{-\beta z_iq\phi(\textbf r)} \tag{1}∇2ϕ(r)=−ϵ4π​i

原创 电化学界面的电化学阻抗-1

An electrochemical double layer (EDL) can be described by Gouy-Chapman-Stern (GCS) model.The EDL is composed of a compact layer between the metal surface and the Helmholtz plane (HP), and a diffuse layer stretching toward the solution bulkThe HP layer ha

原创 Fourier transform 复数形式

Fourier transform 复数形式Z(ω)=∫−∞∞f(t)e−jωtdt Z(\omega)=\int ^{\infty}_{-\infty}f(t)e^{-j\omega t}dtZ(ω)=∫−∞∞​f(t)e−jωtdtwheree−jωt=cos(ωt)−jsin(ωt)e^{-j\omega t}=cos(\omega t)-j sin(\omega t)e−jωt=cos(ωt)−jsin(ωt)ThenZ(ω)=∫−∞∞f(t)[cos(ωt)−jsin(ωt)]dt Z

原创 Convolution Property of the Fourier Transform

Two functions of f(t)f(t)f(t) and g(t)g(t)g(t)The Fourier transform of f(t)f(t)f(t) isF(f(t))=F(ω)=∫−∞∞f(t)exp(−jωt)dt\mathcal{F}(f(t))=F(\omega)=\int ^{\infty} _{-\infty}f(t)exp(-j\omega t)dtF(f(t))=F(ω)=∫−∞∞​f(t)exp(−jωt)dtThe Fourier transform of g(t

原创 Derivative property of the Fourier transform

f(t)f(t)f(t) has the property:The domain is (∞,−∞)(\infty, -\infty)(∞,−∞)The integral of f(t)f(t)f(t) is limited∫−∞∞f(t)dt=A\int ^{\infty} _{-\infty}f(t)dt= A∫−∞∞​f(t)dt=Aandlim⁡t→∞f(t)=0{ \lim_{t \to \infty} f(t)=0}t→∞lim​f(t)=0The Fourier transfor

原创 Inverse Fourier transform

DefinitionF(ω)=∫−∞∞f(t)exp(−jωt)dtF(\omega)=\int ^{\infty} _{-\infty}f(t)exp(-j\omega t)dtF(ω)=∫−∞∞​f(t)exp(−jωt)dtInverse Fourier transformf(t)=F−1(F(ω))=12π∫−∞∞f(ω)exp(jωt)dωf(t)=\mathcal F ^{-1}( F(\omega))=\frac {1}{2\pi}\int ^{\infty} _{-\infty}f(\

原创 2019-2020年催化类专业期刊的影响因子

原创 科学的社会关系

曾经一段时间在思考科学与民主的关系。Science and democracy。正好在机场书店随便翻阅了科学部分，然后然后看到了《极简科学起源课》的一个章节 科学与民主，遂拍了下来，逐字打出来，读之。在希腊的黑暗时代结束之际，一种非常独特的问题名出现在世界面前。这种文明和它之前的迈锡尼威宁截然不同。恢宏的宫殿消失了，被视为半身的君主也消失了。在政治和文化都得到重生的希腊，在也灭有中央集权，没有强...

原创 复杂的拉普拉斯逆变换

Laplace Inverse Transformation (LIT)L(1πte−a2/4t)=e−ass,  a>0\mathcal{L}(\frac{1}{\sqrt{\pi t}}\text{e}^{-a^2/4t})=\frac{\text{e}^{-a\sqrt{s}}}{\sqrt{s}},~~ a>0L(πt​1​e−a2/4t)...

原创 Half-wave potential（半波电位）推导之一

Test conditionNo stirring;Fast kineticsStep to an arbitrary potential suddenlyCottrell-like experimentHalf reaction:(1)O+ne⇌RO+ne{\rightleftharpoons} R \tag{1}O+ne⇌R(1)In the case of Nerstian...

原创 恒电势电流曲线和浓度分布（i-t curve）

Test condition: 1. No stirring; 2. The concentration of reactant at the electrode surface is always 0. Half reaction: O+ne→R(a)(a)O+ne→RO+ne\rightarrow R \tag{a} The equation of the limiting cur...

原创 Laplace transform解传热偏微分方程

Cooling of a Sphere in Contact with a Well-Stirred FluidA homogeneous solid sphere of radius R, initially a t a uniform temperature T1T1T_1, is suddenly immersed at time t=0t=0t=0 in a volume VfVfV...

原创 Reversible Reaction 可逆反应的本质

Reversible Reaction的认知在1803年以前，化学家都认为化学反应式单方向进行的，也就是irreversible reaction。但是就在1803年法国化学家Claude Louis Berthollet先发现了碳酸钠和氯化钙反应生成碳酸钙和氯化钠（Na2CO3+CaCl2→CaCO3+2NaClNa2CO3+CaCl2→CaCO3+2NaClNa_2CO_3+CaCl_2\...

原创 读书方法

《读书方法》这篇文章已经在几年前我买过的一本旧书里看过，感觉张其昀老先生把读书方法和读书的意义阐述的非常清楚和明白，一气呵成，行云流水，字字珠玑。虽然时隔几十年，也甚具实用意义。而今我也经常读这篇文章，常读常新。全篇文字较长，但是并没有影响阅读的负担。读书是学习的一种方法子路曰：“何必读书，然后为学。”这是说读书只是求学方法的一种，求学并非以读书为止境。朱子以即物穷理，释格物致知，物是事物...

原创 Arrhenius 方程的数学表达式背后的物理意义

Arrhenius 方程的数学表达式背后的物理意义当我们谈到一个基元反应方程式，此处重点是基元反应，只有一种过渡态 A+B→C(1)(1)A+B→CA+B \rightarrow C \tag{1} 想要描述在温度TTT的化学反应速率rrr，我们一般用到的表达式是： r=k[cA][cB](2)(2)r=k[cA][cB]r=k[c_A][c_B] \tag{2}如果反应机理不变，改...

原创 碱性溶液中OER动力学分析

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT]i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\beta F(E-...

原创 碱性溶液中HER动力学分析

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT]i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\beta F(E-...

原创 酸性溶液中HER动力学分析

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT]i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\beta F(E-...

原创 应用Bulter-Volmer方程与Langmuir-Hinshelwood 模型分析CO电氧化出峰机理

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT](1)(1)i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\bet...

原创 应用Bulter-Volmer方程与Monte Carlo 模型分析CO电氧化动力学

Kinetics analysis of the CO electrooxidation by Bulter-Volmer equationBulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT](1)(1)i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0...

原创 Butler-Volmer 方程推导交换电流Exchange Current

交换电流Exchange Current 从Bulter-Volmer 方程看出，在平衡状态下，净电流为零，即正逆反应速率相等，电流而相互抵消，这时候的正逆向电流绝对值被称作交换电流Exchange current，i0。 i0=|ia|=|ic|i0=|ia|=|ic|i_0=|i_a|=|i_c|交换电流表达式为 i0=FAk0c∗Oexp−βF(Eeq−E0')RT=FAk0c...

原创 Bulter-Volmer 方程推导 Tafel斜率

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT](1)(1)i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\bet...

原创 Butler-Volmer 方程与Nernst 方程的关系

Bulter-Volmer 方程：i=ic−ia=FAk0[cO(0,t)exp−βF(E−E0′)RT−cR(0,t)exp(1−β)F(E−E0′)RT](1)(1)i=ic−ia=FAk0[cO(0,t)exp⁡−βF(E−E0′)RT−cR(0,t)exp⁡(1−β)F(E−E0′)RT]i=i_c-i_a=FAk^0 \left [c_O(0,t) \exp \dfrac{-\bet...

原创 Butler-Volmer方程推导

巴特勒 - 沃尔默方程推导Derivation of Butler-Volmer Equation阿伦尼乌斯Svante Arrhenius在1889年通过实验经验提出了化学反应速率常数Reaction Constant与活化能Activation Energy 和温度的关系： k=Aexp−EaRT(1)(1)k=Aexp⁡−EaRTk=A\exp\dfrac {-E_a}{...

原创 Nernst 方程推导

电化学势【Electrochemical Potential】定义以及推导在电场中，带电微粒从一种环境变化到另外一种环境的过程中不仅受到化学势的驱使，同时也受到电势的驱使。那么电势如何作用于带电微粒呢？还得从无电场的热力学状态出发，推导出电势对于带电微粒的作用。 在热力学的简单系统中，内能(Internal energy, U)被描述为广度量（extensive variable） 熵(En...

原创 传质方程

传质方程Equation of Continuity for Species A with Constant ρDABρDAB\rho \mathscr{D}_{AB} in terms of wAwAw_A控制方程表达通式： ρDwADt=ρDAB∇2wA+rAρDwADt=ρDAB∇2wA+rA\rho \dfrac {Dw_A}{Dt}=\rho \mathscr{D}_{AB}\...

原创 流体传质方程【Mass Transfer Equation】

流体传质方程Equation of Continuity for Species αα \alpha in terms of jjαjjα \pmb j_\alpha控制方程表达通式： ρDwαDt=−∇⋅jjα+rαρDwαDt=−∇⋅jjα+rα\rho \dfrac {D w_\alpha}{Dt}=-\nabla \cdot \pmb j_\alpha + r_\alpha...

原创 牛顿流体传热方程

牛顿流体传热方程Equation of Energy for Newtonian Fluids with Constant ρρ\rho and kkk控制方程表达通式： ρC^pDTDt=k∇2T+μΦΦvρC^pDTDt=k∇2T+μΦΦv\rho \hat {C}_p \dfrac {DT}{Dt}=k\nabla^2T+\mu \pmb \Phi_v1.直角坐标系(x,y...

原创 流体传热方程【Heat Transfer Equation】

流体传热方程Equation of Energy in terms of qqqq\pmb q控制方程表达通式： ρC^pDTDt=−∇⋅qq−∂lnρ∂lnTDpDt−ττ:∇vvρC^pDTDt=−∇⋅qq−∂lnρ∂lnTDpDt−ττ:∇vv\rho \hat {C}_p \dfrac {DT}{Dt}=-\nabla \cdot \pmb q-\dfrac {\partial ...

原创 牛顿流体耗散方程【Dissipation Function】

牛顿流体耗散方程Dissipation Function for Newtonian Fluids1.直角坐标系(x,y,zx,y,z x, y, z ) 直角坐标系Cartesian coordinates ( x,y,z  x,y,z \textit { x,y,z }): NO. ΦΦv=2[(∂vx∂x)2+(∂...

原创 牛顿流体动量控制方程

流体动量控制方程Equation of Motion for a Newton’s Fluid with Constant ρρ\rho and μμ\mu控制方程通式： ρDvvDt=−∇p+μ∇2vv+ρggρDvvDt=−∇p+μ∇2vv+ρgg\rho \dfrac {D \pmb v}{D t}=- \nabla p + \mu \nabla^2 \pmb v+\rho \...

原创 流体动量控制方程【Motion Equation】

流体动量控制方程The Equation of Motion in terms of ττ\tau控制方程通式： ρDvvDt=−∇p−∇⋅ττ+ρgρDvvDt=−∇p−∇⋅ττ+ρg\rho \dfrac {D \pmb v}{D t}=- \nabla p-\nabla \cdot \pmb \tau+\rho g1.直角坐标系(x,y,zx,y,z x, y, z )...

原创 流体连续性方程【The Equation of Continuity】

流体连续性方程TheEquationofContinuityTheEquationofContinuity\color{Blue} {The\quad Equation \quad of \quad Continuity} 流体连续性方程表达通式 ∂ρ∂t+∇⋅ρvv=0∂ρ∂t+∇⋅ρvv=0\dfrac {\partial \rho}{\partial t}+ {\nabla } \...

原创 菲克扩散(第一)定律【Fick's (First) Law of Binary Diffussion】

菲克扩散(第一)定律Fick’s (First) Law of Binary Diffussion物质扩散通量Diffusion flux表达通式 jjA=−ρDAB∇wAjjA=−ρDAB∇wA\pmb j_A=-\rho \mathscr{D} _{AB} {\nabla } w_A 在单位时间内物质 A 通过垂直于扩散方向的单位截面积的扩散物质流量，kg⋅m−2⋅s−1kg⋅m...

原创 傅里叶热传导定律【Fourier's Law of Heat Conduction】

傅里叶热传导定律Fourier’s Law of Heat Conduction热通量Heat Flux 表达式 qq=−k∇Tqq=−k∇T\pmb q=-k{\nabla }T其中qqqq \pmb q 是单位面积上的热通量， W⋅m−2W·m−2W·m^{−2}, kkk 是材料的热传导系数， W⋅m−1⋅K−1W·m−1·K−1W·m^{−1}·K^{−1}, ...

原创 牛顿黏度定律【Newton's Law of Viscosity】

Newton’s Law of Viscosity先定义 ττ\mathbf{\tau } τ=−μ(∇vv+(∇vv)†)+(23μ−κ)(∇⋅vv)δτ=−μ(∇vv+(∇vv)†)+(23μ−κ)(∇⋅vv)δ \tau=-\mu(\nabla \pmb v+(\nabla \pmb v)^{\dagger} ) +(\frac{2}{3}\mu-\kappa)(\nabla\cd...