Defining Multiphysics Models Manually with User-Defined Couplings
In COMSOL Multiphysics®, you can manually administer couplings between physics interfaces for which no coupling features are available. We refer to this as the manual approach with user-defined couplings. Here, we discuss how to implement the approach with a combination of physics phenomena in your model. The approach requires more effort to implement than the manual with predefined couplings approach and completely contrasts the automatic approach, which provides an almost fully automated implementation.
When to Use the Manual Approach with User-Defined Couplings
When there is no alternative predefined option to define the multiphysics couplings in a model, there is functionality in the software that provides you with the flexibility to manually implement custom multiphysics couplings. This is done through user-defined variables, expressions, equations, or functions to couple different physics interfaces. We refer to this multiphysics modeling approach as the manual with user-defined couplings approach, since it involves manually administering the coupling and interaction between physics interfaces.The Model Builder in COMSOL Multiphysics, with the model tree showing the Pressure condition selected, the Settings window showing the Pressure condition settings, and the Graphics window showing the model geometry.
The Pressure condition in the probe tube microphone tutorial model, which is used for manually coupling the 1D (pipe) and 3D (exterior) domains in the model.
The Model Builder in COMSOL Multiphysics, with the model tree showing the Pressure Acoustics, Transient condition selected, the Settings window showing the Pressure Acoustics, Transient condition settings, and the Graphics window showing the model geometry.
The Pressure Acoustics, Transient interface of the probe tube microphone tutorial model, wherein the dependent variable, p2, is used to manually couple it to the Pipe Acoustics, Transient interface through a Pressure condition (as pictured in the previous screenshot).
This approach is wide and vast with respect to how it can be implemented for defining multiphysics models. There is no standard that dictates how it necessarily should or can be implemented, as it will vary depending on what is being modeled. Some examples include using expressions to make material properties or boundary conditions functions of the dependent variables of another physics interface.
A screenshot of the model tree with the Equation View option under the Pressure Acoustics node selected, with the list of variables shown on the right.
In the Piezoelectric MEMS Speaker tutorial model, a Boundary Load is used to manually couple the Solid Mechanics interface and Pressure Acoustics, Frequency Domain interface. However, this is also available as a predefined multiphysics coupling in the Acoustics Module.
Equation View for the Pressure Acoustics node, from which a variable is used for a boundary condition in another physics interface (as pictured in the previous screenshot) in order to manually couple the two physics interfaces.
Since this approach requires that you know the equation formulation in detail, and since you have to enter equation contributions manually, it is the most error-prone and usually also the most laborious way of defining multiphysics models. This is why we recommend utilizing the fully automatic approach or manual approach with predefined couplings whenever possible. However, even the manual approach is straightforward to use and gives you full control of the model definition. Although there are readily available predefined multiphysics interfaces and couplings, user-defined couplings are often used in teaching to show students how different equations describing different physics phenomena can be coupled manually. The upside is that you have the predefined multiphysics couplings as reference to compare with. But in most cases, the manual, user-defined couplings are used when there is no predefined multiphysics coupling available. Some guidelines and suggestions on how a manual multiphysics coupling can be implemented are discussed below.
Procedure to Implement the Approach
As noted earlier, the manual approach using user-defined couplings can be implemented in a variety of different ways. Although there is no explicit sequence of steps that must be followed, the general procedure for using the approach is as follows:
- Add the physics interfaces
- Define the physics settings (material properties, boundary conditions, initial conditions, etc.)
- Repeat steps 1 and 2 for each subsequent physics interface
- Define the multiphysics couplings
As noted earlier, since this approach is the most prone to errors, it is important to always check your results after computing the model. You'll notice that defining the physics settings is split into two parts: defining each individual physics interface and then defining the multiphysics couplings. Each physics interface defines its dependent variables and a set of derived variables. A dependent variable may be, for example, the temperature variable. Derived variables may be the components of the heat flux vector, which are defined using the temperature variable. The heat transfer equations are solved for the temperature, while the heat flux components may be used to formulate these equations or to generate streamline and arrow plots when the temperature field is obtained in the solution.
You can define multiphysics couplings by creating your own derived variables of the dependent variables. You can also directly type in expressions of dependent variables in any physics setting's edit fields. However, creating your own variables is a good way of administrating multiphysics couplings. For example, you can create a derived variable for the heat source in a heat transfer model that depends on the electric field. You can do this in the Definitions node by adding a Variable. You can then enter the heat source derived variable in the heat source edit field in the heat transfer settings. A derived variable can be an arbitrary expression of any dependent variable and its derivatives. Once a variable is defined, you can use it to describe sources, material properties, fluxes, forces, boundary conditions, initial conditions, or any other input in the physics interfaces. You can simply type in expressions containing dependent variables, derivatives of dependent variables, and derived variables in the input edit fields for any of the the physics settings in the UI.
The derived variables defined and used internally by COMSOL Multiphysics® can be inspected and changed directly in the UI. By clicking Show More Options in the Model Builder toolbar and selecting the Equation View check box in the dialog box, all variables and equation contributions defined by each node in a physics interface become visible and editable.
Enable showing the Equation View nodes in the model tree by clicking the respective check box.
This means that you can directly change or add terms to the already predefined derived variables. For example, assume that you want to add a contribution to an already defined heat source. Then you can locate the corresponding derived variable in the Equation View node settings and add the contribution, which may, for example, be an expression of the electric potential. In the edit fields in Equation View, you can enter any expression of a dependent variable, derivatives of dependent variables, and other derived variables.
Below, we show how you can used derived variables and the Equation View node settings to define multiphysics couplings.
Walkthrough Example: Modeling Resistive Heating and Thermal Expansion
Similar to the other articles discussing the other approaches for defining multiphysics models, let's use the thermal microactuator tutorial model to show how we would define the multiphysics effects manually. Also, similarly to the other articles, it is assumed that all details regarding the model setup, except for the physics, have been completed. A detailed overview of the device and the physics involved in its operation are described in the Learning Center article introducing the fully automatic approach.
Note: For this combination of physics phenomena, a predefined multiphysics interface, the Joule Heating and Thermal Expansion multiphysics interface, is available to be used and is what we would add to our model instead of using this approach. We recommend employing either the fully automatic or manual with predefined couplings approach for defining the multiphysics in your models whenever possible, as the manual with user-defined couplings approach requires a completely manual implementation. The following walkthrough is meant for demonstration purposes only.
Following the procedure we listed above for implementing this approach would result in us first adding and then defining the physics settings for the electrical part of the model, then the thermal part of the model, and finally the structural part of the model. Thus, the physics settings have been defined, omitting all multiphysics effects, for each of the subsequent physics interfaces: Electric Currents, Heat Transfer in Solids, and Solid Mechanics, respectively. Please note that the physics settings can be defined in any order.
The model tree thus far for the thermal microactuator tutorial model, with the multiphysics effects not yet accounted for.
After defining the settings for each physics interface, we need to determine how we should define the couplings between the physics interfaces. Think about how the physics phenomena are interacting. Does one physics phenomena serve as causation of, or contribute to, another physics phenomena? We should allow the multiphysics phenomena occurring itself to naturally dictate where and how we define the multiphysics effect in a model.
Since the electric heating is a heat source, it would be logical to define the effect under the Heat Transfer in Solids interface and use an expression that is a function of the electric field and the electric conductivity. We also know that the heat source may be present everywhere in the model domain in this case. We can therefore select the Physics ribbon tab, while the Heat Transfer in Solids interface is selected, and locate the Heat Source feature under the Domains menu to create the multiphysics coupling.
The Heat Source node.
We want to locate the electric losses and add them as a heat source in the heat transfer equations (temperature equations). So now we need to use an expression that describes this quantity. We could define it manually, since we know that the electric heating is the scalar product of the current density vector and the electric field, . However, in this case, this is already predefined in COMSOL Multiphysics®. We can locate the electric heating term under the Current Conservation node and enable displaying the Equation View node to access all the variables and expressions used by that physics node. There, we'll find the variable ec.Qrh, which describes the electric heating in the domain. We can see how this derived variable is defined. As we expected, it is formed by the scalar product of the current density vector and the electric field vector , which are also defined as derived variables. The dependent variable here is the electric potential, V, and the derived variables for the current density vector and the electric field use this variable in their definition (actually the gradient of the electric potential). In the same way as COMSOL Multiphysics® internally defines its derived variables, we can also use this syntax to create our own derived variables, if needed. We can also edit the derived variables defined by COMSOL Multiphysics® in their corresponding edit field. In this case, we have the benefit of the electric heating term already being defined as a derived variable. We can use this variable as the expression for the value in the Heat Source domain node.
All these steps can be avoided, since the electric heating is predefined and can be added automatically by the multiphysics node. However, in cases when a coupling is not predefined, you would use the strategy outlined above.A screenshot of the Model Builder, with the model tree on the left showing the Equation View node selected under the Current Conservation node, with the Settings window on the right.
The Equation View node for the Current Conservation node. The variable ec.Qrh is used to introduce electric heating as a heat source in the modeling domain, i.e., to couple the equations for the electric potential with the equation for temperature.
The Heat Source domain feature with the updated settings for the value applied.
Next, we need to define the thermal expansion that couples the Solid Mechanics and Heat Transfer in Solids interfaces. Since thermal expansion gives a contribution to the strain tensor, we need to locate the strain tensor in the Solid Mechanics interface and add the expression for thermal strain. In this case, we need the strain tensor formulation in one of the Equation View subnodes to the Structural Mechanics interface. We can search for the strain tensor by using the Find tool by pressing Ctrl-F and entering strain tensor in the search dialog box. The Find Results table is then updated in the window below the Graphics window. Double-clicking on one of the strain tensor components listed in the Find Results table takes us to the correct Equation View node. For the device we are modeling, the relevant Equation View is the subnode to the Linear Elastic Material node. We can add the thermal strain contribution as an expression of temperature, T, in the respective strain tensor component.
Note: You can also couple the physics for this model by adding Thermal Expansion as a subnode to the Linear Elastic Material node. However, we will ignore that such a feature is available for the purposes of this demonstration to implement the manual approach using user-defined couplings.
The dependent variable, T, for the Heat Transfer in Solids interface, which we can incorporate into expressions used for the Solid Mechanics interface.
Enabling the Equation View allows us to access all the variables and the expressions used for the Linear Elastic Material node. Here, we can modify the equations to define any coupling between the strain and temperature.A screenshot of the Equation View subnode to the Linear Elastic Material node selected in the model tree, with the Settings window open on the right.
The Equation View subnode Settings window, wherein the variables and the respective expressions that will be modified are selected in the table and highlighted in blue.
We'll want to modify the strain definition for the variables solid.eXX, solid.eYY, and solid.eZZ using the equation for thermal strain,
wherein α is the material's thermal expansion coefficient, T the temperature field, and T0 is the strain reference temperature. The values for α and T0 are known and defined by the parameters alphaps and T0, respectively. We can add the contribution to the strain tensor components as seen in the screenshot below. Note that this expression for thermal expansion is valid only for relatively small deformations, which we may assume that we have in our case.A screenshot of the Settings window for the Equation View subnode, with warning icons shown for different variables in the list.
The Equation View subnode for the Linear Elastic Material node. The expressions for the variables highlighted with a warning icon have been edited to account for thermal strain.
To see this approach be implemented for modeling this combination of physics, you can watch the video below for a full demonstration.
Tutorial Video: Defining a Multiphysics Model Using the Manual Approach with User-Defined Couplings
Demonstrate the knowledge you have gained from this article by putting it into practice with the follow-up modeling exercises listed below. The directions provided for each modeling exercise contain all the information necessary in order to complete it, but are intentionally generalized to encourage self-guided problem-solving. You can check your implementation with the solution model files provided by comparing them manually or using the Comparison tool to exactly identify the differences.
- Previously, you built the busbar tutorial model using the fully automatic approach and the manual approach with predefined couplings. Now, use the manual approach with user-defined couplings to define the physics for this model. Follow the directions provided to define the physics setup. The geometry is available to download: busbar.mphbin.
- Previously, you built the free convection tutorial model using the fully automatic approach and the manual approach with predefined couplings. Now, use the manual approach with user-defined couplings to define the physics for this model. Follow the directions to define the physics setup. The geometry is available to download: free_convection.mphbin.
- Previously, you built the convective cooling of a busbar tutorial model using the manual approach with predefined couplings. Now, follow the directions provided to define the physics for the model using the manual approach with user-defined couplings. The geometry is available to download: busbar_box.mphbin.