LinFlowMechPDE
LinFlow-Mechanic coupling¶
openCFS provides the coupling between incompressible/compressible viscous fluid with flexible structure by using the Nitsche method.
Coupling conditions¶
In LinFlow-mechanic coupled problems, two transmission conditions along the interfaces have to be fulfilled:
- Dynamic coupling condition is the continuity of traction (force equilibrium) at the interface
$\mathbf \sigma_\mathrm{m} $ and \mathbf \sigma_\mathrm{v} are mechanic and LinFlow stress tensors, respectively and the surface normal vector is \mathbf n = \mathbf {n}_{\rm m} = -\mathbf{n}_{\rm v}
- Kinematic coupling condition is velocity continuity at the interface \begin{equation} \frac{\partial \mathbf u}{\partial t} = \mathbf v \, , \label{eq:veltranscond-mechcoupl} \end{equation} where \mathbf u and \mathbf v are mechanic displacement and LinFlow velocity vector, respectively.
Weak formulation¶
To obtain the coupling PDEs in Nitsche method, the weak form of mechanic PDE and balance of momentum in LinFlow PDE will be combined and the traction continuity condition (\ref{eq:tractranscond-mechcoupl}) will be applied
In this equation, \mathbf{v'}, \mathbf{u’} are fluid velocity test functions, mechanical displacement test functions, \rho_0 and \rho_\mathrm{m} are LinFlow viscous fluid and mechanic density, respectively. The final version of Nitsche coupling formulations will be obtained by adding the penalty term
In this formulation, the penalty term guarantees the velocity continuity, and the term traction guarantees the traction continuity at the interface.
The penalty term is considered dependent on the characteristic length h_{\rm e} (which is mesh-dependent). This characteristic length is the smaller element length between LinFlow or mechanic regions h_{\rm e} = \operatorname{min}(h_{\rm v}, h_{\rm m}).
The penalty factor \beta is defined as
Where \beta_0 is Nitsche factor and p_{\rm m} and p_{\rm v} are the order of basic functions for mechanic and LinFlow fluid domains, respectively.
Assigning the correct value of the Nitsche factor \beta_0 is crucial for solving this formulation: A small value of the Nitsche factor does not meet velocity continuity at the interface. In contrast, a high value of the Nitsche factor makes numerical errors. A Nitsche factor between 10^7-10^{14} is appropriate for LinFlow-mechanic coupling problems.
Analysis Types¶
Same as LinFlow PDE the coupling can be solved in harmonic transient studies.
- Harmonic
not validated yet: Transient
Defining the coupling¶
LinFlow-mechanics coupling benefits form non-conforming interfaces option. Therefore, similar to other non-conforming interface couplings, in simulation xml-file the non-conforming interface option needs to be defined as
<domain geometryType="">
<regionList>
<region name="" material=""/>
</regionList>
<surfRegionList>
<surfRegion name="LinFLow"/>
<surfRegion name="Mechanic"/>
</surfRegionList>
<ncInterfaceList>
<ncInterface name="LinFLow_mech_interface" masterSide="Mechanic" slaveSide="LinFLow"/>
</ncInterfaceList>
</domain>
slaveSide
and the mechanic interface is considered as the masterSide
.
In the direct coupling section the non-conforming interface defined in <ncInterfaceList>
will be entered as
<couplingList>
<direct>
<linFlowMechDirect>
<ncInterfaceList>
<ncInterface name="LinFLow_mech_interface" formulation="Nitsche" nitscheFactor="10e10"/>
</ncInterfaceList>
</linFlowMechDirect>
</direct>
</couplingList>
nitscheFactor
is optional and by default is set to 10^{10}.
Material and postprocessing results¶
In openCFS material-xml file beside the flow material input for the LinFlow PDE, a mechanic material input for the mechanic PDE should also be defined.
The postprocessing results and the boundary conditions for each PDE are described in the "Singlefield" section.