【首先假设任意区域Ω，以及其子区域Ω1和Ω2，且有Ω2=Ω/Ω1。Γ为分割Ω的界面，构造一個正则化函数Φ（通常为光滑的海维赛德函数）隐式表示Γ，使得当经过Γ时，Φ快速地从1变为0，而在子区域Ω1内，则有Φ≈1，子区域Ω2内，有Φ≈0，故而一般将Γ设为Φ的0.5水平集，在计算中，我们希望即使在有小扰动存在的情况也可以保持界面的形状……】

【……】

时间，随着常浩南双手在键盘上的飞速敲击而飞速流逝着。

不过在大约二十分钟后，常浩南就停下了打字的动作。

稍加思索之后，他抽出一张软盘，把才写了个开头的文档复制进去，然后起身径直前往了机房。

实际上，他开发的这种方法在数学上并不复杂，实质就是一种重新初始化过程，通过求解一个粘性项的稳定解来修复水平集方程中的ψ值，以达到维持守恒性的目的。

而之所以过去从未有人涉及到，是因为其在具体的计算方法上难以实现。

比如传统的网格划分方式，就很难适应笛卡尔网格上高速变化的两相流体。

但对于重生后首先就点出了数值计算专精，还全程主导了TORCHMultiphysics软件开发和算法设计的常浩南来说，这反而不是什么问题。

根据他的估计，在引入带自适应网格的有限元离散方法之后，计算效率还是能得到基本保障的。

只是如果想要通过非自主编程，而是图形界面互动的简单形式操作，那估计就只有等到TORCHMultiphysics的下一次大版本更新才行了……

不过么……

“这倒也是个不错的机会。”

正在按照论文中内容敲着代码的常浩南心想到。

TORCHMultiphysics从最初的试用版发布到现在已经过去了差不多3年半时间。

软件嘛，只要你往出卖，代码其实是无论如何都藏不住的。

哪怕常浩南搞了三个门槛不一的版本，但最多也只能在一定程度上拖延这一进程而已。

毕竟，能接触到专业提升版软件的人也有成千上万，基数大了以后，透出去风声也属于意料之中。

再者，同行竞品毕竟也不是吃干饭的。

总之从2000年左右开始，像MaterialsStudio、ANSYSFluent、COMSOLMultiphysics等软件也开始在计算效率和性能上奋起直追，尽管在综合性能上距离TORCHMultiphysics这个先行者尚有距离，但在各自擅长的领域，已经有了一较高下的资格。

尽管后者靠着过去几年的广撒网，以及火炬集团提供的数值计算服务，仍然牢牢占据着市场占有率的头把交椅，但长此以往下去，优势肯定会越来越小。

因此，火炬集团2月份提交上来的工作计划中，就包含一次预计在2001年年底或2002年年初上线的全新大版本更新。

正好可以把笛卡尔网格，以及这一类新的水平集方法给引入进去。

“是时候重新让你们感受到恐惧了……”

常浩南啪地一声按下回车键，然后靠在椅背上，看着屏幕上显示出的计算流程，自言自语道。

然而，就在这个时候，他背后不远处突然传出一个好奇的女声：

“呃……什么恐惧？”