Lattice Boltzmann modeling of dendritic growth in forced and natural convection

AbstractA two-dimensional (2D) coupled model is developed for the simulation of dendritic growth during alloy solidification in the presence of forced and natural convection. Instead of conventional continuum-based Navier–Stokes (NS) solvers, the present model adopts a kinetic-based lattice Boltzmann method (LBM), which describes flow dynamics by the evolution of distribution functions of moving pseudo-particles, for the numerical computations of flow dynamics as well as thermal and solutal transport. The dendritic growth is modeled using a solutal equilibrium approach previously proposed by Zhu and Stefanescu (ZS), in which the evolution of the solid/liquid interface is driven by the difference between the local equilibrium composition and the local actual liquid composition. The local equilibrium composition is calculated from the local temperature and curvature. The local temperature and actual liquid composition, controlled by both diffusion and convection, are obtained by solving the LB equations using the lattice Bhatnagar–Gross–Krook (LBGK) scheme. Detailed model validation is performed by comparing the simulations with analytical predictions, which demonstrates the quantitative capability of the proposed model. Furthermore, the convective dendritic growth features predicted by the present model are compared with those obtained from the Zhu–Stefanescu and Navier–Stokes (ZS–NS) model, in which the fluid flow is calculated using an NS solver. It is found that the evolution of the solid fraction of dendritic growth calculated by both models coincides well. However, the present model has the significant advantages of numerical stability and computational efficiency for the simulation of dendritic growth with melt convection

Effect of Ordering on Spinodal Decomposition of Liquid-Crystal/Polymer Mixtures

Partially phase-separated liquid-crystal/polymer dispersions display highly fibrillar domain morphologies that are dramatically different from the typical structures found in isotropic mixtures. To explain this, we numerically explore the coupling between phase ordering and phase separation kinetics in model two-dimensional fluid mixtures phase separating into a nematic phase, rich in liquid crystal, coexisting with an isotropic phase, rich in polymer. We find that phase ordering can lead to fibrillar networks of the minority polymer-rich phase

Reversible stretching of homopolymers and random heteropolymers

We have analyzed the equilibrium response of chain molecules to stretching. For a homogeneous sequence of monomers, the induced transition from compact globule to extended coil below the $\theta$-temperature is predicted to be sharp. For random sequences, however, the transition may be smoothed by a prevalence of necklace-like structures, in which globular regions and coil regions coexist in a single chain. As we show in the context of a random copolymer, preferential solvation of one monomer type lends stability to such structures. The range of stretching forces over which necklaces are stable is sensitive to chain length as well as sequence statistics.Comment: 14 pages, 4 figure

Rigorous results on superconducting ground states for attractive extended Hubbard models

We show that the exact ground state for a class of extended Hubbard models including bond-charge, exchange, and pair-hopping terms, is the Yang "eta-paired" state for any non-vanishing value of the pair-hopping amplitude, at least when the on-site Coulomb interaction is attractive enough and the remaining physical parameters satisfy a single constraint. The ground state is thus rigorously superconducting. Our result holds on a bipartite lattice in any dimension, at any band filling, and for arbitrary electron hopping.Comment: 12 page

Two-soliton collisions in a near-integrable lattice system

We examine collisions between identical solitons in a weakly perturbed Ablowitz-Ladik (AL) model, augmented by either onsite cubic nonlinearity (which corresponds to the Salerno model, and may be realized as an array of strongly overlapping nonlinear optical waveguides), or a quintic perturbation, or both. Complex dependences of the outcomes of the collisions on the initial phase difference between the solitons and location of the collision point are observed. Large changes of amplitudes and velocities of the colliding solitons are generated by weak perturbations, showing that the elasticity of soliton collisions in the AL model is fragile (for instance, the Salerno's perturbation with the relative strength of 0.08 can give rise to a change of the solitons' amplitudes by a factor exceeding 2). Exact and approximate conservation laws in the perturbed system are examined, with a conclusion that the small perturbations very weakly affect the norm and energy conservation, but completely destroy the conservation of the lattice momentum, which is explained by the absence of the translational symmetry in generic nonintegrable lattice models. Data collected for a very large number of collisions correlate with this conclusion. Asymmetry of the collisions (which is explained by the dependence on the location of the central point of the collision relative to the lattice, and on the phase difference between the solitons) is investigated too, showing that the nonintegrability-induced effects grow almost linearly with the perturbation strength. Different perturbations (cubic and quintic ones) produce virtually identical collision-induced effects, which makes it possible to compensate them, thus finding a special perturbed system with almost elastic soliton collisions.Comment: Phys. Rev. E, in pres

Surface layering of liquids: The role of surface tension

Recent measurements show that the free surfaces of liquid metals and alloys are always layered, regardless of composition and surface tension; a result supported by three decades of simulations and theory. Recent theoretical work claims, however, that at low enough temperatures the free surfaces of all liquids should become layered, unless preempted by bulk freezing. Using x-ray reflectivity and diffuse scattering measurements we show that there is no observable surface-induced layering in water at T=298 K, thus highlighting a fundamental difference between dielectric and metallic liquids. The implications of this result for the question in the title are discussed.Comment: 5 pages, 4 figures, to appear in Phys. Rev. B. 69 (2004

Switching model with two habitats and a predator involving group defence

Switching model with one predator and two prey species is considered. The prey species have the ability of group defence. Therefore, the predator will be attracted towards that habitat where prey are less in number. The stability analysis is carried out for two equilibrium values. The theoretical results are compared with the numerical results for a set of values. The Hopf bifuracation analysis is done to support the stability results

Compressibility and Electronic Structure of MgB2 up to 8 GPa

The lattice parameters of MgB2 up to pressures of 8 GPa were determined using high-resolution x-ray powder diffraction in a diamond anvil cell. The bulk modulus, B0, was determined to be 151 +-5 GPa. Both experimental and first-principles calculations indicate nearly isotropic mechanical behavior under pressure. This small anisotropy is in contrast to the 2 dimensional nature of the boron pi states. The pressure dependence of the density of states at the Fermi level and a reasonable value for the average phonon frequency account within the context of BCS theory for the reduction of Tc under pressure.Comment: REVTeX file. 4 pages, 4 figure

Confront Holographic QCD with Regge Trajectories of vectors and axial-vectors

We derive the general 5-dimension metric structure of the $Dp-Dq$ system in type II superstring theory, and demonstrate the physical meaning of the parameters characterizing the 5-dimension metric structure of the \textit{holographic} QCD model by relating them to the parameters describing Regge trajectories. By matching the spectra of vector mesons $\rho_1$ with deformed $Dp-Dq$ soft-wall model, we find that the spectra of vector mesons $\rho_1$ can be described very well in the soft-wall $D3-Dq$ model, i.e, $AdS_5$ soft-wall model. We then investigate how well the $AdS_5$ soft-wall model can describe the Regge trajectory of axial-vector mesons $a_1$. We find that the constant component of the 5-dimension mass square of axial-vector mesons plays an efficient role to realize the chiral symmetry breaking in the vacuum, and a small negative $z^4$ correction in the 5-dimension mass square is helpful to realize the chiral symmetry restoration in high excitation states.Comment: 9 pages, 3 figure and 3 tables, one section adde

Long-term trends in tropical cyclone tracks around Korea and Japan in late summer and early fall

This study investigates long-term trends in tropical cyclones (TCs) over the extratropical western North Pacific (WNP) over a period of 35 years (1982-2016). The area analyzed extended across 30-45 degrees N and 120-150 degrees E, including the regions of Korea and Japan that were seriously affected by TCs. The northward migration of TCs over the WNP to the mid-latitudes showed a sharp increase in early fall. In addition, the duration of TCs over the WNP that migrated northwards showed an increase, specifically in early to mid-September. Therefore, more recently, TC tracks have been observed to significantly extend into the mid-latitudes. The recent northward extension of TC tracks over the WNP in early fall was observed to be associated with changes in environmental conditions that were favorable for TC activities, including an increase in sea surface temperature (SST), decrease in vertical wind shear, expansion of subtropical highs, strong easterly steering winds, and an increase in relative vorticity. In contrast, northward migrations of TCs to Korea and Japan showed a decline in late August, because of the presence of unfavorable environmental conditions for TC activities. These changes in environmental conditions, such as SST and vertical wind shear, can be partially associated with the Pacific decadal oscillation