We present a numerical study of the spinodal decomposition of a binary fluid undergoing shear flow using the advective Cahn-Hilliard equation, a stiff, nonlinear, parabolic equation characterized by the presence of fourth-order spatial derivatives. Our numerical solution procedure is based on isogeometric analysis, an approximation technique for which basis functions of high-order continuity are employed. These basis functions allow us to directly discretize the advective Cahn-Hilliard equation without resorting to a mixed formulation. We present steady state solutions for rectangular domains in two-dimensions and, for the first time, in three-dimensions. We also present steady state solutions for the two-dimensional Taylor-Couette cell. To enforce periodic boundary conditions in this curved domain, we derive and utilize a new periodic Bézier extraction operator. We present an extensive numerical study showing the effects of shear rate, surface tension, and the geometry of the domain on the phase evolution of the binary fluid. Theoretical and experimental results are compared with our simulations. © 2013 Elsevier Inc.

Isogeometric analysis of the advective Cahn-Hilliard equation: Spinodal decomposition under shear flow

DEDE', LUCA;
2013-01-01

Abstract

We present a numerical study of the spinodal decomposition of a binary fluid undergoing shear flow using the advective Cahn-Hilliard equation, a stiff, nonlinear, parabolic equation characterized by the presence of fourth-order spatial derivatives. Our numerical solution procedure is based on isogeometric analysis, an approximation technique for which basis functions of high-order continuity are employed. These basis functions allow us to directly discretize the advective Cahn-Hilliard equation without resorting to a mixed formulation. We present steady state solutions for rectangular domains in two-dimensions and, for the first time, in three-dimensions. We also present steady state solutions for the two-dimensional Taylor-Couette cell. To enforce periodic boundary conditions in this curved domain, we derive and utilize a new periodic Bézier extraction operator. We present an extensive numerical study showing the effects of shear rate, surface tension, and the geometry of the domain on the phase evolution of the binary fluid. Theoretical and experimental results are compared with our simulations. © 2013 Elsevier Inc.
2013
Bézier extraction; Cahn-Hilliard equation; Isogeometric analysis; Shear flow; Spinodal decomposition; Steady state; Computer Science Applications1707 Computer Vision and Pattern Recognition; Physics and Astronomy (miscellaneous)
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11311/1002679
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 87
  • ???jsp.display-item.citation.isi??? 90
social impact