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.

Measuring the Cosmic Microwave Background : Topics along the analysis pipeline

The first quarter of the 20th century witnessed a rebirth of cosmology, study of our Universe, as a field of scientific research with testable theoretical predictions. The amount of available cosmological data grew slowly from a few galaxy redshift measurements, rotation curves and local light element abundances into the first detection of the cos- mic microwave background (CMB) in 1965. By the turn of the century the amount of data exploded incorporating fields of new, exciting cosmological observables such as lensing, Lyman alpha forests, type Ia supernovae, baryon acoustic oscillations and Sunyaev-Zeldovich regions to name a few. -- CMB, the ubiquitous afterglow of the Big Bang, carries with it a wealth of cosmological information. Unfortunately, that information, delicate intensity variations, turned out hard to extract from the overall temperature. Since the first detection, it took nearly 30 years before first evidence of fluctuations on the microwave background were presented. At present, high precision cosmology is solidly based on precise measurements of the CMB anisotropy making it possible to pinpoint cosmological parameters to one-in-a-hundred level precision. The progress has made it possible to build and test models of the Universe that differ in the way the cosmos evolved some fraction of the first second since the Big Bang. -- This thesis is concerned with the high precision CMB observations. It presents three selected topics along a CMB experiment analysis pipeline. Map-making and residual noise estimation are studied using an approach called destriping. The studied approximate methods are invaluable for the large datasets of any modern CMB experiment and will undoubtedly become even more so when the next generation of experiments reach the operational stage. -- We begin with a brief overview of cosmological observations and describe the general relativistic perturbation theory. Next we discuss the map-making problem of a CMB experiment and the characterization of residual noise present in the maps. In the end, the use of modern cosmological data is presented in the study of an extended cosmological model, the correlated isocurvature fluctuations. Current available data is shown to indicate that future experiments are certainly needed to provide more information on these extra degrees of freedom. Any solid evidence of the isocurvature modes would have a considerable impact due to their power in model selection.

Structure and optical properties of Cd-substituted ZnO (Zn1-xCdxO) nanostructures synthesized by the high-pressure solution route

We report the synthesis of Cd-substituted ZnO nanostructures (Zn1-xCdxO with x up to approximate to 0.09) by the high-pressure solution growth method. The synthesized nanostructures comprise nanocrystals that are both particles (similar to 10-15 nm) and rods which grow along the [002] direction as established by transmission electron microscope (TEM) and x-ray diffraction (XRD) analysis. Rietveld analysis of the XRD data shows a monotonic increase of the unit cell volume with the increase of Cd concentration. The optical absorption, as well as the photoluminescence (PL), shows a red shift on Cd substitution. The line width of the PL spectrum is related to the strain inhomogeneity and it peaks in the region where the CdO phase separates from the Zn1-xCdxO nanostructures. The time-resolved photoemission showed a long-lived (similar to 10 ns) component. We propose that the PL behaviour of the Zn1-xCdxO is dominated by strain in the sample with the red shift of the PL linked to the expansion of the unit cell volume on Cd substitution.

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.

Modified density matrix renormalization group algorithm for the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange: Comparison with field theory at large J(2)/J(1)

A modified density matrix renormalization group (DMRG) algorithm is applied to the zigzag spin-1/2 chain with frustrated antiferromagnetic exchange J(1) and J(2) between first and second neighbors. The modified algorithm yields accurate results up to J(2)/J(1) approximate to 4 for the magnetic gap Delta to the lowest triplet state, the amplitude B of the bond order wave phase, the wavelength lambda of the spiral phase, and the spin correlation length xi. The J(2)/J(1) dependences of Delta, B, lambda, and xi provide multiple comparisons to field theories of the zigzag chain. The twist angle of the spiral phase and the spin structure factor yield additional comparisons between DMRG and field theory. Attention is given to the numerical accuracy required to obtain exponentially small gaps or exponentially long correlations near a quantum phase transition.

Existence and approximations of solutions to some fractional order functional integral equations

In this paper we shall study a fractional order functional integral equation. In the first part of the paper, we proved the existence and uniqueness of mile and global solutions in a Banach space. In the second part of the paper, we used the analytic semigroups theory oflinear operators and the fixed point method to establish the existence, uniqueness and convergence of approximate solutions of the given problem in a separable Hilbert space. We also proved the existence and convergence of Faedo-Galerkin approximate solution to the given problem. Finally, we give an example.

Conformation of the flexible bent helix of Leu1-zervamicin in crystal C and a possible gating action for ion passage

The membrane channel-forming polypeptide, Leu(1)-zervamicin, Ac-Leu-Ile-Gln-Iva-Ile(5)-Thr-Aib-Leu-Aib-Hyp(10) -Gln-Aib-Hyp-Aib-Pro(15)-Phol (Aib: alpha-aminoisobutyric acid; Iva: isovaline; Hyp: 4-hydroxyproline; Phol: phenylalininol) has been analyzed by x-ray diffraction in a third crystal form. Although the bent helix is quite similar to the conformations found in crystals A and B, the amount of bending is more severe with a bending angle approximate to 47 degrees, The water channel formed by the convex polar faces of neighboring helices is larger at the mouth than in crystals A and B, and the water sites have become disordered. The channel is interrupted in the middle by a hydrogen bond between the OH of Hyp(10) and the NH2 of the Gln(11) of a neighboring molecule. The side chain of Gln(11) is wrapped around the helix backbone in an unusual fashion in order that it can augment the polar side of the helix. In the present crystal C there appears to be an additional conformation for the Gln(11) side chain (with approximate to 20% occupancy) that opens the channel for possible ion passage. Structure parameters for C85H140N18O22.xH(2)O.C2H5OH are space group P2(1)2(1)2(1), a = 10.337 (2) Angstrom, b = 28.387 (7) Angstrom, c = 39.864 (11) Angstrom, Z = 4, agreement factor R = 12.99% for 3250 data observed > 3 sigma(F), resolution = 1.2 Angstrom. (C) 1994 John Wiley & Sons, Inc.

Structure and optical properties of Cd-substituted ZnO (Zn1-xCdxO) nanostructures synthesized by the high-pressure solution route

We report the synthesis of Cd-substituted ZnO nanostructures (Zn1-xCdxO with x up to approximate to 0.09) by the high-pressure solution growth method. The synthesized nanostructures comprise nanocrystals that are both particles (similar to 10-15 nm) and rods which grow along the [002] direction as established by transmission electron microscope (TEM) and x-ray diffraction (XRD) analysis. Rietveld analysis of the XRD data shows a monotonic increase of the unit cell volume with the increase of Cd concentration. The optical absorption, as well as the photoluminescence (PL), shows a red shift on Cd substitution. The line width of the PL spectrum is related to the strain inhomogeneity and it peaks in the region where the CdO phase separates from the Zn1-xCdxO nanostructures. The time-resolved photoemission showed a long-lived (similar to 10 ns) component. We propose that the PL behaviour of the Zn1-xCdxO is dominated by strain in the sample with the red shift of the PL linked to the expansion of the unit cell volume on Cd substitution.

Solution of MDPS using simulation-based value iteration

This article proposes a three-timescale simulation based algorithm for solution of infinite horizon Markov Decision Processes (MDPs). We assume a finite state space and discounted cost criterion and adopt the value iteration approach. An approximation of the Dynamic Programming operator T is applied to the value function iterates. This 'approximate' operator is implemented using three timescales, the slowest of which updates the value function iterates. On the middle timescale we perform a gradient search over the feasible action set of each state using Simultaneous Perturbation Stochastic Approximation (SPSA) gradient estimates, thus finding the minimizing action in T. On the fastest timescale, the 'critic' estimates, over which the gradient search is performed, are obtained. A sketch of convergence explaining the dynamics of the algorithm using associated ODEs is also presented. Numerical experiments on rate based flow control on a bottleneck node using a continuous-time queueing model are performed using the proposed algorithm. The results obtained are verified against classical value iteration where the feasible set is suitably discretized. Over such a discretized setting, a variant of the algorithm of [12] is compared and the proposed algorithm is found to converge faster.

On the low temperature deformation mechanism in polycrystalline cadmium

The thermally activated plastic flow of polycrystalline cadmium was investigated by differentialstress creep tests at 86°K and tensile tests in the temperature range 86°–473°K. The activation energy (0.55 eV) at zero effective stress and the activation volume as a function of effective stress were obtained. It is concluded that intersection of glide and forest dislocations becomes rate controlling for low temperature deformation. The approximate stacking-fault width in cadmium is deduced to be “1.5b”.

Effect of static deflection on natural frequency of non-linear spring mass system by direct linearization method

An exact expression for the frequency of a non-linear cubic spring mass system is obtained considering the effect of static deflection. An alternative expression for the approximate frequency is also obtained by the direct linearization procedure; it is shown that this is very accurate as compared with the exact method. This approximate frequency equation is used to explain a “dual behaviour” of the frequency amplitude curves.

On the design of quadratic weirs

This paper considers the problem of the design of the quadratic weir notch, which finds application in the proportionate method of flow measurement in a by-pass, such that the discharge through it is proportional to the square root of the head measured above a certain datum. The weir notch consists of a bottom in the form of a rectangular weir of width 2W and depth a over which a designed curve is fitted. A theorem concerning the flow through compound weirs called the “slope discharge continuity theorem” is discussed and proved. Using this, the problem is reduced to the determination of an exact solution to Volterra's integral equation in Abel's form. It is shown that in the case of a quadratic weir notch, the discharge is proportional to the square root of the head measured above a datum Image a above the crest of the weir. Further, it is observed that the function defining the shape of the weir is rapidly convergent and its value almost approximates to zero at distances of 3a and above from the crest of the weir. This interesting and significant behaviour of the function incidentally provides a very good approximate solution to a particular Fredholm integral equation of the first kind, transforming the notch into a device called a “proportional-orifice”. A new concept of a “notch-orifice” capable of passing a discharge proportional to the square root of the head (above a particular datum) while acting both as a notch, and as an orifice, is given. A typical experiment with one such notch-orifice, having A = 4 in., and W = 6 in., shows a remarkable agreement with the theory and is found to have a constant coefficient of discharge of 0.61 in the ranges of both notch and orifice.

Primordial Perturbations from a Self-interacting Curvaton

Inflation is a period of accelerated expansion in the very early universe, which has the appealing aspect that it can create primordial perturbations via quantum fluctuations. These primordial perturbations have been observed in the cosmic microwave background, and these perturbations also function as the seeds of all large-scale structure in the universe. Curvaton models are simple modifications of the standard inflationary paradigm, where inflation is driven by the energy density of the inflaton, but another field, the curvaton, is responsible for producing the primordial perturbations. The curvaton decays after inflation as ended, where the isocurvature perturbations of the curvaton are converted into adiabatic perturbations. Since the curvaton must decay, it must have some interactions. Additionally realistic curvaton models typically have some self-interactions. In this work we consider self-interacting curvaton models, where the self-interaction is a monomial in the potential, suppressed by the Planck scale, and thus the self-interaction is very weak. Nevertheless, since the self-interaction makes the equations of motion non-linear, it can modify the behaviour of the model very drastically. The most intriguing aspect of this behaviour is that the final properties of the perturbations become highly dependent on the initial values. Departures of Gaussian distribution are important observables of the primordial perturbations. Due to the non-linearity of the self-interacting curvaton model and its sensitivity to initial conditions, it can produce significant non-Gaussianity of the primordial perturbations. In this work we investigate the non-Gaussianity produced by the self-interacting curvaton, and demonstrate that the non-Gaussianity parameters do not obey the analytically derived approximate relations often cited in the literature. Furthermore we also consider a self-interacting curvaton with a mass in the TeV-scale. Motivated by realistic particle physics models such as the Minimally Supersymmetric Standard Model, we demonstrate that a curvaton model within the mass range can be responsible for the observed perturbations if it can decay late enough.

Response spectrum of non-linear spring mass system subjected to transient disturbances

The transient response spectrum of a cubic spring mass system subjected to a step function input is obtained. An approximate method is adopted where non-linear restoring force characteristic is replaced by two linear segments, so that the mean square error between them is a minimum. The effect of viscous damping on the peak response is also discussed for various values of the damping constant and the non-linearity restoring force parameter.

Effect of boundary condition description on convergence of solution in a boundary-value problem

It is well known that the numerical accuracy of a series solution to a boundary-value problem by the direct method depends on the technique of approximate satisfaction of the boundary conditions and on the stage of truncation of the series. On the other hand, it does not appear to be generally recognized that, when the boundary conditions can be described in alternative equivalent forms, the convergence of the solution is significantly affected by the actual form in which they are stated. The importance of the last aspect is studied for three different techniques of computing the deflections of simply supported regular polygonal plates under uniform pressure. It is also shown that it is sometimes possible to modify the technique of analysis to make the accuracy independent of the description of the boundary conditions.