Consistent quaternion interpolation for objective finite element approximation of geometrically exact beam

We explore an isoparametric interpolation of total quaternion for geometrically consistent, strain-objective and path-independent finite element solutions of the geometrically exact beam. This interpolation is a variant of the broader class known as slerp. The equivalence between the proposed interpolation and that of relative rotation is shown without any recourse to local bijection between quaternions and rotations. We show that, for a two-noded beam element, the use of relative rotation is not mandatory for attaining consistency cum objectivity and an appropriate interpolation of total rotation variables is sufficient. The interpolation of total quaternion, which is computationally more efficient than the one based on local rotations, converts nodal rotation vectors to quaternions and interpolates them in a manner consistent with the character of the rotation manifold. This interpolation, unlike the additive interpolation of total rotation, corresponds to a geodesic on the rotation manifold. For beam elements with more than two nodes, however, a consistent extension of the proposed quaternion interpolation is difficult. Alternatively, a quaternion-based procedure involving interpolation of relative rotations is proposed for such higher order elements. We also briefly discuss a strategy for the removal of possible singularity in the interpolation of quaternions, proposed in [I. Romero, The interpolation of rotations and its application to finite element models of geometrically exact rods, Comput. Mech. 34 (2004) 121–133]. The strain-objectivity and path-independence of solutions are justified theoretically and then demonstrated through numerical experiments. This study, being focused only on the interpolation of rotations, uses a standard finite element discretization, as adopted by Simo and Vu-Quoc [J.C. Simo, L. Vu-Quoc, A three-dimensional finite rod model part II: computational aspects, Comput. Methods Appl. Mech. Engrg. 58 (1986) 79–116]. The rotation update is achieved via quaternion multiplication followed by the extraction of the rotation vector. Nodal rotations are stored in terms of rotation vectors and no secondary storages are required.

Bond-order wave phase, spin solitons, and thermodynamics of a frustrated linear spin-1/2 Heisenberg antiferromagnet

The linear spin-1/2 Heisenberg antiferromagnet with exchanges J(1) and J(2) between first and second neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at J(2)/J(1)=0.2411 and is particularly simple at J(2)/J(1)=1/2. The BOW phase has a doubly degenerate singlet ground state, broken inversion symmetry, and a finite-energy gap E-m to the lowest-triplet state. The interval 0.4 < J(2)/J(1) < 1.0 has large E-m and small finite-size corrections. Exact solutions are presented up to N = 28 spins with either periodic or open boundary conditions and for thermodynamics up to N = 18. The elementary excitations of the BOW phase with large E-m are topological spin-1/2 solitons that separate BOWs with opposite phase in a regular array of spins. The molar spin susceptibility chi(M)(T) is exponentially small for T << E-m and increases nearly linearly with T to a broad maximum. J(1) and J(2) spin chains approximate the magnetic properties of the BOW phase of Hubbard-type models and provide a starting point for modeling alkali-tetracyanoquinodimethane salts.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.

Two- and Three-Dimensional Open-Framework Uranium Arsenates: Synthesis, Structure, and Characterization

Hydrothermal reactions between uranium salts and arsenic pentoxide in the presence of two different amines yielded six new uranium arsenate phases exhibiting open-framework structures, ethylenediamine (en): [C2N2H9]-[(UO2)(ASO(4))] I; [C2N2H10][(UO2)F(HASO(4))]2 center dot 4H(2)O, II; [C2N2H9][U2F5(HASO(4))(2)], III; [C2N2H9][UF2(ASO(4))], IV; diethylenetriamine (DETA), [C4N3H16][U2F3(ASO(4))(2)(HAsO4)] V; and [C4N3H16][U2F6(AsO4)(HAsO4)], VI. The structures were determined using single crystal studies, which revealed two- (I, II, V) and three-dimensional (III, IV, VI) structures for the uranium arsenates. The uranium atom, in these compounds, exhibits considerable variations in the coordination (6 to 9) that appears to have some correlation with the synthetic conditions. The water molecules in [C2N2H10][(UO2)F(HAsO4)](2 center dot)4H(2)O, II, could be reversibly removed, and the dehydrated phase, [C2N2H10][(UO2)F(HAsO4)](2), IIa, was also characterized using single crystal studies. The observation of many mineralogical structures in the present compounds suggests that the hydrothermal method could successfully replicate the geothermal conditions. As part of this study, we have observed autunite, Ca[(UO2)(PO4)](2)(H2O)(11), metavauxite, [Fe(H2O)(6)][Al(OH)(H2O)(PO4)](2), finarite, PbCU(SO4)(OH)(2), and tancoite, LiNa2H[Al(PO4)(2)(OH)], structures. The repeated observation of the secondary building unit, SBU-4, in many of the uranium arsenate structures suggests that these are viable building units. Optical studies on the uranium arsenate compound, [C4N3H16][U2F6(AsO4)(HASO(4))), VI, containing uranium in the +4 oxidation state indicates a blue emission through an upconversion process. The compound also exhibits antiferromagnetic behavior.

Graphene-like two-dimensional ionic boron with double dirac cones at ambient condition

Recently, partially ionic boron (γ-B28) has been predicted and observed in pure boron, in bulk phase and controlled by pressure [Nature, 457 (2009) 863]. By using ab initio evolutionary structure search, we report the prediction of ionic boron at a reduced dimension and ambient pressure, namely, the two-dimensional (2D) ionic boron. This 2D boron structure consists of graphene-like plane and B2 atom pairs, with the P6/mmm space group and 6 atoms in the unit cell, and has lower energy than the previously reported α-sheet structure and its analogues. Its dynamical and thermal stability are confirmed by the phonon-spectrum and ab initio molecular dynamics simulation. In addition, this phase exhibits double Dirac cones with massless Dirac fermions due to the significant charge transfer between the graphene-like plane and B2 pair that enhances the energetic stability of the P6/mmm boron. A Fermi velocity (vf) as high as 2.3 x 106 m/s, which is even higher than that of graphene (0.82 x 106 m/s), is predicted for the P6/mmm boron. The present work is the first report of the 2D ionic boron at atmospheric pressure. The unique electronic structure renders the 2D ionic boron a promising 2D material for applications in nanoelectronics.

Uniqueness of Walkup's 9-vertex 3-dimensional Klein bottle

Via a computer search, Altshuler and Steinberg found that there are 1296+1 combinatorial 3-manifolds on nine vertices, of which only one is non-sphere. This exceptional 3-manifold View the MathML source triangulates the twisted S2-bundle over S1. It was first constructed by Walkup. In this paper, we present a computer-free proof of the uniqueness of this non-sphere combinatorial 3-manifold. As opposed to the computer-generated proof, ours does not require wading through all the 9-vertex 3-spheres. As a preliminary result, we also show that any 9-vertex combinatorial 3-manifold is equivalent by proper bistellar moves to a 9-vertex neighbourly 3-manifold.

The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Dimensional Reduction Near the Deconfinement Transition

When ordinary nuclear matter is heated to a high temperature of ~ 10^12 K, it undergoes a deconfinement transition to a new phase, strongly interacting quark-gluon plasma. While the color charged fundamental constituents of the nuclei, the quarks and gluons, are at low temperatures permanently confined inside color neutral hadrons, in the plasma the color degrees of freedom become dominant over nuclear, rather than merely nucleonic, volumes. Quantum Chromodynamics (QCD) is the accepted theory of the strong interactions, and confines quarks and gluons inside hadrons. The theory was formulated in early seventies, but deriving first principles predictions from it still remains a challenge, and novel methods of studying it are needed. One such method is dimensional reduction, in which the high temperature dynamics of static observables of the full four-dimensional theory are described using a simpler three-dimensional effective theory, having only the static modes of the various fields as its degrees of freedom. A perturbatively constructed effective theory is known to provide a good description of the plasma at high temperatures, where asymptotic freedom makes the gauge coupling small. In addition to this, numerical lattice simulations have, however, shown that the perturbatively constructed theory gives a surprisingly good description of the plasma all the way down to temperatures a few times the transition temperature. Near the critical temperature, the effective theory, however, ceases to give a valid description of the physics, since it fails to respect the approximate center symmetry of the full theory. The symmetry plays a key role in the dynamics near the phase transition, and thus one expects that the regime of validity of the dimensionally reduced theories can be significantly extended towards the deconfinement transition by incorporating the center symmetry in them. In the introductory part of the thesis, the status of dimensionally reduced effective theories of high temperature QCD is reviewed, placing emphasis on the phase structure of the theories. In the first research paper included in the thesis, the non-perturbative input required in computing the g^6 term in the weak coupling expansion of the pressure of QCD is computed in the effective theory framework at an arbitrary number of colors. The two last papers on the other hand focus on the construction of the center-symmetric effective theories, and subsequently the first non-perturbative studies of these theories are presented. Non-perturbative lattice simulations of a center-symmetric effective theory for SU(2) Yang-Mills theory show --- in sharp contrast to the perturbative setup --- that the effective theory accommodates a phase transition in the correct universality class of the full theory. This transition is seen to take place at a value of the effective theory coupling constant that is consistent with the full theory coupling at the critical temperature.

Dancing into the second phase of creative industries

As we enter the second phase of creative industries there is a shift away from the early 1990s ideology of the arts as a creative content provider for the wealth generating ‘knowledge’ economy to an expanded rhetoric encompassing ‘cultural capital’ and its symbolic value. A renewed focus on culture is examined through a regional scan of creative industries in which social engineering of the arts occurs through policy imperatives driven by ‘profit oriented conceptualisations of culture’ (Hornidge 2011, p. 263) In the push for artists to become ‘culturpreneurs’ a trend has emerged where demand for ‘embedded creatives’ (Cunningham 2013) sees an exodus from arts-based employment through use of transferable skills into areas outside the arts. For those that stay, within the performing arts in particular, employment remains project-based, sporadic, underpaid, self-initiated and often self-financed, requiring adaptive career paths. Artist entrepreneurs must balance creation and performance of their art with increasing amounts of time spent on branding, compliance, fundraising and the logistical and commercial requirements of operating in a CI paradigm. The artists’ key challenge thus becomes one of aligning core creative and aesthetic values with market and business considerations. There is also the perceived threat posed by the ‘prosumer’ phenomenon (Bruns 2008), in which digital on-line products are created and produced by those formerly seen as consumers of art or audiences for art. Despite negative aspects to this scenario, a recent study (Steiner & Schneider 2013) reveals that artists are happier and more satisfied than other workers within and outside the creative industries. A lively hybridisation of creative practice is occurring through mobile and interactive technologies with dynamic connections to social media. Continued growth in arts festivals attracts participation in international and transdisciplinary collaborations, whilst cross-sectoral partnerships provide artists with opportunities beyond a socio-cultural setting into business, health, science and education. This is occurring alongside a renewed engagement with place through the rise of cultural precincts in ‘creative cities’ (Florida 2008, Landry 2000), providing revitalised spaces for artists to gather and work. Finally, a reconsideration of the specialist attributes and transferable skills that artists bring to the creative industries suggests ways to dance through both the challenges and opportunities occasioned by the current complexities of arts’ practices.

Quantum phase transition in capacitively coupled double quantum dots

We investigate two equivalent, capacitively coupled semiconducting quantum dots, each coupled to its own lead, in a regime where there are two electrons on the double dot. With increasing interdot coupling, a rich range of behavior is uncovered: first a crossover from spin- to charge-Kondo physics, via an intermediate SU(4) state with entangled spin and charge degrees of freedom, followed by a quantum phase transition of Kosterlitz-Thouless type to a non-Fermi-liquid "charge-ordered" phase with finite residual entropy and anomalous transport properties. Physical arguments and numerical renormalization group methods are employed to obtain a detailed understanding of the problem.

Thermal stresses in a finite hollow cylinder due to an axisymmetric temperature field at the end surface

A three-dimensional exact solution for determining the thermal stresses in a finite hollow cylinder subject to a steady state axisymmetric temperature field over one of its end surfaces has been given. Numerical results for a hollow cylinder, having lenght to outer diameter ratio equal to one and inner to outer diameter ratio equal to 0.75, subjected to a symmetric temperature variation over the end surfaces of the cylinder have been given.

Thermal stresses in a finite solid cylinder due to an axisymmetric temperature field at the end surface.

A three-dimensional rigorous solution for determining thermal stresses in a finite solid cylinder due to a steady state axisymmetric temperature field over one of its end surfaces is given. Numerical results for a solid cylinder having a length to diameter ratio equal to one and subjected to a symmetric temperature variation over half the radius of the cylinder at the end surfaces are included. These results have been compared with the results of the approximate solution given by W. Nowacki.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

The Line Spectral Frequency Model of a Finite-Length Sequence

The line spectral frequency (LSF) of a causal finite length sequence is a frequency at which the spectrum of the sequence annihilates or the magnitude spectrum has a spectral null. A causal finite-length sequencewith (L + 1) samples having exactly L-LSFs, is referred as an Annihilating (AH) sequence. Using some spectral properties of finite-length sequences, and some model parameters, we develop spectral decomposition structures, which are used to translate any finite-length sequence to an equivalent set of AH-sequences defined by LSFs and some complex constants. This alternate representation format of any finite-length sequence is referred as its LSF-Model. For a finite-length sequence, one can obtain multiple LSF-Models by varying the model parameters. The LSF-Model, in time domain can be used to synthesize any arbitrary causal finite-length sequence in terms of its characteristic AH-sequences. In the frequency domain, the LSF-Model can be used to obtain the spectral samples of the sequence as a linear combination of spectra of its characteristic AH-sequences. We also summarize the utility of the LSF-Model in practical discrete signal processing systems.

Effective models of two-flavor QCD: finite $\mu$ and $m_q$-dependence

We study effective models of chiral fields and Polyakov loop expected to describe the dynamics responsible for the phase structure of two-flavor QCD at finite temperature and density. We consider chiral sector described either using linear sigma model or Nambu-Jona-Lasinio model and study the phase diagram and determine the location of the critical point as a function of the explicit chiral symmetry breaking (i.e. the bare quark mass $m_q$). We also discuss the possible emergence of the quarkyonic phase in this model.