Fem for Schrödinger equation with Rashba effect

A finite element method for treating two-dimensional electron systems with Rashba spin–orbit interaction is developed. The Rashba spin–orbit interaction removes spin degeneracy, so that each spin contributes to the conductance differently. By accounting for the connection between a system and leads, this method yields the conductance of a nanoscale quantum device for each spin state. As an example, this calculation method is applied to a model of a quantum point contact. The results of this calculation indicate conductance quantization and a large spin polarization. We discuss the estimated accuracies of this calculation.     

An evolutionary game based particle swarm optimization algorithm

Particle swarm optimization (PSO) is an evolutionary algorithm used extensively. This paper presented a new particle swarm optimizer based on evolutionary game (EGPSO). We map particles’ finding optimal solution in PSO algorithm to players’ pursuing maximum utility by choosing strategies in evolutionary games, using replicator dynamics to model the behavior of particles. And in order to overcome premature convergence a multi-start technique was introduced. Experimental results show that EGPSO can overcome premature convergence and has great performance of convergence property over traditional PSO.     

Coupled model of hygro-thermal behavior of concrete during fire

We present a nonlinear mathematical model for numerical analysis of the behavior of concrete subject to transient heating according to the standard ISO fire curve. The mathematical model consists of the balance equations (conservations laws) with boundary and initial conditions, constitutive laws and material data for concrete at high temperatures. The numerical algorithm based on finite element method for the numerical solution of the energy equation and finite difference method for the mass balance equations is presented. Distributions of temperature, saturation of water and water vapor pressure are demonstrated.     

A RELAXATION RESULT OF FUNCTIONALS IN THE SPACE SBD（Ω）

In this article,the authors obtain an integral representation for the relaxation of the functional
F（x,u,Ω）：={∫^f（x,u（x）,εu（x））dx Ω if u∈W^1,1（Ω,R^N）, ＋∞ otherwise, in the space of functions of bounded deformation,with respect to L^1-convergence.Here Eu represents the absolutely continuous part of the symmetrized distributional derivative Eu.f（x,p,ξ）satisfying weak convexity assumption.     

An explicit solution of the large deformation of a cantilever beam under point load at the free tip

The large deformation of a cantilever beam under point load at the free tip is investigated by an analytic method, namely the homotopy analysis method (HAM). The explicit analytic formulas for the rotation angle at the free tip are given, which provide a convenient and straightforward approach to calculate the vertical and horizontal displacements of a cantilever beam with large deformation. These explicit formulas are valid for most practical problems, thus providing a useful reference for engineering applications. The corresponding Mathematica code is given in the Appendix.     

Formulation and analysis of two energy-consistent methods for nonlinear elastodynamic frictional contact problems

Energy-conserving algorithms are necessary to solve nonlinear elastodynamic problems in order to recover long term time integration accuracy and stability. Furthermore, some physical phenomena (such as friction) can generate dissipation; then in this work, we present and analyse two energy-consistent algorithms for hyperelastodynamic frictional contact problems which are characterised by a conserving behaviour for frictionless impacts but also by an admissible frictional dissipation phenomenon. The first approach permits one to enforce, respectively, the Kuhn–Tucker and persistency conditions during each time step by combining an adapted continuation of the Newton method and a Lagrangean formulation. In addition the second method which is based on the work in [P. Hauret, P. Le Tallec, Energy-controlling time integration methods for nonlinear elastodynamics and low-velocity impact, Comput. Methods Appl. Mech. Eng. 195 (2006) 4890–4916] represents a specific penalisation of the unilateral contact conditions. Some numerical simulations are presented to underscore the conservative or dissipative behaviour of the proposed methods.     

Mixed Formulations of Bending Problems for Homogeneous Elastic Plates and Beams

Mixed formulations of bending problems for homogeneous plates (beams) are proposed, whose essence is that the deformation of a plate (beam) near its fixed boundary is described by the three74qmm240056/xxlarge8208.gif" alt="dash" align="MIDDLE" BORDER="0">dimensional elasticity equations, and the remaining part by the conventional equations of plate (beam) bending. At the interface between these regions, the solutions of these equations are joined. The mixed formulation allows one to describe the three74qmm240056/xxlarge8208.gif" alt="dash" align="MIDDLE" BORDER="0">dimensional stress state in the neighborhood of the fixed boundaries of plates (beams) and take into account the complex nature of the fixing conditions. Finite74qmm240056/xxlarge8208.gif" alt="dash" align="MIDDLE" BORDER="0">element implementation is more efficient for the mixed formulations of plate (beam) bending problems than for the well74qmm240056/xxlarge8208.gif" alt="dash" align="MIDDLE" BORDER="0">known three74qmm240056/xxlarge8208.gif" alt="dash" align="MIDDLE" BORDER="0">dimensional formulations.     

Competition Between Radial Expansion and Thickening in the Enlargement of an Intracranial Saccular Aneurysm

Rupture of intracranial saccular aneurysms is the leading cause of spontaneous subarachnoid hemorrhage, which results in significant morbidity and mortality. Although many have suggested that saccular aneurysms enlarge and rupture due to mechanical instabilities, our recent nonlinear analyses suggest that at least certain classes of aneurysms do not exhibit a quasi-static limit point instability or dynamic instabilities in response to periodic loading. Based on an increased understanding of the ubiquitous role of growth and remodeling within the vasculature and recent histopathological data on saccular aneurysms, it is hypothesized that a stress-mediated regulation of collagen turnover causes their enlargement. There is a need, however, for a theoretical framework to explore this and competing hypotheses. In this paper, we present a 2-D constrained mixture model for growth and remodeling of an ellipsoidally shaped saccular aneurysm and numerically simulate enlargement and changes in material symmetry in the aneurysmal wall. Results suggest that ellipsoidal aneurysms tend toward spherical shapes, and a competition between radial expansion and wall thickening plays a critical role in determining the stability of an enlarging lesion.     

Site-specific ϕ- and ψ-torsion angle determination in a uniformly/extensively 13C- and 15N-labeled peptide

A solid-state rotational-echo double resonance (REDOR) NMR method was introduced to identify the ?- and ψ-torsion angle from a ¹H–¹⁵N or ¹H–¹³C′ spin system of alanine-like residues in a selectively, uniformly, or extensively ¹⁵N-/¹³C-labeled peptide. When a C_α(i) or a ¹⁵N peak is site-specifically obtainable in the NMR spectrum of a uniformly ¹⁵N/¹³C-labeled sample system, the ψ- or ?-torsion angle specified by the conformational structure of peptide geometry involving ¹⁵N(i)–¹H_αi–¹⁵N(i + 1) or ¹³C′(i − 1)–¹H^Ni–¹³C′(i) spin system can be identified based on ¹³C_α- or ¹⁵N-detected ¹H_α–¹⁵N or ¹H_N–¹³C REDOR experiment. This method will conveniently be utilized to identify major secondary motifs, such as α-helix, β-sheet, and β-turn, from a uniformly ¹⁵N-/¹³C-labled peptide sample system. When tested on a ¹³C-/¹⁵N-labeled model system of a three amino acid peptide Gly–[U–¹³C, ¹⁵N]Ala–[U–¹³C, ¹⁵N]Leu, the ψ-angle of alanine obtained experimentally, ψ = −40 ± 30°, agreed reasonably well with the X-ray determined angle, ψ = −39°.     

Cross sections for (d–t) neutron interaction with germanium isotopes

The cross sections for (n,x)$(n,x)$ reactions with Ge isotopes were measured at (d–t) neutron energies around 14 MeV with the activation technique using metal discs of natural composition. Calculations of detector efficiency, incident neutron spectrum and correction factors were performed with the Monte Carlo technique (MCNP4C code). Cross sections data are presented for ⁷⁰Ge(n,2n$n,2n$)⁶⁹Ge, ⁷⁴Ge(n,α$n,\alpha$)^71mZn, ⁷⁶Ge(n,2n$n,2n$)^75(m + g)Ge, ⁷⁰Ge(n,p$n,p$)⁷⁰Ga and ⁷²Ge(n,2n$n,2n$)^71gGe reactions. The cross section results for ⁷²Ge(n,2n$n,2n$)^71gGe reaction were reported for the first time. Some other cross sections were obtained with higher precision, including the ⁷⁰Ge(n,p$n,p$)⁷⁰Ga reaction. Theoretical calculations of excitation functions were performed with the TALYS-1.0 code and compared with the experimental cross section values. Data were included in the EXFOR database.