Aaron Dettmann
COMSOL Employee
Please login with a confirmed email address before reporting spam
Posted:
1 year ago
Sep 25, 2023, 4:04 a.m. EDT
Updated:
1 year ago
Sep 25, 2023, 4:12 a.m. EDT
Hello Y Ko, if you are only interested in the effect of the fluid pressure on the enclosing structure, your case may be analyzed in such a way, that the fluid domain does not need to be modeled explicitly. This modeling technique is shown in the Hyperelastic Seal model:
There are essentially three steps involved:
- Compute the enclosing volume in the undeformed and deformed states.
- Formulate an equation of state for the fluid, p = f(Delta_V, T).
- Apply the computed pressure as a follower load.
The first step involves converting a volume integral into a surface integral using Gauss' theorem (https://en.wikipedia.org/wiki/Divergence_theorem). Assuming that your study is geometrically nonlinear (which is necessary to model the pressure load as a true follower load), you can directly make use of the material and spatial frame coordinates. For 3D, for instance, the undeformed and deformed volumes can be defined as
V_undeformed = -(1/3)*intop1(X*nX + Y*nY + Z*nZ)
V_deformed = -(1/3)*intop2(x*nx + y*ny + z*nz)
where intop1() is an integration operator defined on the enclosing boundaries using the Material frame setting. Similarly, intop2() is an integration operator using the Spatial frame setting. The variables nX, nx etc. are the normal directions in the material and spatial frames, respectively. If the normals point from the solid domain into the fluid domain, there will be a negative sign, otherwise it will be positive. In 2D, the factor 1/3 is replaced by 1/2, and the out-of-plane thickness needs to be taken into account.
Once you have defined the initial and current volumes, you can formulate the pressure law, which in your case will depend on the volume change and the temperature. To do this, you can add a new variable definition for your pressure variable. Depending on the physics interface, add a Boundary Load or Face Load node, and choose the load type Pressure to apply the pressure you defined.
In the upcoming version of COMSOL Multiphysics there will be built-in functionality in the Solid Mechanics interface for this type of simplified fluid-structure analysis.
Hello Y Ko, if you are only interested in the effect of the fluid pressure on the enclosing structure, your case may be analyzed in such a way, that the fluid domain does not need to be modeled explicitly. This modeling technique is shown in the *Hyperelastic Seal* model:
* [https://www.comsol.com/model/hyperelastic-seal-206](https://www.comsol.com/model/hyperelastic-seal-206)
* [https://www.comsol.com/blogs/computing-controlling-volume-cavity/](https://www.comsol.com/blogs/computing-controlling-volume-cavity/)
There are essentially three steps involved:
1. Compute the enclosing volume in the undeformed and deformed states.
2. Formulate an equation of state for the fluid, p = f(Delta\_V, T).
3. Apply the computed pressure as a follower load.
The first step involves converting a volume integral into a surface integral using Gauss' theorem ([https://en.wikipedia.org/wiki/Divergence_theorem](https://en.wikipedia.org/wiki/Divergence_theorem)). Assuming that your study is geometrically nonlinear (which is necessary to model the pressure load as a true follower load), you can directly make use of the material and spatial frame coordinates. For 3D, for instance, the undeformed and deformed volumes can be defined as
V\_undeformed = -(1/3)\*intop1(X\*nX + Y\*nY + Z\*nZ)
V\_deformed = -(1/3)\*intop2(x\*nx + y\*ny + z\*nz)
where intop1() is an integration operator defined on the enclosing boundaries using the Material frame setting. Similarly, intop2() is an integration operator using the Spatial frame setting. The variables nX, nx etc. are the normal directions in the material and spatial frames, respectively. If the normals point from the solid domain into the fluid domain, there will be a negative sign, otherwise it will be positive. In 2D, the factor 1/3 is replaced by 1/2, and the out-of-plane thickness needs to be taken into account.
Once you have defined the initial and current volumes, you can formulate the pressure law, which in your case will depend on the volume change and the temperature. To do this, you can add a new variable definition for your pressure variable. Depending on the physics interface, add a Boundary Load or Face Load node, and choose the load type Pressure to apply the pressure you defined.
In the upcoming version of COMSOL Multiphysics there will be built-in functionality in the Solid Mechanics interface for this type of simplified fluid-structure analysis.