In dry spinning solvent evaporates out of the spun jets and leads to thinning and solidification of the produced fibers. The production process is significantly driven by the interaction of the fibers with the surrounding airflow. Faced with industrial applications producing up to several hundred fibers simultaneously, the direct numerical simulation of the three-dimensional multiphase, multiscale problem is computationally extremely demanding and thus in general not possible. In this paper, we hence propose a dimensionally reduced, efficiently evaluable fiber model that enables the realization of fiber-air interactions in a two-way coupling with airflow computations. For viscous dry spinning of an uni-axial two-phase flow, we deduce one-dimensional equations for fiber velocity and stress from cross-sectional averaging and combine them with two-dimensional advection-diffusion equations for polymer mass fraction and temperature revealing the radial effects that are observably present in experiments. For the numerical treatment of the resulting parametric boundary value problem composed of one-dimensional ordinary differential equations and two-dimensional partial differential equations we develop an iterative coupling algorithm. Thereby, the solution of the advection-diffusion equations is implicitly given in terms of Green's functions and leads for the surface values to Volterra integral equations of second kind with singular kernel, which we can solve very efficiently by the product integration method. For the ordinary differential equations a suitable collocation-continuation procedure is presented. Compared with the referential solution of a three-dimensional setting, the numerical results are very convincing. They provide a good approximation while drastically reducing the computational time. This efficiency allows the two-way coupled simulation of industrial dry spinning in airflow.
A novel numerical scheme including time and spatial discretization is offered for coupled Cahn-Hilliard and Navier-Stokes governing equation sys-tem in this paper. Variable densities and viscosities are considered in the nu-merical scheme. By introducing an intermediate velocity in both Cahn-Hilliard equation and momentum equation, the scheme can keep discrete energy law. A decouple approach based on pressure stabilization is implemented to solve the Navier-Stokes part, while the stabilization or convex splitting method is adopt-ed for the Cahn-Hilliard part. This novel scheme is totally decoupled, linear, unconditionally energy stable for incompressible two-phase flow diffuse inter-face model. Numerical results demonstrate the validation, accuracy, robustness and discrete energy law of the proposed scheme in this paper.
In the present study, we have studied the electro-hydrodynamic of a physical system where a Newtonian dielectric liquid column or droplet suspended in another Newtonian dielectric liquid medium in presence of a simple shear flow. Taking both the phases as leaky dielectric and perfect dielectric in to consideration, we have performed 2D numerical solution for capturing the essential features of droplet deformation in between the parallel plate configuration. For a perfect dielectric system, this study shows that the deformation characteristic follows a monotonic as well as non-monotonic variation with domain confinement depending on the values of electrical permittivity ratio of the droplet and the surrounding fluid. For a leaky-dielectric system, presence of small conductivity further alters the deformation characteristic and it is happened that, at low electric field strength, the deformation increases with confinement monotonically. On contrary, deformation parameter shows non-monotonic variation with the domain confinement at higher electric field strength. Furthermore, in confined domain, the transient evolution of the deformation parameter is also markedly altered by the electric field strength in terms of steady state value of the deformation parameter and steady state time. Finally, the present analysis shows that the domain confinement significantly augment the deformation parameter in presence of electric field that leads to possible droplet break up phenomenon. From the present study, it is worthy to mention that domain confinement can be used to modulate the droplet morphology that has potential applications in modern-days droplet-based micro-fluidic devices.