Stability analysis of core-annular flow and neutral stability wave number

A general transition criterion is proposed in order to locate the core-annular flow pattern in horizontal and vertical oil-water flows. It is based on a rigorous one-dimensional two-fluid model of liquid-liquid two-phase flow and considers the existence of critical interfacial wave numbers related to a non-negligible interfacial tension term to which the linear stability theory still applies. The viscous laminar-laminar flow problem is fully resolved and turbulence effects on the stability are analyzed through experimentally obtained shape factors. The proposed general transition criterion includes in its formulation the inviscid Kelvin-Helmholtz`s discriminator. If a theoretical maximum wavelength is considered as a necessary condition for stability, a stability criterion in terms of the Eotvos number is achieved. Effects of interfacial tension, viscosity ratio, density difference, and shape factors on the stability of core-annular flow are analyzed in detail. The more complete modeling allowed for the analysis of the neutral-stability wave number and the results strongly suggest that the interfacial tension term plays an indispensable role in the correct prediction of the stable region of core-annular flow pattern. The incorporation of a theoretical minimum wavelength into the transition model produced significantly better results. The criterion predictions were compared with recent data from the literature and the agreement is encouraging. (C) 2007 American Institute of Chemical Engineers.

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

Optimal design of periodic functionally graded composites with prescribed properties

The computational design of a composite where the properties of its constituents change gradually within a unit cell can be successfully achieved by means of a material design method that combines topology optimization with homogenization. This is an iterative numerical method, which leads to changes in the composite material unit cell until desired properties (or performance) are obtained. Such method has been applied to several types of materials in the last few years. In this work, the objective is to extend the material design method to obtain functionally graded material architectures, i.e. materials that are graded at the local level (e.g. microstructural level). Consistent with this goal, a continuum distribution of the design variable inside the finite element domain is considered to represent a fully continuous material variation during the design process. Thus the topology optimization naturally leads to a smoothly graded material system. To illustrate the theoretical and numerical approaches, numerical examples are provided. The homogenization method is verified by considering one-dimensional material gradation profiles for which analytical solutions for the effective elastic properties are available. The verification of the homogenization method is extended to two dimensions considering a trigonometric material gradation, and a material variation with discontinuous derivatives. These are also used as benchmark examples to verify the optimization method for functionally graded material cell design. Finally the influence of material gradation on extreme materials is investigated, which includes materials with near-zero shear modulus, and materials with negative Poisson`s ratio.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Optical study of a light diaphragm rupture process in an expansion tube

Rupture of a light cellophane diaphragm in an expansion tube has been studied by an optical method. The influence of the light diaphragm on test flow generation has long been recognised, however the diaphragm rupture mechanism is less well known. It has been previously postulated that the diaphragm ruptures around its periphery due to the dynamic pressure loading of the shock wave, with the diaphragm material at some stage being removed from the flow to allow the shock to accelerate to the measured speeds downstream. The images obtained in this series of experiments are the first to show the mechanism of diaphragm rupture and mass removal in an expansion tube. A light diaphragm was impulsively loaded via a shock wave and a series of images was recorded holographically throughout the rupture process, showing gradual destruction of the diaphragm. Features such as the diaphragm material, the interface between gases, and a reflected shock were clearly visualised. Both qualitative and quantitative aspects of the rupture dynamics were derived from the images and compared with existing one-dimensional theory.

An open-boundary integrable model of three coupled XY spin chains

The integrable open-boundary conditions for the model of three coupled one-dimensional XY spin chains are considered in the framework of the quantum inverse scattering method. The diagonal boundary K-matrices are found and a class of integrable boundary terms is determined. The boundary model Hamiltonian is solved by using the coordinate space Bethe ansatz technique and Bethe ansatz equations are derived. (C) 1998 Elsevier Science B.V.

Characterization of the sequential non-covalent and covalent interactions of the antitumour antibiotic hedamycin with double stranded DNA by NMR spectroscopy

Hedamycin, a member of the pluramycin class of antitumour antibiotics, consists of a planar anthrapyrantrione chromophore to which is attached two aminosugar rings at one end and a bisepoxide-containing sidechain at the other end, Binding to double-stranded DNA is known to involve both reversible and non-reversible modes of interaction. As a part of studies directed towards elucidating the structural basis for the observed 5'-pyGT-3' sequence selectivity of hedamycin, we conducted one-dimensional NMR titration experiments at low temperature using the hexadeoxyribonucleotide duplexes d(CACGTG)(2) and d(CGTACG)(2). Spectral changes which occurred during these titrations are consistent with hedamycin initially forming a reversible complex in slow exchange on the NMR timescale and binding through intercalation of the chromophore. Monitoring of this reversible complex over a period of hours revealed a second type of spectral change which corresponds with formation of a non-reversible complex. Copyright (C) 1999 John Wiley & Sons, Ltd.

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.

Phonon anomalies due to strong-electronic correlations in layered organic metals

We show how the coupling between the phonons and electrons in a strongly correlated metal can result in phonon frequencies that have a nonmonotonic temperature dependence. Dynamical mean-field theory is used to study the Hubbard-Holstein model that describes the kappa-(BEDT-TTF)(2)X [where BEDT-TTF is bis-(ethylenedithia-tetrathiafulvalene)] family of superconducting molecular crystals. The crossover with increasing temperature from a Fermi liquid to a bad metal produces phonon anomalies that are relevant to recent Raman scattering and acoustic experiments.

Production efficiency and financing of public health: an analysis of small municipalities in the state of Sao Paulo - Brazil

This paper measured the variations in performance of small municipalities in the State of Sao Paulo, Brazil, regarding the technical efficiency in the use of public funds in public primary health care actions concerning the funding profile, in a scenario of fiscal federalism. Technical efficiency is one of the parameters of evaluation of public sector performance and was measured by means of Data Envelopment Analysis (DEA). The correlation analysis of DEA score was used to verify possible associations between technical efficiency and the funding profile of expenses with health care. The results showed that 6.41% of the municipalities were considered efficient. They also showed that the level of municipality dependence to inter-governmental general purpose grants and the national health funding specific purpose grants have negative correlation with efficiency scores.

Teleportation in a non-inertial frame

In this work, we describe the process of teleportation between Alice in an inertial frame, and Rob who is in uniform acceleration with respect to Alice. The fidelity of the teleportation is reduced due to Davies-Unruh radiation in Rob's frame. In so far as teleportation is a measure of entanglement, our results suggest that quantum entanglement is degraded in non-inertial frames. We discuss this reduction in fidelity for both bosonic and fermionic resources.

Algebraic Bethe ansatz for the anisotropic supersymmetric U model

We present an algebraic Bethe ansatz for the anisotropic supersymmetric U model for correlated electrons on the unrestricted 4(L)-dimensional electronic Hilbert space x(n=l)(L)C(4)(where L is the lattice length). The supersymmetry algebra of the local Hamiltonian is the quantum superalgebra U-q[gl(2\1)] and the model contains two symmetry-preserving free real parameters; the quantization parameter q and the Hubbard interaction parameter U. The parameter U arises from the one-parameter family of inequivalent typical four-dimensional irreps of U-q[gl(2\1)]. Eigenstates of the model are determined by the algebraic Bethe ansatz on a one-dimensional periodic lattice.

Discrete and continuous models for dry masonry columns

The dynamic response of dry masonry columns can be approximated with finite-difference equations. Continuum models follow by replacing the difference quotients of the discrete model by corresponding differential expressions. The mathematically simplest of these models is a one-dimensional Cosserat theory. Within the presented homogenization context, the Cosserat theory is obtained by making ad hoc assumptions regarding the relative importance of certain terms in the differential expansions. The quality of approximation of the various theories is tested by comparison of the dispersion relations for bending waves with the dispersion relation of the discrete theory. All theories coincide with differences of less than 1% for wave-length-block-height (L/h) ratios bigger than 2 pi. The theory based on systematic differential approximation remains accurate up to L/h = 3 and then diverges rapidly. The Cosserat model becomes increasingly inaccurate for L/h < 2 pi. However, in contrast to the systematic approximation, the wave speed remains finite. In conclusion, considering its relative simplicity, the Cosserat model appears to be the natural starting point for the development of continuum models for blocky structures.

Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.