Numerical Simulation of Liquid Resin Impregnation in a Real Digitized Reinforcement during Liquid Composites Molding Processes

Abstract:

Liquid Composite Molding (LCM) processes, used for producing high-quality, complex composite parts, rely on the uniform infiltration of liquid resin into fibrous fabrics. These fabrics possess a dual-scale structure: highly porous inter-tow spaces surrounding denser fiber tows. The resulting disparity in permeability and flow rates is a primary cause of defect formation, such as voids. To minimize these defects, accurate simulation, incorporating the critical influence of capillary pressure on resin infiltration within the fiber tows, is essential. This work presents a robust numerical model developed to simulate the two-phase resin flow and impregnation dynamics within a digitized, real plain-weave E-glass reinforcement obtained via X-ray micro-computed tomography (CT). The simulation utilizes the open-source multiscale multiphase solver, hybridPorousInterFoam, which employs a Darcy-Brinkman approach, transitioning between Darcy's law in porous regions and Navier-Stokes in free space. A key methodological enhancement involved modifying the advection algorithm using the isoAdvector scheme to mitigate numerical instabilities caused by the high viscosity ratio between the resin and air. Capillary effects at the mesoscale are incorporated through multiscale parameters, specifically the drag and surface tension forces. The key findings demonstrate that the modified solver successfully handles the fluid-fluid interface advection for high viscosity ratios. A parametric study highlighted the significant effect of capillary pressure on multiphase flow within the dual-scale porous media. The numerical results for flow front advancement showed very good agreement when compared against dedicated experimental validation data, confirming the model's high predictive accuracy and its potential for optimizing LCM injection conditions.

You might also be interested in these eBooks

Info:

Periodical:

Materials Science Forum (Volume 1182)

Pages:

103-112

Citation:

Online since:

April 2026

Export:

Share:

Citation:

* - Corresponding Author

[1] Facciotto, S., Simacek, P., Advani, S.G., Pickett, A., & Middendorf, P. (2023). Modeling formation and evolution of voids in unsaturated dual scale preforms in Resin Transfer Molding processes. Composites Part A: Applied Science and Manufacturing, 173, 107675.

DOI: 10.1016/j.compositesa.2023.107675

Google Scholar

[2] Park, C. H., & Woo, L. (2011). Modeling void formation and unsaturated flow in liquid composite molding processes: A survey and review. Journal of Reinforced Plastics and Composites, 30(12), 957–977.

DOI: 10.1177/0731684411411338

Google Scholar

[3] LeBel, F., Ruiz, E., & Trochu, F. (2019). Void content analysis and processing issues to minimize defects in liquid composite molding. Polymer Composites, 40(1), 109-120.

DOI: 10.1002/pc.24609

Google Scholar

[4] Gascón, L., García, J.A., LeBel, F., Ruiz, E., & Trochu, F. (2016). A two-phase flow model to simulate mold filling and saturation in Resin Transfer Molding. International Journal of Material Forming, 9(2), 229–239.

DOI: 10.1007/s12289-015-1225-z

Google Scholar

[5] Raeini, A.Q., Blunt, M.J., & Bijeljic, B. (2012). Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method. Journal of Computational Physics, 231(17), 5653-5668.

DOI: 10.1016/j.jcp.2012.04.011

Google Scholar

[6] Adhikari, D., Lisegaard, J.J., Pierce, R.S., Mikkelsen, L.P., Hattel, J. H., & Mohanty, S. (2025). Meso-scale permeability prediction from the multi-block 3D reconstruction of fiber reinforced polymer composite with experimental validation. Composites Science and Technology, 111290.

DOI: 10.1016/j.compscitech.2025.111290

Google Scholar

[7] Carrillo, F. J., Bourg, I. C., & Soulaine, C. (2020). Multiphase flow modeling in multiscale porous media: An open-source micro-continuum approach. Journal of Computational Physics: X, 8, 100073.

DOI: 10.1016/j.jcpx.2020.100073

Google Scholar

[8] Tong, A. Y., & Wang, Z. (2007). A numerical method for capillarity-dominant free surface flows. Journal of Computational Physics, 221(2), 506–523.

DOI: 10.1016/j.jcp.2006.06.034

Google Scholar

[9] Renardy, Y., & Renardy, M. (2002). PROST: A parabolic reconstruction of surface tension for the volume–of–fluid method. Journal of Computational Physics, 183(2), 400–421.

DOI: 10.1006/jcph.2002.7190

Google Scholar

[10] Inguva, V. (2020). On methods to reduce spurious currents within VOF solver frameworks. Part 1: A review of the static bubble/droplet. Chemical Product and Process Modeling, 1

DOI: 10.1515/cppm-2020-0052

Google Scholar

[11] Deshpande, S. S., Anumolu, L., & Trujillo, M. F. (2012). Evaluating the performance of the two-phase flow solver interFoam. Computational Science & Discovery, 5(1), 014016.

DOI: 10.1088/1749-4699/5/1/014016

Google Scholar

[12] Roenby, J., Larsen, B. E., Bredmose, H., & Jasak, H. (2017). A new volume-of-fluid method in OpenFOAM. In Marine vi: Proceedings of the vi International Conference on Computational Methods in Marine Engineering (p.266–277). https://upcommons.upc.edu/handle/2117/331052

DOI: 10.23967/marine.2025.120

Google Scholar

[13] Roenby, J., Bredmose, H., & Jasak, H. (2019). IsoAdvector: Geometric VOF on general meshes. In OpenFOAM® (p.281–296). Springer International Publishing

DOI: 10.1007/978-3-319-60846-4_21

Google Scholar

[14] isoAdvector.(2025). isoAdvector/isoAdvector[C++]. GitHub. https://github.com/isoAdvector/isoAdvector (Original work published 2016) (Accessed June 24, 2025)

DOI: 10.3795/ksme-b.2019.43.11.739

Google Scholar

[15] Wijaya, W., Ali, M. A., Umer, R., Khan, K. A., Kelly, P. A., & Bickerton, S. (2019). An automatic methodology to CT-scans of 2D woven textile fabrics to structured finite element and voxel meshes. Composites Part A: Applied Science and Manufacturing, 125, 105561.

DOI: 10.1016/j.compositesa.2019.105561

Google Scholar

[16] Pucci, M. F., Liotier, P-J., & Drapier, S. (2015). Capillary wicking in a fibrous reinforcement - orthotropic issues to determine the capillary pressure components. Composites Part A: Applied Science and Manufacturing, 77, 133–141.

DOI: 10.1016/j.compositesa.2015.05.031

Google Scholar