浸没式光滑点插值法在流固耦合算法中的应用与优化研究

浸没式光滑点插值法在流固耦合算法中的应用与优化研究一、引言1.1研究背景与意义在众多科学与工程领域中，流固耦合现象广泛存在，且发挥着关键作用。航空航天领域里，飞机飞行时机翼与气流间的相互作用即为典型的流固耦合情况。气流作用于机翼，会使其产生变形，而机翼的变形又反过来影响气流的流动状态。这种相互作用对飞机的飞行性能，如稳定性、操纵性等有着重大影响，若不能准确处理，可能导致飞行事故。在船舶与海洋工程方面，船舶在海浪中航行，船身与海水之间存在着复杂的流固耦合。海浪的冲击力会使船身结构承受巨大载荷，引发船身振动、变形，严重时甚至威胁船舶的安全航行；同时，船身的运动也会改变周围海水的流动特性，影响船舶的航行阻力和推进效率。在生物医学工程中，人体血管内血液流动与血管壁的相互作用也涉及流固耦合。血液的流动会对血管壁产生压力和剪切力，长期作用可能导致血管壁的损伤和病变，如动脉粥样硬化等；而血管壁的弹性和变形又会影响血液的流动分布，对人体的血液循环系统产生重要影响。由此可见，流固耦合问题的研究对于保障工程安全、优化设计以及深入理解自然现象和生物生理过程都具有极其重要的意义。随着工程技术的不断进步，对复杂系统中流固耦合问题的研究提出了更高要求。传统的流固耦合求解方法在面对复杂几何形状、大变形以及多尺度问题时，常常面临诸多挑战，如计算精度不足、计算效率低下等。浸没式光滑点插值法（ImmersedSmoothedPointInterpolationMethod，IS-PIM）作为一种新兴的数值方法，为解决流固耦合问题开辟了新途径。它结合了光滑点插值法（SmoothedPointInterpolationMethod，S-PIM）和浸没边界法（ImmersedBoundaryMethod，IBM）的优点，能够有效处理复杂边界条件和大变形问题，具有较高的计算精度和效率。在解决流固耦合问题时，IS-PIM通过引入虚拟流体，将固体区域嵌入到流体计算域中，采用分区域求解的技术，能够灵活地处理固体与流体之间的相互作用。其基于光滑点插值法的梯度光滑技术，可软化固体模型刚度，基于容易剖分的线性背景网格，能改善固体求解精度。采用不同的光滑域构建方式还可得到不同的固体求解器，从而在不同程度上提高计算精度。因此，深入研究基于浸没式光滑点插值法的流固耦合算法，对于提高流固耦合问题的求解能力，推动相关工程领域的发展具有重要的理论和实际应用价值。1.2国内外研究现状流固耦合算法的研究一直是国际学术界和工程领域的热点话题。在国外，早在20世纪初，随着航空航天技术的兴起，流固耦合问题开始受到关注，当时主要聚焦于飞行器的气弹稳定性问题，多采用理论分析和简化模型进行研究。随着计算机技术和数值方法的迅猛发展，20世纪80-90年代，计算流体力学（CFD）和计算结构力学（CSM）迅速崛起，为流固耦合问题的数值模拟带来了重大突破，出现了结合CFD和CSM的耦合算法。进入21世纪，高性能计算技术的进步使得大规模复杂流固耦合问题的求解成为可能，应用领域也从航空航天逐步拓展到土木、水利、能源、生物医学等多个领域。在流固耦合算法的具体研究方面，国外学者在不同的数值方法上取得了诸多成果。有限元法（FEM）作为一种经典的数值方法，在流固耦合领域应用广泛。比如，在船舶与海洋工程中，利用有限元法对船舶结构在波浪载荷作用下的响应进行分析，但传统有限元法在处理固体模型时存在刚度过硬的问题，导致低阶单元求解精度较低。为了克服这一缺陷，光滑点插值法（S-PIM）应运而生。S-PIM得益于梯度光滑技术能软化固体模型刚度，基于容易剖分的线性背景网格能改善固体求解精度。采用不同的光滑域构建方式可以得到不同的固体求解器，从而在不同程度上提高计算精度。浸没式光滑点插值法（IS-PIM）是在S-PIM基础上发展而来的，它基于浸没类方法的思想，采用分区域求解的技术，通过引入虚拟流体来求解流固耦合系统。国外有学者运用IS-PIM对生物流体问题进行研究，成功模拟了血液在血管中的流动与血管壁的相互作用。在国内，流固耦合算法的研究起步相对较晚，但发展迅速。近年来，随着国家对航空航天、海洋工程等领域的大力支持，国内学者在流固耦合算法研究方面取得了显著成果。在航空领域，研究人员针对飞机机翼颤振问题，运用各种流固耦合算法进行数值模拟和分析，通过改进算法提高了颤振预测的准确性。在海洋工程中，对海洋平台在海浪作用下的流固耦合响应开展了深入研究，采用不同的数值方法建立了海洋平台的流固耦合模型，分析了平台的受力和变形情况。对于浸没式光滑点插值法，国内学者也进行了大量研究。例如，针对原始IS-PIM中流固耦合力计算的缺陷，提出了一种将粘性力基于流体欧拉网格进行计算的流固耦合力计算形式。通过对典型流固耦合数值算例进行计算，证明了该方法具有更高的准确性。同时，针对原始IS-PIM存在的流固边界不准确、“新鲜点”等问题，分别采用虚拟点修正算法、质量源/沉算法以及尖锐界面方法等对其进行修正，修正后的方法在计算精度和稳定性方面都有了明显提升。总体而言，国内外在流固耦合算法研究方面已经取得了丰硕的成果，但随着工程技术的不断发展，对算法的精度、效率和适用范围提出了更高的要求。浸没式光滑点插值法作为一种新兴的流固耦合算法，虽然在理论和应用方面取得了一定进展，但仍存在一些问题需要进一步研究和解决，具有广阔的研究空间。1.3研究内容与方法1.3.1研究内容本论文围绕基于浸没式光滑点插值法的流固耦合算法展开深入研究，主要内容包括以下几个方面：浸没式光滑点插值法的算法原理与理论基础研究：详细剖析浸没式光滑点插值法的基本原理，包括光滑点插值法的梯度光滑技术、不同光滑域构建方式及其对固体求解精度的影响，以及浸没边界法引入虚拟流体、分区域求解的技术细节。深入研究该方法所涉及的流体动力学和固体力学的基本方程，如Navier-Stokes方程、弹性力学方程等，以及它们在流固耦合界面的耦合条件，为后续算法的改进和应用奠定坚实的理论基础。基于浸没式光滑点插值法的流固耦合问题分析：运用浸没式光滑点插值法对典型的流固耦合问题进行数值模拟和分析，如圆柱绕流与弹性支撑结构的耦合问题、管道内流体流动与管壁的相互作用问题等。通过模拟，深入研究流固耦合过程中流体与固体的相互作用机制，包括流体对固体的作用力、固体变形对流体流动的影响等。分析不同参数（如流体的雷诺数、固体的弹性模量等）对耦合结果的影响规律，为实际工程应用提供理论指导。浸没式光滑点插值法的算法改进与优化：针对现有浸没式光滑点插值法存在的问题，如流固边界不准确、“新鲜点”问题以及流固耦合力计算的缺陷等，提出相应的改进措施。采用虚拟点修正算法、质量源/沉算法以及尖锐界面方法等对原始算法进行修正，提高算法对固体边界的捕捉能力和计算精度。提出新的流固耦合力计算形式，将粘性力基于流体欧拉网格进行计算，以更准确地计算固体所受的流固耦合力。通过数值算例对比分析改进前后算法的性能，验证改进措施的有效性。基于改进算法的流固耦合问题应用研究：将改进后的浸没式光滑点插值法应用于实际工程中的流固耦合问题，如航空发动机叶片的气弹稳定性分析、船舶在海浪中的运动响应模拟等。结合具体工程问题，建立相应的流固耦合模型，进行数值模拟和分析。根据模拟结果，为工程设计提供优化建议，如调整结构参数、改进材料性能等，以提高工程结构在流固耦合环境下的性能和可靠性。1.3.2研究方法为实现上述研究内容，本论文拟采用以下研究方法：理论分析方法：对流体动力学、固体力学以及浸没式光滑点插值法的相关理论进行深入研究，推导流固耦合问题的基本方程和耦合条件。分析现有算法存在的问题及其理论根源，为算法的改进提供理论依据。通过理论分析，明确流固耦合过程中各物理量之间的关系，揭示流固耦合现象的本质。数值模拟方法：利用数值计算软件，基于浸没式光滑点插值法建立流固耦合模型，对各种流固耦合问题进行数值模拟。通过设置不同的参数和边界条件，模拟不同工况下的流固耦合现象，获取丰富的数值结果。对数值结果进行分析和处理，研究流固耦合的规律和特性，验证理论分析的正确性。对比分析方法：将改进后的浸没式光滑点插值法与传统的流固耦合算法以及原始的浸没式光滑点插值法进行对比。对比不同算法在计算精度、计算效率、稳定性等方面的性能差异，评估改进算法的优势和不足。通过对比分析，明确改进算法的适用范围和应用前景，为算法的进一步优化提供参考。二、流固耦合及浸没式光滑点插值法理论基础2.1流固耦合基本理论2.1.1流固耦合的定义与分类流固耦合是流体力学与固体力学交叉而生成的一门力学分支，主要研究变形固体在流场作用下的各种行为以及固体位形对流场的影响，其核心在于两相介质之间的相互作用。当流体作用于固体时，会使固体产生变形或运动；而固体的变形或运动又会反过来改变流体的运动状态，这种相互作用在不同条件下会产生各种各样的流固耦合现象。例如，飞机机翼在气流作用下产生的颤振，气流的作用力使机翼发生振动变形，而机翼的振动变形又改变了周围气流的流动特性，进而影响飞机的飞行性能。按照耦合机理，流固耦合问题大致可分为两类。第一类问题的耦合作用仅发生在流体与固体的相交界面上，在方程上的耦合是通过两相耦合面上的平衡及协调关系来引入的，像气动弹性、水动弹性等问题都属于这一类。在这类问题中，场间不相互重叠与渗透，耦合效应主要通过界面力起作用，计算时只需满足耦合界面力平衡和界面相容条件即可。例如桥梁的风致振动，风作为流体作用于桥梁结构，桥梁在风力作用下产生振动变形，其变形又影响风的流动状态，而这种相互作用主要集中在桥梁与风接触的表面，即耦合界面。第二类问题中，流体域与固体域部分或全部重叠在一起，难以清晰地分开，描述物理现象的方程，尤其是本构方程需要针对具体的物理现象来建立，其耦合效应通过建立与不同单相介质的本构方程等微分方程来体现。从耦合程度的角度来看，流固耦合又可分为强耦合和弱耦合。强耦合问题中，流体和固体之间的相互作用非常强烈，二者的物理场紧密关联，需要同时考虑流体和固体的所有控制方程，通过将流场和结构场的控制方程耦合到同一方程矩阵中进行求解。这种方法理论上较为先进，适用于处理大固体变形、生物隔膜运动等复杂情况，但在实际应用中，由于需要将计算流体动力学和计算固体力学技术深度融合，同步求解的收敛难度较大，耗时也较长，目前主要应用于一些简单的热－结构耦合和电磁－结构耦合问题的模拟分析，在流体－结构耦合的实际工程应用中还存在一定困难。例如，在模拟生物心脏中血液流动与心肌组织的相互作用时，由于心肌的大变形和血液流动的复杂性，属于强耦合问题，求解难度较大。弱耦合问题中，流体和固体之间的相互作用相对较弱，可分别求解流体和固体的控制方程，然后通过流固耦合交界面进行数据传递。这种方法对计算机性能的要求相对较低，能够用于求解实际的大规模问题，目前在商业软件的流固耦合分析中基本都采用这种分离解法。例如，在分析飞机机翼在气流中的受力情况时，若机翼变形相对较小，对气流的影响不太显著，就可以采用弱耦合的分离解法，先计算气流对机翼的作用力，再将该力施加到机翼结构上计算其变形。2.1.2流固耦合的数学模型流固耦合数学模型的建立基于流体力学和固体力学的基本方程，以及它们在耦合界面上的相互作用条件。在流体力学中，常用的控制方程包括连续性方程、动量方程（Navier-Stokes方程，简称N-S方程）和能量方程。连续性方程描述了流体在流动过程中的质量守恒，其数学表达式为：\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0其中，\rho是流体密度，t是时间，\vec{v}是流体速度矢量。动量方程（N-S方程）体现了流体的动量守恒，对于不可压缩粘性流体，其形式为：\rho(\frac{\partial\vec{v}}{\partialt}+\vec{v}\cdot\nabla\vec{v})=-\nablap+\mu\nabla^{2}\vec{v}+\vec{f}这里，p是流体压力，\mu是流体动力粘度，\vec{f}是作用在流体上的体积力。能量方程用于描述流体的能量守恒，在考虑热传导和粘性耗散的情况下，其表达式较为复杂，对于一些不涉及明显热效应的流固耦合问题，可能暂不考虑能量方程。在固体力学中，基本方程主要基于弹性力学理论。对于线性弹性小变形问题，常用的方程包括平衡方程、几何方程和物理方程。平衡方程表示固体微元体上的力平衡，在笛卡尔坐标系下的表达式为：\frac{\partial\sigma_{ij}}{\partialx_{j}}+f_{i}=0其中，\sigma_{ij}是应力张量分量，x_{j}是坐标分量，f_{i}是单位体积的体力分量。几何方程描述了固体的应变与位移之间的关系，对于小变形情况，有：\varepsilon_{ij}=\frac{1}{2}(\frac{\partialu_{i}}{\partialx_{j}}+\frac{\partialu_{j}}{\partialx_{i}})这里，\varepsilon_{ij}是应变张量分量，u_{i}和u_{j}分别是位移矢量在x_{i}和x_{j}方向上的分量。物理方程则建立了应力与应变之间的关系，对于各向同性线性弹性材料，服从胡克定律，其表达式为：\sigma_{ij}=\lambda\varepsilon_{kk}\delta_{ij}+2\mu\varepsilon_{ij}其中，\lambda和\mu是拉梅常数，\varepsilon_{kk}是体积应变，\delta_{ij}是克罗内克符号。在流固耦合问题中，流体域和固体域通过耦合界面相互作用，因此需要满足一定的耦合界面条件。主要包括以下几个方面：一是力平衡条件，即作用在耦合界面上的流体应力和固体应力相等，在法向和切向都应满足该条件，可表示为\vec{\sigma}_{f}\cdot\vec{n}=\vec{\sigma}_{s}\cdot\vec{n}和\vec{\tau}_{f}\cdot\vec{t}=\vec{\tau}_{s}\cdot\vec{t}，其中\vec{\sigma}_{f}和\vec{\sigma}_{s}分别是流体和固体的应力矢量，\vec{n}是耦合界面的法向单位矢量，\vec{\tau}_{f}和\vec{\tau}_{s}是切应力矢量，\vec{t}是切向单位矢量；二是位移和速度连续条件，即耦合界面上流体和固体的位移和速度相等，可表示为\vec{u}_{f}=\vec{u}_{s}和\vec{v}_{f}=\vec{v}_{s}，其中\vec{u}_{f}和\vec{u}_{s}分别是流体和固体的位移矢量，\vec{v}_{f}和\vec{v}_{s}是速度矢量。基于上述流体力学和固体力学的基本方程以及耦合界面条件，可建立流固耦合的数学模型。通常采用数值方法，如有限元法、有限差分法、光滑点插值法等对该数学模型进行离散求解。在浸没式光滑点插值法中，通过引入虚拟流体，将固体区域嵌入到流体计算域中，基于光滑点插值法对固体部分进行离散，采用合适的流体求解器对流体部分进行求解，然后根据耦合界面条件实现流体和固体之间的信息传递和耦合求解。2.2浸没式光滑点插值法原理2.2.1光滑点插值法（S-PIM）光滑点插值法（S-PIM）是一种基于点的数值方法，其核心在于通过梯度光滑技术来软化固体模型刚度。在传统的有限元法中，固体模型的刚度往往相对较大，这会导致在使用低阶单元进行求解时精度较低。而S-PIM通过在每个节点周围构建一个光滑域，利用该光滑域内节点的信息进行插值计算，从而有效地软化了模型刚度。以二维问题为例，假设在求解域内有一系列离散的节点，对于每个节点，选取其周围一定范围内的节点构成光滑域。在光滑域内，基于拉格朗日插值法构建插值函数。例如，对于一个三角形光滑域，可利用三角形三个顶点的函数值来构建线性插值函数。通过对插值函数求梯度，得到该光滑域内的梯度近似值。将各个光滑域的梯度进行组合，就可以得到整个求解域的近似梯度。这种通过光滑域进行梯度计算的方式，与传统有限元法直接基于单元进行梯度计算不同，它能更灵活地处理节点分布，从而软化固体模型刚度。S-PIM基于容易剖分的线性背景网格来改善固体求解精度。线性背景网格的剖分相对简单，能够快速地对复杂形状的求解域进行离散。在这些线性背景网格上布置节点，利用节点间的插值关系，可以更准确地逼近固体的真实力学行为。与传统的有限元网格相比，线性背景网格在处理复杂几何形状时具有更高的适应性，不需要像有限元网格那样进行复杂的网格划分和网格质量控制。而且，不同的光滑域构建方式可以得到不同的固体求解器。例如，除了三角形光滑域外，还可以采用四边形光滑域、圆形光滑域等。不同形状的光滑域在计算精度和计算效率上可能会有所差异。圆形光滑域在处理各向同性问题时可能具有更好的精度，因为它在各个方向上对节点信息的利用更加均匀；而四边形光滑域在处理矩形或近似矩形的求解域时，可能在计算效率上更具优势，因为其与线性背景网格的匹配度更好。通过选择合适的光滑域构建方式，可以在不同程度上提高固体求解的精度和效率。2.2.2浸没式光滑点插值法（IS-PIM）的基本思想浸没式光滑点插值法（IS-PIM）基于浸没类方法的思想，采用分区域求解的技术，通过引入虚拟流体来求解流固耦合系统。在处理流固耦合问题时，传统方法常常面临复杂边界条件和大变形问题带来的挑战，而IS-PIM通过独特的分区域求解策略有效地克服了这些问题。其基本思路是将流固耦合系统划分为流体区域和固体区域，在流体区域采用适合流体计算的方法，在固体区域则采用光滑点插值法进行离散和求解。为了实现固体区域与流体区域的耦合，IS-PIM引入了虚拟流体。虚拟流体被定义在固体区域内，它与真实流体具有相似的物理性质，如密度、粘度等。通过虚拟流体，将固体区域嵌入到流体计算域中，使得整个计算域在形式上成为一个统一的流体域。这样，在求解过程中，可以使用统一的流体求解器对整个计算域进行求解。在模拟圆柱绕流与弹性支撑结构的流固耦合问题时，将弹性支撑结构所在区域视为固体区域，在该区域内定义虚拟流体。将包含固体区域和流体区域的整个计算域作为一个大的流体计算域，采用有限体积法等流体求解方法对该计算域进行离散和求解。在固体区域，利用光滑点插值法对固体的力学行为进行描述。通过在固体区域内布置离散节点，构建光滑域，基于光滑点插值法计算固体的位移、应力等物理量。在流固耦合界面，IS-PIM满足一定的耦合条件。根据力平衡条件，在耦合界面上，流体对固体的作用力与固体对流体的反作用力大小相等、方向相反。从动量守恒的角度来看，流体的动量变化会通过耦合界面传递给固体，使固体产生相应的运动和变形；反之，固体的运动和变形也会影响流体的动量分布。从能量守恒的角度分析，流固耦合过程中，流体和固体之间会发生能量的传递和转换。在圆柱绕流问题中，流体的动能会传递给弹性支撑结构，使其发生振动，结构的振动又会将部分能量反馈给流体，影响流体的流动状态。通过满足这些耦合条件，实现了流体和固体之间的信息传递和相互作用，从而准确地模拟流固耦合现象。2.2.3IS-PIM的求解流程IS-PIM的求解流程包括控制方程离散、求解、耦合条件施加及最终求解等多个关键步骤。在控制方程离散阶段，对于流体部分，基于Navier-Stokes方程，采用有限体积法等数值方法将其离散到欧拉网格上。以二维不可压缩粘性流体为例，将计算域划分为一系列的矩形或三角形网格单元。对于每个网格单元，根据连续性方程和动量方程，将其积分形式离散为代数方程。在离散连续性方程时，通过对网格单元边界上的质量通量进行积分，得到关于网格单元内流体密度变化的离散方程；在离散动量方程时，对网格单元边界上的压力和粘性应力进行积分，得到关于网格单元内流体速度变化的离散方程。对于固体部分，基于光滑点插值法，将弹性力学方程离散到拉格朗日网格上。在固体区域内布置离散节点，以每个节点为中心构建光滑域。对于每个光滑域，根据弹性力学的平衡方程、几何方程和物理方程，利用光滑点插值法构建插值函数，将方程离散为关于节点位移的代数方程。在求解阶段，先对离散后的流体控制方程进行求解，得到流体的速度和压力分布。采用合适的迭代算法，如SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）等，对离散后的代数方程组进行迭代求解。在每次迭代中，先根据上一次迭代得到的压力场计算速度场，然后根据速度场的变化更新压力场，直到速度场和压力场收敛。接着对离散后的固体控制方程进行求解，得到固体的位移和应力分布。同样采用迭代算法，如牛顿-拉夫逊迭代法等，对关于节点位移的代数方程组进行求解。在迭代过程中，根据当前的位移猜测值计算应力，然后根据应力与位移的关系更新位移猜测值，直到位移和应力收敛。在耦合条件施加阶段，将求解得到的流体速度和压力传递到流固耦合界面，作为固体的边界条件；同时，将求解得到的固体位移和速度传递到流固耦合界面，作为流体的边界条件。在圆柱绕流与弹性支撑结构的耦合问题中，将流体在耦合界面上的压力和粘性力作为弹性支撑结构的外力，施加到固体控制方程中；将弹性支撑结构在耦合界面上的位移和速度作为流体的边界位移和边界速度，施加到流体控制方程中。通过这种方式，实现了流体和固体之间的信息传递和相互作用。对耦合后的方程进行最终求解，得到整个流固耦合系统的解。在每次时间步长内，重复上述离散、求解和耦合条件施加的过程，逐步推进时间，直到达到设定的时间终点。在求解过程中，需要不断地检查计算结果的收敛性和稳定性。如果计算结果不收敛，需要调整迭代算法的参数，如松弛因子等；如果计算结果不稳定，可能需要调整网格尺寸、时间步长等计算参数，以确保计算的准确性和可靠性。三、基于IS-PIM的流固耦合算法分析3.1算法实现过程3.1.1流体求解器在基于浸没式光滑点插值法（IS-PIM）的流固耦合算法中，流体求解器的选择至关重要，它直接影响到对流体运动的模拟精度和效率。以半隐式特征线分离算法（CBS）为例，其在求解不可压缩粘性流体问题时展现出独特的优势。CBS算法的核心在于对不可压缩粘性流体控制方程的有效处理。不可压缩粘性流体的控制方程主要包括连续性方程和动量方程（Navier-Stokes方程）。连续性方程体现了流体在流动过程中的质量守恒特性。对于二维不可压缩粘性流体，在笛卡尔坐标系下，连续性方程可表示为：\frac{\partialu}{\partialx}+\frac{\partialv}{\partialy}=0其中，u和v分别是流体速度在x和y方向上的分量。动量方程则描述了流体的动量守恒，其表达式为：\begin{cases}\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy})=-\frac{\partialp}{\partialx}+\mu(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}})+f_{x}\\\rho(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy})=-\frac{\partialp}{\partialy}+\mu(\frac{\partial^{2}v}{\partialx^{2}}+\frac{\partial^{2}v}{\partialy^{2}})+f_{y}\end{cases}这里，\rho是流体密度，p是流体压力，\mu是流体动力粘度，f_{x}和f_{y}分别是作用在流体上的体积力在x和y方向上的分量。半隐式CBS算法的具体求解步骤如下：在第一步，先将动量方程中的对流项和扩散项显式处理，而压力项则隐式处理。通过这种方式，得到一个关于速度的预测值。在二维情况下，假设当前时间步为n，下一时间步为n+1，对于速度预测值\widetilde{u}^{n+1}和\widetilde{v}^{n+1}，可通过如下公式计算：\begin{cases}\frac{\rho(\widetilde{u}^{n+1}-u^{n})}{\Deltat}=-\rho(u^{n}\frac{\partialu^{n}}{\partialx}+v^{n}\frac{\partialu^{n}}{\partialy})+\mu(\frac{\partial^{2}u^{n}}{\partialx^{2}}+\frac{\partial^{2}u^{n}}{\partialy^{2}})+f_{x}^{n}-\frac{\partialp^{n}}{\partialx}\\\frac{\rho(\widetilde{v}^{n+1}-v^{n})}{\Deltat}=-\rho(u^{n}\frac{\partialv^{n}}{\partialx}+v^{n}\frac{\partialv^{n}}{\partialy})+\mu(\frac{\partial^{2}v^{n}}{\partialx^{2}}+\frac{\partial^{2}v^{n}}{\partialy^{2}})+f_{y}^{n}-\frac{\partialp^{n}}{\partialy}\end{cases}其中，\Deltat是时间步长。在第二步，利用连续性方程来修正速度预测值，以满足质量守恒条件。由于连续性方程要求\frac{\partial\widetilde{u}^{n+1}}{\partialx}+\frac{\partial\widetilde{v}^{n+1}}{\partialy}=0，但通过第一步计算得到的速度预测值可能并不满足这一条件，因此需要引入压力修正项。设压力修正值为p^{'}，速度修正值为u^{'}和v^{'}，则有：\begin{cases}u^{n+1}=\widetilde{u}^{n+1}-\frac{\Deltat}{\rho}\frac{\partialp^{'}}{\partialx}\\v^{n+1}=\widetilde{v}^{n+1}-\frac{\Deltat}{\rho}\frac{\partialp^{'}}{\partialy}\end{cases}将上述修正后的速度代入连续性方程，可得到一个关于压力修正值p^{'}的泊松方程：\frac{\partial^{2}p^{'}}{\partialx^{2}}+\frac{\partial^{2}p^{'}}{\partialy^{2}}=\frac{\rho}{\Deltat}(\frac{\partial\widetilde{u}^{n+1}}{\partialx}+\frac{\partial\widetilde{v}^{n+1}}{\partialy})通过求解该泊松方程，得到压力修正值p^{'}，进而得到满足连续性方程的速度u^{n+1}和v^{n+1}。在第三步，根据得到的速度和压力，更新压力场。压力更新公式为：p^{n+1}=p^{n}+\rho\frac{\Deltat}{\Deltax}(\frac{\partialu^{n+1}}{\partialx}+\frac{\partialv^{n+1}}{\partialy})通过以上三个步骤的循环迭代，不断更新速度和压力场，直至满足收敛条件，从而得到不可压缩粘性流体在不同时刻的速度和压力分布。3.1.2固体求解器光滑点插值方法（S-PIM）作为IS-PIM中的固体求解器，其工作原理基于独特的节点插值和光滑域构建技术。在固体力学分析中，传统的有限元法存在一些局限性，如对复杂几何形状的适应性较差、低阶单元求解精度不足等。而S-PIM通过创新的方式有效地克服了这些问题。S-PIM基于节点构建光滑域，以实现对固体力学行为的准确描述。在求解域内，离散分布着一系列节点。以二维问题为例，对于每个节点，选取其周围一定范围内的节点构成光滑域。常见的光滑域形状有三角形、四边形等。以三角形光滑域为例，在该光滑域内，基于拉格朗日插值法构建插值函数。假设三角形光滑域的三个顶点分别为i、j、k，其对应的函数值分别为u_{i}、u_{j}、u_{k}，则该光滑域内任意一点的函数值u可通过如下线性插值公式计算：u=N_{i}u_{i}+N_{j}u_{j}+N_{k}u_{k}其中，N_{i}、N_{j}、N_{k}是基于三角形面积坐标定义的插值基函数。通过对插值函数求梯度，可得到该光滑域内的梯度近似值。例如，对于位移函数u的梯度\nablau，在三角形光滑域内的近似值可表示为：\nablau=\nablaN_{i}u_{i}+\nablaN_{j}u_{j}+\nablaN_{k}u_{k}将各个光滑域的梯度进行组合，就可以得到整个求解域的近似梯度。这种基于光滑域的梯度计算方式，与传统有限元法直接基于单元进行梯度计算不同，它能更灵活地处理节点分布，从而软化固体模型刚度。在处理复杂形状的固体结构时，传统有限元法可能需要划分大量的小单元来逼近结构形状，这会导致计算量大幅增加且单元质量难以保证。而S-PIM通过构建光滑域，能够在较少的节点下准确地描述结构的力学行为，提高了计算效率和精度。S-PIM基于容易剖分的线性背景网格，进一步改善了固体求解精度。线性背景网格的剖分相对简单，能够快速地对复杂形状的求解域进行离散。在这些线性背景网格上布置节点，利用节点间的插值关系，可以更准确地逼近固体的真实力学行为。与传统的有限元网格相比，线性背景网格在处理复杂几何形状时具有更高的适应性，不需要像有限元网格那样进行复杂的网格划分和网格质量控制。不同的光滑域构建方式会对固体求解器的性能产生影响。除了三角形光滑域外，采用四边形光滑域时，由于其在处理矩形或近似矩形的求解域时与线性背景网格的匹配度更好，可能在计算效率上更具优势。而圆形光滑域在处理各向同性问题时，由于其在各个方向上对节点信息的利用更加均匀，可能具有更好的精度。通过选择合适的光滑域构建方式，可以在不同程度上提高固体求解的精度和效率。3.1.3流固耦合条件的施加在IS-PIM中，流固耦合条件的施加是实现准确模拟流固耦合现象的关键环节，主要通过虚拟流体来实现。虚拟流体的引入是IS-PIM的核心技术之一，它为流固耦合问题的求解提供了一种有效的途径。在流固耦合问题中，固体区域被嵌入到流体计算域中，为了实现固体与流体之间的相互作用，在固体区域内定义虚拟流体。虚拟流体与真实流体具有相似的物理性质，如密度、粘度等。通过虚拟流体，将固体区域视为流体计算域的一部分，使得在整个计算域上可以采用统一的数值方法进行求解。在模拟圆柱绕流与弹性支撑结构的流固耦合问题时，将弹性支撑结构所在区域定义为虚拟流体区域。在该区域内，虚拟流体的密度和粘度等参数可以根据固体的材料特性进行合理设置。通过虚拟流体施加流固耦合速度条件是实现流固耦合的重要步骤。在流固耦合界面，需要满足流体和固体的速度相等这一条件。在IS-PIM中，通过对虚拟流体和真实流体的速度进行匹配来实现这一条件。具体来说，在每个时间步，将固体表面节点的速度传递给与之相邻的虚拟流体节点，使得虚拟流体节点的速度与固体表面节点的速度一致。在模拟弹性梁在流体中振动的问题时，在每个时间步，根据固体求解器计算得到的弹性梁表面节点的速度，更新与之相邻的虚拟流体节点的速度。这样，在流固耦合界面，虚拟流体和真实流体的速度能够保持连续，从而满足流固耦合速度条件。计算流固耦合力是流固耦合条件施加的另一个关键方面。在原始的IS-PIM中，流固耦合力的计算是通过虚拟流体域的拉格朗日网格进行的，但这种方式存在一定的缺陷，它忽略了作用在固体边界上的粘性力大小，导致计算出的固体所受流固耦合力的大小不准确。针对这一问题，改进后的方法将粘性力基于流体欧拉网格进行计算。具体做法是，先根据流体求解器得到的流体速度和压力分布，在流体欧拉网格上计算作用在固体边界上的粘性力。根据Navier-Stokes方程，粘性力与流体速度的梯度相关。通过对流体速度在欧拉网格上进行差分计算，可以得到速度梯度，进而计算出粘性力。然后，将计算得到的粘性力与通过虚拟流体域拉格朗日网格计算得到的惯性力等其他力分量相结合，得到准确的流固耦合力。在模拟低雷诺数下的圆柱绕流问题时，改进后的方法能够更准确地计算圆柱所受的流固耦合力，提高了模拟结果的准确性。通过准确计算流固耦合力，并将其作为固体的外力施加到固体求解器中，实现了流体和固体之间的力的传递和相互作用，从而完成了流固耦合条件的施加。3.2算法优势与特点3.2.1处理复杂边界的能力在处理复杂边界问题上，浸没式光滑点插值法（IS-PIM）展现出了卓越的能力，这是其相较于传统流固耦合算法的显著优势之一。传统的有限元法（FEM）在面对复杂几何形状的固体边界时，往往需要花费大量的时间和精力进行网格划分。在对航空发动机叶片进行流固耦合分析时，叶片的形状复杂，具有扭曲、变截面等特点。使用FEM进行网格划分时，为了准确地描述叶片的几何形状，需要生成大量的小单元，这不仅增加了网格划分的难度和时间，还可能导致网格质量下降，影响计算精度。而且，在固体发生大变形的情况下，传统FEM的网格会发生严重扭曲，甚至出现网格畸变，使得计算无法继续进行。相比之下，IS-PIM基于浸没类方法的思想，采用分区域求解的技术，通过引入虚拟流体来求解流固耦合系统。在处理复杂边界时，它不需要像传统方法那样对固体边界进行精确的网格划分。在模拟具有复杂外形的船舶在海浪中的流固耦合问题时，IS-PIM只需将船舶的固体区域嵌入到流体计算域中，通过虚拟流体来处理固体与流体之间的相互作用。在固体区域内定义虚拟流体，将其视为流体计算域的一部分，这样就避免了对船舶复杂外形进行繁琐的网格划分。同时，由于虚拟流体的引入，使得在整个计算域上可以采用统一的数值方法进行求解，大大简化了计算过程。IS-PIM还能有效处理固体的大变形问题。当固体发生大变形时，传统方法的网格会随着固体的变形而发生严重扭曲，导致计算精度下降甚至计算失败。而IS-PIM采用拉格朗日网格来描述固体的运动和变形，拉格朗日网格能够跟随固体一起运动和变形，从而避免了网格畸变的问题。在模拟弹性梁在流体中发生大变形的情况时，IS-PIM通过拉格朗日网格能够准确地捕捉到弹性梁的变形过程，保证了计算的准确性。通过虚拟点修正算法、质量源/沉算法以及尖锐界面方法等对原始IS-PIM存在的流固边界不准确、“新鲜点”等问题进行修正，进一步提高了其对固体边界的捕捉能力和计算精度。这些改进措施使得IS-PIM在处理复杂边界和大变形问题时更加稳定和可靠。3.2.2计算效率的提升在计算效率方面，IS-PIM通过独特的算法设计和优化，实现了显著的提升，为大规模流固耦合问题的高效求解提供了有力支持。从离散化方式来看，IS-PIM对流体和固体分别采用了不同且高效的离散化方法。对于流体部分，基于Navier-Stokes方程，采用有限体积法等数值方法将其离散到欧拉网格上。以二维不可压缩粘性流体为例，将计算域划分为一系列的矩形或三角形网格单元。这种基于规则网格的离散方式，在计算流体的速度和压力分布时，具有较高的计算效率。有限体积法在计算流体通量时，通过对网格单元边界上的物理量进行积分，能够有效地利用网格的规则性，减少计算量。在计算流体的动量方程时，有限体积法可以快速地计算出每个网格单元上的动量变化，从而高效地更新流体的速度场。对于固体部分，基于光滑点插值法（S-PIM），将弹性力学方程离散到拉格朗日网格上。S-PIM基于容易剖分的线性背景网格，能够快速地对复杂形状的固体求解域进行离散。线性背景网格的剖分相对简单，不需要像传统有限元网格那样进行复杂的网格划分和网格质量控制。在处理复杂形状的固体结构时，S-PIM可以在较少的节点下准确地描述结构的力学行为，减少了计算所需的节点数量，从而提高了计算效率。而且，S-PIM通过在每个节点周围构建光滑域，利用该光滑域内节点的信息进行插值计算，有效地软化了固体模型刚度，使得在求解固体力学方程时，迭代收敛速度更快。在求解过程中，IS-PIM采用了分区域求解的技术。先对离散后的流体控制方程进行求解，得到流体的速度和压力分布。采用合适的迭代算法，如SIMPLE算法（Semi-ImplicitMethodforPressure-LinkedEquations）等，对离散后的代数方程组进行迭代求解。在每次迭代中，先根据上一次迭代得到的压力场计算速度场，然后根据速度场的变化更新压力场，直到速度场和压力场收敛。接着对离散后的固体控制方程进行求解，得到固体的位移和应力分布。同样采用迭代算法，如牛顿-拉夫逊迭代法等，对关于节点位移的代数方程组进行求解。这种分区域求解的方式，避免了同时求解整个流固耦合系统带来的巨大计算量，使得计算过程更加高效。IS-PIM还通过改进流固耦合力的计算方式来提高计算效率。在原始的IS-PIM中，流固耦合力的计算存在缺陷，忽略了作用在固体边界上的粘性力大小，导致计算出的固体所受流固耦合力的大小不准确。改进后的方法将粘性力基于流体欧拉网格进行计算，通过准确计算流固耦合力，并将其作为固体的外力施加到固体求解器中，实现了流体和固体之间的力的传递和相互作用。这种改进不仅提高了计算精度，还减少了不必要的计算步骤，从而提高了计算效率。在模拟低雷诺数下的圆柱绕流问题时，改进后的方法能够更准确地计算圆柱所受的流固耦合力，同时计算时间也有所缩短。3.2.3计算精度的提高IS-PIM在计算精度方面表现出色，通过多种技术手段有效地提升了对流固耦合问题的模拟准确性。从固体求解器的角度来看，光滑点插值方法（S-PIM）作为IS-PIM中的固体求解器，通过独特的节点插值和光滑域构建技术，显著提高了固体力学行为的模拟精度。在固体力学分析中，传统的有限元法存在一些局限性，如对复杂几何形状的适应性较差、低阶单元求解精度不足等。而S-PIM基于节点构建光滑域，以实现对固体力学行为的准确描述。在求解域内，离散分布着一系列节点。以二维问题为例，对于每个节点，选取其周围一定范围内的节点构成光滑域。常见的光滑域形状有三角形、四边形等。以三角形光滑域为例，在该光滑域内，基于拉格朗日插值法构建插值函数。通过对插值函数求梯度，可得到该光滑域内的梯度近似值。将各个光滑域的梯度进行组合，就可以得到整个求解域的近似梯度。这种基于光滑域的梯度计算方式，与传统有限元法直接基于单元进行梯度计算不同，它能更灵活地处理节点分布，从而软化固体模型刚度，提高了固体求解的精度。在处理复杂形状的固体结构时，S-PIM能够在较少的节点下准确地描述结构的力学行为，相比传统有限元法，在相同节点数量下，S-PIM能够提供更准确的位移和应力计算结果。在流固耦合条件的施加方面，IS-PIM通过虚拟流体实现了准确的流固耦合。虚拟流体的引入是IS-PIM的核心技术之一，它为流固耦合问题的求解提供了一种有效的途径。在流固耦合问题中，固体区域被嵌入到流体计算域中，为了实现固体与流体之间的相互作用，在固体区域内定义虚拟流体。虚拟流体与真实流体具有相似的物理性质，如密度、粘度等。通过虚拟流体，将固体区域视为流体计算域的一部分，使得在整个计算域上可以采用统一的数值方法进行求解。在模拟圆柱绕流与弹性支撑结构的流固耦合问题时，通过虚拟流体，准确地实现了流体和固体之间的速度和力的传递。在流固耦合界面，通过将固体表面节点的速度传递给与之相邻的虚拟流体节点，使得虚拟流体节点的速度与固体表面节点的速度一致，从而满足流固耦合速度条件。在计算流固耦合力时，改进后的方法将粘性力基于流体欧拉网格进行计算，避免了原始方法中忽略粘性力导致的计算误差，提高了流固耦合力的计算精度。IS-PIM还针对原始算法中存在的问题进行了修正，进一步提高了计算精度。原始IS-PIM存在流固边界不准确问题，即流体域并不能准确捕捉到固体域的边界，进而导致计算流固耦合问题时的流固边界不准确，从而影响计算精度。通过采用虚拟点修正算法（Ghost-nodetechnique）等方法对其进行修正，有效地解决了这一问题。虚拟点修正算法通过在非结构网格中标记节点，对局部速度进行修正，并利用质量源/沉算法移除多余流量，从而准确地捕捉到固体边界的作用，提高了计算精度。原始IS-PIM存在“新鲜点（freshnodes）”问题，即流体节点从固体域内部到固体域外部过程中，会存在速度以及压力振荡情况。通过采用尖锐界面方法（Sharp-interfacemethod）等对其进行修正，有效地消除了速度和压力振荡，保证了计算结果的稳定性和准确性。3.3存在的问题与挑战尽管浸没式光滑点插值法（IS-PIM）在流固耦合问题求解中展现出诸多优势，但原始算法仍存在一些亟待解决的问题，这些问题在一定程度上限制了其应用范围和计算精度。原始IS-PIM在流固耦合力求解方面存在不准确的问题。在采用原始IS-PIM求解流固耦合问题时，流固耦合力的计算是通过虚拟流体域的拉格朗日网格进行的。这种求解方式忽略了作用在固体边界上的粘性力大小，导致计算出的固体所受流固耦合力的大小不准确。在固体刚度较大或者流体雷诺数较小时，这种问题更加明显。在模拟低雷诺数下的圆柱绕流问题时，由于粘性力在流固相互作用中起着重要作用，而原始IS-PIM忽略了边界粘性力，使得计算得到的圆柱所受流固耦合力与实际情况存在较大偏差，进而影响了对圆柱运动和周围流场的准确模拟。原始IS-PIM存在流固边界不准确的问题。作为典型的浸没类方法，原始IS-PIM中流体域并不能准确捕捉到固体域的边界，进而导致计算流固耦合问题时的流固边界不准确，从而影响计算精度。在模拟弹性梁在流体中的流固耦合问题时，由于流体域对弹性梁边界的捕捉不准确，使得在流固耦合界面上，流体与固体之间的速度和力的传递存在误差，最终影响了整个流固耦合系统的计算结果。在对流固边界不准确问题进行修正后，还需要保证修正后的算法在物理层面满足连续性的要求，这增加了算法改进的复杂性。原始IS-PIM还存在“新鲜点（freshnodes）”问题。当流体节点从固体域内部移动到固体域外部过程中，会出现速度以及压力振荡情况。在模拟圆盘落水的流固耦合问题时，圆盘进入流体的过程中，流体节点从原本位于固体圆盘内部逐渐移动到外部，此时这些“新鲜点”处的速度和压力会出现剧烈振荡，使得计算结果不稳定，难以准确反映实际的流固耦合现象。这种振荡不仅影响了计算精度，还可能导致计算过程的不收敛，增加了计算的难度和时间成本。四、IS-PIM算法的改进与优化4.1针对流固耦合力计算问题的改进4.1.1传统计算方法的缺陷分析在原始的浸没式光滑点插值法（IS-PIM）中，流固耦合力的计算是通过虚拟流体域的拉格朗日网格进行的。这种计算方式在一定程度上简化了计算过程，但却存在明显的缺陷，尤其是在处理固体边界上的粘性力时。粘性力在流固耦合问题中起着重要作用，它反映了流体与固体表面之间的摩擦和相互作用。在低雷诺数流动情况下，粘性力对固体的运动和受力有着显著影响。然而，原始IS-PIM的计算方法却忽略了作用在固体边界上的粘性力大小。从理论角度分析，根据流体力学的基本原理，粘性力与流体速度的梯度密切相关。在流固耦合界面，流体速度的变化会产生速度梯度，从而导致粘性力的产生。而通过虚拟流体域拉格朗日网格计算流固耦合力时，没有充分考虑这种速度梯度的影响，使得计算出的固体所受流固耦合力中缺失了粘性力这一重要分量。在模拟低雷诺数下的圆柱绕流问题时，由于粘性力在流固相互作用中占据主导地位，忽略粘性力会导致计算得到的圆柱所受流固耦合力与实际情况存在较大偏差。这种偏差会进一步影响对圆柱运动轨迹和周围流场的准确模拟，使得模拟结果无法真实反映实际的流固耦合现象。当圆柱在低雷诺数流体中运动时，粘性力会对圆柱产生一个阻碍其运动的阻力，同时还会影响圆柱周围的流场分布。如果忽略粘性力，计算得到的圆柱所受阻力会偏小，流场分布也会与实际情况不同，从而导致对整个流固耦合系统的理解和分析出现偏差。在固体刚度较大的情况下，这种忽略粘性力的问题也会更加突出。当固体刚度较大时，固体对流体的反作用相对较小，流体的粘性力对固体的影响相对更加显著。此时，忽略粘性力会使得计算出的流固耦合力无法准确反映实际的力学作用，进而影响对整个流固耦合系统的分析和预测。在模拟刚性结构在粘性流体中的受力情况时，如果忽略粘性力，可能会低估结构所受的阻力，从而在工程设计中导致结构的安全性和可靠性受到影响。4.1.2改进的流固耦合力计算方法针对原始IS-PIM中流固耦合力计算的缺陷，提出一种将粘性力基于流体欧拉网格进行计算的改进方法。该方法的原理基于流体力学中粘性力的计算原理以及欧拉网格和拉格朗日网格的特点。在流体力学中，粘性力与流体速度的梯度相关。对于牛顿流体，粘性力的表达式可以通过牛顿内摩擦定律推导得到。在二维情况下，粘性力在x和y方向上的分量分别为：\tau_{xy}=\mu(\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx})\tau_{yx}=\mu(\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx})其中，\tau_{xy}和\tau_{yx}分别是粘性力在x和y方向上的分量，\mu是流体动力粘度，u和v分别是流体速度在x和y方向上的分量。基于流体欧拉网格进行粘性力计算，首先根据流体求解器得到的流体速度和压力分布，在流体欧拉网格上对速度进行差分计算，从而得到速度梯度。以中心差分格式为例，对于二维笛卡尔坐标系下的速度分量u，在x方向上的一阶导数（速度梯度）可以近似表示为：\frac{\partialu}{\partialx}\approx\frac{u_{i+1,j}-u_{i-1,j}}{2\Deltax}其中，u_{i+1,j}和u_{i-1,j}分别是相邻网格点在x方向上的速度值，\Deltax是x方向上的网格间距。同理，可以计算出y方向上的速度梯度。得到速度梯度后，根据上述粘性力的表达式，就可以计算出作用在固体边界上的粘性力。将计算得到的粘性力与通过虚拟流体域拉格朗日网格计算得到的惯性力等其他力分量相结合，得到准确的流固耦合力。在模拟圆柱绕流问题时，通过流体欧拉网格计算得到的粘性力与通过虚拟流体域拉格朗日网格计算得到的惯性力相加，得到圆柱所受的总流固耦合力。这种改进方法具有多方面的优势。从计算精度角度来看，通过准确计算粘性力，弥补了原始方法中忽略粘性力导致的计算误差，使得计算出的流固耦合力更加准确，从而提高了整个流固耦合模拟的精度。在模拟低雷诺数下的圆柱绕流问题时，改进后的方法能够更准确地计算圆柱所受的流固耦合力，模拟结果与实际情况更加吻合。从物理意义角度分析，该方法更符合流固耦合的物理本质。在实际的流固耦合过程中，粘性力是流固相互作用的重要组成部分，改进后的方法将粘性力纳入计算，能够更真实地反映流固之间的力学作用。通过基于流体欧拉网格计算粘性力，充分利用了欧拉网格在描述流体连续介质特性方面的优势，使得粘性力的计算更加合理和准确。4.2解决流固边界不准确问题4.2.1虚拟点修正算法（Ghost-nodealgorithm）虚拟点修正算法作为解决原始浸没式光滑点插值法（IS-PIM）中流固边界不准确问题的有效手段，通过一系列精细的操作来实现对固体边界的准确捕捉。该算法首先在非结构网格中对节点进行标记，将节点分为内部节点、边界节点和虚拟节点。内部节点位于流体域内部，其周围均为流体；边界节点位于流固耦合界面，一侧为流体，另一侧为固体；虚拟节点则位于固体域内部，是为了修正流固边界而引入的特殊节点。在模拟弹性梁在流体中的流固耦合问题时，通过对非结构网格中的节点进行标记，能够清晰地区分不同类型的节点，为后续的计算提供基础。对于边界节点，虚拟点修正算法采用局部速度修正的方法。由于在原始IS-PIM中，流体域不能准确捕捉固体域边界，导致边界节点的速度计算存在误差。通过局部速度修正，利用边界节点周围的流体信息和固体信息，对边界节点的速度进行调整。具体来说，根据边界节点与周围流体节点和固体节点的相对位置关系，采用合适的插值方法，如线性插值或样条插值，对边界节点的速度进行修正。在一个二维的流固耦合模型中，对于位于弹性梁边界的节点，通过对其周围流体节点和固体节点的速度进行线性插值，得到更准确的边界节点速度。这样可以使边界节点的速度更准确地反映流固耦合界面的真实情况，从而提高流固边界的准确性。虚拟点修正算法还利用质量源/沉算法来移除多余流量，保证算法在物理层面的连续性。在修正流固边界的过程中，可能会引入多余的流量，这会破坏流体的质量守恒。质量源/沉算法通过在虚拟节点处设置合适的质量源或沉，来调整流体的质量分布。如果在某个虚拟节点处出现了多余的流量，将该虚拟节点设置为质量沉，吸收多余的流量；反之，如果流量不足，将其设置为质量源，补充流量。通过这种方式，确保了整个计算域内的质量守恒，保证了修正后的算法在物理上的合理性和连续性。在模拟圆柱绕流问题时，通过质量源/沉算法，有效地移除了由于边界修正而产生的多余流量，使得计算结果更加准确和稳定。4.2.2质量源/沉算法（Masssource/sinkalgorithm）质量源/沉算法在解决流固边界不准确问题中起着关键作用，它通过巧妙的流量调整机制，保证了算法在物理层面的连续性，弥补了原始IS-PIM在这方面的不足。从物理原理上看，质量源/沉算法基于流体的质量守恒定律。在流固耦合问题中，当对流固边界进行修正时，如采用虚拟点修正算法等方法，可能会导致局部区域的流量发生变化，从而破坏质量守恒。质量源/沉算法通过在特定节点（如虚拟节点）处设置质量源或沉，来调整流体的质量分布，使整个计算域内的质量保持守恒。在一个简单的二维流固耦合模型中，假设在流固边界修正过程中，某个区域出现了流量的增加，这意味着该区域的质量有多余。质量源/沉算法会在该区域对应的虚拟节点处设置质量沉，将多余的质量吸收掉，以保证该区域的质量与周围区域的质量平衡。反之，如果某个区域的流量减少，即质量不足，会在相应的虚拟节点处设置质量源，向该区域补充质量。在实际应用中，质量源/沉算法与虚拟点修正算法等配合使用，能够有效地解决流固边界不准确问题。在采用虚拟点修正算法对非结构网格中的节点进行标记和局部速度修正后，利用质量源/沉算法对流量进行调整。在模拟弹性板在流体中的流固耦合问题时，首先通过虚拟点修正算法对弹性板边界的节点进行处理，修正其速度。然后，由于边界修正可能会导致局部流量的变化，采用质量源/沉算法对这些变化进行调整。通过监测计算域内各个节点的流量情况，确定哪些节点需要设置质量源或沉。根据流量的多余或不足程度，精确计算质量源或沉的强度，以确保质量守恒。这样，通过质量源/沉算法与虚拟点修正算法的协同作用，既准确地捕捉了固体边界，又保证了算法在物理层面的连续性，提高了流固耦合计算的精度和稳定性。4.3“新鲜点（freshnodes）”问题的处理针对原始IS-PIM中存在的“新鲜点（freshnodes）”问题，即流体节点从固体域内部到固体域外部过程中出现的速度以及压力振荡情况，采用尖锐界面方法进行处理。尖锐界面方法的核心在于对流体节点跨越固体边界时的物理量变化进行精确处理，以消除速度和压力的振荡，确保计算结果的稳定性和准确性。当流体节点从固体域内部移动到外部时，在流固边界处，由于流体的物理特性发生突变，传统的IS-PIM方法无法准确地描述这种变化，从而导致速度和压力的振荡。尖锐界面方法通过在流固边界上设置特殊的处理机制，来准确地捕捉这种物理量的突变。在数学处理上，尖锐界面方法基于流体的连续性方程和动量方程，对跨越边界的流体节点的速度和压力进行重新计算。根据连续性方程，在流固边界处，流体的质量通量应该保持连续。通过在边界上设置合适的质量通量条件，可以确保流体节点在跨越边界时，质量的变化是连续的，从而避免因质量不连续导致的速度和压力振荡。对于动量方程，在边界上考虑流体与固体之间的相互作用力，通过精确计算这些力对流体节点动量的影响，来准确地更新流体节点的速度和压力。从物理原理上看，尖锐界面方法类似于在流固边界上设置了一个过渡层，该过渡层能够平滑地连接流体域和固体域，使得流体节点在跨越边界时，物理量能够连续变化。在模拟圆盘落水的流固耦合问题时，当流体节点从圆盘内部移动到外部时，尖锐界面方法通过在圆盘表面设置过渡层，对流体节点的速度和压力进行调整。在过渡层内，根据流体与圆盘之间的相互作用，逐渐改变流体的速度和压力，使其从固体内部的状态平稳地过渡到流体外部的状态。通过这种方式，有效地消除了“新鲜点”处的速度和压力振荡，保证了计算结果的稳定性和准确性。尖锐界面方法还能够准确地捕捉流固边界的位置和形状变化，对于处理固体的大变形和复杂边界问题具有重要意义。在固体发生大变形时，流固边界的形状会不断变化，尖锐界面方法能够及时地跟踪这种变化，并对边界处的流体节点进行相应的处理，从而保证了流固耦合计算的精度。五、改进后算法的数值算例验证5.1二维圆柱绕流算例5.1.1模型建立与参数设置为了验证改进后浸没式光滑点插值法（IS-PIM）算法的有效性，建立二维圆柱绕流模型。在笛卡尔坐标系下，定义一个矩形计算域，其长为L_x=20D，宽为L_y=10D，其中D为圆柱的直径，取D=1m。圆柱位于计算域中心，圆心坐标为(10D,5D)。在流体参数设置方面，假设流体为不可压缩粘性流体，密度\rho=1.225kg/m^3，动力粘度\mu=1.7894\times10^{-5}Pa\cdots。通过改变来流速度U_0来调整雷诺数Re，本次算例中，重点研究Re=100和Re=200两种工况，对应Re=100时，来流速度U_0=1.45m/s；Re=200时，来流速度U_0=2.9m/s。对于计算域的边界条件，入口采用速度入口边界条件，给定来流速度U_0，方向沿x轴正方向；出口采用压力出口边界条件，表压设为0Pa；圆柱表面采用无滑移边界条件，即流体在圆柱表面的速度与圆柱表面的速度相同（此处圆柱静止，速度为0）；计算域的上下边界采用对称边界条件。在数值计算中，时间步长\Deltat的选择对计算结果的准确性和稳定性有重要影响。根据Courant-Friedrichs-Lewy（CFL）条件，\Deltat需满足一定的限制，以确保数值计算的稳定性。在本次算例中，经过多次调试和验证，对于Re=100的工况，选取\Deltat=0.001s；对于Re=200的工况，选取\Deltat=0.0005s。空间离散方面，采用有限体积法对流体控制方程进行离散，对计算域进行结构化网格划分，在圆柱周围进行网格加密，以提高对圆柱绕流复杂流场的分辨率。对于固体部分，基于光滑点插值法（S-PIM），在圆柱区域布置离散节点，构建三角形光滑域，进行固体力学方程的离散和求解。5.1.2计算结果与分析通过改进后的IS-PIM算法对二维圆柱绕流进行数值模拟，得到了不同雷诺数下的流场分布、圆柱所受的阻力系数和升力系数等结果，并与原始IS-PIM算法以及相关实验数据进行对比分析。在流场分布方面，对于Re=100的工况，改进后的IS-PIM算法模拟得到的流场中，圆柱后方形成了明显的对称涡街结构。从流线图可以清晰地看到，流体在绕过圆柱时，在圆柱的两侧边界层发生分离，形成了一对稳定的对称漩涡，且漩涡的脱落呈现出周期性。与原始IS-PIM算法模拟结果相比，改进后的算法得到的涡街结构更加清晰和稳定，漩涡的位置和形状与相关实验观测结果更为接近。原始IS-PIM算法由于流固边界不准确以及流固耦合力计算的缺陷，导致在模拟涡街时，漩涡的位置和大小存在一定偏差，涡街的稳定性也较差。对于Re=200的工况，流场变得更加复杂，圆柱后方的涡街呈现出不对称性，且漩涡的脱落频率增加。改进后的IS-PIM算法能够准确地捕捉到这种不对称涡街的形成和发展过程，流场中的速度矢量分布和压力分布也更加合理。而原始IS-PIM算法在模拟高雷诺数下的复杂流场时，误差更为明显，无法准确地描述涡街的不对称性和漩涡的脱落特性。在阻力系数和升力系数方面，将改进后的IS-PIM算法计算得到的结果与原始算法以及实验数据进行对比。以Re=100为例，经过长时间的数值模拟计算，得到改进后算法的阻力系数平均值约为Cd=1.32，升力系数的幅值约为Cl_{max}=0.35，升力系数的变化呈现出明显的周期性。原始IS-PIM算法计算得到的阻力系数平均值约为Cd=1.45，与改进后算法相比，偏差较大。通过查阅相关实验数据，在Re=100时，实验测得的阻力系数平均值约为1.30-1.35之间，改进后的IS-PIM算法计算结果与实验值更为接近，验证了改进算法在阻力系数计算方面的准确性。在升力系数方面，实验测得的升力系数幅值约为0.32-0.38之间，改进后的算法计算结果也在这个合理范围内，而原始IS-PIM算法计算得到的升力系数幅值与实验值偏差较大，无法准确反映实际情况。对于Re=200的工况，改进后算法计算得到的阻力系数平均值约为Cd=1.25，升力系数幅值约为Cl_{max}=0.42，同样与实验数据和相关文献结果更为吻合，而原始算法的计算结果偏差较大。从阻力系数和升力系数随时间的变化曲线来看，改进后的IS-PIM算法得到的曲线更加平滑和稳定，波动较小，说明改进后的算法在计算流固耦合力时更加准确，能够更真实地反映圆柱在绕流过程中的受力情况。而原始IS-PIM算法得到的曲线存在较大的波动，这是由于其在计算流固耦合力时忽略了粘性力等因素，导致计算结果不稳定。通过对不同雷诺数下的二维圆柱绕流算例的计算结果分析，充分验证了改进后的IS-PIM算法在提高计算精度和稳定性方面的有效性，能够更准确地模拟流固耦合问题中的复杂流场和流固相互作用。5.2圆盘落水算例5.2.1模型建立与参数设置为了进一步验证改进后算法的有效性，针对圆盘落水这一复杂的流固耦合问题建立数值模型。在笛卡尔坐标系下构建一个矩形计算域，长为L_x=20D，宽为L_y=15D，其中D为圆盘的直径，设定D=0.5m。圆盘初始位置位于计算域顶部，圆心坐标为(10D,14.5D)，初始速度为v_0=5m/s，方向垂直向下。在流体参数设置方面，假设流体为不可压缩粘性水，密度\rho_f=1000kg/m^3，动力粘度\mu=1.0\times10^{-3}Pa\cdots。圆盘的密度设为\rho_s=2000kg/m^3，弹性模量E=2\times10^7Pa，泊松比\nu=0.3。对于计算域的边界条件，入口采用速度入口边界条件，给定圆盘的初始速度；出口采用压力出口边界条件，表压设为0Pa；计算域的左右边界采用对称边界条件；圆盘表面采用无滑移边界条件，即流体在圆盘表面的速度与圆盘表面的速度相同。在数值计算中，时间步长\Deltat根据Courant-Friedrichs-Lewy（CFL）条件进行选取，经过多次调试，最终确定\Deltat=0.0005s。空间离散方面，采用有限体积法对流体控制方程进行离散，对计算域进行结构化网格划分，在圆盘周围进行网格加密，以提高对圆盘落水过程中复杂流场的分辨率。对于固体部分，基于光滑点插值法（S-PIM），在圆盘区域布置离散节点，构建三角形光滑域，进行固体力学方程的离散和求解。5.2.2计算结果与分析利用改进后的IS-PIM算法对圆盘落水过程进行数值模拟，得到了圆盘的运动轨迹、受力情况以及流场分布等结果，并与原始IS-PIM算法的结果进行对比分析。在圆盘的运动轨迹方面，改进后的IS-PIM算法能够准确地捕捉到圆盘在水中的运动过程。从模拟结果可以看出，圆盘在入水瞬间，由于受到流体的阻力作用，速度迅速减小。随着圆盘的下沉，流体对圆盘的作用力逐渐发生变化，圆盘开始出现一定的摆动。在圆盘下沉到一定深度后，由于流体的浮力和阻力的共同作用，圆盘的速度逐渐趋于稳定。与原始IS-PIM算法模拟结果相比，改进后的算法得到的圆盘运动轨迹更加平滑和准确。原始IS-PIM算法由于存在“新鲜点”问题，导致在圆盘入水初期，速度和压力振荡明显，使得圆盘的运动轨迹出现较大偏差。在圆盘的受力情况方面，改进后的算法能够更准确地计算出圆盘所受的流固耦合力。圆盘在落水过程中，主要受到流体的浮力、阻力和粘性力的作用。改进后的算法通过将粘性力基于流体欧拉网格进行计算，弥补了原始算法中忽略粘性力的缺陷，使得计算出的流固耦合力更加准确。在圆盘入水的初期，流体的阻力和粘性力对圆盘的运动起着主导作用，改进后的算法能够准确地反映出这些力的变化情况。随着圆盘的下沉，浮力逐渐增大，改进后的算法也能够准确地计算出浮力的大小及其对圆盘运动的影响。通过对比改进前后算法计算得到的圆盘所受合力随时间的变化曲线，可以明显看出改进后的算法得到的曲线更加平滑和稳定，波动较小，说明改进后的算法在计算流固耦合力时更加准确，能够更真实地反映圆盘在落水过程中的受力情况。在流场分布方面，改进后的IS-PIM算法模拟得到的流场中，圆盘周围的流场结构更加清晰和合理。在圆盘入水时，流体在圆盘周围形成了复杂的流动结构，包括边界层分离、漩涡的形成和发展等。改进后的算法能够准确地捕捉到这些流动现象，并且能够清晰地展示出漩涡的位置、大小和演化过程。与原始IS-PIM算法相比，原始算法由于流固边界不准确等问题，导致在模拟流场时，漩涡的位置和大小存在一定偏差，流场的稳定性也较差。通过对不同时刻的流场流线图和速度矢量图的分析，可以发现改进后的算法得到的流场分布更加符合实际物理现象，能够为进一步研究圆盘落水过程中的流固耦合机制提供更准确的依据。5.3流道中的稳态梁算例5.3.1模型建立与参数设置为了进一步验证改进后算法在处理复杂流固耦合问题时对结构变形和应力分布计算的准确性，建立流道中的稳态梁模型。在笛卡尔坐标系下构建一个