Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Discrete approximation to the global spectrum of the tangent operator for flow past a circular cylinder

This paper deals with the calculation of the discrete approximation to the full spectrum for the tangent operator for the stability problem of the symmetric flow past a circular cylinder. It is also concerned with the localization of the Hopf bifurcation in laminar flow past a cylinder, when the stationary solution loses stability and often becomes periodic in time. The main problem is to determine the critical Reynolds number for which a pair of eigenvalues crosses the imaginary axis. We thus present a divergence-free method, based on a decoupling of the vector of velocities in the saddle-point system from the vector of pressures, allowing the computation of eigenvalues, from which we can deduce the fundamental frequency of the time-periodic solution. The calculation showed that stability is lost through a symmetry-breaking Hopf bifurcation and that the critical Reynolds number is in agreement with the value presented in reported computations. (c) 2007 IMACS. Published by Elsevier B.V. All rights reserved.

Extensions of discrete triangular distributions and boundary bias in kernel estimation for discrete functions

Asymmetric discrete triangular distributions are introduced in order to extend the symmetric ones serving for discrete associated kernels in the nonparametric estimation for discrete functions. The extension from one to two orders around the mode provides a large family of discrete distributions having a finite support. Establishing a bridge between Dirac and discrete uniform distributions, some different shapes are also obtained and their properties are investigated. In particular, the mean and variance are pointed out. Applications to discrete kernel estimators are given with a solution to a boundary bias problem. (C) 2010 Elsevier B.V. All rights reserved.

Entanglement in a simple quantum phase transition

What entanglement is present in naturally occurring physical systems at thermal equilibrium? Most such systems are intractable and it is desirable to study simple but realistic systems that can be solved. An example of such a system is the one-dimensional infinite-lattice anisotropic XY model. This model is exactly solvable using the Jordan-Wigner transform, and it is possible to calculate the two-site reduced density matrix for all pairs of sites. Using the two-site density matrix, the entanglement of formation between any two sites is calculated for all parameter values and temperatures. We also study the entanglement in the transverse Ising model, a special case of the XY model, which exhibits a quantum phase transition. It is found that the next-nearest-neighbor entanglement (though not the nearest-neighbor entanglement) is a maximum at the critical point. Furthermore, we show that the critical point in the transverse Ising model corresponds to a transition in the behavior of the entanglement between a single site and the remainder of the lattice.

Two-level atom in a squeezed vacuum with finite bandwidth

The resonance fluorescence of a two-level atom driven by a coherent laser field and damped by a finite bandwidth squeezed vacuum is analysed. We extend the Yeoman and Barnett technique to a non-zero detuning of the driving field from the atomic resonance and discuss the role of squeezing bandwidth and the detuning in the level shifts, widths and intensities of the spectral lines. The approach is valid for arbitrary values of the Rabi frequency and detuning but for the squeezing bandwidths larger than the natural linewidth in order to satisfy the Markoff approximation. The narrowing of the spectral lines is interpreted in terms of the quadrature-noise spectrum. We find that, depending on the Rabi frequency, detuning and the squeezing phase, different factors contribute to the line narrowing. For a strong resonant driving field there is no squeezing in the emitted field and the fluorescence spectrum exactly reveals the noise spectrum. In this case the narrowing of the spectral lines arises from the noise reduction in the input squeezed vacuum. For a weak or detuned driving field the fluorescence exhibits a large squeezing and, as a consequence, the spectral lines have narrowed linewidths. Moreover, the fluorescence spectrum can be asymmetric about the central frequency despite the symmetrical distribution of the noise. The asymmetry arises from the absorption of photons by the squeezed vacuum which reduces the spontaneous emission. For an appropriate choice of the detuning some of the spectral lines can vanish despite that there is no population trapping. Again this process can be interpreted as arising from the absorption of photons by the squeezed vacuum. When the absorption is large it may compensate the spontaneous emission resulting in the vanishing of the fluorescence lines.

Three-dimensional solution structure of alpha-conotoxin MII by NMR spectroscopy: Effects of solution environment on helicity

alpha-Conotoxin MII, a 16-residue polypeptide from the venom of the piscivorous cone snail Conus magus, is a potent and highly specific blocker of mammalian neuronal nicotinic acetylcholine receptors composed of alpha 3 beta 2 subunits. The role of this receptor type in the modulation of neurotransmitter release and its relevance to the problems of addiction and psychosis emphasize the importance of a structural understanding of the mode of interaction of MII with the alpha 3 beta 2 interface. Here we describe the three-dimensional solution structure of MIT determined using 2D H-1 NMR spectroscopy. Structural restraints consisting of 376 interproton distances inferred from NOEs and 12 dihedral restraints derived from spin-spin coupling constants were used as input for simulated annealing calculations and energy minimization in the program X-PLOR. The final set of 20 structures is exceptionally well-defined with mean pairwise rms differences over the whole molecule of 0.07 Angstrom for the backbone atoms and 0.34 Angstrom for all heavy atoms. MII adopts a compact structure incorporating a central segment of alpha-helix and beta-turns at the N- and C-termini. The molecule is stabilized by two disulfide bonds, which provide cross-links between the N-terminus and both the middle and C-terminus of the structure. The susceptibility of the structure to conformational change was examined using several different solvent conditions. While the global fold of MII remains the same, the structure is stabilized in a more hydrophobic environment provided by the addition of acetonitrile or trifluoroethanol to the aqueous solution. The distribution of amino acid side chains in MII creates distinct hydrophobic and polar patches on its surface that may be important for the specific interaction with the alpha 3 beta 2 neuronal nAChR. A comparison of the structure of MII with other neuronal-specific alpha-conotoxins provides insights into their mode of interaction with these receptors.

Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Quantum oscillations in quasi-one-dimensional metals with spin-density-wave ground states

We consider the magnetoresistance oscillation phenomena in the Bechgaard salts (TMTSF)(2)X, where X = ClO4, PF6, and AsF6 in pulsed magnetic fields to 51 T. Of particular importance is the observation of a new magnetoresistance oscillation for X = ClO4 in its quenched state. In the absence of any Fermi-surface reconstruction due to anion order at low temperatures, all three materials exhibit nonmonotonic temperature dependence of the oscillation amplitude in the spin-density-wave (SDW) state. We discuss a model where, below a characteristic temperature T* within the SDW state, a magnetic breakdown gap opens. [S0163-1829(99)00904-2].

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Incorporating phase retardation effects into radiofrequency resonator models: application to magnetic resonance microscopy

A method is presented for including path propagation effects into models of radiofrequency resonators for use in magnetic resonance imaging. The method is based on the use of Helmholtz retarded potentials and extends our previous work on current density models of resonators based on novel inverse finite Hilbert transform solutions to the requisite integral equations. Radiofrequency phase retardation effects are most pronounced at high field strengths (frequencies) as are static field perturbations due to the magnetic materials in the resonators themselves. Both of these effects are investigated and a novel resonator structure presented for use in magnetic resonance microscopy.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.

Synthesis and structure elucidation of kappa-casein (45-87) - A two dimensional nuclear magnetic resonance study

We have shown that 44 amino acid residues N-terminal segment of kappa-casein exhibits considerable a-helical structure. This prompted us to investigate the structures of the remaining segments of kappa-casein. Thus, in this study the chemical synthesis and structure elucidation of the peptide 45-87 amino acid residues of kappa-casein is reported. The peptide was assembled using solid phase peptide synthesis methodology on pam resin, cleaved via HF, freeze dried and, after purification, characterised by mass spectrometry (observed m/z 4929; calculated mit 4929.83). The amino acid sequence of the peptide is: CKPVALINNQFLPYPYYAKPAAVRSPAQILQWQVLSNTVPAKA Its structure elucidation has been carried out using circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques. CD spectrum of the peptide shows it to be a random structure in water but in 30% trifluoroethanol the peptide exhibits considerable structure. The 1D and 2D NMR spectra corroborated the results of CD. The structure elucidation of the peptide using TOCSY and NOESY NMR techniques will be discussed.

Scope for further analytical solutions for constant flux infiltration into a semi-infinite soil profile or redistribution in a finite soil profile

[1] We attempt to generate new solutions for the moisture content form of the one-dimensional Richards' [1931] equation using the Lisle [1992] equivalence mapping. This mapping is used as no more general set of transformations exists for mapping the one-dimensional Richards' equation into itself. Starting from a given solution, the mapping has the potential to generate an infinite number of new solutions for a series of nonlinear diffusivity and hydraulic conductivity functions. We first seek new analytical solutions satisfying Richards' equation subject to a constant flux surface boundary condition for a semi-infinite dry soil, starting with the Burgers model. The first iteration produces an existing solution, while subsequent iterations are shown to endlessly reproduce this same solution. Next, we briefly consider the problem of redistribution in a finite-length soil. In this case, Lisle's equivalence mapping is generalized to account for arbitrary initial conditions. As was the case for infiltration, however, it is found that new analytical solutions are not generated using the equivalence mapping, although existing solutions are recovered.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.