DOE OSTI2015
This report details the study of PMDI structural foam and the development of a computational model to represent it. This polyurethane has a fast catalyst, such that foaming and polymerization occur simultaneously. The foam is over-packed to twice or more of its free rise density to reach the density of interest. Our approach is to combine model development closely with experiments to discover underlying physical phenomena, to parameterize models, and to validate the models once they have been developed. The model must be able to represent the expansion, filling of molds, curing, and final foam properties. PMDI is a chemically blown foam, where carbon dioxide is produced via the reaction of water and isocyanate. The isocyanate also reacts with polyol in a competing reaction, which produces the polymer. A new kinetic model is developed and implemented, which follows a simplified mathematical formalism that decouples these two reactions. The model predicts the polymerization reaction via condensation chemistry, where vitrification and glass transition temperature evolution must be included to correctly predict the temperature dependence of the kinetics. The foam gas generation kinetics are determined by tracking the molar concentration of both water and carbon dioxide. Understanding the thermal history and loads on the foam due to exothermicity and oven curing is very important to the results, since the kinetics and material properties are all sensitive to temperature. The conservation equations, including the equations of motion, an energy balance, and three rate equations are solved via a stabilized finite element method. We assume generalized-Newtonian rheology that is dependent on the extent of cure, gas fraction, and temperature. The conservation equations are combined with a level set method to determine the location of the free surface over time. Results from the model are compared to experimental flow visualization data and post-test X-ray computed tomography (CT) data for the density. Several geometries are investigated including a mock encapsulation part, two configurations of a mock structural part, and a bar geometry to specifically test the density model. We have found that the model predicts both average density and filling profiles well. However, it under predicts density gradients, especially in the gravity direction. Further model improvements are also discussed.