Periodic forcing in viscous fingering of a nematic liquid crystal

We study viscous fingering of an air-nematic interface in a radial Hele-Shaw cell when periodically switching on and off an electric field, which reorients the nematic and thus changes its viscosity, as well as the surface tension and its anisotropy (mainly enforced by a single groove in the cell). We observe undulations at the sides of the fingers which correlate with the switching frequency and with tip oscillations which give maximal velocity to smallest curvatures. These lateral undulations appear to be decoupled from spontaneous (noise-induced) side branching. We conclude that the lateral undulations are generated by successive relaxations between two limiting finger widths. The change between these two selected pattern scales is mainly due to the change in the anisotropy. This scenario is confirmed by numerical simulations in the channel geometry, using a phase-field model for anisotropic viscous fingering.Comment: completely rewritten version, more clear exposition of results (14 pages in Revtex + 7 eps figures

A stabilized finite element predictor–corrector scheme for the incompressible Navier–Stokes equations using a nodal‐based implementation

A finite element model to solve the incompressible Navier–Stokes equations based on the stabilization with orthogonal subscales, a predictor–corrector scheme to segregate the pressure and a nodal based implementation is presented in this paper. The stabilization consists of adding a least‐squares form of the component orthogonal to the finite element space of the convective and pressure gradient terms, which allows to deal with convection‐dominated flows and to use equal velocity–pressure interpolation. The pressure segregation is inspired in fractional step schemes, although the converged solution corresponds to that of a monolithic time integration. Finally, the nodal‐based implementation is based on an a priori calculation of the integrals appearing in the formulation and then the construction of the matrix and right‐hand side vector of the final algebraic system to be solved. After appropriate approximations, this matrix and this vector can be constructed directly for each nodal point, without the need to loop over the elements and thus making the calculations much faster. Some issues related to this implementation for fractional step and our predictor–corrector scheme, which is the main contribution of this paper, are discusse

Towards a quantitative phase-field model of two-phase solidification

We construct a diffuse-interface model of two-phase solidification that quantitatively reproduces the classic free boundary problem on solid-liquid interfaces in the thin-interface limit. Convergence tests and comparisons with boundary integral simulations of eutectic growth show good accuracy for steady-state lamellae, but the results for limit cycles depend on the interface thickness through the trijunction behavior. This raises the fundamental issue of diffuse multiple-junction dynamics.Comment: 4 pages, 2 figures. Better final discussion. 1 reference adde

Quantitative Phase Field Model of Alloy Solidification

We present a detailed derivation and thin interface analysis of a phase-field model that can accurately simulate microstructural pattern formation for low-speed directional solidification of a dilute binary alloy. This advance with respect to previous phase-field models is achieved by the addition of a phenomenological "antitrapping" solute current in the mass conservation relation [A. Karma, Phys. Rev. Lett 87, 115701 (2001)]. This antitrapping current counterbalances the physical, albeit artificially large, solute trapping effect generated when a mesoscopic interface thickness is used to simulate the interface evolution on experimental length and time scales. Furthermore, it provides additional freedom in the model to suppress other spurious effects that scale with this thickness when the diffusivity is unequal in solid and liquid [R. F. Almgren, SIAM J. Appl. Math 59, 2086 (1999)], which include surface diffusion and a curvature correction to the Stefan condition. This freedom can also be exploited to make the kinetic undercooling of the interface arbitrarily small even for mesoscopic values of both the interface thickness and the phase-field relaxation time, as for the solidification of pure melts [A. Karma and W.-J. Rappel, Phys. Rev. E 53, R3017 (1996)]. The performance of the model is demonstrated by calculating accurately for the first time within a phase-field approach the Mullins-Sekerka stability spectrum of a planar interface and nonlinear cellular shapes for realistic alloy parameters and growth conditions.Comment: 51 pages RevTeX, 5 figures; expanded introduction and discussion; one table and one reference added; various small correction

Phase-field simulations of viscous fingering in shear-thinning fluids

A phase-field model for the Hele-Shaw flow of non-Newtonian fluids is developed. It extends a previous model for Newtonian fluids to a wide range of shear-dependent fluids. The model is applied to perform simulations of viscous fingering in shear- thinning fluids, and it is found to be capable of describing the complete crossover from the Newtonian regime at low shear rate to the strongly shear-thinning regime at high shear rate. The width selection of a single steady-state finger is studied in detail for a 2-plateaux shear-thinning law (Carreau law) in both its weakly and strongly shear-thinning limits, and the results are related to previous analyses. In the strongly shear-thinning regime a rescaling is found for power-law (Ostwald-de-Waehle) fluids that allows for a direct comparison between simulations and experiments without any adjustable parameters, and good agreement is obtained

β-decay half-lives and β-delayed neutron emission probabilities for several isotopes of Au, Hg, Tl, Pb, and Bi, beyond N = 126

Background: There have been measurements on roughly 230 nuclei that are β-delayed neutron emitters. They range from 8 He up to 150La. Apart from 210Tl, with a branching ratio of only 0.007%, no other neutron emitter has been measured beyond A = 150. Therefore, new data are needed, particularly in the region of heavy nuclei around N = 126, in order to guide theoretical models and help understand the formation of the third r-process peak at A ∼ 195. Purpose: To measure both β-decay half-lives and neutron branching ratios of several neutron-rich Au, Hg, Tl, Pb, and Bi isotopes beyond N = 126. Method: Ions of interest were produced by fragmentation of a 238U beam, selected and identified via the GSI-FRS fragment separator. A stack of segmented silicon detectors (SIMBA) was used to measure ion implants and β decays. An array of 30 3 He tubes embedded in a polyethylene matrix (BELEN) was used to detect neutrons with high efficiency and selectivity. A self-triggered digital system is employed to acquire data and to enable time correlations. The latter were analyzed with an analytical model and results for the half-lives and neutron-branching ratios were derived by using the binned maximum-likelihood method. Results: Twenty new β-decay half-lives are reported for 204−206Au, 208–211Hg, 211–216Tl, 215–218Pb, and 218–220Bi, nine of them for the first time. Neutron emission probabilities are reported for 210,211Hg and 211–216Tl. Conclusions: The new β-decay half-lives are in good agreement with previous measurements on nuclei in this region. The measured neutron emission probabilities are comparable to or smaller than values predicted by global models such as relativistic Hartree Bogoliubov plus the relativistic quasi-particle random phase approximation (RHB + RQRPA).Spanish Ministerio de Economía y Competitividad-FPA2011- 28770-C03-03, FPA2008-04972-C03-3, AIC-D2011-0705, FPA2011-24553, FPA2008-6419, FPA2010-17142, FPA2014-52823-C2-1-P, FPA2014- 52823-C2-2-P, and CPAN CSD-2007-00042 (Ingenio2010)Program Severo Ochoa-SEV-2014-0398German Helmholtz Association (Young Investigators)-VH-NG 627 (LISA-Lifetime Spectroscopy for Astrophysics)Nuclear Astrophysics Virtual Institute-VH-VI-417German Bundesministerium für Bildung und Forschung-06MT7178 / 05P12WOFNFSpanish Nuclear Security Council (CSN)-Catedra ArgosUK Science & Technology Facilities Council (STFC)-ST/F012012/

FPLUME-1.0: An integrated volcanic plume model accounting for ash aggregation

Eruption source parameters (ESP) characterizing volcanic eruption plumes are crucial inputs for atmospheric tephra dispersal models, used for hazard assessment and risk mitigation.We present FPLUME-1.0, a steady-state 1-D (one-dimensional) cross-section-averaged eruption column model based on the buoyant plume theory (BPT). The model accounts for plume bending by wind, entrainment of ambient moisture, effects of water phase changes, particle fallout and re-entrainment, a new parameterization for the air entrainment coefficients and a model for wet aggregation of ash particles in the presence of liquid water or ice. In the occurrence of wet aggregation, the model predicts an effective grain size distribution depleted in fines with respect to that erupted at the vent. Given a wind profile, the model can be used to determine the column height from the eruption mass flow rate or vice versa. The ultimate goal is to improve ash cloud dispersal forecasts by better constraining the ESP (column height, eruption rate and vertical distribution of mass) and the effective particle grain size distribution resulting from eventual wet aggregation within the plume. As test cases we apply the model to the eruptive phase-B of the 4 April 1982 El Chichón volcano eruption (México) and the 6 May 2010 Eyjafjallajökull eruption phase (Iceland). The modular structure of the code facilitates the implementation in the future code versions of more quantitative ash aggregation parameterization as further observations and experiment data will be available for better constraining ash aggregation processes

Radiation Induced Membrane Changes and Programmed Cell Death: Possible Interrelationships

A short review of the evidence that lymphocyte membranes are a target for the initiation of irradiation induced programmed cell death (PCD) is given. It is assumed that for lymphocytes PCD represents an essential physiological mechanism in order to prevent degeneration of the biological system involved. Initiation of PCD can be obtained by a pharmacological activation as well as with irradiation. In both cases, protein kinase-C (PKC) is involved in the signal transduction from the cellular membrane to the nucleus where, by means of a metabolically active process, DNA fragmentation is induced. It is hypothesized that processes connected to lipid peroxidation in the cell membrane constitute a primary effect of irradiation induced PCD, where membrane fluidization or a compensatory process aimed to the maintenance of membrane fluidity (membrane homeoviscosity hypothesis) are likely to be involved