# Influence of Quasiperiodic Gravitational Modulation on Convective Instability of Liquid-Liquid Polymerization Front

Authors: Saadia Assiyad; Karam Allali; Mohamed Belhaq
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

The influence of quasiperiodic gravitational modulation on convective instability of polymerization front with liquid monomer and liquid polymer is studied. The model includes the heat
equation, the concentration equation, and the Navier-Stokes equations under the Boussinesq approximation. The linear stability analysis of the problem is carried out and the interface problem is derived. Using numerical simulations, the convective instability threshold is determined and the boundary of the convective instability is obtained for different amplitudes and frequencies ratio.


---

## Body

## 1. Introduction

Frontal polymerization is the process of polymer production in propagating reaction fronts [1–4]. In the absence of vibration, the influence of convective instability on polymerization front when the monomer is liquid and the polymer is solid was studied in [5], while the case of liquid polymer was considered in [6]. The influence of periodic gravitational modulation on the convective instability in the case of liquid-solid polymerization front was studied [7] and it was shown that the propagation of polymerization reaction front is strongly affected by the amplitude and the frequency of vibrations. In particular, the polymerization front can be stable or unstable depending on the values of vibration parameters. The influence of periodic vibrations on convective instability of reaction front was also studied in the case of liquids [8] and it was concluded that, for small vibration amplitudes, the reaction front remains stable and it loses its stability for sufficiently large amplitude of vibrations.Recent works were devoted to the influence of a quasiperiodic (QP) gravitational modulation on the convective instability of reaction front. For instance, the influence of the QP gravitational modulation on reaction front was examined in the case of porous media described by the Darcy equation [9]. On the other hand, the case of liquid-solid polymerization front was considered in [10] using the Navier-Stokes equations instead of the Darcy equation and it was revealed that both the amplitudes and the frequencies ratio influence the stability domain of the polymerization front. Specifically, for appropriate values of vibration amplitudes and increasing values of the frequencies ratio, a stabilizing effect is observed. The effect of the wave number on the reaction front was also examined showing that increasing the wave number widens the stability domain.The present work studies the effect of QP gravitational modulation on the convective instability of the polymerization front, but in the case of liquid-liquid frontal polymerization. This case is different from the previous one [10] in the sense that in [10] the equation of motion is considered only after the reaction zone since the polymer is in the solid phase. In the present work instead, the equation of motion is considered before and after the reaction zone because both the monomer and the polymer are liquids.The next section presents the frontal polymerization model. Section3 develops the perturbation analysis, while the interface problem is examined in Section 4. The linear stability analysis is discussed in Section 5 and numerical investigations are carried out in Section 6. The last section concludes the work.


## 2. Frontal Polymerization Model

We consider a polymerization process with a liquid reactant and a liquid product by assuming that the reaction front propagates in the direction opposite to the direction of gravity. In this case, the model of the frontal polymerization is given by the system of equations:(1)
∂
T
∂
t
+
v
·
∇
T
=
κ
Δ
T
+
q
W
,
∂
α
∂
t
+
v
·
∇
α
=
W
,
∂
v
∂
t
+
v
·
∇
v
=
-
1
ρ
∇
p
+
ν
Δ
v
+
g
1
+
b
t
β
T
-
T
0
γ
,
div
⁡
v
=
0
,with the following boundary conditions:(2)
z
⟶
+
∞
:
T
=
T
i
,
α
=
0
,
v
=
0
,
(3)
z
⟶
-
∞
:
T
=
T
b
,
α
=
1
,
v
=
0
.The gradient, divergence, and Laplace operators are defined by(4)
∇
v
=
∂
v
∂
x
,
∂
v
∂
y
,
∂
v
∂
z
,
div
⁡
v
→
=
∂
v
1
∂
x
+
∂
v
2
∂
y
+
∂
v
3
∂
z
,
Δ
v
=
∂
2
v
∂
2
x
+
∂
2
v
∂
2
y
+
∂
2
v
∂
2
z
,where the variables (
x, y, and z
) are the spatial coordinates such that -
∞
<
x
,
y
,
z
<
+
∞, T is the temperature, α is the concentration of the reaction product, v is the velocity, p is the pressure, κ is the coefficient of thermal diffusivity, q is the adiabatic temperature heat release, ρ is an average value of the density, ν is the coefficient of kinematic viscosity, γ is the unit vector in the upward direction, β is the coefficient of thermal expansion, g is the gravitational acceleration, and b
(
t
) is the QP acceleration acting on the fluid which is given by b
(
t
)
=
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
(
σ
2
t
) such that λ
1, λ
2 are the amplitudes and σ
1, σ
2 are the incommensurate frequencies of the QP gravitational modulation. The quantity T
0 is a mean value of temperature, T
i is the initial temperature, and T
b
=
T
i
+
q is the temperature of the reacted mixture. We consider one-step reaction of zero order where the reaction rate is defined as follows:(5)
W
=
k
T
ϕ
α
,
ϕ
α
=
1
i
f
α
<
1
0
i
f
α
=
1
.The temperature dependence of the reaction rate is given by the Arrhenius law [11] (6)
k
T
=
k
0
exp
⁡
-
E
R
0
T
,where k
0 is the preexponential factor, E is the activation energy, and R
0 is the universal gas constant. We assume that the two liquids are incompressible and the diffusivity coefficient is very small comparing to the thermal diffusivity coefficient such that the diffusivity will be neglected in the concentration equation.To obtain the dimensionless model, we introduce the dimensionless spatial variables as(7)
x
1
=
x
c
1
κ
,
y
1
=
y
c
1
κ
,
z
1
=
z
c
1
κ
,
t
1
=
t
c
1
2
κ
,
p
1
=
p
c
1
2
ρ
,
c
1
=
c
2
,
v
1
=
v
c
1
,
θ
=
T
-
T
b
q
,
c
2
=
2
k
0
κ
R
0
T
b
2
q
E
exp
⁡
-
E
R
0
T
b
,where c defines the stationary reaction front velocity and can be calculated asymptotically for large Zeldovich number [12]. For simplicity, we keep the same notation for the other variables and pressure. System (1) with the two boundary conditions (2)-(3) can be written in the form(8)
∂
θ
∂
t
+
v
·
∇
θ
=
Δ
θ
+
Z
exp
⁡
θ
Z
-
1
+
δ
θ
ϕ
α
,
(9)
∂
α
∂
t
+
v
·
∇
α
=
Z
exp
⁡
θ
Z
-
1
+
δ
θ
ϕ
α
,
(10)
∂
v
∂
t
+
v
·
∇
v
=
-
∇
p
+
P
Δ
v
+
P
R
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
+
θ
0
γ
,
(11)
div
⁡
v
=
0
,with the following conditions at infinity:(12)
z
⟶
+
∞
:
θ
=
-
1
,
α
=
0
,
v
=
0
,
z
⟶
-
∞
:
θ
=
0
,
α
=
1
,
v
=
0 and P
=
ν
/
κ is the Prandtl number, R
=
g
β
q
κ
2
/
(
ν
c
3
) is the Rayleigh number, Z
=
q
E
/
R
0
T
b
2 is the Zeldovich number, and δ and θ
0 are given, respectively, by δ
=
R
0
T
b
/
E and θ
0
=
(
T
b
-
T
0
)
/
q and μ
=
2
κ
σ
/
c
2.Next, we perform the linear stability analysis to tackle the interface problem.


## 3. Approximation of Infinitely Narrow Reaction Zone

To study the interface problem analytically, we use a singular perturbation analysis where the reaction zone is infinitely narrow and the reaction term is neglected outside the zone [13]. In this way, the problem can be reduced to an interface problem.To perform a formal asymptotic analysis,ϵ
=
Z
-
1 is considered as a small parameter. The new independent variable is given by z
1
=
z
-
ζ
(
x
,
y
,
t
) where ζ
(
x
,
y
,
t
) defines the location of the reaction zone. Introducing new functions θ
1, α
1, v
1, and p
1 as (13)
θ
x
,
y
,
z
,
t
=
θ
1
x
,
y
,
z
1
,
t
,
α
x
,
y
,
z
,
t
=
α
1
x
,
y
,
z
1
,
t
,
v
x
,
y
,
z
,
t
=
v
1
x
,
y
,
z
1
,
t
,
p
x
,
y
,
z
,
t
=
p
1
x
,
y
,
z
1
,
t
. The system of (8)–(11) can be written in the form (index 1 for the new function is omitted)(14)
∂
θ
∂
t
-
∂
θ
∂
z
1
∂
ζ
∂
t
+
v
·
∇
~
θ
=
Δ
~
θ
+
Z
exp
⁡
θ
Z
-
1
+
δ
θ
ϕ
α
,
∂
α
∂
t
-
∂
α
∂
z
1
∂
ζ
∂
t
+
v
·
∇
~
α
=
Z
exp
⁡
θ
Z
-
1
+
δ
θ
ϕ
α
,
∂
v
∂
t
-
∂
v
∂
z
1
∂
ζ
∂
t
+
v
·
∇
~
v
=
-
∇
~
p
+
P
Δ
~
v
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
+
θ
0
γ
,
∂
v
x
∂
x
-
∂
v
x
∂
z
1
∂
ζ
∂
x
+
∂
v
y
∂
y
-
∂
v
y
∂
z
1
∂
ζ
∂
y
+
∂
v
z
∂
z
1
=
0
,where Δ
~, ∇
~, and Q are given by (15)
Δ
~
=
∂
2
∂
x
2
+
∂
2
∂
y
2
+
∂
2
∂
z
1
2
-
2
∂
2
∂
x
∂
z
1
∂
ζ
∂
x
-
2
∂
2
∂
y
∂
z
1
∂
ζ
∂
y
+
∂
2
∂
z
1
2
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
∂
∂
z
1
∂
2
ζ
∂
x
2
+
∂
2
ζ
∂
y
2
,
∇
~
=
∂
∂
x
-
∂
∂
z
1
∂
ζ
∂
x
,
∂
∂
y
-
∂
∂
z
1
∂
ζ
∂
y
,
∂
∂
z
1
,
Q
=
P
R
.To approximate the jump conditions and resolve the interface problem, we use the matched asymptotic expansion by seeking the outer solution of problem (14) in the form (16)
θ
=
θ
0
+
ϵ
θ
1
+
⋯
,
α
=
α
0
+
ϵ
α
1
+
⋯
,
v
=
v
0
+
ϵ
v
1
+
⋯
,
p
=
p
0
+
ϵ
p
1
+
⋯
.Introducing the stretched coordinate η
=
z
1
ϵ
-
1 where ϵ
=
Z
-
1, the inner solution, can be approximated in the following form:(17)
θ
=
ϵ
θ
~
1
+
⋯
,
α
=
α
~
0
+
ϵ
α
~
1
+
⋯
,
v
=
v
~
0
+
ϵ
v
~
1
+
⋯
,
p
=
p
~
0
+
ϵ
p
~
1
+
⋯
,
ζ
=
ζ
0
+
ϵ
ζ
1
+
⋯
.Substituting the inner and outer solutions in (14) leads toOrder
ϵ
-
2:(18)
P
1
+
∂
ζ
0
∂
x
2
+
∂
ζ
0
∂
y
2
∂
2
v
~
0
∂
η
2
=
0
.
Order
ϵ
-
1:(19)
1
+
∂
ζ
0
∂
x
2
+
∂
ζ
0
∂
y
2
∂
2
θ
~
1
∂
η
2
+
exp
⁡
θ
~
1
1
+
δ
θ
~
1
ϕ
α
~
0
=
0
,
(20)
-
∂
α
~
0
∂
η
∂
ζ
0
∂
t
-
∂
α
~
0
∂
η
v
~
0
x
∂
ζ
0
∂
x
+
v
~
0
y
∂
ζ
0
∂
y
-
v
~
0
z
=
exp
⁡
θ
~
1
1
+
δ
θ
~
1
ϕ
α
~
0
,
(21)
-
∂
v
~
0
x
∂
η
∂
ζ
0
∂
t
-
v
~
0
x
∂
v
~
0
∂
η
∂
ζ
0
∂
x
-
v
~
0
y
∂
v
~
0
∂
η
∂
ζ
0
∂
y
+
v
~
0
z
∂
v
~
0
∂
η
=
t
0
∂
p
~
0
∂
η
+
P
A
∂
2
v
~
1
∂
η
2
+
F
0
∂
v
~
0
∂
η
,
(22)
-
∂
v
~
0
x
∂
η
∂
ζ
0
∂
x
-
∂
v
~
0
y
∂
η
∂
ζ
0
∂
y
+
∂
v
~
0
z
∂
η
=
0
.
Order
ϵ
0: (23)
∂
v
~
0
∂
t
-
∂
v
~
1
∂
η
∂
ζ
0
∂
t
-
∂
v
~
0
∂
η
∂
ζ
1
∂
t
+
v
~
0
x
∂
v
~
0
∂
x
-
∂
v
~
1
∂
η
∂
ζ
0
∂
x
-
∂
v
~
0
∂
η
∂
ζ
1
∂
x
+
v
~
1
x
∂
v
~
0
∂
η
∂
ζ
0
∂
x
+
v
~
1
y
∂
v
~
0
∂
η
∂
ζ
0
∂
y
+
v
~
0
y
∂
v
~
0
∂
y
-
∂
v
~
1
∂
η
∂
ζ
0
∂
y
-
∂
v
~
0
∂
η
∂
ζ
1
∂
y
+
v
~
0
z
∂
v
~
1
∂
η
+
v
~
1
z
∂
v
~
0
∂
η
=
-
∇
0
p
~
0
+
t
1
∂
p
~
0
∂
η
+
t
0
∂
p
~
1
∂
η
+
P
A
∂
2
v
~
2
∂
η
2
+
t
3
∂
2
v
~
1
∂
η
2
+
F
0
∂
v
~
1
∂
η
+
t
4
∂
2
v
~
0
∂
η
2
+
F
1
∂
v
~
0
∂
η
+
Δ
1
v
~
0
+
Q
1
+
λ
1
s
i
n
σ
1
t
+
λ
2
s
i
n
σ
2
t
γ
θ
0
,
(24)
∂
v
~
0
x
∂
x
-
∂
v
~
1
x
∂
η
∂
ζ
0
∂
x
-
∂
v
~
0
x
∂
η
∂
ζ
1
∂
x
+
∂
v
~
0
y
∂
y
-
∂
v
~
1
y
∂
η
∂
ζ
0
∂
y
-
∂
v
~
0
y
∂
η
∂
ζ
1
∂
y
+
∂
v
~
1
z
∂
η
=
0
.
Order
ϵ
1:(25)
∂
v
~
1
∂
t
-
∂
v
~
2
∂
η
∂
ζ
0
∂
t
+
∂
v
~
1
∂
η
∂
ζ
1
∂
t
+
∂
v
~
0
∂
η
∂
ζ
2
∂
t
+
v
~
0
x
∂
v
~
1
∂
x
-
∂
v
~
1
∂
η
∂
ζ
1
∂
x
-
∂
v
~
0
∂
η
∂
ζ
2
∂
x
-
∂
v
~
2
∂
η
∂
ζ
0
∂
x
+
v
~
1
x
∂
v
~
0
∂
x
-
∂
v
~
0
∂
η
∂
ζ
1
∂
x
-
∂
v
~
1
∂
η
∂
ζ
0
∂
x
-
v
~
2
x
∂
v
~
0
∂
η
∂
ζ
0
∂
x
+
v
~
0
y
∂
v
~
1
∂
y
-
∂
v
~
1
∂
η
∂
ζ
1
∂
y
-
∂
v
~
0
∂
η
∂
ζ
2
∂
y
-
∂
v
~
2
∂
η
∂
ζ
0
∂
y
+
v
~
1
y
∂
v
~
0
∂
y
-
∂
v
~
0
∂
η
∂
ζ
1
∂
y
-
∂
v
~
1
∂
η
∂
ζ
0
∂
y
-
v
~
2
y
∂
v
~
0
∂
η
∂
ζ
0
∂
y
+
v
~
0
z
∂
v
~
2
∂
η
+
v
~
1
z
∂
v
~
1
∂
η
+
v
~
2
z
∂
v
~
0
∂
η
=
t
0
∂
p
~
2
∂
η
+
t
1
∂
p
~
1
∂
η
-
∇
0
p
~
1
+
t
2
∂
p
~
0
∂
η
+
P
A
∂
2
v
~
3
∂
η
2
+
t
3
∂
2
v
~
2
∂
η
2
+
F
0
∂
v
~
2
∂
η
+
F
1
∂
v
~
1
∂
η
+
Δ
1
v
~
1
+
F
2
∂
v
~
0
∂
η
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
0
γ
,
(26)
∂
v
~
1
x
∂
x
-
∂
v
~
2
x
∂
η
∂
ζ
0
∂
x
-
∂
v
~
1
x
∂
η
∂
ζ
1
∂
x
-
∂
v
~
0
x
∂
η
∂
ζ
2
∂
x
+
∂
v
~
1
y
∂
y
-
∂
v
~
2
y
∂
η
∂
ζ
0
∂
y
-
∂
v
~
1
y
∂
η
∂
ζ
1
∂
y
-
∂
v
~
0
y
∂
η
∂
ζ
2
∂
y
+
∂
v
~
2
z
∂
η
=
0
.
Order
ϵ
2:(27)
∂
v
~
2
x
∂
x
-
∂
v
~
3
x
∂
η
∂
ζ
0
∂
x
-
∂
v
~
2
x
∂
η
∂
ζ
1
∂
x
-
∂
v
~
1
x
∂
η
∂
ζ
2
∂
x
-
∂
v
~
0
x
∂
η
∂
ζ
3
∂
x
+
∂
v
~
2
y
∂
y
-
∂
v
~
3
y
∂
η
∂
ζ
0
∂
y
-
∂
v
~
2
y
∂
η
∂
ζ
1
∂
y
-
∂
v
~
1
y
∂
η
∂
ζ
2
∂
y
-
∂
v
~
0
y
∂
η
∂
ζ
3
∂
y
+
∂
v
~
3
z
∂
η
=
0
,where (28)
Δ
1
=
∂
2
∂
x
2
+
∂
2
∂
y
2
,
∇
0
=
∂
∂
x
,
∂
∂
y
,
0
,
F
i
=
-
2
∂
ζ
i
∂
x
∂
∂
x
+
∂
ζ
i
∂
y
∂
∂
y
-
∂
2
ζ
i
∂
x
2
+
∂
2
ζ
i
∂
y
2
I
,
i
=
0,1
,
2
,
t
0
=
∂
ζ
0
∂
x
,
∂
ζ
0
∂
y
,
-
1
,
t
i
=
∂
ζ
i
∂
x
,
∂
ζ
i
∂
y
,
-
1
i
=
1,2
,
t
3
=
2
∂
ζ
1
∂
x
∂
ζ
0
∂
x
+
∂
ζ
0
∂
y
∂
ζ
1
∂
y
,
t
4
=
∂
ζ
1
∂
x
2
+
∂
ζ
1
∂
y
2
+
∂
ζ
0
∂
x
∂
ζ
2
∂
x
+
∂
ζ
0
∂
y
∂
ζ
2
∂
y
,
A
=
1
+
∂
ζ
0
∂
x
2
+
∂
ζ
0
∂
y
2and I is the identity operator. The matching conditions as η
→
±
∞ are given by(29)
v
~
0
~
v
0
z
1
=
±
0
,
(30)
v
~
1
~
∂
v
0
∂
z
1
z
1
=
±
0
η
+
v
1
z
1
=
±
0
,
(31)
v
~
2
~
1
2
∂
2
v
0
∂
z
1
2
z
1
=
±
0
η
2
+
∂
v
1
∂
z
1
z
1
=
±
0
η
+
v
2
z
1
=
±
0
,
(32)
v
~
3
~
1
6
∂
3
v
0
∂
z
1
3
z
1
=
±
0
η
3
+
1
2
∂
2
v
1
∂
z
1
2
z
1
=
±
0
η
2
+
∂
v
2
∂
z
1
z
1
=
±
0
η
+
v
3
z
1
=
±
0
.As η
→
+
∞
(33)
θ
~
1
≈
θ
1
z
1
=
+
0
+
∂
θ
0
∂
z
1
z
1
=
+
0
η
,
α
~
0
⟶
0
.As η
→
-
∞
(34)
θ
~
1
~
θ
1
z
1
=
-
0
,
α
~
0
⟶
1
.From (18), we obtain (35)
∂
2
v
~
0
∂
η
2
=
0
.Consequently, we can conclude thatv
~
0
(
η
) is a linear function of η and identically constant since the velocity is bounded. Equation (29) becomes (36)
v
0
z
1
=
+
0
=
v
0
z
1
=
-
0
,
(37)
∂
v
~
0
∂
η
=
0
.From (36)-(37), we deduce that the first term in the expression of the velocity v
0 is continuous at the front. By substituting (37) into (21), we find (38)
P
A
∂
2
v
~
1
∂
η
2
+
t
0
∂
p
~
0
∂
η
=
0
,
(39)
∂
v
~
0
x
∂
x
-
∂
v
~
1
x
∂
η
∂
ζ
0
∂
x
+
∂
v
~
0
y
∂
y
-
∂
v
~
1
y
∂
η
∂
ζ
0
∂
y
+
∂
v
~
1
z
∂
η
=
0
.Differentiating (39) with respect to η, one obtains (40)
∂
2
v
~
1
x
∂
η
2
∂
ζ
0
∂
x
+
∂
2
v
~
1
y
∂
η
2
∂
ζ
0
∂
y
-
∂
2
v
~
1
z
∂
η
2
=
0
.As a result, (38) is a vectorial equation with three components. We multiply the first component by ∂
ζ
0
/
∂
x, the second by ∂
ζ
0
/
∂
y, and the third by -
1 and adding we have (41)
P
A
∂
2
v
~
1
x
∂
η
2
∂
ζ
0
∂
x
+
∂
2
v
~
1
y
∂
η
2
∂
ζ
0
∂
y
-
∂
2
v
~
1
z
∂
η
2
+
A
∂
p
~
0
∂
η
=
0
.Considering the two equations (38) and (40), we obtain (42)
∂
p
~
0
∂
η
=
0
,
(43)
∂
2
v
~
1
∂
η
2
=
0
.From the boundedness of the velocity and (30), the continuity of the first order term of the outer expansion and of the first derivative of the zero-order term, we can write(44)
∂
v
0
∂
z
1
z
1
=
+
0
=
∂
v
0
∂
z
1
z
1
=
-
0
,
v
1
z
1
=
+
0
=
v
1
z
1
=
-
0
.Differentiating (23) once and (26) twice with respect to η and using the three equations (37), (42), and (43), we get (45)
t
0
∂
2
p
~
1
∂
η
2
+
P
A
∂
3
v
~
2
∂
η
3
=
0
,
(46)
∂
3
v
~
2
x
∂
η
3
∂
ζ
0
∂
x
+
∂
3
v
~
2
y
∂
η
3
∂
ζ
0
∂
y
-
∂
3
v
~
2
z
∂
η
3
=
0
.As above, we multiply the three components of (45), respectively, by ∂
ζ
0
/
∂
x, ∂
ζ
0
/
∂
y, and -
1 and adding, we get (47)
P
A
∂
3
v
~
2
x
∂
η
3
∂
ζ
0
∂
x
+
∂
3
v
~
2
y
∂
η
3
∂
ζ
0
∂
y
-
∂
3
v
~
2
z
∂
η
3
+
A
∂
2
p
~
1
∂
η
2
=
0
.From the previous equation and using (45)-(46), we have (48)
∂
2
p
~
1
∂
η
2
=
0
,
∂
3
v
~
2
∂
η
3
=
0
.Knowing that the velocity is bounded and considering (31), we obtain the jump conditions in the following forms: (49)
∂
2
v
0
∂
z
1
2
z
1
=
+
0
=
∂
2
v
0
∂
z
1
2
z
1
=
-
0
,
∂
v
1
∂
z
1
z
1
=
+
0
=
∂
v
1
∂
z
1
z
1
=
-
0
,
v
2
z
1
=
+
0
=
v
2
z
1
=
-
0
.Differentiating (25) twice and (27) three times with respect to η and taking into account (37), (42), (48), we obtain (50)
t
0
∂
3
p
~
2
∂
η
3
+
P
A
∂
4
v
~
3
∂
η
4
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
γ
∂
2
θ
~
1
∂
η
2
=
0
,
∂
4
v
~
3
x
∂
η
4
∂
ζ
0
∂
x
+
∂
4
v
~
3
y
∂
η
4
∂
ζ
0
∂
y
-
∂
4
v
~
3
z
∂
η
4
=
0
,
(51)
P
A
∂
4
v
~
3
x
∂
η
4
∂
ζ
0
∂
x
+
∂
4
v
~
3
y
∂
η
4
∂
ζ
0
∂
y
-
∂
4
v
~
3
z
∂
η
4
+
A
∂
3
p
~
2
∂
η
3
-
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
∂
2
θ
~
1
∂
η
2
=
0
.From the last equation and using (50), we can write (52)
A
∂
3
p
~
2
∂
η
3
-
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
∂
2
θ
~
1
∂
η
2
=
0
,
(53)
γ
0
∂
2
θ
~
1
∂
η
2
=
∂
4
v
~
3
∂
η
4
,where γ
0 is defined by (54)
γ
0
=
-
∂
ζ
0
∂
x
R
A
2
,
-
∂
ζ
0
∂
y
R
A
2
,
R
A
2
-
R
A
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
.Now, integrating (53) with respect to η using the system (32)–(34) lead to (55)
∂
3
v
0
∂
z
1
3
z
1
=
-
0
-
∂
3
v
0
∂
z
1
3
z
1
=
+
0
=
-
γ
0
∂
θ
0
∂
z
1
z
1
=
+
0
.As a first conclusion, the velocity jump conditions across the reaction front are given by (36), (44), (49), and (55).From (20), (22), (37), we conclude that α
~
0 is a monotonic function satisfying 0
<
α
~
0
<
1. Therefore, the reaction is of zero order and ϕ
(
α
~
0
)
≡
1. Multiplying (19) by ∂
θ
~
1
/
∂
η and integrating the result yield (56)
∂
θ
~
1
∂
η
2
+
∞
-
∂
θ
~
1
∂
η
2
-
∞
=
2
A
-
1
∫
-
∞
θ
1
∣
z
1
=
-
0
exp
⁡
τ
1
+
τ
δ
d
τ
.Subtracting (19) from (20) and integrating the result one obtains (57)
∂
θ
~
1
∂
η
+
∞
-
∂
θ
~
1
∂
η
-
∞
=
-
A
-
1
∂
ζ
0
∂
t
+
s
,where s is given by (58)
s
=
v
~
0
x
∂
ζ
0
∂
x
+
v
~
0
y
∂
ζ
0
∂
y
-
v
~
0
z
.As a second conclusion, the temperature jump conditions across the reaction front are given by the two last equations (56)-(57). By using the matching conditions above and truncating the expansion as (59)
θ
≈
θ
0
,
θ
1
z
1
=
-
0
≈
Z
θ
z
1
=
+
0
,
ζ
≈
ζ
0
,
v
≈
v
0
,the jump conditions can be written as(60)
∂
θ
∂
z
1
2
+
0
-
∂
θ
∂
z
1
2
-
0
=
2
Z
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
∫
-
∞
θ
1
z
1
=
-
0
exp
⁡
τ
Z
-
1
+
τ
δ
d
τ
,
∂
θ
∂
z
1
z
1
=
+
0
-
∂
θ
∂
z
1
z
1
=
-
0
=
-
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
∂
ζ
∂
t
+
v
x
∂
ζ
∂
x
+
v
y
∂
ζ
∂
y
-
v
z
z
1
=
+
0
,
v
z
z
1
=
+
0
=
v
z
z
1
=
-
0
,
∂
v
z
∂
z
1
z
1
=
+
0
=
∂
v
z
∂
z
1
z
1
=
-
0
,
∂
2
v
z
∂
z
1
2
z
1
=
+
0
=
∂
v
z
∂
z
1
2
z
1
=
-
0
,
∂
3
v
z
∂
z
1
3
z
1
=
-
0
-
∂
3
v
z
∂
z
1
3
z
1
=
+
0
=
-
R
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
-
1
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
∂
θ
∂
z
1
z
1
=
+
0
.


## 4. The Interface Problem and Perturbation

To study the propagation of polymerization front with a liquid reactant and liquid product, the equation of motion has been considered before and after the reaction zone. This changes the jump conditions and influences the stability conditions of the frontal polymerization process. Notice that in the case of liquid-solid polymerization front, the equation of motion is considered after the reaction zone. The original system (2), (3), (8), (9) leads to the following interface problem:In the liquid monomerz
>
ζ
(61)
∂
θ
∂
t
+
v
·
∇
θ
=
Δ
θ
,
(62)
α
=
0
,
(63)
∂
v
∂
t
+
v
·
∇
v
=
-
∇
p
+
P
Δ
v
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
+
θ
0
γ
,
(64)
div
⁡
v
=
0
.
In the liquid polymerz
<
ζ
(65)
∂
θ
∂
t
+
v
·
∇
θ
=
Δ
θ
,
(66)
α
=
1
,
(67)
∂
v
∂
t
+
v
·
∇
v
=
-
∇
p
+
P
Δ
v
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
+
θ
0
γ
,
(68)
div
⁡
v
=
0
.
At the interfacez
=
ζ
(69)
θ
ζ
-
0
=
θ
ζ
+
0
,
∂
θ
∂
z
ζ
=
-
0
-
∂
θ
∂
z
ζ
=
+
0
=
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
∂
ζ
∂
t
+
v
x
∂
ζ
∂
x
+
v
y
∂
ζ
∂
y
-
v
z
ζ
,
∂
θ
∂
z
2
ζ
-
0
-
∂
θ
∂
z
2
ζ
1
+
0
=
-
2
Z
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
∫
-
∞
θ
ζ
exp
⁡
τ
Z
-
1
+
δ
τ
d
τ
,
v
z
ζ
=
-
0
=
v
z
ζ
=
+
0
,
∂
v
z
∂
z
ζ
=
-
0
=
∂
v
z
∂
z
ζ
=
+
0
,
∂
2
v
z
∂
z
2
ζ
=
-
0
=
∂
2
v
z
∂
z
2
ζ
=
+
0
,
∂
3
v
z
∂
z
3
ζ
=
-
0
-
∂
3
v
z
∂
z
3
ζ
=
+
0
=
-
R
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
1
+
∂
ζ
∂
x
2
+
∂
ζ
∂
y
2
-
1
-
1
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
∂
θ
∂
z
ζ
=
+
0
,with the conditions at infinity (70)
z
=
-
∞
:
θ
=
0
,
v
=
0
,
z
=
+
∞
:
θ
=
-
1
,
v
=
0
.


## 5. Stability Analysis

The interface problem analyzed in Section4 has a travelling wave solution in the following form: (71)
θ
x
,
y
,
z
,
t
,
α
x
,
y
,
z
,
t
,
v
=
θ
s
z
-
u
t
,
α
s
z
-
u
t
,
0
,
θ
s
z
-
u
t
,
α
s
z
-
u
t
=
0,1
,
z
2
<
0
e
x
p
-
u
z
2
-
1,0
,
z
2
>
0
,where z
2
=
z
-
u
t and u is the speed of the stationary reaction front. This solution is a basic and stationary solution of interface problem (61)–(70) written in the moving coordinates. Then, (61), (63), (65), and (67) can be replaced by (72)
∂
θ
∂
t
+
v
·
∇
θ
=
Δ
θ
+
u
∂
θ
∂
z
2
,
∂
v
∂
t
+
v
∇
v
=
-
∇
p
+
P
Δ
v
+
u
∂
θ
∂
z
2
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
+
θ
0
γ
.The other equations remain unchanged. To study the stability of reaction front, we seek the solution of the problem as follows:(73)
θ
=
θ
s
+
θ
~
,
p
=
p
s
+
p
~
,
v
=
v
s
+
v
~
,where θ
~, p
~, and v
~ are, respectively, small perturbation of temperature, pressure, and velocity.Subtituting (73) into (64), (72), one obtains the first-order terms as follows:In the liquid monomerz
2
>
ξ
(74)
∂
θ
~
∂
t
=
Δ
θ
~
+
u
∂
θ
~
∂
z
2
-
v
~
z
θ
s
′
,
(75)
∂
v
~
∂
t
=
-
∇
p
~
+
P
Δ
v
~
+
u
∂
θ
~
∂
z
2
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
~
γ
,
(76)
div
⁡
v
~
=
0
.
In the liquid polymerz
2
<
ξ
(77)
∂
θ
~
∂
t
=
Δ
θ
~
+
u
∂
θ
~
∂
z
2
,
(78)
∂
v
~
∂
t
=
-
∇
p
~
+
P
Δ
v
~
+
u
∂
v
~
∂
z
2
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
~
γ
,
(79)
div
⁡
v
~
=
0
,where ξ is defined by ξ
=
ζ
-
u
t.We note(80)
θ
~
,
v
~
z
=
θ
^
1
,
v
^
z
1
f
o
r
z
2
<
ξ
θ
^
2
,
v
^
z
2
f
o
r
z
2
<
ξand consider the perturbation as follows: (81)
θ
^
i
=
θ
i
z
2
,
t
exp
⁡
j
k
1
x
+
k
2
y
,
v
^
z
i
=
v
z
i
z
2
,
t
exp
⁡
j
k
1
x
+
k
2
y
,
ξ
=
ϵ
1
t
exp
⁡
j
k
1
x
+
k
2
y
,where k
i
,
i
=
1,2, and ϵ
1 denote, respectively, the wave numbers (in x and y directions) and the amplitude of the perturbation and j
2
=
-
1. Then, jump conditions (69) are linearized by taking into account the fact that (82)
θ
ξ
=
±
0
=
θ
s
±
0
+
ξ
θ
s
′
±
0
+
θ
~
±
0
,
∂
θ
∂
z
2
ξ
=
±
0
=
θ
s
′
±
0
+
ξ
θ
s
′
′
±
0
+
∂
θ
~
∂
z
2
ξ
=
±
0
.Then the higher-order terms are given as (83)
θ
^
2
z
2
=
0
-
θ
^
1
z
2
=
0
=
u
ξ
,
(84)
∂
θ
^
2
∂
z
2
z
2
=
0
-
∂
θ
^
1
∂
z
2
z
2
=
0
=
-
u
2
ξ
-
∂
ξ
∂
t
+
v
~
z
z
2
=
0
,
(85)
u
2
ξ
+
∂
θ
^
2
∂
z
2
z
2
=
0
=
-
Z
u
θ
^
1
z
2
=
0
,
(86)
v
^
2
z
z
2
=
0
=
v
^
1
z
z
2
=
0
,
(87)
∂
v
^
z
2
∂
z
2
z
2
=
0
=
∂
v
^
1
z
∂
z
2
z
2
=
0
,
(88)
∂
2
v
^
2
z
∂
z
2
2
z
2
=
0
=
∂
2
v
^
1
z
∂
z
2
2
z
2
=
0
,
(89)
∂
3
v
^
2
z
∂
z
2
3
z
2
=
0
=
∂
3
v
^
1
z
∂
z
2
3
z
2
=
0
.Substituting (81) into (83)–(89) leads to (90)
θ
2
0
,
t
-
θ
1
0
,
t
=
u
ϵ
1
t
,
θ
2
′
0
,
t
-
θ
1
′
0
,
t
=
-
ϵ
1
t
u
2
-
ϵ
1
′
t
+
v
1
0
,
t
,
ϵ
1
t
u
2
+
θ
2
′
0
,
t
=
-
Z
u
θ
1
0
,
t
,
v
2
i
0
,
t
=
v
1
i
0
,
t
i
=
0,1
,
2,3
,where (91)
v
i
=
∂
i
v
∂
v
i
,
θ
′
=
∂
θ
i
∂
z
,
ϵ
1
′
t
=
d
ϵ
1
t
d
t
.After simplification, (75) and (78) become (92)
∂
∂
t
Δ
v
~
z
-
u
∂
∂
z
2
Δ
v
~
z
=
P
Δ
Δ
v
~
z
+
Q
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
∂
2
∂
x
2
+
∂
2
∂
y
2
θ
~
.Substituting (71), (87), (88) into (74), (77), (92) and introducing the variable w
=
v
′
′
-
k
2
v, where k
=
k
1
2
+
k
2
2, one obtains two systems of equations as follows:z
2
<
ξ:(93)
∂
w
∂
t
-
u
w
′
-
P
w
′
′
-
k
2
w
=
-
Q
k
2
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
,
w
=
v
′
′
-
k
2
v
,
∂
θ
∂
t
-
θ
′
′
-
u
θ
′
+
k
2
θ
=
0
.
z
2
>
ξ:(94)
∂
w
∂
t
-
u
w
′
-
P
w
′
′
-
k
2
w
=
-
Q
k
2
1
+
λ
1
sin
⁡
σ
1
t
+
λ
2
sin
⁡
σ
2
t
θ
,
w
=
v
′
′
-
k
2
v
,
∂
θ
∂
t
-
θ
′
′
-
u
θ
′
+
k
2
θ
=
u
exp
⁡
-
u
z
2
v
,with the following boundary conditions: (95)
θ
2
0
,
t
-
θ
1
0
,
t
=
u
ϵ
1
t
,
(96)
θ
2
′
0
,
t
-
θ
1
′
0
,
t
=
-
u
2
ε
1
t
-
ε
1
′
t
+
v
1
0
,
t
,
(97)
θ
2
′
0
,
t
+
Z
u
θ
1
0
,
t
=
-
u
2
ε
1
t
,
(98)
v
1
i
0
,
t
=
v
2
i
0
,
t
,
i
=
0,1
,
2,3
.Combining the three boundary conditions (95)–(97) yields (99)
θ
1
′
0
,
t
+
Z
u
θ
1
0
,
t
=
1
u
θ
2
0
,
t
-
θ
1
0
,
t
t
′
-
v
1
0
,
t
.


## 6. Numerical Results

The numerical simulations of the problem are performed using the finite-difference approximation with implicit scheme. The algorithm is given in Appendix.Evolution of maximum temperature versus time provides the onset of stability of the polymerization front. The convective instability occurs when a jump from bounded to unbounded values of maximum temperature is achieved.Figure1 shows the critical Rayleigh number as a function of the frequency σ
1 for different frequencies ratio and for σ
2
=
0. In the absence of gravitational modulation (σ
1
=
σ
2
=
0), we find the same result (R
c
≈
27) as in [8, 14]. It can be observed that the frequencies ratio of the QP gravitational modulation has a significant influence on the convective instability of the reaction front. As indicated, increasing the frequencies ratio, the stability domain becomes larger. In particular, this can be seen in the range of the frequency σ
1 belonging to the interval (
0,10
). For large values of the frequency σ
1, the frequencies ratio σ
2
/
σ
1 has no effect on the convective instability. It is also worthy to notice that the QP gravitational modulation has a stabilizing effect for sufficiently large frequencies ratio (σ
2
=
37
σ
1) comparing to the periodic modulation case (σ
2
=
0).Figure 1
Critical Rayleigh number as a function ofσ
1 for different frequencies ratio and for σ
2
=
0, λ
1
=
λ
2
=
5, P
=
0.5, Z
=
8, and k
=
0.7.The critical Rayleigh number as a function of the amplitudeλ
2 is shown in Figure 2 for different values of λ
1. We can conclude that, for a fixed value of λ
2 and increased values of λ
1, the polymerization front loses its stability monotonically revealing that the stability of the reaction front can be controlled by acting on the amplitudes of the QP modulation.Figure 2
Critical Rayleigh number as a function ofλ
2 for different values of λ
1, σ
2
=
2
σ
1, σ
1
=
5, P
=
0.5, Z
=
8, and k
=
0.7.Figure3 illustrates the critical Rayleigh number as a function of the amplitude λ
2 for different frequencies ratio. It can be observed that increasing the frequencies ratio the stability domain becomes large for certain values of λ
2. For large values of λ
2 (λ
2
≥
100), the effect of the frequencies ratio on the convective instability boundary is insignificant.Figure 3
Critical Rayleigh number as a function ofλ
2 for different frequencies ratio, k
=
0.7, λ
1
=
5, P
=
0.5, Z
=
8, and σ
1
=
5.Figure4 shows for a fixed frequencies ratio, σ
2
/
σ
1
=
7, the influence of the wave number on the critical Rayleigh number versus λ
2. The plots indicate that for large values of the wave number the reaction front gains stability.Figure 4
Critical Rayleigh number as a function ofλ
2 for different values of number wave, σ
2
=
7
σ
1, σ
1
=
10, P
=
0.5, Z
=
8, and λ
1
=
5.The influence ofλ
2 on the stability of the reaction front is given in Figure 5 for two different frequencies ratio. From this figure, one can conclude that an increase of the frequencies ratio has a stabilizing effect.Figure 5
Critical Rayleigh number as a function ofλ
2 for two different frequencies ratios, k
=
0.7, λ
1
=
5, P
=
0.5, Z
=
8, and σ
1
=
10.


## 7. Conclusion

In this work, we have studied the influence of the QP vibrations on the convective instability of liquid-liquid polymerization front. We have used the model which includes the heat equation, the concentration equation, and the Navier-Stokes equations under Boussinesq approximation. The Zeldovich Frank-Kamenetskii method has been used assuming that the reaction occurs in a narrow zone. To obtain the convective instability threshold using the linear stability analysis, the reduced system of equations has been discretized and resolved using the finite-difference method with implicit scheme. The results have shown that, for fixed values of amplitudes, an increase of the frequencies ratio stabilizes the reaction front, and for large values of the frequencyσ
1 the critical Rayleigh number tends to the unmodulated critical value. In addition, it is observed that for a fixed frequencies ratio and for a given amplitude λ
2 an increase of the amplitude λ
1 destabilizes the reaction front. Also, the reaction front becomes more stable by increasing the wave number.



---

*Source: 484562-2015-12-20.xml*