Minimize Your Model File Size with Storing Solution Techniques

A COMSOL Multiphysics® simulation typically includes one or more field quantities in its output. Depending on the number of field quantities, the geometry’s complexity, and the mesh density required for valid results, simulations can include millions of degrees of freedom (DOFs). Oftentimes, storing one or more scalar quantities or the results on a small geometry part is sufficient. Here, we explore tools for storing selected output quantities and minimizing model file sizes and the time required to display this data. Also, a preference setting is available for storing model files in a more compact way.

Editor’s note: This blog post was updated on March 4, 2026, to reflect updated modeling functionality.

Optimizing the Size of Saved Model Files

On the Save page in the Preferences window, there is a setting for optimizing the saving of model files for speed or for file size. The default value for the Optimize MPH files for list is Speed, which uses uncompressed files that are faster to save. If you instead choose File size, the model is saved using compressed files, which can reduce the file size significantly.

The Preferences window with the Optimize MPH file for setting set to File size.

Two Options for Storing Relevant Simulation Output

COMSOL Multiphysics provides two ways for you to include only selected parts of a solution in your output. The first option is to define one or more selections, including the points, boundaries, or domains of interest, so you can restrict the output from the study to only include the fields in the parts of the geometry that those selections define. This method is straightforward and suitable if you are only interested in the simulation output within a certain part of the geometry, where you can evaluate and plot results as well as access the fields and derived quantities as usual.

Note that if you store the solution for some boundaries or points, only the solution fields (the dependent variables) are available using this method. This means that derivatives and quantities that include derivatives, such as stresses and fluxes, are not available as they require that the solution is also available in the domain. If you are interested in storing some derived values such as stresses on just a few boundaries or points of interest, you will want to use our second option.

The second option is to add an ODE and DAE interface (global or distributed) to define a new dependent variable, to which you can transfer the quantity of interest. The quantity of interest could be a global scalar quantity, such as the average or maximum of some field, or it could be the field value along some boundary. (In the latter case, using a selection provides an easier way of achieving the same result.) This method is useful in situations where the output of interest is a global scalar value that you transfer to a single degree of freedom as a variable within a simple algebraic equation. As mentioned above, it is also useful if the quantity of interest is a derived value based on derivatives, such as stresses and fluxes.

Creating a Selection to Store Selected Parts of a Solution

To only store the solution in a selected part of a geometry, create a named selection in the model component. To do so, right-click Definitions and choose a selection from the Selections submenu. An Explicit selection is simple to use if you want the selection to include some specific geometric entities — for example, domains, boundaries, or points. You can give the selection node a descriptive label like Surrounding Air or Substrate Contact. You can use several selection nodes to represent different parts of the geometry and combine them, either when determining what to store in the output or by creating another selection node that implements a Boolean operation, such as a Union or an Intersection selection node. Use one or more of these created selections in the study step settings to define the parts of the geometry for which the solution will be stored.

Example: A Selection for Storing the Deflection on the Top Surfaces of a Mechanical Design

In a solid mechanics simulation, the value of the deflection (displacement) on one or more specific surfaces or points of interest can provide valuable information, which may be sufficient for you as the output from a simulation. The following plot, for instance, indicates the selected top surfaces on the model geometry for a feeder clamp within the Graphics window.

Selected surfaces for a solid mechanics simulation.

To store only the solution on those surfaces (boundaries), create an Explicit selection with the input entities set to the four top surfaces within the geometry where you want to store the solution. If a surface or point is not part of the created geometry, you can add extra Curve or Point nodes in the geometry sequence to divide a boundary or add a point (and a mesh node) at the desired location. The settings for the Explicit selection node appear as shown below.

Explicit selection node settings.

Here, the label has been changed to Deflection Surfaces to reflect what the selection contains.

You can now modify what to store in the output. In the Settings window of the Study step node, locate and expand the Store in Output section. This section contains settings for storing parts of the solution as output. First choose Selection from the drop-down menu in the Output column. You can then click the Add button under Selections to open a list with available selections. In this case, choose the Deflection Surfaces (Boundary) selection and click OK. To completely exclude dependent variables from storing them in a model, choose None from the drop-down menus in the Output column.

Settings for storing the deflection for the selected surfaces in the geometry.

You can now modify what to store in the output. In the Settings window of the Study step node, locate and expand the Values of Dependent Variables section. This section contains Store fields in output settings, where you can choose For selections from the drop-down menu. You can then click the Add button to open a list with available selections. In this case, choose the Deflection Surface selection and click OK.

Settings for storing the deflection for the selected surfaces in the geometry.

Now you can run your simulation. The only output will be the solution on the top surfaces, which you can visualize using a Surface plot.

The displacement in the y direction on the selected surfaces in the geometry. For the rest of the geometry, there is no solution output.

For the solution at points, you can postprocess the output using a Point Evaluation node under Derived Values and display the total displacement, for example. You can also use a Point Graph node in a 1D Plot Group to plot the displacement versus time in a time-dependent simulation or a parameter value in a parametric sweep.

Creating a Variable of Interest to Store in the Output

If you are interested in a scalar quantity, you only need a single degree of freedom (variable) that represents that quantity in the output. You can create such a variable as a simple algebraic equation using a global equation defined via the Global ODEs and DAEs interface. This interface and similar physics interfaces for defining ODEs and DAEs on domains or at points, for instance, are available under Mathematics > ODE and DAE Interfaces in the Add Physics window and in the Model Wizard. In the Settings window for the Global Equations node, define the name of the variable and the simple equation that sets it to the scalar quantity that you want to include in your simulation output.

Scalar coupling operators are important features in COMSOL Multiphysics to utilize for this purpose as well. With these powerful tools, you have the ability to create scalar quantities that are made globally available within the model.

Example: A Variable that Represents the Average Temperature

Say that you are primarily interested in the average temperature within the model geometry. To extract its value, add an Average Coupling Operator (aveop1, for example), defining it to be valid in the entire geometry (all domains) or in the domains of interest. You can then use this operator in a global algebraic equation to make it available in the simulation output as a scalar variable. If you want to store the maximum temperature instead, you can use a Maximum Coupling Operator.

The equation that the Global ODEs and DAEs interface solves for is

In this case, however, the equation should only set the variable equal to the average temperature. As such, it is sufficient to enter aveop1(T)-avtemp if you have called the variable in the algebraic equation avtemp. Remember that the Global ODEs and DAEs interface sets the expression to zero to form the equation, thus solving the simple equation avtemp = aveop1(T). Here, T represents the temperature field in the geometry, which is the dependent variable that you solve for but don’t want to include in the output. The following screenshot shows the Settings window for a Global Equations node in the Global ODEs and DAEs interface for this case.

Global Equations node settings for creating a variable to store.

Notice the Units section within the Settings window shown above. To avoid unit inconsistencies and to evaluate the created variable with the correct units, you should adjust the respective units for the dependent variable and for the source term.

Your final step is to set up the simulation so just the new scalar variable that you have defined is stored in the output. The following plot shows the average temperature in a geometry for a time-dependent simulation as a Global plot defined in a 1D Plot Group.

The mean (average) temperature versus time, computed as a single scalar output from a time-dependent simulation.

Controlling What Is Stored in the Output

In the study that you want to perform, click the study step node for the heat transfer simulation and, in its settings, open the Store in Output section. Then select None from the drop-down menu in the Output column for the Heat Transfer in Solids interface (or another heat transfer interface in your model). The Global ODEs and DAEs interface has no settings for its only degree of freedom.

Turning off the temperature field in the Store in Output section in the settings for the time-dependent study step.

When you have computed the solution, the scalar quantity will be available as a variable that you can display in a table using a Global Evaluation node. You can also plot a Global graph plot to show the deflection versus time or the mean temperature versus some parametric solutions. Since there are no field solutions, any plots that show such variables, or quantities using such variables, will appear with values of zero.

Example: An Algebraic Equation for the Equivalent Stress on the Top Surfaces

Let’s revisit our feeder clamp model to show how to use an algebraic equation to store the equivalent stress (von Mises stress) values on the top surfaces. To do so, add a Boundary ODEs and DAEs interface, available under Mathematics > ODE and DAE Interfaces in the Add Physics window and in the Model Wizard.

In the Settings window for the main Boundary ODEs and DAEs interface, use the same Deflection Boundaries selection for the top surfaces. Make sure to specify the units so that they are compatible with a stress quantity. You can do so by choosing Stress tensor as the quantity for the new dependent variable and for the source, which is the von Mises stress from the main Solid Mechanics interface.

The settings for the Boundary ODEs and DAEs interface, set up for a mechanical stress quantity.

In the Distributed ODE subnode, you define the algebraic equation that sets the new dependent variable bndstress (defined on the top surfaces only) equal to the predefined variable for the von Mises stress, solid.mises, from the Solid Mechanics interface: bndstress-nojac(solid.mises).

The nojac() operator is required, since we do not want this equation to contribute to the Jacobian (the system matrix) for the full model. You enter this expression in the source term and set all other coefficients to zero to solve the equation 0 = bndstress-nojac(solid.mises) (that is, bndstress = nojac(solid.mises)).

The settings in the Distributed ODE node for defining an algebraic equation for the boundary stress.

Before solving, select None from the drop-down menu for Solid Mechanics in the Output column under Store in Output for the Stationary study step. You can then compute the solution and only store the stress values on the top surfaces, which you can plot using a Surface plot.

A plot of the von Mises stress, stored only for the top surfaces.

Optimize Your Simulations with Storing Solution Techniques

As we have shown here today, with just a few simple steps, you can set up a simulation where only part of the fields computed for, or even just one or a few scalar quantities of interest, can be stored in the output. These modeling techniques can significantly reduce your model file size and the time it takes to calculate and display these scalar quantities of interest. This is especially true if you plan to make simulations that include large parametric sweeps or long, detailed time-dependent studies.

To learn about other tools in COMSOL Multiphysics that are designed to optimize your simulation workflow, browse other posts within our General category on the COMSOL Blog. If you have further questions relating to the modeling techniques presented here, please feel free to contact us.