Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userDirect Numerical Simulation of Electrokinetic Instability and Transition to Chaotic Motion
118
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A new type of instability - electrokinetic instability - and an unusual transition to chaotic motion near a charge-selective surface was studied by numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. Two kinds of initial conditions were considered: (a) white noise initial conditions to mimic room disturbances and subsequent natural evolution of the solution; (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution; two-dimensional steady electroconvective vortices (stationary point in a proper phase space); unsteady vortices aperiodically changing their parameters (homoclinic contour); periodic motion (limit cycle); and chaotic motion. The transition to chaotic motion did not include Hopf bifurcation. Numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt-current characteristics.
rate research
Read More
A direct numerical simulation of the three-dimensional elektrokinetic instability near a charge selective surface (electric membrane, electrode, or system of micro-/nanochannels) is carried out and analyzed. A special finite-difference method was used for the space discretization along with a semi-implicit $3frac{1}{3}$-step Runge-Kutta scheme for the integration in time. The calculations employed parallel computing. Three characteristic patterns, which correspond to the overlimiting currents, are observed: (a) two-dimensional electroconvective rolls, (b) three-dimensional regular hexagonal structures, and (c) three-dimensional structures of spatiotemporal chaos, which are a combination of unsteady hexagons, quadrangles and triangles. The transition from (b) to (c) is accompanied by the generation of interacting two-dimensional solitary pulses.
The influence of the texture of a hydrophobic surface on the electro-osmotic slip of the second kind and the electrokinetic instability near charge-selective surfaces (permselective membranes, electrodes, or systems of micro- and nanochannels) is investigated theoretically using a simple model based on the Rubinstein-Zaltzman approach. A simple formula is derived to evaluate the decrease in the instability threshold due to hydrophobicity. The study is complemented by numerical investigations both of linear and nonlinear instabilities near a hydrophobic membrane surface. Theory predicts a significant enhancement of the ion flux to the surface and shows a good qualitative agreement with the available experimental data.
Rod bundle flows are prevalent in nuclear engineering for both light water reactors (LWR) and advanced reactor concepts. Unlike canonical channel flow, the flow in rod bundles presents some unique characteristics, notably due to the inhomogeneous cross section which can present different local conditions of turbulence as well as localized effects characteristic of external flows. Despite the ubiquity of rod bundle flows and the decades of experimental and numerical knowledge acquired in this field, there are no publicly available direct numerical simulations (DNS) of the flow in multiple pin rod bundles with heat transfer. A multiple pin DNS study is of great value as it would allow for assessment of the reliability of various turbulence models in the presence of heat transfer, as well as allow for a deeper understanding of the flow physics. We present work towards DNS of the flow in a square 5x5 rod bundle representative of LWR fuel. We consider standard configurations as well as configurations where the central pin is replaced with a guide thimble. We perform simulations in STAR-CCM+ to design the numerical DNS, which is to be conducted using the open source spectral element code Nek5000. Large Eddy Simulations are also performed in Nek5000 to confirm that the resolution requirements are adequate. We compare results from STAR-CCM+ and Nek5000, which show very good agreement in the wide gaps with larger discrepancies in the narrow gaps. In particular, evidence of a gap vortex street is seen in the edge subchannels in LES but is not predicted by STAR-CCM+.
Results of direct numerical simulation of isotropic turbulence of surface gravity waves in the framework of Hamiltonian equations are presented. For the first time simultaneous formation of both direct and inverse cascades was observed in the framework of primordial dynamical equations. At the same time, strong long waves background was developed. It was shown, that obtained Kolmogorov spectra are very sensitive to the presence of this condensate. Such situation has to be typical for experimental wave tanks, flumes, and small lakes.
A comprehensive and detailed overview of the flow topology over a cambered NACA 65(1)-412 airfoil at Re = 20,000 is presented for angles of attack ranging from 0{deg} to 10{deg} using high-order direct numerical simulations. It is shown that instabilities bifurcate the flow and cause it to change at a critical angle of attack from laminar separation without reattachment over a laminar separation bubble at the trailing edge to a bubble at the leading edge. The transition of the flow regimes is governed by the Karman vortex shedding of the pressure side boundary layer at the trailing edge, Kelvin-Helmholtz instabilities within the separated shear layer on the suction side, as well as three-dimensional instabilities of elliptic flow within the vortex cores and hyperbolic flow in the shear layer between subsequent Karman vortices. As the suction side shear layer transitions and reattaches, the interaction of the two and three-dimensional instabilities results in three-dimensional tubular structures and large-scale turbulent puffs. The formation and shifting of the laminar separation bubble defines the far-wake topology several chord-lengths behind the airfoil and is accompanied by a sudden increase of the lift force and decrease in the drag that underscores the sensitive nature of low-Reynolds number airfoil aerodynamics. Lift and drag polars are presented for direct numerical simulations, wind tunnel experiments, and simplified numerical procedures where incorrect prediction of the force coefficients is caused by the failure to correctly model the low-pressure region at the trailing edge that is caused by the time-dependent generation of the Karman vortices.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
