The central idea behind a numerical approach for parameter tolerance design is the concept that multiple numerical process simulations can be used in accordance with the Taguchi orthogonal experiment arrays to perform parameter sensitivity analysis and for tolerance design using a quality loss function. The quality loss function is analytically defined using a novel material modeling technique based upon the second law of thermodynamics.