4. One-way versus two-way coupling
In thermal–structural multiphysics analysis, the coupling strategy determines how the thermal and mechanical fields exchange information during the solution. The two dominant approaches are one-way (sequential) coupling and two-way (fully coupled) coupling. A third variant, two-way sequential coupling, iterates between the two physics solvers until convergence but still uses separate solver kernels. Choosing among these strategies is one of the most consequential decisions in a multiphysics workflow: it affects accuracy, computational cost, convergence robustness, and even whether a solution can be obtained at all.
What Each Strategy Actually Does
In one-way coupling, the thermal solution is computed first and treated as a fixed input to the structural analysis. The temperature field produces thermal strains, and those strains drive the mechanical response. The structural deformation is assumed not to influence the thermal field. This is valid when the mechanical response is small relative to the thermal driving force, or when the thermal boundary conditions are insensitive to deformation.
In two-way coupling, the thermal and structural fields influence each other. The temperature field drives deformation, and the deformation (or the stress state) feeds back into the thermal problem. This feedback can arise through several mechanisms: contact conductance that depends on contact pressure, convective boundary conditions that change with surface geometry, frictional heating that depends on contact pressure and sliding velocity, or temperature-dependent material properties that are themselves affected by stress. Two-way coupling can be implemented either as a monolithic solution (both fields assembled into a single system of equations and solved simultaneously) or as a staggered/partitioned scheme (separate solvers exchanging data at each iteration or time step).
The Decision Framework
The choice is not simply "more coupling is better." A one-way analysis that captures the dominant physics is often the correct engineering choice. The decision should be driven by three questions:
- Does the mechanical state measurably change the thermal boundary conditions or material properties? If contact pressure changes gap conductance by more than a few percent, or if deformation alters convective area significantly, two-way coupling is warranted.
- Is the feedback strong enough to change the engineering conclusion? A small feedback that shifts peak stress by 2% may be irrelevant; a feedback that changes the sign of a stress or triggers contact opening is decisive.
- Can the coupled system be solved robustly and affordably? Monolithic coupling is more robust for strongly coupled problems but requires a single solver that handles both fields and can be memory-intensive. Staggered coupling is more flexible and modular but may converge slowly or diverge for strong coupling.
Practical Example: Brake Disc Thermal Analysis
Consider a ventilated brake disc subjected to a single emergency stop. The thermal problem involves frictional heat generation at the pad–disc interface, conduction into the disc, and convection to the surrounding air. The structural problem involves thermal expansion, contact pressure distribution, and stress. The coupling mechanisms are:
- Contact pressure affects the real contact area and thus the heat partition between pad and disc.
- Thermal expansion changes the disc thickness and the contact pressure distribution.
- Temperature-dependent yield strength and elastic modulus alter the stress response.
For a first-pass analysis, a one-way approach is often used: apply a heat flux distribution based on nominal contact pressure, solve the transient thermal problem, then map the temperature field to a static structural model. This is fast and provides a reasonable estimate of peak temperature and thermal stress. However, if the analysis reveals that the contact pressure redistributes significantly (e.g., the pad edges lift off due to thermal crowning), the heat flux distribution assumed in the one-way analysis is no longer valid. At that point, a two-way staggered approach is needed: at each time step, solve the thermal problem, update the contact pressure and heat partition, solve the structural problem, and iterate until the temperature and displacement changes fall below a tolerance.
Monolithic versus Staggered Two-Way Coupling
When two-way coupling is required, the next decision is monolithic versus staggered. Monolithic coupling assembles the thermal and structural degrees of freedom into a single Jacobian and solves them simultaneously. This is the most robust approach for strongly coupled problems, such as those involving contact conductance with high sensitivity to gap size, or phase-change problems where latent heat release is coupled to stress. The downside is implementation complexity and memory cost: the coupled system is larger, and the solver must handle the different scaling of thermal and mechanical variables.
Staggered coupling solves each field separately and exchanges loads and boundary conditions at the interface. It is easier to implement, allows the use of specialized solvers for each physics, and is often more efficient when the coupling is moderate. The convergence of staggered schemes depends on the strength of the coupling: weak coupling converges in a few iterations, while strong coupling may require under-relaxation or Aitken acceleration to avoid divergence.
Implementation Notes and Common Pitfalls
When implementing a staggered two-way coupling, the order of operations matters. A common scheme is:
- Solve the thermal problem at time step n using the structural state from step n−1.
- Map the temperature field to the structural mesh.
- Solve the structural problem at step n.
- Check convergence: compare the change in interface temperature and displacement between iterations. If the change exceeds a tolerance (e.g., 0.1% of the maximum temperature or 0.5% of the maximum displacement), return to step 1 with the updated structural state.
A frequent pitfall is inconsistent mapping between non-matching meshes. If the thermal and structural meshes differ, interpolation errors can introduce artificial oscillations or energy loss. Use conservative interpolation schemes where possible, and verify that the total heat flux is preserved across the interface.
Another pitfall is neglecting the time-step size. In transient problems, the thermal and structural time scales may differ by orders of magnitude. A time step that resolves the thermal transient may be unnecessarily small for the structural problem, leading to excessive runtimes. Conversely, a time step that is too large for the thermal problem can cause the staggered iteration to diverge. A good practice is to perform a time-step sensitivity study: halve the time step and confirm that the peak temperature and stress change by less than a few percent.
Summary
One-way coupling is the appropriate starting point for most thermal–structural problems. It is fast, robust, and often sufficient when the mechanical response is small or the thermal boundary conditions are insensitive to deformation. Two-way coupling is necessary when the feedback between fields is non-negligible. Within two-way coupling, staggered schemes are preferred for moderate coupling due to their modularity and lower memory cost, while monolithic schemes are reserved for strongly coupled problems where staggered iteration fails to converge. The decision should be revisited as the analysis progresses: a one-way baseline can reveal whether the feedback is large enough to justify the additional cost of two-way coupling.