A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Temperature (over) compensation in an oscillatory surface reaction

Biological rhythms are regulated by homeostatic mechanisms that assure that physiological clocks function reliably independent of temperature changes in the environment. Temperature compensation, the independence of the oscillatory period on temperature, is known to play a central role in many biological rhythms, but it is rather rare in chemical oscillators. We study the influence of temperature on the oscillatory dynamics during the catalytic oxidation of formic acid on a polycrystalline platinum electrode. The experiments are performed at five temperatures from 5 to 25 degrees C, and the oscillations are studied under galvanostatic control. Under oscillatory conditions, only non-Arrhenius behavior is observed. Overcompensation with temperature coefficient (q(10), defined as the ratio between the rate constants at temperature T + 10 degrees C and at T) < I is found in most cases, except that temperature compensation with q(10) approximate to I predominates at high applied currents. The behavior of the period and the amplitude result from a complex interplay between temperature and applied current or, equivalently, the distance from thermodynamic equilibrium. High, positive apparent activation energies were obtained under voltammetric, nonoscillatory conditions, which implies that the non-Arrhenius behavior observed under oscillatory conditions results from the interplay among reaction steps rather than, from a weak temperature dependence of the individual steps.