Linear stability analysis of time-dependent algorithms with spectral element methods for the simulation of viscoelastic flows

2003 ◽  
Vol 115 (2-3) ◽  
pp. 157-190 ◽  
Author(s):  
Nicolas Fiétier ◽  
Michel O. Deville
Author(s):  
Elia Merzari ◽  
Paul Fischer ◽  
W. David Pointer

Buoyancy-driven systems are subject to several types of flow instabilities. To evaluate the performance of such systems it is becoming increasingly crucial to be able to predict the stability of a given base flow configuration. Traditional Modal Linear stability Analysis requires the solution of very large eigenvalue systems for three-dimensional flows, which make this problem difficult to tackle. An alternative to modal Linear stability Analysis is the use of adjoint solvers [1] in combination with a power iteration [2]. Such methodology allows for the identification of an optimal disturbance or forcing and has been recently used to evaluate the stability of several isothermal flow systems [2]. In this paper we examine the extension of the methodology to non-isothermal flows driven by buoyancy. The contribution of buoyancy in the momentum equation is modeled through the Boussinesq approximation. The method is implemented in the spectral element code Nek5000. The test case is the flow is a two-dimensional cavity with differential heating and conductive walls and the natural circulation flow in a toroidal thermosiphon.


2019 ◽  
Author(s):  
Edwin Khoo ◽  
Hongbo Zhao ◽  
Martin Z. Bazant

We study the linear stability analysis of time-dependent electrodeposition in a charged random porous medium, whose pore surface charges can generally be of any sign, that is flanked by a pair of planar metal electrodes. Discretization of the linear stability problem results in a generalized eigenvalue problem for the dispersion relation that is solved numerically, which agrees well with the analytical approximation obtained from a boundary layer analysis valid at high wavenumbers. Under galvanostatic conditions in which an overlimiting current is applied, in the classical case of zero pore surface charges, the voltage and electric field at the cathode diverge when the bulk electrolyte concentration there vanishes at Sand's time. The same phenomenon happens for positive surface charges but at a time earlier than Sand's time. In contrast, negative surface charges allow the electrochemical system to sustain an overlimiting current via surface conduction past Sand's time, keeping the voltage and electric field bounded. Therefore, at Sand's time, negative surface charges greatly reduce the electrode surface instabilities while zero and positive surface charges magnify them. We compare theoretical predictions for overall electrode surface stabilization from the linear stability analysis with published experimental data for copper electrodeposition in cellulose nitrate membranes and demonstrate good agreement between theory and experiment. We also use the linear stability analysis as a tool to analyze how the crystal grain size changes with duty cycle during pulse electroplating.


Sign in / Sign up

Export Citation Format

Share Document