Thanks for your thoughts. I have tried to reduce the overall process as much as possible already above, it is as follows:
Step 1: Plasma treatment of surfaces
Step 2: UV fix of adhesive used to bond two surfaces
Step 3: Final curing of the adhesive
Response: Die shear strength in shear test
I am considering omitting step 1 in case I can get a response in the desired range without it as it is costly. I can measure some output of step 1, e.g. surface energy, but only whether or not it is above a threshold, not the exact value (ink test). I understand you are suggesting to treat surface energy as a covariate factor instead of having the plasma process parameters in the model itself, which is a good hint, thanks.
There are more steps involved in the overall process, I have reduced them to 3 to keep it short here. I suspect strong interactions between some of the steps, e.g. UV and final curing, adhesive layer thickness and UV curing. Unfortunately, a meaningful measurement of the response can only be conducted after the final curing. Of course, other measurements could be conducted after the individual steps, e.g. polymer chain length, degree of polymerization, but I lack the equipment to do so, which is why we kind of try to condense everything into one response measurement at the end of the process, which is of course less than ideal.
I am not entirely sure what direct sampling refers to in this context.