[ripudaman]

I have a created a custom fixedGradient Boundary condition based on the tractionDisplacement boundary condition in OF 2.3.x. In this boundary condition I am calculating the traction distribution along a particular boundary patch. In my model I have two overlapping patches (baffles) that I obtain by using snappyHexMesh and snapping a 2D surface in the blockMesh grid. The traction boundary condition on the two overlapping patches causes the opening of the patches as shown in the figure below.(Here the traction condition has a uniform distribution)

The boundary condition on the two faces is equal and opposite. Thus I want to calculate the traction distribution for one patch and just apply the equal and opposite traction distribution on the opposite (corresponding) patch. Since my current implementation is in the boundary condition file, the solution of the traction distribution along both the patches is being calculated and my simulation is wasting precious computation time. I want to avoid that and calculate only once. Do you have an solution to this problem?