Shell Model Study of Local and Global Energy Barriers in KDP

We perform a study of the energetics of KH2PO4 (KDP) by using a shell model (SM) which was constructed by adjusting the interaction parameters to ab initio calculations, and was fitted to reproduce phonons, polarization-inversion energies and structural properties. We calculate the energy profiles by performing global displacements and local distortions following the ferroelectric (FE) mode pattern in clusters of different sizes embedded in a paraelectric (PE) phase matrix. These properties are expected to be relevant to the PE-FE phase transition. The obtained SM results are compared to corresponding ab initio (AI) data. The global instabilities are found in good agreement for both KDP and DKDP. We also find qualitative good agreement in the KDP structure and even quantitative agreement in the expanded DKDP structure for the local distortions. The SM results reproduce well different trends like increasing instabilities as the cluster sizes grows, as the heavier atoms are included, and as the volume is increased, in accordance with the corresponding data from AI calculations.

Nonmesonic weak decay of Λ-hypernuclei within independent-particle shell-model

After a short introduction to the nonmesonic weak decay (NMWD) ΛN→nN of Λ-hypernuclei we discuss the long-standing puzzle on the ratio Γn/Γp, and some recent experimental evidences that signalized towards its final solution. Two versions of the Independent-Particle-Shell-Model (IPSM) are employed to account for the nuclear structure of the final residual nuclei. They are: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into account, and (b) IPSM-b, where the highly excited hole states are considered to be quasi-stationary and are described by Breit-Wigner distributions, whose widths are estimated from the experimental data. We evaluate the coincidence spectra in Λ 4He, Λ 5He, Λ 12C, Λ 16O, and Λ 28Si, as a function of the sum of kinetic energies EnN=En+EN for N=n, p. The recent Brookhaven National Laboratory experiment E788 on Λ 4He, is interpreted within the IPSM. We found that the shapes of all the spectra are basically tailored by the kinematics of the corresponding phase space, depending very weakly on the dynamics, which is gauged here by the one-meson-exchange- potential. In spite of the straightforwardness of the approach a good agreement with data is achieved. This might be an indication that the final-state- interactions and the two-nucleon induced processes are not very important in the decay of this hypernucleus. We have also found that the π+K exchange potential with soft vertex-form-factor cutoffs (Λπ≈0. 7GeV, ΛK≈0.9GeV), is able to account simultaneously for the available experimental data related to Γp and Γn for Λ 4H, and Λ 5He. © 2010 American Institute of Physics.

Soft tissue human thigh and buttock finite element model to simulate vehicle seat cushion indentation

Accurate modelling of automotive occupant posture is strongly related to the mechanical interaction between human body soft tissue and flexible seat components. This paper presents a finite-element study simulating the deflection of seat cushion foam and supportive seat structures, as well as human buttock and thigh soft tissue when seated. The thigh-buttock surface shell model was based on 95th percentile male subject scan data and made of two layers, covering thin to moderate thigh and buttock proportions. To replicate the effects of skin and fat, the neoprene rubber layer was modelled as a hyperelastic material with viscoelastic behaviour. The analytical seat model is based on a Ford production seat. The result of the finite-element indentation simulation is compared to a previous simulation of an indentation with a hard shell human model of equal geometry, and to the physical indentation result. We conclude that SAE composite buttock form and human-seat indentation of a suspended seat cushion can be validly simulated.

Shell effect and capture cross sections in the synthesis of superheavy nuclei

The shell effect is included in the improved isospin dependent quantum molecular dynamics model in which the shell correction energy of the system is calculated by using the deformed two-center shell model. A switch function is introduced to connect the shell correction energy of the projectile and the target with that of the compound nucleus during the dynamical fusion process. It is found that the calculated capture cross sections reproduce the experimental data quantitatively at the energy near the Coulomb barrier. The capture cross sections for reaction (35) (80) Br + (82) (208) Pb -> (117) (288) X are also calculated and discussed.

Structure of upper-g(9/2)-shell nuclei and shape effect in the Ag-94 isomeric states

Using a shell model which is capable of describing the spectra of upper g(9/2)-shell nuclei close to the N = Z line, we study the structure of two isomeric states 7(+) and 21(+) in the odd-odd N = Z nucleus Ag-94. It is found that both isomeric states exhibit a large collectivity. The 7(+) state is oblately deformed, and is suggested to be a shape isomer in nature. The 21(+) state becomes isomeric because of level inversion of the 19(+) and 21(+) states due to core excitations across the N = Z = 50 shell gap. Calculation of spectroscopic quadrupole moment indicates clearly an enhancement in these states due to the core excitations. However, the present shell model calculation that produces the 19(+)-21(+) level inversion cannot accept the large-deformation picture of Mukha et al.

Classification of states in S0(8) proton-neutron pairing model

lsoscalar (T = 0) plus isovector (T = 1) pairing Hamiltonian in LS-coupling. which is important for heavy N = Z nuclei, is solvable in terms of a SO(8) Lie algebra for three special values of the mixing parameter that measures the competition between the T = 0 aid T = 1 pairing. The SO(8) algebra is generated, amongst others, by the S = 1, T = 0 and S = 0, T = 1 pair creation and annihilation operators and corresponding to the three values of the mixing parameter, there are three chains of subalgebras: SO(8) superset of SOST (6) superset of SOS(3) circle times SOT(3), SO(8) superset of [SOS(5) superset of SOS(3)] circle times SOT(3) and SO(8) superset of [SOT(5) superset of SOT(3)] circle times SOS(3). Shell model Lie algebras, with only particle number conserving generators, that are complementary to these three chains of subalgebras are identified and they are used in the classification of states for a given number of nucleons. The classification problem is solved explicitly tor states with SO(8) seniority nu = 0, 1, 2, 3 and 4. Using them, hand structures in isospin space are identified for states with nu = 0, 1, 2 and 3. (c) 2005 Elsevier B.V. All rights reserved.

Finite Element Modeling of Shell Shape in the Freshwater Turtle Pseudemys concinna Reveals a Trade-Off between Mechanical Strength and Hydrodynamic Efficiency

Aquatic species can experience different selective pressures on morphology in different flow regimes. Species inhabiting lotic regimes often adapt to these conditions by evolving low-drag (i.e., streamlined) morphologies that reduce the likelihood of dislodgment or displacement. However, hydrodynamic factors are not the only selective pressures influencing organismal morphology and shapes well suited to flow conditions may compromise performance in other roles. We investigated the possibility of morphological trade-offs in the turtle Pseudemys concinna. Individuals living in lotic environments have flatter, more streamlined shells than those living in lentic environments; however, this flatter shape may also make the shells less capable of resisting predator-induced loads. We tested the idea that ‘‘lotic’’ shell shapes are weaker than ‘‘lentic’’ shell shapes, concomitantly examining effects of sex. Geometric morphometric data were used to transform an existing finite element shell model into a series of models corresponding to the shapes of individual turtles. Models were assigned identical material properties and loaded under identical conditions, and the stresses produced by a series of eight loads were extracted to describe the strength of the shells. ‘‘Lotic’’ shell shapes produced significantly higher stresses than ‘‘lentic’’ shell shapes, indicating that the former is weaker than the latter. Females had significantly stronger shell shapes than males, although these differences were less consistent than differences between flow regimes. We conclude that, despite the potential for many-to-one mapping of shell shape onto strength, P. concinna experiences a trade-off in shell shape between hydrodynamic and mechanical performance. This trade-off may be evident in many other turtle species or any other aquatic species that also depend on a shell for defense. However, evolution of body size may provide an avenue of escape from this trade-off in some cases, as changes in size can drastically affect mechanical performance while having little effect on hydrodynamic performance.

Out-of-plane deformation and failure of masonry walls with various forms of reinforcement

Out-of-plane behaviour of mortared and mortarless masonry walls with various forms of reinforcement, including unreinforced masonry as a base case is examined using a layered shell element based explicit finite element modelling method. Wall systems containing internal reinforcement, external surface reinforcement and intermittently laced reinforced concrete members and unreinforced masonry panels are considered. Masonry is modelled as a layer with macroscopic orthotropic properties; external reinforcing render, grout and reinforcing bars are modelled as distinct layers of the shell element. Predictions from the layered shell model have been validated using several out-of-plane experimental datasets reported in the literature. The model is used to examine the effectiveness of two retrofitting schemes for an unreinforced masonry wall.

An explicit finite element modelling method for masonry walls under out-of-plane loading

An explicit finite element modelling method is formulated using a layered shell element to examine the behaviour of masonry walls subject to out-of-plane loading. Masonry is modelled as a homogenised material with distinct directional properties that are calibrated from datasets of a “C” shaped wall tested under pressure loading applied to its web. The predictions of the layered shell model have been validated using several out-of-plane experimental datasets reported in the literature. Profound influence of support conditions, aspect ratio, pre-compression and opening to the strength and ductility of masonry walls is exhibited from the sensitivity analyses performed using the model.

Threshold photoproduction of charged pions on14N

Pion photoproduction processes14Ngs(gamma, pgr +)14C and14Ngs(gamma, pgr –)14O have been studied in the threshold region. These processes provide an excellent tool to study the corrections to soft pion theorems and Kroll-Ruderman limit as applied to nuclear processes. The agreement with the available experimental data for these processes is better with the empirical wave functions while the shell-model wave functions predict a much higher value. Detailed experimental studies of these reactions at threshold, it is shown, are expected to lead to a better understanding of the shell-model inputs and radial distributions in the 1p state. We thank Dr. S.C.K. Nair for a helpful discussion during the initial stages of this work. One of us (MVN) thanks Dr. J.M. Laget for sending some unpublished data on pion photoproduction. He is also thankful to Dr. J. Pasupathy and Dr. R. Rajaraman for their interest and encouragement.

Raman spectrum of rubidium iodide

Raman spectrum of rubidium iodide has been recorded for the first time using the resonance radiation of mercury (λ 2537 ) as the exciter. The frequencies of the 24p limiting modes (p = 2, the number of non-equivalent atoms in the unit cell), postulated by Raman in 1943, which correspond to the frequencies from the critical points Γ, L and X, have been worked out using the shell model of Cochran, taking into account the nearest and the next-nearest neighbour short-range interactions and the polarization of both the ions. The observed Raman lines have been assigned to the overtones and the combinations of the phonon branches from Γ, L and X.

The plastic collapse and energy absorption capacity of egg-box panels

The plastic collapse response of aluminium egg-box panels subjected to out-of-plane compression has been measured and modelled. It is observed that the collapse strength and energy absorption are sensitive to the level of in-plane constraint, with collapse dictated either by plastic buckling or by a travelling plastic knuckle mechanism. Drop weight tests have been performed at speeds of up to 6 m s-1, and an elevation in strength with impact velocity is noted. A 3D finite element shell model is needed in order to reproduce the observed behaviours. Additional calculations using an axisymmetric finite element model give the correct collapse modes but are less accurate than the more sophisticated 3D model. The finite element simulations suggest that the observed velocity dependence of strength is primarily due to strain-rate sensitivity of the aluminium sheet, with material inertia playing a negligible role. Finally, it is shown that the energy absorption capacity of the egg-box material is comparable to that of metallic foams. © 2003 Elsevier Ltd. All rights reserved.

Surface-waves, stability, and scattering for a lined duct with flow

We consider a straight cylindrical duct with a steady subsonic axial flow and a reacting boundary (e.g. an acoustic lining). The wave modes are separated into ordinary acoustic duct modes, and surface modes confined to a small neighbourhood of the boundary. Many researchers have used a mass-spring-damper boundary model, for which one surface mode has previously been identified as a convective instability; however, we show the stability analysis used in such cases to be questionable. We investigate instead the stability of the surface modes using the Briggs-Bers criterion for a Flügge thin-shell boundary model. For modest frequencies and wavenumbers the thin-shell has an impedance which is effectively that of a mass-spring-damper, although for the large wavenumbers needed for the stability analysis the thin-shell and mass-spring-damper impedances diverge, owing to the thin shell's bending stiffness. The thin shell model may therefore be viewed as a regularization of the mass-spring-damper model which accounts for nonlocally-reacting effects. We find all modes to be stable for realistic thin-shell parameters, while absolute instabilities are demonstrated for extremely thin boundary thicknesses. The limit of vanishing bending stiffness is found to be a singular limit, yielding absolute instabilities of arbitrarily large temporal growth rate. We propose that the problems with previous stability analyses are due to the neglect of something akin to bending stiffness in the boundary model. Our conclusion is that the surface mode previously identified as a convective instability may well be stable in reality. Finally, inspired by Rienstra's recent analysis, we investigate the scattering of an acoustic mode as it encounters a sudden change from a hard-wall to a thin-shell boundary, using a Wiener-Hopf technique. The thin-shell is considered to be clamped to the hard-wall. The acoustic mode is found to scatter into transmitted and reflected acoustic modes, and surface modes strongly linked to the solid waves in the boundary, although no longitudinal or transverse waves within the boundary are excited. Examples are provided that demonstrate total transmission, total reflection, and a combination of the two. This thin-shell scattering problem is preferable to the mass-spring-damper scattering problem presented by Rienstra, since the thin-shell problem is fully determined and does not need to appeal to a Kutta-like condition or the inclusion of an instability in order to avoid a surface-streamline cusp at the boundary change.