3D Long-Wave 3D Long-Wave Oscillatory Patterns Oscillatory Patterns in Thermocapillary in Thermocapillary Convection with Convection with Soret Effect Soret Effect A. Nepomnyashchy, A. Oron A. Nepomnyashchy, A. Oron Technion, Haifa, Israel, Technion, Haifa, Israel, and and S. Shklyaev S. Shklyaev , , Technion, Haifa, Israel, Technion, Haifa, Israel, Perm State University, Russia Perm State University, Russia
33
Embed
3D Long-Wave Oscillatory Patterns in Thermocapillary Convection with Soret Effect
3D Long-Wave Oscillatory Patterns in Thermocapillary Convection with Soret Effect. A. Nepomnyashchy, A. Oron Technion, Haifa, Israel, and S. Shklyaev , Technion, Haifa, Israel, Perm State University, Russia. This work is supported by the Israel Science Foundation - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
3D Long-Wave Oscillatory 3D Long-Wave Oscillatory Patterns in Thermocapillary Patterns in Thermocapillary
Convection with Soret Convection with Soret EffectEffect
A. Nepomnyashchy, A. OronA. Nepomnyashchy, A. Oron
Technion, Haifa, Israel,Technion, Haifa, Israel,
and and S. ShklyaevS. Shklyaev,,
Technion, Haifa, Israel,Technion, Haifa, Israel,
Perm State University, RussiaPerm State University, Russia
22
This work is supported by the Israel This work is supported by the Israel Science FoundationScience Foundation
I am grateful to Isaac Newton I am grateful to Isaac Newton Institute for the invitation and for the Institute for the invitation and for the
Double-diffusive Marangoni convection:• J.L. Castillo and M.G. Velarde, JFM (1982);
• C.L. McTaggart, JFM (1983).
Linear stability problem with Soret effect:• C.F. Chen, C.C. Chen, Phys. Fluids (1994);
• J.R.L. Skarda, D.Jackmin, and F.E. McCaughan, JFM (1998).
55
Nonlinear analysis of long-wave Nonlinear analysis of long-wave perturbationsperturbations
Marangoni convection in pure liquids:• E. Knobloch, Physica D (1990);
• A.A. Golovin, A.A. Nepomnyashchy,nd L.M. Pismen, Physica D (1995);
Marangoni convection in solutions:• L. Braverman, A. Oron, J. Eng. Math. (1997);
• A. Oron and A.A. Nepomnyashchy, Phys. Rev. E (2004).
Oscillatory mode in Rayleigh-Benard convection• L.M. Pismen, Phys. Rev. A (1988).
66
Basic assumptionsBasic assumptions Gravity is negligible;Gravity is negligible; Free surface is nondeformable;Free surface is nondeformable; Surface tension linearly depends on both Surface tension linearly depends on both
the temperature and the concentration:the temperature and the concentration: 0 0 0 ;T CT T C C
D C T j
0.zkT q T T
The heat flux is fixed at the rigid plate;The heat flux is fixed at the rigid plate; The Newton law of cooling governs the The Newton law of cooling governs the heat transfer at the free surface:heat transfer at the free surface:
Soret effect plays an important role:Soret effect plays an important role:
77
Governing equationsGoverning equations
1 2
2
1 2 2
,
,
,
div 0.
P pt
TP T T
t
CS L C C T
t
vv v v
v
v
v
88
Boundary conditionsBoundary conditions
2
1: 0, 0, 0,
0.
z z
z
z w T BT C BT
M T C
u
At the rigid wall:
0 : 0, 1, ;z zz T C v
At the interface:
,zw v u eHere
is the differential operator in plane x-y 2
99
Dimensionless parametersDimensionless parameters
The Prandtl numberThe Prandtl number
The Schmidt numberThe Schmidt number
The Soret numberThe Soret number
The Marangoni numberThe Marangoni number
The Biot numberThe Biot number
The Lewis numberThe Lewis number
SD
C
T
2TAHM
P
PL
S
qHB
k
1010
Basic stateBasic state
0 0
0
0
0, ,
1,
.
p const
BT z
BC z const
v
There exist the equilibrium state corresponding to the linear temperature and concentration distribution:
1111
Equation for perturbationsEquation for perturbations
1 2
1 2
2
1 2 2
0
0 : 0;
1: 0, 0, 0,
0
t z zz
t z z zz
t z zz
t z zz zz
z
z z
z z
z
P w
w P w ww w w
P w w
S L w
w
z w
z w B B
M
u u u u u u
u
u
u
u
u
u
are the perturbations of the pressure, the temperature and the concentration, respectively; here and below 2
1212
Previous resultsPrevious results
Linear stability problem was studied;Linear stability problem was studied; Monotonous mode was found and weakly Monotonous mode was found and weakly
nonlinear analysis was performed;nonlinear analysis was performed; Oscillatory mode was revealed;Oscillatory mode was revealed; The set of amplitude equations to study The set of amplitude equations to study
2D oscillatory convective motion was 2D oscillatory convective motion was obtained.obtained.
Linear and nonlinear stability analysis of above conductive state with respect to long-wave perturbations was carried out by A.Oron and A.Nepomnyashchy (PRE, 2004):
1313
Multi-scale expansion for the Multi-scale expansion for the analysis of long wave perturbationsanalysis of long wave perturbations
, ,X x Y y Z z
2,W w U u
Rescaled coordinates:
“Slow” times :
Rescaled components of the velocity:
2 4,T t t
1414
Multi-scale expansion for the Multi-scale expansion for the analysis of long wave perturbationsanalysis of long wave perturbations
20 2
0 2 0 22 2
0 22
0 2 0 22 2
,
, ,
,
,
M M M
W W W
U U U
Expansion with respect to
Small Biot number:4B
1515
The zeroth order solutionThe zeroth order solution
0 0
00
0 0 2 20 0
, , , , , , ,
3, ,
21 1
3 2 , 1 ,4 4
F X Y T G X Y T
M h h F G
M Z Z h W M Z Z h
U
X
Z
1616
The second orderThe second order
2 20 0
1 2 1 20 0
1 ,
1 1
T
T
PF m F m G
SG m L F m L G
The solvability conditions:
The plane wave solution:
0
exp . .,
exp . .,1
h A i i T c c
mF A i i T c c G F h
i P
k R
k R
1717
The second orderThe second order
The dispersion relation:
2 22
2
2
1, ,
1
11
L L Lk P
LL L
00
148 1M L
m
Critical Marangoni number:
The solution of the second order:
2 2
2 2
, , , , , , ,
, , , , , , ,
z X Y T Q X Y T
z X Y T R X Y T
1818
The fourth orderThe fourth orderThe solvability conditions:
2 2 220 0
4 1002
2 222 20 0
2 2
2 20 0
2
1
2 160
2 div10 10
48 312div , ,
35 35
T
mm mF Q Q R F h
P P P Pm
m LP
m mP h h h
P Pm m
h F h J FP P
1919
2 2 2200
4 1002
2 222 20 0
2 2
2 210 0
2
2 160
2 div10 1048 312
div , ,35 35
T
mL mG R L m Q R h
P P Pm
m LP
m mS P h h h
P Pm m
h F L G h J GP P
0
1
3 2 1 3 1 3 ,
6 1 ,
, .X Y Y X
P m h L F LG
P h F L G
J f g f g f g
2 2 2 .Y X X Yh h h h
2020
2 20 0
1 2 1 20 0
1 ,
1 1
T
T
PF m F m G
SG m L F m L G
2121
Linear stability analysisLinear stability analysisOron, Nepomnyashchy, PRE, 2004