How to Compute Distances Between Objects in COMSOL Multiphysics®

March 2, 2017

Have you ever modeled deforming objects in the COMSOL Multiphysics® software and wanted to know the distance between them? In today’s blog post, we will look at how to compute distances between objects using methods for determining the closest distance field. We’ll also find out how to use the distance field as a part of a multiphysics model.

Computing Distances with the General Extrusion Component Coupling

Let’s consider two objects labeled A and B, shown below. The three distances that we want to compute are:

  1. Distance to object A as a distance field. In this case, we calculate the distance and direction from all of the points surrounding and inside object A to the closest point on its boundary (dA).
  2. Distance from every point on the boundary of object B to the closest point on object A (dAB).
  3. Endpoints of the line that is the shortest distance between objects A and B (dAB,min).

A schematic showing two objects and equations for the distances between them.
Two objects, A and B, and the distances that we want to compute.

We can compute all of these various distances using a combination of the General Extrusion and Minimum component couplings in COMSOL Multiphysics. Let’s first look at how to use the General Extrusion component coupling. We name the operator A_b and define its Source Selection to be the boundaries of object A. Within the Advanced section, we use the Mesh search method of Closest point. These settings are shown below. All other settings for this operator can be left at their defaults.

A screenshot showing the settings for the General Extrusion Component Coupling, used to compute distances between objects in COMSOL Multiphysics®.
The settings for the General Extrusion component coupling used to compute the closest point distance. Note that the Mesh search method is set to Closest point.

We use this operator within the definition of a variable called d_A, defined as:

sqrt((x-A_b(x))^2+(y-A_b(y))^2)

This variable is defined over the domains where we want to compute the distance field; in this case, just the surrounding domain. We can also compute the negative of the gradient of this distance field, -\nabla d_A(\mathbf{x}). This gives us the components of a vector field that points toward the closest point on the boundary of A. We can use the Differentiation operator d(d_A,x) and d(d_A,y) to take the spatial derivatives, as shown in the screenshot below.

A screen capture showing the variable definitions.
The variable definitions.

We can use these variables anywhere that we want. For example, we can plot the distance field or make material properties dependent upon distance. The image below plots the contours of the distance and the direction vectors. Note that the distance is computed even in the region behind object B. We clearly get quite a bit of information here, but there is a substantial computational cost, since the shortest distance is computed at every point in the surrounding domain. There are also times when we don’t need all of this information and just want the distances between objects.

A plot of the distance field and shortest direction between the two objects.
The distance field (contour lines) and shortest direction to the boundary of object A (arrows) in the domain surrounding the two objects.

Calculating Clearances Between Objects

Let’s make things a bit easier and only concern ourselves with the distance between two objects and not the direction. We use the same General Extrusion component coupling, but only need to define a variable on the boundary of object B to compute the distance.

Image showing the variable definition for the distance between the two objects.
The variable defining the distance between the objects.

While this is the same distance function we used before, we don’t need a mesh in the intermediate space. We don’t even need a mesh over domains A and B; there just needs to be a mesh on the boundaries of the objects. This approach takes much less time, but it gives us only the shortest distance from object A to every point on the boundary of object B. We cannot recover the direction vector. We can also flip all of these definitions around and compute the shortest distance from object B to every point on the boundary of object A. These distances, shown in the plot below, are available along the boundaries of the objects.

A plot of the closest distances between points on the boundaries of the two objects.
The distance from every point on the boundary of object B to the closest point on object A and vice versa.

Computing the Shortest Line Between Objects Using the Minimum Component Coupling

Now, let’s find the line that describes the shortest distance between the two objects. In the previous section, we saw that we can compute two variables, d_AB and d_BA, which describe the shortest distances between A and B and vice versa. We now want to find the minimum distance between the boundaries of these domains. Thus, we set up two different Minimum component couplings: one for the boundary of object A and another for object B. We call these operators minA and minB, as shown in the screenshot below.

A screenshot of the settings for the Minimum component coupling in COMSOL Multiphysics®.
The definition of the Minimum component coupling over the boundary of object A.

We then call these Minimum component couplings to extract the minimum distance. We can also provide a second argument to the Minimum component coupling to find the coordinates at which the distance is at a minimum. For example, by defining the variable A_x as the expression minA(d_BA,x), it takes on the value of the x-coordinate at which d_BA is at a minimum over the boundary of A.

A graphic showing the Settings window for the coordinate definitions for the line distances between domains.
The definitions for the coordinates of the shortest line segment between two domains.

We can call the variables defining these coordinates anywhere we want. For example, we can use the Cut Line feature to show the shortest line segment connecting the two objects, as seen in the following image. If we have a meshed domain and a solution between the two objects, then we can plot the fields just along the shortest line between the two.

Two images: a screenshot of the settings for the Cut Line feature and an overlay showing a plot of the shortest distances between the two objects.
The Cut Line feature, used to determine the shortest line between objects.

Limitations and Uses of Calculating Distances in COMSOL Multiphysics®

These techniques for determining distances can be used in any model. Although the examples presented here are in 2D, they can all be generalized to 3D as well. However, computing the 3D distance field does take a relatively long time, whereas calculating distances between boundaries and clearances is less intensive.

Computing the distance field around nonsmooth shapes also requires a bit more care. As shown in the figure below, the distance field around reentrant corners are nonsmooth, hence the direction vector will be undefined along those lines that are equidistant from two different parts of the boundary. Resolving this nonsmoothness of the distance field requires a finer mesh.

Four images showing results plots of the distance fields around and inside an object.
The distance field around and inside an object with reentrant corners on a coarse mesh (left) and a more refined mesh (right). The smoothness of the distance field is mesh dependent in such cases.

Once we have computed this distance field on an appropriately fine mesh, we treat it like any other variable in our model. For example, we can make material properties a function of distance from a surface. The image below shows such a representative material distribution.

A graphic showing a representative material distribution that is a function of distance to the surface.
A representative material distribution that is a function of distance to the surface.

It is also possible to use the distance function to help visualize our results. Suppose we are only interested in the part of the solution that is within a specific distance of the surface. In this case, we can use the Filter subfeature when making a volume plot. We then enter a logical expression to only display the results that are within a certain distance of the object’s surface, an example of which is shown below.

Two images: a screenshot of the settings for the Filter subfeature with an overlay of the results plot for surface distances.
Using the distance function to plot only the solution within 5 mm of the surface.

Concluding Thoughts on Computing Distances Between Objects in COMSOL Multiphysics®

We have demonstrated how to compute a distance field to a boundary within a model, the distances between boundaries, and the shortest line segment between two boundaries. This approach also works to calculate distance fields from edges and points in 3D models. The computed distances can be used anywhere within the setup, physics definitions, and results evaluation of a model. We’ve shared a couple of examples here, but now it’s your turn. We would love to hear what you come up with!

Further Resources


Comments (9)

Leave a Comment
Log In | Registration
Loading...
Amr Al Abed
Amr Al Abed
April 9, 2017

Thanks Walter for an excellent post and tutorial.

So far, I’ve been using COMSOL’s Wall Distance physics to do these calculations.

Just wondering what are the main differences between the two approaches (extrusion variables vs. WD physics) in terms of pros and cons?

Do you have any recommendations when it’s preferable to use one approach over the other?

Best regards
Amr

David Ashirov
David Ashirov
January 7, 2022

Another difference besides what Dr. Frei already mentioned is that the Wall Distance solution can be used only once computed in a specialized study step. General extrusion based distance variables can be used immediately after the definition in the same component or even in a different component without storing any solution.

Walter Frei
Walter Frei
April 10, 2017

Hello Amr, Excellent question. The Wall Distance physics solves a form of the Eikonal equation (https://en.wikipedia.org/wiki/Eikonal_equation) and we use it, for example, in turbulent fluid flow problems where the distance to the wall is needed. (see: https://www.comsol.com/blogs/which-turbulence-model-should-choose-cfd-application/)
If you want to compute the distance field within an entire domain, the Wall Distance physics has the advantage of being faster and using less memory than this approach. However I don’t think that you can use the Wall Distance interface to simply find the distances, or minimum distance, between boundaries. There may be other differences, of course, but these spring easily to mind.

Ali Kadivar
Ali Kadivar
April 11, 2018

Many Thanks Walter for an excellent and useful lemma you provided. it is perfect to calculate distance between two boundaries but how about inside a contour.
Let’s imagine we have a time variable closed contour between two point for instance “sigmaxx>A”
between two electrodes while contour is in a fluid and going larger and changing shape.

we want to find the width and minimum width i.e. distance between two boraders surrunding the material with sigmaxx>A and acts as contour wall.

do you have any idea about it?
it will be appreciated if you can furnish me up with a solution.

Best regards
Ali

Caty Fairclough
Caty Fairclough
April 26, 2018

Hi Ali,

Thanks for your comment!

For your question, I suggest contacting our Support team.

Online Support Center: https://www.comsol.com/support
Email: support@comsol.com

ZHEN LI
ZHEN LI
August 7, 2019

Dear Author,

Thank your paper for broadening my horizon. However, I have a question that how to set to obtain the Figures 4 and 6 that this paper shows. Thank you for your time.

Best wishes
Li

Nick Devlin
Nick Devlin
November 23, 2019

Hi – thanks for this article. Is it possible to automatically measure the maximum x, y or z dimension of a domain in a similar manner? (The geometry is created in Rhino and imported into Comsol) For example, I need to complete heat transfer calculations where the thermal conductivity of the material applied to a domain is dependent on the dimensions of the domain in one or two directions and the aspect ratio between the two.

Amirhossein Omidvar
Amirhossein Omidvar
March 2, 2021

Thanks, Walter for the tutorial!
I’m trying to find how the gap between two membranes in 3D changes while they are under certain loads. I applied the same method to find the “Clearances” between them-including defining the general extrusion component and a variable with x,y,z parameters. I wonder how I can plot the results?

Jahangier Ahmad
Jahangier Ahmad
November 12, 2021

Hi Walter
Can you please share the .mph file of this model (Computing distances between objects in COMSOL)

EXPLORE COMSOL BLOG