journal article·2026-04-17Open access

Numerical Simulation of Co-Continuous Morphologies in PEO/PS Polymer Blends

Seungjae Lee, Yongho Choi, Junseok Kim

Applied Sciences· Volume 16· Issue 8

Abstract

This paper investigates co-continuous structures in immiscible polymer blends through three-dimensional (3D) computational calculations based on a multiphase phase-field equation for fluid flow. The mathematical model describes phase separation with the Cahn–Hilliard (CH) equation and fluid motion with the incompressible Navier–Stokes (NS) equations. Both polymers are treated as Newtonian viscous fluids, and the model includes surface tension, viscosity, and volume fraction effects. A semi-implicit finite difference method (FDM) solves the CH equation, and a projection method maintains the incompressibility of the flow field. Multigrid techniques solve the nonlinear systems efficiently. In addition, a connectivity-based detection algorithm determines whether a phase forms a connected structure that reaches all boundaries of the numerical domain. The numerical results show that the morphology changes from a droplet–matrix structure to a co-continuous structure as the volume fraction increases. The interfacial area per unit volume reaches a local maximum near the transition between these two regimes.
Source: crossref · Quality: baseline

Research context

OpenAlex topics: Block Copolymer Self-Assembly · Solidification and crystal growth phenomena · Fluid Dynamics and Thin Films

SciVX discovery tags are terminology-based. OpenAlex topics are source-provided classifications.

Affiliations

Department of Mathematics, Korea University, Seoul 02841, Republic of Korea

Department of Computer & Information Engineering, Daegu University, Gyeongsan-si 38453, Republic of Korea

Funding

National Research Foundation of Korea · Ministry of Education

References (41)

Show source-provided references
1. Polymer alloys and blends
2. Shamsuri, A.A., and Jamil, S.N.A.M. (2021). Application of quaternary ammonium compounds as compatibilizers for polymer blends and polymer composites—A concise review. Appl. Sci., 11.10.3390/app11073167 ↗
3. López-Martínez, E.I., Zaragoza-Contreras, E.A., Vega-Rios, A., and Flores-Gallardo, S.G. (2024). Effect of the compounding method on the development of high-performance binary and ternary blends based on PPE. Appl. Sci., 14.10.3390/app142210264 ↗
4. Phase field modeling of fiber-based thermal diffusion and phase transitions in the fused deposition modeling process10.1016/j.cnsns.2025.109071 ↗
5. Formation of co-continuous structures in melt-mixed immiscible polymer blends10.1081/MC-120018022 ↗
6. Blends of polydimethylsiloxane-based polyurethane and poly (propylene glycol)-based polyurethane with co-continuous structures: Morphology evolution, synergistic effects and application in strain sensors10.1016/j.polymer.2024.127615 ↗
7. Preparation of self-healing thermoplastic polysiloxane–polyurea/polyether–polyurea elastomer blends with a co-continuous microphase structure and in-depth research on their synergistic effects10.1021/acsami.4c12019 ↗
8. A dual-energy physics-informed multi-material topology optimization method within the phase-field framework10.1016/j.cma.2025.118338 ↗
9. On the unconditionally stable phase field model of ternary components system considering liquid-solid phase transition10.1016/j.jcp.2025.114404 ↗
10. A novel ensemble Kalman filter based data assimilation method with an adaptive strategy for dendritic crystal growth10.1016/j.jcp.2024.113711 ↗
11. Efficient prediction of phase-field crystal dynamics via β-variational autoencoders and time-series transformers on coupled physical fields10.1016/j.camwa.2025.12.024 ↗
12. Comparison of methods for the detection of cocontinuity in poly(ethylene oxide)/polystyrene blends10.1002/pen.20064 ↗
13. Cocontinuous polymer blends: Influence of viscosity and elasticity ratios of the constituent polymers on phase inversion10.1016/S0032-3861(01)00100-8 ↗
14. Evolution of polymer blend morphology during compounding in an internal mixer10.1016/S0032-3861(99)00022-1 ↗
15. Role of rheology in morphology development and advanced processing of thermoplastic polymer materials: A review10.1021/acsomega.3c03310 ↗
16. Image analysis for interfacial area and cocontinuity detection in polymer blends10.1016/S0032-3861(02)00282-3 ↗
17. Flexible and sustainable PLA/PBAT-g-GMA nanocomposites based on carbon nanotubes with potential for electrostatic control10.1002/pat.70369 ↗
18. Detection of compatibilizer efficiency in conductive polymer blends by impedance spectroscopy: Microstructure evaluation approach10.1002/pen.26633 ↗
19. Liu, C., Yu, C., Sang, G., Xu, P., and Ding, Y. (2019). Improvement in EMI shielding properties of silicone rubber/POE blends containing ILs modified with carbon black and MWCNTs. Appl. Sci., 9.10.3390/app9091774 ↗
20. Alhamidi, A., Anis, A., Al-Zahrani, S.M., Bashir, Z., and Alrashed, M.M. (2022). Conductive plastics from Al platelets in a PBT-PET polyester blend having co-continuous morphology. Polymers, 14.10.3390/polym14061092 ↗
21. Light-colored conductive nanorod-mediated co-continuous polypropylene/poly (butylene succinate) blends: Regulatory mechanism, morphology evolution, electrical, and optical properties10.1002/app.54747 ↗
22. Jariyasakoolroj, P., Kumsang, P., Phattarateera, S., and Kerddonfag, N. (2024). Enhanced impact resistance, oxygen barrier, and thermal dimensional stability of biaxially processed miscible poly (lactic acid)/poly (butylene succinate) thin films. Polymers, 16.10.3390/polym16213033 ↗
23. Co-continuous polymer systems: A numerical investigation10.1016/j.commatsci.2014.10.039 ↗
24. Phase-field modeling and simulation of high-strength porous polymer fabrication via immiscible polymer blending for bone implants10.1115/1.4068548 ↗
25. Continuum-scale modelling of polymer blends using the Cahn–Hilliard equation: Transport and thermodynamics10.1039/D1SM00272D ↗
26. Phase-field modeling and numerical computation for the hydrodynamics coupled diblock copolymer in complex domains10.1016/j.cnsns.2025.109381 ↗
27. An effective numerical method for solving ternary Cahn–Hilliard-type Nakazawa–Ohta phase-field model on three-dimensional surfaces10.1016/j.cnsns.2025.109427 ↗
28. A stable numerical framework for long-time dynamic crack analysis10.1016/j.ijsolstr.2024.112768 ↗
29. Dynamic analysis of bi-material interfacial cracks by the high-order GFDM with an enhanced Krylov deferred correction technique10.1016/j.ijsolstr.2025.113451 ↗
30. On quadratic and curvature-dependent variable mobilities for ternary phase-field fluid simulations with matching densities10.1016/j.apm.2026.116937 ↗
31. Ottino, J.M. (1989). The Kinematics of Mixing: Stretching, Chaos, and Transport, Cambridge University Press.
32. A quantitative image analysis method for the determination of cocontinuity in polymer blends10.1016/0032-3861(95)90929-V ↗
33. Kim, J. (2002). Modeling and Simulation of Multi-Component, Multi-Phase Fluid Flows. [Ph.D. Thesis, University of Minnesota].
34. On the numerical approximation of a phase-field volume reconstruction model: Linear and energy-stable leap-frog finite difference scheme10.1016/j.cnsns.2025.109104 ↗
35. Three-dimensional narrow volume reconstruction method with unconditional stability based on a phase-field lagrange multiplier approach10.1016/j.camwa.2025.11.009 ↗
36. Trottenberg, U., Oosterlee, C.W., and Schüller, A. (2001). Multigrid, Academic Press.
37. A three-dimensional multi-phase-field vesicles model and its practical finite difference solver10.1016/j.cpc.2026.110053 ↗
38. A numerical method for the incompressible Navier–Stokes equations based on an approximate projection10.1137/S1064827593244213 ↗
39. Consistently energy-stable decoupled method with second-order accuracy and lower density bounds for the incompressible fluid flows with variable density10.1016/j.cnsns.2025.109303 ↗
40. Energy-stable decoupled numerical approximation with practical correction technique for the binary phase field Darcy fluid system10.1016/j.cma.2026.118791 ↗
41. An unsplit, higher order Godunov method for scalar conservation laws in multiple dimensions10.1016/0021-9991(88)90065-4 ↗