The present research work is focussed on establishing the complex nonlinear input-output relations of a furan resinbased molding sand system. Further, a set of input parameters, which will result in optimized mold properties, is determined. Grain fineness number, setting time, percentage of resin, and hardener are considered as process variables. Mold properties, such as green compression strength, shear strength, mold hardness, gas evolution, permeability, and collapsibility are treated as the process outputs. Nonlinear input-output relations have been developed and statistical analysis has been carried out by utilizing design of experiments, central composite design. Surface plots are developed to study and analyze the input-output relations. The input parameters that will result in best molding conditions and improve casting quality characteristics are determined by utilizing desirability function approach and multiple particle swarm optimization-based crowding distance (MOPSO-CD) techniques. The optimum value for the process variables namely grain fineness number, furan resin, hardener, and setting time are found to be equal to 55, 1.85, 1.2, and 60, respectively. The quality characteristics of the castings namely yield strength, ultimate tensile strength, hardness, density, and secondary dendrite arm spacing are found to improve by 14.03%, 15.08%, 14.14%, 12%, 2.22%, and 12.24%, respectively for the castings made in optimized molding sand conditions.