新放入的球体对原有球体的影响

新放入的球体对原有球体的影响，本质上是二者周围流场通过流体的非线性耦合

。这种相互作用主要通过改变自由液面的形状和流体速度场分布来实现，进而改变原有球体表面的压力场和剪应力场。要定量分析这一过程，我们需要将系统作为一个整体来建立数学模型。

考虑两个球体浸没或漂浮于不可压缩、无粘且运动无旋的理想流体中。第一个球体记为A，第二个新放入的球体记为B。整个系统的速度势函数Phi满足拉普拉斯方程

：

nabla^2 Phi = 0

该方程在整个流体域内成立。为了求解这个偏微分方程，需要给定边界条件。在自由水面z=eta(x,y,t)处，同时满足运动学与动力学边界条件。运动学条件描述自由液面的演化，即水面上的流体微团始终停留在水面上：

partial eta / partial t + partial Phi / partial x * partial eta / partial x + partial Phi / partial y * partial eta / partial y = partial Phi / partial z

动力学条件则要求水面上的压力等于大气压，通过伯努利方程

可表达为：

partial Phi / partial t + 1/2 (nabla Phi)^2 + g eta = 0

在球体A和B的瞬时可动湿表面S_A(t)和S_B(t)上，流体法向速度必须等于该点球体表面的法向运动速度。若球体作为刚体运动，其表面任一点的速度由质心平动速度U和绕质心转动速度Omega决定，则边界条件为：

partial Phi / partial n = (U_A + Omega_A times r_A) dot n_A on S_A(t)

partial Phi / partial n = (U_B + Omega_B times r_B) dot n_B on S_B(t)

这里的n是物面的单位法向量，r是从球心指向表面点的位置矢量。在无穷远处，流体的扰动速度衰减为零。

问题的核心复杂性在于，球体B的运动是未知的，且其运动方程必须与流场方程耦合求解。球体B的运动由牛顿-欧拉方程

描述：

m_B dU_B/dt = m_B g + F_{fluid,B}

I_B dOmega_B/dt = T_{fluid,B}

其中，流体动力F_{fluid,B}和力矩T_{fluid,B}是通过对球体B表面积分流体应力得到的。应力不仅包含压力p，还包括由速度梯度引起的粘性应力（如果考虑粘性）。压力场p通过伯努利方程与速度势Phi关联：p = -rho (partial Phi/partial t + 1/2 (nabla Phi)^2 + gz)。可以看到，球体B的运动改变了流场Phi，而流场的改变又反过来影响其自身的受力，这是一个典型的流固耦合问题。

球体B的引入对流场最直接的影响是产生了扰动势。总速度势可以分解为未受扰动时的势（如均匀流、波浪势或由球体A单独产生的势）与由球体B产生的扰动势之和：Phi = Phi_A + Phi_B。这个扰动势Phi_B会改变球体A周围的局部流速和压力分布。具体而言，由于球体B的存在，流体需要绕流经过两个物体，导致球体A背向球体B一侧的流体速度可能降低，而面向球体B一侧的速度可能升高或产生复杂的涡旋结构（尤其在粘性流体中）。根据伯努利原理，速度的变化直接导致表面压力分布的变化，从而改变了作用在球体A上的净力和力矩。

为了获得解析解，通常采用摄动法或势流理论中的镜像法。假设球体B的引入对流场是一个小扰动，可以将Phi_B视为小量，进而在平衡态附近线性化控制方程和边界条件。对于小振幅运动，我们可以将自由表面条件线性化并在平均自由表面z=0上表达。同时，球体表面的边界条件也需要在平均位置处展开。

对于球体几何，一个强大的分析工具是球谐函数展开

。将扰动势Phi_B在球坐标系下展开为球谐函数的级数和，其系数由边界条件确定。当两个球体相距较远时，可以采用多极子展开法。每个球体对远处流场的影响近似于一个位于球心的奇点（如源、汇或偶极子）。那么，球体B的存在对球体A的影响，就相当于用一个位于B球心的等效奇点产生的流场去修正A周围的流动。通过这种方法，可以将复杂的相互作用问题简化为求解一系列代数方程，从而得到作用在球体A上的附加质量力、阻尼力和恢复力的变化量。

在更一般的粘性流体情形下，控制方程变为不可压缩Navier-Stokes方程：

rho (partial u/partial t + u dot nabla u) = -nabla p + mu nabla^2 u + rho g

nabla dot u = 0

此时，问题变得更加复杂。球体B的引入会改变整个流场的涡量分布。新球体表面会形成新的边界层，其脱落的位置和形态会受到原有球体A尾流的影响。同时，球体A的尾流也会因球体B的存在而发生偏转或破碎。这种强非线性作用通常无法解析求解，必须借助计算流体力学方法。在实际数值模拟中，常采用任意拉格朗日-欧拉方法

来处理移动的球体边界，结合VOF或相场法捕捉自由液面的剧烈变化。通过这些方法，可以精确地计算出球体A所受的流体合力，该合力等于其表面应力张量的积分：

F_A = integral_{S_A} (-p I + tau) dot n dS

其中tau是粘性应力张量。这个合力决定了球体A后续的运动轨迹和姿态变化。最终，新球体的放入对原有球体的影响，正是通过这种复杂的流体动力学相互作用，改变了其表面的应力积分，进而使其偏离了原有的运动状态。