High--Resolution 3D Simulations of Relativistic Jets

We have performed high-resolution 3D simulations of relativistic jets with beam flow Lorentz factors up to 7, a spatial resolution of 8 cells per beam radius, and for up to 75 normalized time units to study the morphology and dynamics of 3D relativistic jets. Our simulations show that the coherent fast backflows found in axisymmetric models are not present in 3D models. We further find that when the jet is exposed to non-axisymmetric perturbations, (i) it does not display the strong perturbations found for 3D classical hydrodynamic and MHD jets (at least during the period of time covered by our simulations), and (ii) it does propagate according to the 1D estimate. Small 3D effects in the relativistic beam give rise to a lumpy distribution of apparent speeds like that observed in M87. The beam is surrounded by a boundary layer of high specific internal energy. The properties of this layer are briefly discussed.Comment: 15 pages, 4 figures. Accepted to be publish in the ApJ Letters. Tar+gzip documen

Postbuckling behaviour of beams with discrete nonlinear restraints

A beam with nonlinearly ‐ elastic lateral restraints attached at discrete points along its span is investigated via analytical and numerical methods. Previous results for the critical moment and the deflected shape based on an eigenvalue analysis of a similar beam with linearly ‐ elastic restraints are discussed, along with a validation of these results against an equivalent finite element model and results from numerical continuation. A beam with nonlinearly ‐ elastic restraints is then analysed with treatments for both quadratic and cubic restraint force–displacement relationships being provided. After formulation of the potential energy functionals, the governing differential equations of the system are derived via the calculus of variations and appropriate boundary conditions are applied. The equations are then solved using the numerical continuation software AUTO ‐ 07p for a standard I ‐ section beam. The variation in elastic critical buckling moment with the linear component of the restraint stiffness is tracked via a two ‐ parameter numerical continuation, allowing determination of the stiffness values at which the critical buckling modes changes qualitatively. Using these stiffness values, subsequent analyses are conducted to examine the influence of the nonlinear component of the restraint stiffness, from which post ‐ buckling equilibrium paths and deformation modes are extracted. The results of these analyses are then compared with an equivalent Rayleigh–Ritz formulation whereby the displacement components are represented by Fourier series. Equilibrium equations are derived by minimizing the potential energy functional with respect to the amplitudes of the constituent harmonics of the Fourier series. The amplitudes are solved for in the post ‐ buckling range by AUTO ‐ O7p and equilibrium paths are produced and compared to the equivalent solutions of the differential equations, with good agreement observed

The evolution-dominated hydrodynamic model and the pseudorapidity distributions in high energy physics

By taking into account the effects of leading particles, we discuss the pseudorapidity distributions of the charged particles produced in high energy heavy ion collisions in the context of evolution-dominated hydrodynamic model. The leading particles are supposed to have a Gaussian rapidity distribution normalized to the number of participants. A comparison is made between the theoretical results and the experimental measurements performed by BRAHMS and PHOBOS Collaboration at BNL-RHIC in Au-Au and Cu-Cu collisions at sqrt(s_NN) =200 GeV and by ALICE Collaboration at CERN-LHC in Pb-Pb collisions at sqrt(s_NN) =2.76 TeV.Comment: 17 pages,4 figures, 2 table

Improving the viability of mental models held by novice programmers

Recent research has found that many novice programmers often hold non-viable mental models of basic programming concepts such as assignment and object reference. This paper proposes a constructivist-based teaching model, integrating a cognitive conflict strategy with program visualization, with the aim of improving novice programmers’ mental models. The results of a preliminary empirical study suggest that, for the relatively straightforward concept of assignment, tight integration of program visualization with a cognitive conflict event that highlights a student’s inappropriate understanding can help improve students’ non-viable mental models. 14 out of 18 participants who held non-viable mental models of the assignment process successfully changed their model to be viable as a result of the proposed teaching model

Origin of sawtooth domain walls in ferroelectrics

Domains and domain walls are among the key factors that determine the performance of ferroelectric materials. In recent years, a unique type of domain walls, i.e., the sawtooth-shaped domain walls, has been observed in BiFeO$_{3}$ and PbTiO$_{3}$. Here, we build a minimal model to reveal the origin of these sawtooth-shaped domain walls. Incorporating this model into Monte-Carlo simulations shows that (i) the competition between the long-range Coulomb interaction (due to bound charges) and short-range interaction (due to opposite dipoles) is responsible for the formation of these peculiar domain walls and (ii) their relative strength is critical in determining the periodicity of these sawtooth-shaped domain walls. Necessary conditions to form such domain walls are also discussed

Frequency response in short thermocouple wires

Theoretical expressions are derived for the steady state frequency response of a thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for a nonuniform wire with unequal material properties and wire diameters across the junction. The amplitude ratio at low frequency omega approaches 0 agrees with the results of Scadron and Warshawsky (1952) for a steady state temperature distribution. Moreover, the frequency response for a nonuniform wire in the limit of infinite length l approaches infinity is shown to reduce to a simple expression that is analogous to the classic first order solution for a thermocouple wire with uniform properties. Theoretical expressions are also derived for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire and a two material wire with unequal material properties across the junction. For the case of a one material supported wire, an exact solution is derived which compares favorably with an approximate expression that only matches temperatures at the support junction. Moreover, for the case of a two material supported wire, an analytical expression is derived that closely correlates numerical results. Experimental measurements are made for the steady state frequency response of a supported thermocouple wire. In particular, the effects of axial heat conduction are demonstrated for both a supported one material wire (type K) and a two material wire (type T) with unequal material properties across the junction. The data for the amplitude ratio and phase angle are correlated to within 10 pct. with the theoretical predictions of Forney and Fralick (1991). This is accomplished by choosing a natural frequency omega sub n for the wire data to correlate the first order response at large gas temperature frequencies. It is found that a large bead size, however, will increase the amplitude ratio at low frequencies but decrease the natural frequency of the wire. The phase angle data are also distorted for imperfect junctions

Gluon GPDs and Exclusive Photoproduction of a Quarkonium in Forward Region

Forward photoproduction of $J/\psi$ can be used to extract Generalized Parton Distributions(GPD's) of gluons. We analyze the process at twist-3 level and study relevant classifications of twist-3 gluon GPD's. At leading power or twist-2 level the produced $J/\psi$ is transversely polarized. We find that at twist-3 the produced $J/\psi$ is longitudinally polarized. Our study shows that in high energy limit the twist-3 amplitude is only suppressed by the inverse power of the heavy quark mass relatively to the twist-2 amplitude. This indicates that the power correction to the cross-section of unpolarized $J/\psi$ can have a sizeable effect. We have also derived the amplitude of the production of $h_c$ at twist-3, but the result contains end-point singularities. The production of other quarkonia has been briefly discussed.Comment: Discussions of results are adde