IP Library Granted Patent US 12,547,798
Granted Patent B2
US 12,547,798 · App. 17/574,036 · Granted Feb 10, 2026

Machine learning approach to piecewise linear interface construction

Inventors: Markus Bussmann (Mississauga, CA); Mohammadmehdi Ataei (Toronto, CA)
Assignee: Autodesk, Inc.
G06F30/23G06F30/28G06N3/08G06F30/27G06F2111/10G06T17/20
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,547,798
App. No.
17/574,036
Granted
Feb 10, 2026
Kind
B2
Abstract

Methods, systems, and apparatus, including medium-encoded computer program products, for a machine learning approach for piecewise linear interface construction include: obtaining a cell fraction for a mesh cell in a mesh, a normal vector perpendicular to a linear interface for the mesh cell, and geometry information of the mesh cell, wherein the geometry information includes at least two data values; normalizing the geometry information of the mesh cell to reduce the at least two data values to at least one data value; obtaining an output value from a machine learning algorithm, previously trained using normalized geometry information of multiple different cells of a same mesh type as the mesh cell; determining the linear interface for the mesh cell based on the output value, the normal vector and the geometry information; using the linear interface for the mesh cell in numerical simulation processing of the mesh.

Claims (60)

1 . A method comprising:

obtaining, in a computer program, properties of interface cells from a set of mesh cells in a computer model;

providing, by the computer program, the properties of the interface cells to a machine learning algorithm, wherein the properties of the interface cells comprise a normal vector and a cell fraction for each of the interface cells, wherein the machine learning algorithm was previously trained using multiple different cells of a same mesh type as the interface cells from the set of mesh cells;

determining, by the computer program, a linear interface, which divides the interface cells, using output values from the machine learning algorithm; and

using, by the computer program, the linear interface during processing of the computer model.

2 . The method of claim 1 , wherein the same mesh type is a regular mesh cell type, and the properties of the interface cells consist of the normal vector and the cell fraction for each of the interface cells.

3 . The method of claim 1 , wherein the same mesh type is an irregular mesh cell type, and the properties of the interface cells comprise at least one data value representing normalized geometry information for each of the interface cells.

4 . The method of claim 3 , comprising: normalizing, by the computer program, initial geometry information for each of the interface cells to reduce a number of data values representing the normalized geometry information in the properties of the interface cells.

5 . The method of claim 4 , wherein each of the interface cells has three or more vertices, and

wherein the normalizing comprises: mapping the initial geometry information of each interface cell to a new coordinate system; and using coordinate information of one or more vertices of the normalized geometry information for the interface cell in the new coordinate system, wherein a number of the one or more vertices in the new coordinate system is less than a number of the three or more vertices; and

wherein the determining comprises: obtaining final output values by remapping the output values back to a coordinate system for the interface cells; and determining the linear interface using the final output values.

6 . The method of claim 5 , wherein the same mesh type is triangular, each interface cell is a triangular interface cell having three vertices, and the one or more vertices of the normalized geometry information for the interface cell is one vertex.

7 . The method of claim 1 , further comprising training the machine learning algorithm, wherein training the machine learning algorithm comprises:

generating multiple different cells of the same mesh type as the interface cells as training inputs;

obtaining, for each of the training inputs, a respective target output for the training input; and

training the machine learning algorithm on the training inputs and the corresponding target outputs.

8 . The method of claim 1 , wherein using the linear interface during the processing of the computer model comprises using the linear interface for the interface cells in numerical simulation processing of the computer model for computational fluid dynamics simulation, plastic injection molding simulation, or topology or geometry optimization.

9 . The method of claim 1 , wherein using the linear interface during the processing of the computer model comprises using the linear interface for the interface cells to extract a polyhedral mesh surface from a three-dimensional iso-surface defined in the computer model for use in mesh-based calculation or visualization.

10 . A system comprising:

a first computer communicatively coupled with a network, the first computer comprising a first non-transitory storage medium having instructions of a computer program stored thereon; and

a second computer communicatively coupled with the network, the second computer comprising a second non-transitory storage medium having instructions of the computer program stored thereon;

wherein the computer program is configured to cause the first and second computers to

obtain properties of interface cells from a set of mesh cells in a computer model,

provide the properties of the interface cells to a machine learning algorithm, wherein the properties of the interface cells comprise a normal vector and a cell fraction for each of the interface cells, wherein the machine learning algorithm was previously trained using multiple different cells of a same mesh type as the interface cells from the set of mesh cells,

determine a linear interface, which divides the interface cells, using output values from the machine learning algorithm, and

use the linear interface during processing of the computer model.

11 . The system of claim 10 , wherein the same mesh type is a regular mesh cell type, and the properties of the interface cells consist of the normal vector and the cell fraction for each of the interface cells.

12 . The system of claim 10 , wherein the same mesh type is an irregular mesh cell type, and the properties of the interface cells comprise at least one data value representing normalized geometry information for each of the interface cells.

13 . The system of claim 12 , wherein the computer program is configured to cause the first and second computers to normalize initial geometry information for each of the interface cells to reduce a number of data values representing the normalized geometry information in the properties of the interface cells.

14 . The system of claim 13 , wherein each of the interface cells has three or more vertices, and

wherein the computer program is configured to cause the first and second computers to

map the initial geometry information of each interface cell to a new coordinate system,

use coordinate information of one or more vertices of the normalized geometry information for the interface cell in the new coordinate system, wherein a number of the one or more vertices in the new coordinate system is less than a number of the three or more vertices,

obtain final output values by remapping the output values back to a coordinate system for the interface cells, and

determine the linear interface using the final output values.

15 . The system of claim 14 , wherein the same mesh type is triangular, each interface cell is a triangular interface cell having three vertices, and the one or more vertices of the normalized geometry information for the interface cell is one vertex.

16 . The system of claim 10 , wherein the computer program or another computer program is configured to cause the first and second computers to train the machine learning algorithm, comprising:

generate multiple different cells of the same mesh type as the interface cells as training inputs,

obtain, for each of the training inputs, a respective target output for the training input, and

train the machine learning algorithm on the training inputs and the corresponding target outputs.

17 . The system of claim 10 , wherein the computer program is configured to cause the first and second computers to use the linear interface for the interface cells in numerical simulation processing of the computer model for computational fluid dynamics simulation, plastic injection molding simulation, or topology or geometry optimization.

18 . The system of claim 10 , wherein the computer program is configured to cause the first and second computers to use the linear interface for the interface cells to extract a polyhedral mesh surface from a three-dimensional iso-surface defined in the computer model for use in mesh-based calculation or visualization.

19 . A non-transitory computer-readable medium encoding instructions operable to cause data processing apparatus to perform operations comprising:

obtaining, in a computer program, properties of interface cells from a set of mesh cells in a computer model;

providing, by the computer program, the properties of the interface cells to a machine learning algorithm, wherein the properties of the interface cells comprise a normal vector and a cell fraction for each of the interface cells, wherein the machine learning algorithm was previously trained using multiple different cells of a same mesh type as the interface cells from the set of mesh cells;

determining, by the computer program, a linear interface, which divides the interface cells, using output values from the machine learning algorithm; and

using, by the computer program, the linear interface during processing of the computer model.

20 . The non-transitory computer-readable medium of claim 19 , wherein the same mesh type is an irregular mesh cell type, and the properties of the interface cells comprise at least one data value representing normalized geometry information for each of the interface cells.

21 . The non-transitory computer-readable medium of claim 20 , the operations comprising: normalizing, by the computer program, initial geometry information for each of the interface cells to reduce a number of data values representing the normalized geometry information in the properties of the interface cells.

22 . The non-transitory computer-readable medium of claim 21 , wherein each of the interface cells has three or more vertices, and

wherein the normalizing comprises: mapping the initial geometry information of each interface cell to a new coordinate system; and using coordinate information of one or more vertices of the normalized geometry information for the interface cell in the new coordinate system, wherein a number of the one or more vertices in the new coordinate system is less than a number of the three or more vertices; and

wherein the determining comprises: obtaining final output values by remapping the output values back to a coordinate system for the interface cells; and determining the linear interface using the final output values.

23 . The non-transitory computer-readable medium of claim 22 , wherein the same mesh type is triangular, each interface cell is a triangular interface cell having three vertices, and the one or more vertices of the normalized geometry information for the interface cell is one vertex.

24 . The non-transitory computer-readable medium of claim 19 , wherein the same mesh type is a regular mesh cell type, and the properties of the interface cells consist of the normal vector and the cell fraction for each of the interface cells.

25 . The non-transitory computer-readable medium of claim 19 , the operations further comprising training the machine learning algorithm, wherein training the machine learning algorithm comprises:

generating multiple different cells of the same mesh type as the interface cells as training inputs;

obtaining, for each of the training inputs, a respective target output for the training input; and

training the machine learning algorithm on the training inputs and the corresponding target outputs.

26 . The non-transitory computer-readable medium of claim 19 , wherein using the linear interface during the processing of the computer model comprises using the linear interface for the interface cells in numerical simulation processing of the computer model for computational fluid dynamics simulation, plastic injection molding simulation, or topology or geometry optimization.

27 . The non-transitory computer-readable medium of claim 19 , wherein using the linear interface during the processing of the computer model comprises using the linear interface for the interface cells to extract a polyhedral mesh surface from a three-dimensional iso-surface defined in the computer model for use in mesh-based calculation or visualization.

Assignments (2)
CHANGE OF ADDRESS FOR ASSIGNEE Recorded Aug 19, 2022
From: AUTODESK, INC.
To: AUTODESK, INC.
Reel/Frame 061572/0061 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Feb 3, 2022
From: BUSSMANN, MARKUS; ATAEI, MOHAMMADMEHDI
To: AUTODESK, INC.
Reel/Frame 058883/0212 →
Continuity (2)
Continuation 16593890 · Oct 4, 2019
Related Publication 20220138379A1 · May 5, 2022
References Cited (46)
US 10643384B2 · Brettle · 2020 [cited by examiner]
US 20110238394A1 · Shin et al. · 2011 [cited by applicant]
US 20180319087A1 · Eom · 2018 [cited by applicant]
US 20210103638A1 · Bussman · 2021 [cited by applicant]
CN 102129517 · 2011 [cited by applicant]
CN 104750943 · 2017 [cited by applicant]
KR 100568562 · 2006 [cited by applicant]
Patel, H.V., et al. “Computing interface curvature from volume fractions: A machine learning approach” Computers & Fluids, vol. 193, 104263 (Sep. 2019) (Year: 2019). [cited by examiner]
Scheufler, H. & Roenby, J. “Accurate and efficient surface reconstruction from volume fraction data on general meshes” arXiv:1801.05382v2 (Jan. 2019) (Year: 2019). [cited by examiner]
Maric, T., et al. “voFoam—A geometrical Volume of Fluid algorithm on arbitrary unstructured meshes with local dynamic adaptive mesh refinement using OpenFOAM” arXiv:1305.3417v1 (2013) (Year: 2013). [cited by examiner]
Van Sint Annaland, M., et al. “Numerical simulation of gas bubbles behaviour using a three-dimensional volume of fluid method” Chemical Engineering Science, vol. 60, pp. 2999-3011 (2005) (Year: 2005). [cited by examiner]
A. Alizadehdakhel, M. Rahimi, J. Sanjari, A. A. Alsairafi, CFD and artificial neural network modeling of two-phase flow pressure drop, International Communications in Heat and Mass Transfer, 36 (2009) 850-856. [cited by applicant]
B. Y. Choi, M. Bussmann, A piecewise linear approach to vol. tracking a triple point, International Journal for Numerical Methods in Fluids (2007) 53:1005-1018. [cited by applicant]
C. B. Ivey, P. Moin, Accurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes, Center for Turbulence Research, Annual Research Briefs (2014), pp. 53-68. [cited by applicant]
D. Gueyffier, J. Li, A. Nadim, R. Scardovelli, S. Zaleski, Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows, Journal of Computational Physics 152 (2)(1999) 423-456. [cited by applicant]
D. L. Youngs, Time-dependent multi-material flow with large fluid distortion, Numerical Methods in Fluid Dynamics (1982) 273-486. [cited by applicant]
D. P. Kingma, J. Ba, Adam: A Method for Stochastic Optimization, arXiv e-prints (2017) arXiv:1412.6980. [cited by applicant]
E. Aulisa, S. Manservisi, R. Scardovelli, S. Zaleski, A geometrical area-preserving volume-of-fluid advection method, Journal of Computational Physics 192 (1) (2003) 355-364. [cited by applicant]
E. S. C. Fan, M. Bussmann, Piecewise linear volume tracking in spherical coordinates, Applied Mathematical Modelling 37 (2013) 3077-3092. [cited by applicant]
G. Agbaglah, S. Delaux, D. Fuster, J. Hoepffner, C. Josserand, S. Popinet, P. Ray, R. Scardovelli, S. Zaleski, Parallel simulation of multiphase flows using octree adaptivity and the volume-of-fluid method, Comptes Rend… [cited by applicant]
J. E. Pilliod Jr., E. G. Puckett, Second-order accurate volume-of-fluid algorithms for tracking material interfaces, Journal of Computational Physics 199 (2004) 465-502. [cited by applicant]
J. Ling, A. Kurzawski, J. Templeton, Reynolds averaged turbulence modelling using deep neural networks with embedded invariance, Journal of Fluid Mechanics 807 (2016) 155-166. [cited by applicant]
J. López, C. Zanzi, P. Gómez, F. Faura, J. Hernández, A new volume of fluid method in three dimensions—Part II: Piecewise-planar interface reconstruction with cubic-Bézier fit, International Journal for Numerical Method… [cited by applicant]
J. López, J. Hernández, P. Gómez, F. Faura, A new volume conservation enforcement method for PLIC reconstruction in general convex grids, Journal of Computational Physics 316 (2016) 338-359. [cited by applicant]
J. López, J. Hernández, P. Gómez, F. Faura, VOFTools—A software package of calculation tools for volume of fluid methods using general convex grids, Computer Physics Communications 223 (2018) 45-54. [cited by applicant]
J. Mencinger, I. Žun, A PLIC-VOF method suited for adaptive moving grids, Journal of Computational Physics, 230 (2011) 644-663. [cited by applicant]
K. Kleefsman, G. Fekken, A. Veldman, B. Iwanowski, B. Buchner, A volume-of-fluid based simulation method for wave impact problems, Journal of Computational Physics 206 (1) (2005) 363-393. [cited by applicant]
K. Wu, C. Wang, L. Zhang, Li Hao, Moving boundary problem arising in explosion and impact: An algorithm for Youngs' interface reconstruction, International Journal of Nonlinear Sciences & Numerical Simulation, 10(8), 10… [cited by applicant]
M. Ataei, V. Shaayegan, F. Costa, S. Han, C. B. Park, M. Bussmann, Simulation of polymer foaming using the lattice Boltzmann method, arXiv e-prints (2019) arXiv:1908.04262v1. [cited by applicant]
M. Huang, L. Wu, B. Chen, A piecewise linear interface-capturing volume-of-fluid method based on unstructured grids, Numerical Heat Transfer, Part B: Fundamentals 61 (5) (2012) 412-437. [cited by applicant]
M. Renardy, Y. Renardy, J. Li, Numerical simulation of moving contact line problems using a volume-of-fluid method, Journal of Computational Physics 171 (1) (2001) 243-263. [cited by applicant]
M. Rudman, Volume-tracking methods for interfacial flow calculations, International Journal for Numerical Methods in Fluids (1997) 24:671-691. [cited by applicant]
M. Skarysz, A. Garmory, M. Dianat, An iterative interface reconstruction method for PLIC in general convex grids as part of a coupled level set volume of fluid solver, Journal of Computational Physics 368 (2018) 254-276. [cited by applicant]
M. Sussman, E. G. Puckett, A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flows, Journal of Computational Physics 162, 301-337 (2000). [cited by applicant]
R. Scardovelli, S. Zaleski, Analytical relations connecting linear interfaces and volume fractions in rectangular grids, Journal of Computational Physics 164 (1) (2000) 228-237. [cited by applicant]
R. Scardovelli, S. Zaleski, Direct numerical simulation of free-surface and interfacial flow, Annual Review of Fluid Mechanics 31 (1) (1999) 567-603. [cited by applicant]
R. V. Garimella, V. Dyadechko, B. K. Swartz, M. J. Shashkov, Interface reconstruction in multi-fluid, multi-phase flow simulations, (2005) In: Hanks B.W. (eds) Proceedings of the 14th International Meshing Roundtable. S… [cited by applicant]
S. Popinet, An accurate adaptive solver for surface-tension-driven interfacial flows, Journal of Computational Physics 228 (16) (2009) 5838-5866. [cited by applicant]
S. W. Welch, J. Wilson, A volume of fluid based method for fluid flows with phase change, Journal of Computational Physics 160 (2) (2000) 662-682. [cited by applicant]
T. Milcent, A. Lemoine, Moment-of-fluid analytic reconstruction on a 3D rectangular hexahedrons, Journal of Computational Physics (2019) pp. 1-33. [cited by applicant]
W. E. Lorensen, H. E. Cline, Marching cubes: A high resolution 3D surface construction algorithm, SIGGRAPH Computer Graphics 21 (4) (1987) 163-169. [cited by applicant]
W. J. Rider, D. B. Kothe, Reconstructing volume tracking, Journal of Computational Physics 141 (2) 1998) 112-152. [cited by applicant]
X. Guo, W. Li, F. Iorio, Convolutional neural networks for steady flow approximation, Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (2016) 481-490. [cited by applicant]
X. Yang, A. J. James, Analytic relations for reconstructing piecewise linear interfaces in triangular and tetrahedral grids, Journal of Computational Physics 214 (1) (2006) 41-54. [cited by applicant]
Y. Qi, J. Lu, R. Scardovelli, S. Zaleski, G. Tryggvason, Computing curvature for volume of fluid methods sing machine learning, Journal of Computational Physics 377 (2019) 155-161. [cited by applicant]
Y. Renardy, M. Renardy, PROST: A parabolic reconstruction of surface tension for the volume-of-fluid method, Journal of Computational Physics 183, 400-421 (2002). [cited by applicant]