一种沿程加密的水滴释放方法、收集率获取方法和介质-9479威尼斯

专利名称：一种沿程加密的水滴释放方法、收集率获取方法和介质

专利类型：实用新型专利

专利申请号：cn202210797065.8

专利申请（专利权）人：四川大学
权利人地址：四川省成都市一环路南一段24号

专利发明（设计）人：易贤,周志宏,曾涛,陈宇,田晓宝

专利摘要：本申请实施例公开了一种沿程加密的水滴释放方法、收集率获取方法和介质，涉及水滴撞击特性研究技术领域。该方法包括：在初始截面释放第一水滴；根据第一水滴的运动轨迹，得到第一水滴在加密截面上的第一水滴初始信息；根据第一水滴初始信息对加密截面上的水滴进行加密，并在加密截面中的各网格处释放加密后的第二水滴。根据上述方法中释放的第二水滴的运动轨迹，得到物面收集到的第二水滴的数量和位置；根据物面收集到的第二水滴的数量和位置，以及加密截面释放的第二水滴的数量和位置，采用粒子统计法得到物面的水滴收集率。通过上述方法，可以解决在获取水滴收集率时，存在的计算量大、计算步骤复杂、对于复杂构型的三维物体适应性较差等问题。

主权利要求：
1.一种基于沿程加密水滴释放方式的水滴收集率获取方法，其特征在于，所述方法包括：s110.在初始截面释放第一水滴；
s120.根据所述第一水滴的运动轨迹，得到所述第一水滴穿过加密截面时在所述加密截面上的第一水滴初始信息，其中，所述加密截面位于所述初始截面与物面之间,所述第一水滴初始信息包括所述第一水滴在所述加密截面上的位置和速度，所述第一水滴初始信息还包括所述加密截面上的各网格中穿过的所述第一水滴的数量；
s130.根据所述第一水滴初始信息对所述加密截面上的水滴进行加密，并根据所述加密截面上划分的网格，在所述加密截面中的各网格处释放加密后的第二水滴；
s140.根据释放的所述第二水滴的运动轨迹，得到所述物面收集到的所述第二水滴的数量和位置；
s150.根据所述物面收集到的所述第二水滴的数量和位置，以及所述加密截面释放的第二水滴的数量和位置，采用粒子统计法得到所述物面的水滴收集率；
所述步骤s120包括：根据所述第一水滴在所述加密截面上的位置，得到所述加密截面上的各网格中穿过的所述第一水滴的数量；
所述步骤s130包括：根据所述加密截面上的各网格中穿过的所述第一水滴的数量，得到对所述加密截面中的水滴进行加密时所述加密截面上的各网格中的第二水滴的数量。
2.根据权利要求1所述的基于沿程加密水滴释放方式的水滴收集率获取方法，其特征在于，所述步骤s130包括：将所述加密截面上的各网格中穿过的所述第一水滴的数量乘加密系数，作为对所述加密截面中的水滴进行加密时所述加密截面上的各网格中的第二水滴的数量。
3.根据权利要求1所述的基于沿程加密水滴释放方式的水滴收集率获取方法，其特征在于，所述步骤s150包括：根据所述第一水滴穿过所述加密截面时在所述加密截面上的位置，将包含所述位置的网格的总和标记为所述加密截面的释放区域；
根据所述物面收集到的所述第二水滴的位置，得到所述物面的收集区域；
根据所述物面收集到的所述第二水滴的数量和所述物面的收集区域，以及所述加密截面释放的第二水滴的数量和所述加密截面的释放区域，采用粒子统计法得到所述物面的水滴收集率。
4.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有程序代码，所述程序代码可被处理器调用执行所述权利要求1 3任一项所述方法。
~ 说明书 : 一种沿程加密的水滴释放方法、收集率获取方法和介质技术领域[0001] 本申请涉及水滴撞击特性研究技术领域，更具体地，涉及一种沿程加密的水滴释放方法、收集率获取方法和介质。背景技术[0002] 飞机穿越云层时，过冷水滴撞击到机体后，可能发生相变并导致结冰现象。结冰会改变飞机的外形与绕流流场，破坏气动性能，降低操纵性与稳定性，威胁飞行安全，严重时导致空难事故。水滴收集率可以定量表征撞击到物体上的水滴，是研究水滴撞击特性的重要参数。因此，水滴收集率可以用于研究飞机结冰特性，还可以用于飞机防除冰系统设计。[0003] 目前，可以采用数值模拟方法研究水滴收集率，其中，拉格朗日方法是研究水滴收集率的常用方法之一。利用拉格朗日方法，获得远场平面释放的水滴的运动轨迹，从而得到水滴撞击物面的情况。之后采用面积比法或粒子统计法得到水滴收集率。[0004] 粒子统计法根据远场平面释放的水滴的数量和位置，得到远场平面释放的水滴的浓度，根据物面收集的水滴的数量和位置，得到物面收集的水滴的浓度，再根据远场平面释放的水滴的浓度与物面收集的水滴的浓度，得到物面的水滴收集率。[0005] 现有技术在利用拉格朗日方法结合粒子统计法，在远场平面释放水滴并获取水滴收集率时，为保证物面的水滴收集率的获取精度，需要在远场平面处加密水滴，并且，从远场平面释放加密后的水滴，到加密后的水滴撞击到物面的过程中，需要持续获取加密后的水滴的运动轨迹，计算量很大。[0006] 面积比法根据根据远场平面释放的水滴的运动轨迹，得到远场平面释放区域的面积与物面撞击区域的面积，再根据远场平面释放区域的面积与物面撞击区域的面积，得到物面的水滴收集率。[0007] 现有技术在利用拉格朗日方法结合面积比法计算水滴收集率时，存在计算步骤复杂、适应性差等问题。如面积比法需要判断加密水滴撞击物面的撞击极限，在撞击区域内，还需要判断加密水滴分布形状，以及，需要在物面使用插值方法计算加密水滴的分布，计算步骤复杂。此外，面积比法适用于简单三维外形，对于复杂构型的三维物体适应性较差。[0008] 因此，现有技术在获取水滴收集率时，存在计算量大、计算步骤复杂、对于复杂构型的三维物体适应性较差等问题。发明内容[0009] 本申请提出了一种沿程加密的水滴释放方法，在初始截面释放第一水滴，第一水滴的数量较少，因此计算量较小，根据第一水滴的运动轨迹，得到所述第一水滴穿过初始截面与物面之间的加密截面时，在加密截面上的第一水滴初始信息，第一水滴初始信息包括第一水滴在加密截面上的位置和速度；根据第一水滴初始信息对加密截面上的水滴进行加密，第一水滴初始信息用于获得加密后的水滴的运动轨迹，根据加密截面上划分的网格，在加密截面中的各网格处释放加密后的第二水滴。基于在离物面较远的初始截面处空气流场变化并不剧烈的特点，在初始截面减少第一水滴释放的数量，在离物面较近距离的加密截面处再释放加密后的第二水滴，后续获取水滴收集率时计算量更小。本申请还提出了一种基于沿程加密水滴释放方式的水滴收集率获取方法，根据上述方法中的第二水滴的运动轨迹，得到物面收集到的所述第二水滴的数量和位置；并根据物面收集到的所述第二水滴的数量和位置，以及加密截面释放的第二水滴的数量和位置，采用粒子统计法得到所述物面的水滴收集率。由于加密截面离物面比初始截面离物面更近，因此在加密截面释放加密后的水滴，比在初始截面释放加密后的水滴，需要计算的水滴轨迹更短，计算量更小，且本申请的方法无需像面积比法一样计算加密水滴的撞击极限、加密水滴分布形状，以及无需在物面使用插值方法计算加密水滴的分布。如此，可以有效解决现有技术在获取水滴收集率时，存在的计算量大、计算步骤复杂、对于复杂构型的三维物体适应性较差等问题。[0010] 第一方面，本申请实施例提供了一种沿程加密的水滴释放方法，该方法包括：在初始截面释放第一水滴；根据所述第一水滴的运动轨迹，得到所述第一水滴穿过加密截面时在所述加密截面上的第一水滴初始信息，其中，所述加密截面位于所述初始截面与物面之间,所述第一水滴初始信息包括所述第一水滴在所述加密截面上的位置和速度；根据所述第一水滴初始信息对所述加密截面上的水滴进行加密，并根据所述加密截面上划分的网格，在所述加密截面中的各网格处释放加密后的第二水滴。[0011] 第二方面，本申请实施例提供了一种基于沿程加密水滴释放方式的水滴收集率获取方法，该方法包括：根据上述方法中释放的所述第二水滴的运动轨迹，得到所述物面收集到的所述第二水滴的数量和位置；根据所述物面收集到的所述第二水滴的数量和位置，以及所述加密截面释放的第二水滴的数量和位置，采用粒子统计法得到所述物面的水滴收集率。[0012] 第三方面，本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质中存储有程序代码，所述程序代码可被处理器调用执行上述方法。[0013] 综上所述，本申请至少具有如下技术效果：[0014] 1.本申请通过在初始截面释放第一水滴，且第一水滴的数量较少，计算量也较小，根据第一水滴的运动轨迹，得到第一水滴初始信息，并根据第一水滴初始信息在加密截面上加密得到第二水滴，在离物面更近的加密截面释放第二水滴，为获取物面的水滴收集率建立较为简单的前提条件。[0015] 2.本申请通过根据加密截面上的各网格中穿过的第一水滴的数量，得到对加密截面中的水滴进行加密时加密截面上的各网格中的第二水滴的数量，确保获取的水滴收集率的精度。[0016] 3.本申请利用加密截面上的释放的第二水滴的运动轨迹获取物面的水滴收集率，由于加密截面离物面比初始截面离物面更近，因此在加密截面释放加密后的水滴，比在初始截面释放加密后的水滴，需要计算的水滴轨迹更短，计算量更小。[0017] 4.使用本申请提供的沿程加密的水滴释放方法，以及基于沿程加密水滴释放方式的水滴收集率获取方法获取水滴收集率，无需判断加密水滴撞击物面的撞击极限和撞击区域内的加密水滴分布形状，也无需在物面使用插值方法计算加密水滴的分布，计算步骤更简单。同时，使用本申请提供的沿程加密的水滴释放方法，以及基于沿程加密水滴释放方式的水滴收集率获取方法获取水滴收集率，可以适用于复杂构型的三维物体的水滴收集率计算，适应性更好。[0018] 因此，本申请提供的方案可以有效解决现有技术在获取水滴收集率时，存在的计算量大、计算步骤复杂、对于复杂构型的三维物体适应性较差等问题。附图说明[0019] 为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。[0020] 图1示出了本申请实施例1提供的远场平面释放水滴的示意图；[0021] 图2示出了本申请实施例1提供的一种沿程加密的水滴释放方法的流程示意图；[0022] 图3示出了本申请实施例1提供的初始截面的示意图；[0023] 图4示出了本申请实施例1提供的第一水滴运动轨迹的示意图；[0024] 图5示出了本申请实施例1提供第一水滴在初始截面的速度和第一水滴在加密截面的速度的示意图；[0025] 图6示出了本申请实施例1提供的第二水滴运动轨迹的示意图；[0026] 图7示出了本申请实施例2提供的加密截面释放的第二水滴的浓度的示意图；[0027] 图8示出了本申请实施例3提供的用于保存或者携带实现本申请实施例的沿程加密的水滴释放方法和基于沿程加密水滴释放方式的水滴收集率获取方法的程序代码的存储单元。具体实施方式[0028] 为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。[0029] 目前，拉格朗日方法是研究水滴收集率的常用方法之一。可以利用拉格朗日方法结合粒子统计法或面积比法获取物面的水滴收集率。如图1所示，粒子统计法根据远场平面释放的水滴的数量和位置，得到远场平面释放的水滴的浓度 ，根据物面收集的水滴的数量和位置，得到物面收集的水滴的浓度 ，再根据远场平面释放的水滴的浓度 与物面收集的水滴的浓度 ，得到物面的水滴收集率。现有技术在利用拉格朗日方法结合粒子统计法，在远场平面释放水滴并获取水滴收集率时，为保证物面的水滴收集率的获取精度，需要在远场平面处加密水滴，并且，从远场平面释放加密后的水滴，到加密后的水滴撞击到物面的过程中，需要持续获取加密后的水滴的运动轨迹，计算量很大。而在利用拉格朗日方法结合面积比法获取水滴收集率时，存在计算步骤复杂、适应性差等问题。如面积比法需要判断加密水滴撞击物面的撞击极限，在撞击区域内，还需要判断加密水滴分布形状，以及，需要在物面使用插值方法计算加密水滴的分布，计算步骤复杂。此外，面积比法适用于简单三维外形，对于复杂构型的三维物体适应性较差。[0030] 因此，为了解决上述缺陷，本申请实施例提供了一种沿程加密的水滴释放方法，该方法包括：在初始截面释放第一水滴；根据所述第一水滴的运动轨迹，得到所述第一水滴穿过加密截面时在所述加密截面上的第一水滴初始信息，其中，所述加密截面位于所述初始截面与物面之间,所述第一水滴初始信息包括所述第一水滴在所述加密截面上的位置和速度；根据所述第一水滴初始信息对所述加密截面上的水滴进行加密，并根据所述加密截面上划分的网格，在所述加密截面中的各网格处释放加密后的第二水滴。在采用粒子统计法计算水滴收集率时，由于加密截面离物面比初始截面离物面更近，因此在加密截面释放加密后的水滴，比在初始截面释放加密后的水滴，需要计算的水滴轨迹更短，计算量更小，且本申请的方法无需像面积比法一样计算加密水滴的撞击极限、加密水滴分布形状，以及无需在物面使用插值方法计算加密水滴的分布。从而有效解决现有技术在获取水滴收集率时，存在的计算量大、计算步骤复杂、对于复杂构型的三维物体适应性较差等问题。[0031] 下面对本申请所涉及到的沿程加密的水滴释放方法进行介绍。[0032] 实施例1[0033] 请参照图2，图2为本申请实施例1提供的一种沿程加密的水滴释放方法的流程示意图。本实施例中，该沿程加密的水滴释放方法可以包括以下步骤：[0034] 步骤s110：在初始截面释放第一水滴。[0035] 在本申请实施例中，初始截面中释放第一水滴的区域可以覆盖物面待检测部分的全部区域。具体地，如图3所示，该区域可以是较大的矩形区域，如图3中示出的标记为初始截面的矩形区域，其中释放的第一水滴可以完全覆盖物面待检测部分的全部区域；该区域也可以是如图3中示出的初始截面矩形框中的阴影区域，其中释放的第一水滴可以刚好覆盖物面待检测部分的全部区域。[0036] 本申请提供的沿程加密的水滴释放方法，以及基于沿程加密水滴释放方式的水滴收集率获取方法，适用于飞机上水滴收集率的计算，也适用于风力发电机等其他设备上的水滴收集率的计算。[0037] 以获取飞行中的飞机上的水滴收集率为例进行说明，飞机在飞行时，其附近的空气流场变化较为剧烈，水滴的运动轨迹也较为复杂，但是在离飞机较远的初始截面处，空气流场的变化并不剧烈，水滴的运动轨迹较为简单。[0038] 步骤s120：根据所述第一水滴的运动轨迹，得到所述第一水滴穿过加密截面时在所述加密截面上的第一水滴初始信息，其中，所述加密截面位于所述初始截面与物面之间,所述第一水滴初始信息包括所述第一水滴在所述加密截面上的位置和速度。[0039] 在本申请实施例中，第一水滴的运动轨迹可以根据该第一水滴的受力情况得到。具体地，根据第一水滴受到的重力、浮力和阻力，可以得到：,其中， 是空气流场的密度， 是第一水滴的密度， 是第一水滴的体积， 是第一水滴的位移， 是空气流场的速度， 是第一水滴的速度， 是第一水滴迎风面积，是重力加速度， 是第一水滴阻力系数，且， 是相对雷诺数， ，且 ， 是第一水滴当量直径， 是空气粘性系数。根据上述公式可以得到： ，其中， 为第一水滴的加速度。根据第一水滴的加速度，可以得到第一水滴的速度和位移，从而得到第一水滴的运动轨迹。[0040] 在本申请实施例中，加密截面位于所述初始截面与物面之间,具体地，加密截面与初始截面之间的距离，可以根据经验得到，例如，加密截面与初始截面之间的距离为物面与初始截面之间距离的k倍，且 ，在满足计算精度和较小的计算量的条件下，k可以取0.8，k也可以取其他值，本申请对k的大小不做限制。[0041] 如图4所示，在 点处释放一个第一水滴，以 为起点的虚线表示该第一水滴的运动轨迹，根据该第一水滴的运动轨迹，可以得到该第一水滴穿过加密截面时在加密截面上的位置 点，还可以得到该第一水滴运动到 点时的速度。[0042] 在 点处释放一个第一水滴，以 为起点的虚线表示该第一水滴的运动轨迹，根据该第一水滴的运动轨迹，可以得到该第一水滴穿过加密截面时在加密截面上的位置点，还可以得到该第一水滴运动到 点时的速度。[0043] 如此对初始截面中的每个第一水滴获取运动轨迹，并得到每个第一水滴穿过加密截面时在加密截面上的位置和速度。[0044] 在示例性实施例中，第一水滴初始信息还包括所述加密截面上的各网格中穿过的所述第一水滴的数量，步骤s120还包括子步骤s121。[0045] 子步骤s121：根据所述第一水滴在所述加密截面上的位置，得到所述加密截面上的各网格中穿过的所述第一水滴的数量。[0046] 在本申请实施例中，加密截面上的网格的个数、大小以及划分方式可以根据所需精度进行划分，本申请对此不做限制。网格的形状可以是正方形、长方形、平行四边形或其他形状。[0047] 作为一种可选实施方式，本申请以网格是长方形为例进行说明。[0048] 如图4所示，将空气流场的运动方向作为x轴正方向，初始截面与x轴的交点为原点，将y轴方向和z轴方向作为加密截面上的网格的边的方向。图4示出了加密截面上的25个网格，具体的，网格的排列方式为5行5列。[0049] 在初始截面上的 点处释放一个第一水滴，该第一水滴穿过加密截面时在加密截面上的位置为 点， 点位于第2行第2列的网格内，因此第2行第2列的网格中穿过一个第一水滴。[0050] 在初始截面上的 点处释放一个第一水滴，该第一水滴穿过加密截面时在加密截面上的位置为 点， 点位于第2行第3列的网格内，因此第2行第3列的网格中穿过一个第一水滴。[0051] 如此对每个第一水滴在加密截面上的位置进行判断，并得到加密截面上的每个网格中穿过的第一水滴的数量。[0052] 步骤s130：根据所述第一水滴初始信息对所述加密截面上的水滴进行加密，并根据所述加密截面上划分的网格，在所述加密截面中的各网格处释放加密后的第二水滴。[0053] 在本申请实施例中，步骤s130中的网格的个数、大小、划分方式以及形状可以参考步骤s121中的网格的个数、大小、划分方式以及形状，本申请对此不再赘述。[0054] 在本申请实施例中，对加密截面上的水滴进行加密，可以是采用均匀分布的方式对水滴进行加密，也可以是采用随机分布的方式对水滴进行加密。作为一种可选实施方式，基于计算的便捷性，以及真实空气流场的随机性，本申请可以采用基于蒙特卡洛方法的随机释放方式对水滴进行加密。[0055] 在本申请实施例中，根据所述第一水滴初始信息对所述加密截面上的水滴进行加密，可以指：根据第一水滴在加密截面上的位置，在该第一水滴附近对之后要释放的第二水滴进行加密，并且将该第一水滴在加密截面上的速度作为附近的加密后的第二水滴的初始速度。[0056] 具体地，如图5所示，在初始截面 释放第一水滴，且第一水滴的初始速度为 ，该第一水滴撞击到物面的过程中，会穿过加密截面 ，根据第一水滴的加速度，可以得到第一水滴在穿过加密截面 时的速度为 。根据第一水滴在加密截面 上的位置，在该第一水滴附近对之后要释放的第二水滴进行加密，并且附近的加密后的第二水滴的初始速度为 。[0057] 在本申请实施例中，根据第一水滴在加密截面上的位置，还可以利用插值法得到该位置上的空气流场的速度。具体地，预测点 处空气流场的速度可以用反距离权重插值函数得到： ，其中， 为样本点 处的速度预测值，m为用于预测的空气流场的样本点总数量， 为各样本点的权重系数，且 ，，p为指数值， 为预测点 与各样本点 之间的距离。[0058] 在本申请实施例中，根据第一水滴在加密截面上的位置附近的第二水滴的初始速度，以及该位置附近的空气流场的速度，可以得到第二水滴的加速度。根据第二水滴的加速度，可以得到第二水滴的速度和位移，从而得到第二水滴的运动轨迹，为后续获取水滴收集率建立基础。[0059] 在示例性实施例中，第一水滴初始信息还包括所述加密截面上的各网格中穿过的所述第一水滴的数量，步骤s130还包括子步骤s131。[0060] 子步骤s131：根据所述加密截面上的各网格中穿过的所述第一水滴的数量，得到对所述加密截面中的水滴进行加密时所述加密截面上的各网格中的第二水滴的数量。[0061] 在本申请实施例中，加密截面上的网格中的第二水滴的数量，可以是该网格中穿过的第一水滴的数量乘加密系数，也可以是该网格中穿过的第一水滴的数量的预设次方，还可以是以该网格中穿过的第一水滴的数量为自变量的预设函数值。[0062] 作为一种可选实施方式，将所述加密截面上的各网格中穿过的所述第一水滴的数量乘加密系数，作为对所述加密截面中的水滴进行加密时所述加密截面上的各网格中的第二水滴的数量。[0063] 具体地，加密系数可以是10。如图4所示，第2行第2列的网格中穿过1个第一水滴，第2行第2列的网格中释放的第二水滴的数量为 个，即10个。若第2行第2列的网格中穿过3个第一水滴，则第2行第2列的网格中释放的第二水滴的数量为 个，即30个。加密系数可以根据水滴收集率的计算精度确定，本申请对加密系数的取值不做限制。[0064] 作为另一种可选实施方式，将所述加密截面上的各网格中穿过的所述第一水滴的数量的预设次方，作为对所述加密截面中的水滴进行加密时所述加密截面上的各网格中的第二水滴的数量。[0065] 具体地，预设次方可以是3。如图4所示，第2行第2列的网格中穿过1个第一水滴，第2行第2列的网格中释放的第二水滴的数量为 个，即1个。若第2行第2列的网格中穿过3个第一水滴，则第2行第2列的网格中释放的第二水滴的数量为 个，即27个。预设次方可以根据水滴收集率的计算精度确定，本申请对预设次方的取值不做限制。[0066] 如图6所示，图6示出了根据本申请提供的沿程加密的水滴释放方法得到的第一水滴的运动轨迹和第二水滴的运动轨迹。[0067] 本申请通过根据加密截面上的各网格中穿过的第一水滴的数量，得到对加密截面中的水滴进行加密时加密截面上的各网格中的第二水滴的数量，确保获取的水滴收集率的精度。[0068] 本申请通过在初始截面释放第一水滴，且第一水滴的数量较少，计算量也较小，根据第一水滴的运动轨迹，得到第一水滴初始信息，并根据第一水滴初始信息在加密截面上加密得到第二水滴，在离物面更近的加密截面释放第二水滴，为获取物面的水滴收集率建立较为简单的前提条件。[0069] 实施例2[0070] 本申请实施例2提供了一种基于沿程加密水滴释放方式的水滴收集率获取方法。在实施例1的基础上，实施例2还包括步骤s140和步骤s150。[0071] 步骤s140：根据实施例1所述方法中释放的所述第二水滴的运动轨迹，得到所述物面收集到的所述第二水滴的数量和位置。[0072] 在本申请实施例中，第二水滴的运动轨迹的获取方式可以参照第一水滴的运动轨迹的获取方式，本申请对此不再赘述。[0073] 步骤s150：根据所述物面收集到的所述第二水滴的数量和位置，以及所述加密截面释放的第二水滴的数量和位置，采用粒子统计法得到所述物面的水滴收集率。[0074] 作为一种可选实施方式，步骤s150还可以包括子步骤s151至子步骤s153。[0075] 子步骤s151：根据实施例1所述的方法中所述第一水滴穿过所述加密截面时在所述加密截面上的位置，将包含所述位置的网格标记为所述加密截面的释放区域。[0076] 子步骤s152：根据所述物面收集到的所述第二水滴的位置，得到所述物面的收集区域。[0077] 在本申请实施例中，可以在物面上划分网格，并根据第二水滴撞击到物面上的位置，得到物面上收集到第二水滴的网格，并将这些网格的总和标记为物面的收集区域。[0078] 子步骤s153：根据所述物面收集到的所述第二水滴的数量和所述物面的收集区域，以及所述加密截面释放的第二水滴的数量和所述加密截面的释放区域，采用粒子统计法得到所述物面的水滴收集率。[0079] 在本申请实施例中，采用粒子统计法得到所述物面的水滴收集率，可以是：如图7所示，左侧的矩形为加密截面，右侧的不规则形状为物面，加密截面中的黑色点为释放的第二水滴，物面中的黑色点为收集的第二水滴。物面上的水滴收集率为 ， 为加密截面释放的第二水滴的浓度， 为物面收集的第二水滴的浓度，且 ，，因此可以得到： ，其中， 为加密截面释放的第二水滴的个数， 为第二水滴的质量， 为加密截面的释放区域的面积， 为物面收集的第二水滴的个数， 为物面的收集区域的面积。[0080] 在获取水滴收集率时，需要获取 个第二水滴的运动轨迹，再判断 个第二水滴中会撞击到物面上的第二水滴的个数 ，以及该 个第二水滴撞击到物面上的位置。[0081] 现有技术中，从初始截面释放加密后的水滴，到加密后的水滴撞击到物面的过程中，需要持续获取加密后的水滴的运动轨迹，例如，在图4中，需要持续获取水滴从 点到物面间的运动轨迹。[0082] 根据本申请提供的基于沿程加密水滴释放方式的水滴收集率获取方法，需要获取的水滴运动轨迹，从初始截面到物面间的运动轨迹，缩减到加密截面到物面间的运动轨迹，例如，在图4中，只需要获取水滴从 点到物面间的运动轨迹。[0083] 本申请利用加密截面上的释放的第二水滴的运动轨迹获取物面的水滴收集率，由于加密截面离物面比初始截面离物面更近，因此在加密截面释放加密后的水滴，比在初始截面释放加密后的水滴，需要计算的水滴轨迹更短，计算量更小。[0084] 使用本申请提供的沿程加密的水滴释放方法，以及基于沿程加密水滴释放方式的水滴收集率获取方法获取水滴收集率，无需判断加密水滴撞击物面的撞击极限和撞击区域内的加密水滴分布形状，也无需在物面使用插值方法计算加密水滴的分布，计算步骤更简单。同时，使用本申请提供的沿程加密的水滴释放方法，以及基于沿程加密水滴释放方式的水滴收集率获取方法获取水滴收集率，可以适用于复杂构型的三维物体的水滴收集率计算，适应性更好。[0085] 实施例3[0086] 请参考图8，图8示出了本申请实施例3提供的一种计算机可读存储介质的结构框图。该计算机可读存储介质800中存储有程序代码，所述程序代码可被处理器调用执行上述方法实施例中所描述的方法。[0087] 计算机可读存储介质800可以是诸如闪存、eeprom（电可擦除可编程只读存储器）、eprom（可擦除可编程只读存储器）、硬盘或者rom之类的电子存储器。可选地，计算机可读存储介质800包括非易失性计算机可读存储介质（non‑transitorycomputer‑readablestoragemedium）。计算机可读存储介质800具有执行上述方法中的任何方法步骤的程序代码810的存储空间。这些程序代码可以从一个或者多个计算机程序产品中读取或者写入到这一个或者多个计算机程序产品中。程序代码810可以例如以适当形式进行压缩。[0088] 最后应说明的是：以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不驱使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。

专利地区：四川

专利申请日期：2022-07-08

专利公开日期：2024-07-09

专利公告号：cn115098825b