Consistent dirichlet boundary conditions for numerical solution of moving boundary problems

We consider the imposition of Dirichlet boundary conditions in the finite element mod-elling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass

Scale-invariant moving finite elements for nonlinear partial differential equations in two dimensions

A scale-invariant moving finite element method is proposed for the adaptive solution of nonlinear partial differential equations. The mesh movement is based on a finite element discretisation of a scale-invariant conservation principle incorporating a monitor function, while the time discretisation of the resulting system of ordinary differential equations is carried out using a scale-invariant time-stepping which yields uniform local accuracy in time. The accuracy and reliability of the algorithm are successfully tested against exact self-similar solutions where available, and otherwise against a state-of-the-art h-refinement scheme for solutions of a two-dimensional porous medium equation problem with a moving boundary. The monitor functions used are the dependent variable and a monitor related to the surface area of the solution manifold

A moving mesh finite element algorithm for the adaptive solution of time-dependent partial differential equations with moving boundaries

A moving mesh finite element algorithm is proposed for the adaptive solution of nonlinear diffusion equations with moving boundaries in one and two dimensions. The moving mesh equations are based upon conserving a local proportion, within each patch of finite elements, of the total “mass” that is present in the projected initial data. The accuracy of the algorithm is carefully assessed through quantitative comparison with known similarity solutions, and its robustness is tested on more general problems. Applications are shown to a variety of problems involving time-dependent partial differential equations with moving boundaries. Problems which conserve mass, such as the porous medium equation and a fourth order nonlinear diffusion problem, can be treated by a simplified form of the method, while problems which do not conserve mass require the full theory

Influence of fast interstellar gas flow on dynamics of dust grains

The orbital evolution of a dust particle under the action of a fast interstellar gas flow is investigated. The secular time derivatives of Keplerian orbital elements and the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle's orbit are derived. The secular time derivatives of the semi-major axis, eccentricity, and of the radial, transversal, and normal components of the gas flow velocity vector at the pericentre of the particle's orbit constitute a system of equations that determines the evolution of the particle's orbit in space with respect to the gas flow velocity vector. This system of differential equations can be easily solved analytically. From the solution of the system we found the evolution of the Keplerian orbital elements in the special case when the orbital elements are determined with respect to a plane perpendicular to the gas flow velocity vector. Transformation of the Keplerian orbital elements determined for this special case into orbital elements determined with respect to an arbitrary oriented plane is presented. The orbital elements of the dust particle change periodically with a constant oscillation period or remain constant. Planar, perpendicular and stationary solutions are discussed. The applicability of this solution in the Solar system is also investigated. We consider icy particles with radii from 1 to 10 micrometers. The presented solution is valid for these particles in orbits with semi-major axes from 200 to 3000 AU and eccentricities smaller than 0.8, approximately. The oscillation periods for these orbits range from 10^5 to 2 x 10^6 years, approximately.Comment: 22 pages, 3 figures; Accepted for publication in Celestial Mechanics and Dynamical Astronom

Origin and Evolution of Saturn's Ring System

The origin and long-term evolution of Saturn's rings is still an unsolved problem in modern planetary science. In this chapter we review the current state of our knowledge on this long-standing question for the main rings (A, Cassini Division, B, C), the F Ring, and the diffuse rings (E and G). During the Voyager era, models of evolutionary processes affecting the rings on long time scales (erosion, viscous spreading, accretion, ballistic transport, etc.) had suggested that Saturn's rings are not older than 100 My. In addition, Saturn's large system of diffuse rings has been thought to be the result of material loss from one or more of Saturn's satellites. In the Cassini era, high spatial and spectral resolution data have allowed progress to be made on some of these questions. Discoveries such as the ''propellers'' in the A ring, the shape of ring-embedded moonlets, the clumps in the F Ring, and Enceladus' plume provide new constraints on evolutionary processes in Saturn's rings. At the same time, advances in numerical simulations over the last 20 years have opened the way to realistic models of the rings's fine scale structure, and progress in our understanding of the formation of the Solar System provides a better-defined historical context in which to understand ring formation. All these elements have important implications for the origin and long-term evolution of Saturn's rings. They strengthen the idea that Saturn's rings are very dynamical and rapidly evolving, while new arguments suggest that the rings could be older than previously believed, provided that they are regularly renewed. Key evolutionary processes, timescales and possible scenarios for the rings's origin are reviewed in the light of tComment: Chapter 17 of the book ''Saturn After Cassini-Huygens'' Saturn from Cassini-Huygens, Dougherty, M.K.; Esposito, L.W.; Krimigis, S.M. (Ed.) (2009) 537-57

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Jet size dependence of single jet suppression in lead-lead collisions at sqrt(s(NN)) = 2.76 TeV with the ATLAS detector at the LHC

Measurements of inclusive jet suppression in heavy ion collisions at the LHC provide direct sensitivity to the physics of jet quenching. In a sample of lead-lead collisions at sqrt(s) = 2.76 TeV corresponding to an integrated luminosity of approximately 7 inverse microbarns, ATLAS has measured jets with a calorimeter over the pseudorapidity interval |eta| < 2.1 and over the transverse momentum range 38 < pT < 210 GeV. Jets were reconstructed using the anti-kt algorithm with values for the distance parameter that determines the nominal jet radius of R = 0.2, 0.3, 0.4 and 0.5. The centrality dependence of the jet yield is characterized by the jet "central-to-peripheral ratio," Rcp. Jet production is found to be suppressed by approximately a factor of two in the 10% most central collisions relative to peripheral collisions. Rcp varies smoothly with centrality as characterized by the number of participating nucleons. The observed suppression is only weakly dependent on jet radius and transverse momentum. These results provide the first direct measurement of inclusive jet suppression in heavy ion collisions and complement previous measurements of dijet transverse energy imbalance at the LHC.Comment: 15 pages plus author list (30 pages total), 8 figures, 2 tables, submitted to Physics Letters B. All figures including auxiliary figures are available at http://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/HION-2011-02

Measurement of the polarisation of W bosons produced with large transverse momentum in pp collisions at sqrt(s) = 7 TeV with the ATLAS experiment

This paper describes an analysis of the angular distribution of W->enu and W->munu decays, using data from pp collisions at sqrt(s) = 7 TeV recorded with the ATLAS detector at the LHC in 2010, corresponding to an integrated luminosity of about 35 pb^-1. Using the decay lepton transverse momentum and the missing transverse energy, the W decay angular distribution projected onto the transverse plane is obtained and analysed in terms of helicity fractions f0, fL and fR over two ranges of W transverse momentum (ptw): 35 < ptw < 50 GeV and ptw > 50 GeV. Good agreement is found with theoretical predictions. For ptw > 50 GeV, the values of f0 and fL-fR, averaged over charge and lepton flavour, are measured to be : f0 = 0.127 +/- 0.030 +/- 0.108 and fL-fR = 0.252 +/- 0.017 +/- 0.030, where the first uncertainties are statistical, and the second include all systematic effects.Comment: 19 pages plus author list (34 pages total), 9 figures, 11 tables, revised author list, matches European Journal of Physics C versio

Observation of a new chi_b state in radiative transitions to Upsilon(1S) and Upsilon(2S) at ATLAS

The chi_b(nP) quarkonium states are produced in proton-proton collisions at the Large Hadron Collider (LHC) at sqrt(s) = 7 TeV and recorded by the ATLAS detector. Using a data sample corresponding to an integrated luminosity of 4.4 fb^-1, these states are reconstructed through their radiative decays to Upsilon(1S,2S) with Upsilon->mu+mu-. In addition to the mass peaks corresponding to the decay modes chi_b(1P,2P)->Upsilon(1S)gamma, a new structure centered at a mass of 10.530+/-0.005 (stat.)+/-0.009 (syst.) GeV is also observed, in both the Upsilon(1S)gamma and Upsilon(2S)gamma decay modes. This is interpreted as the chi_b(3P) system.Comment: 5 pages plus author list (18 pages total), 2 figures, 1 table, corrected author list, matches final version in Physical Review Letter