• Particle physics has an ambitious and broad experimental programme for the coming decades. This programme requires large investments in detector hardware, either to build new facilities and experiments, or to upgrade existing ones. Similarly, it requires commensurate investment in the R&D of software to acquire, manage, process, and analyse the shear amounts of data to be recorded. In planning for the HL-LHC in particular, it is critical that all of the collaborating stakeholders agree on the software goals and priorities, and that the efforts complement each other. In this spirit, this white paper describes the R&D activities required to prepare for this software upgrade.
• ### On the Stability of Continuous-Discontinuous Galerkin Methods for Advection-Diffusion-Reaction Problems(1211.0907)

Nov. 5, 2012 math.NA
We consider a finite element method which couples the continuous Galerkin method away from internal and boundary layers with a discontinuous Galerkin method in the vicinity of layers. We prove that this consistent method is stable in the streamline diffusion norm if the convection field flows non-characteristically from the region of the continuous Galerkin to the region of the discontinuous Galerkin method. The stability properties of the coupled method are illustrated by numerical experiments.
• ### A Bayesian technique for improving the sensitivity of the atmospheric neutrino L/E analysis(1208.2899)

Aug. 17, 2012 hep-ex
This paper outlines a method for improving the precision of atmospheric neutrino oscillation measurements. One experimental signature for these oscillations is an observed deficit in the rate of $\nu_{\mu}$ charged-current interactions with an oscillatory dependence on $L_{\nu}/E_{\nu}$, where $L_{\nu}$ is the neutrino propagation distance, and $E_{\nu}$ is the neutrino energy. For contained-vertex atmospheric neutrino interactions, the $L_{\nu}/E_{\nu}$ resolution varies significantly from event to event. The precision of the oscillation measurement can be improved by incorporating information on $L_{\nu}/E_{\nu}$ resolution into the oscillation analysis. In the analysis presented here, a Bayesian technique is used to estimate the $L_{\nu}/E_{\nu}$ resolution of observed atmospheric neutrinos on an event-by-event basis. By separating the events into bins of $L_{\nu}/E_{\nu}$ resolution in the oscillation analysis, a significant improvement in oscillation sensitivity can be achieved.
• ### On Local Super-Penalization of Interior Penalty Discontinuous Galerkin Methods(1205.5672)

May 25, 2012 math.NA, cs.NA
We prove in an abstract setting that standard (continuous) Galerkin finite element approximations are the limit of interior penalty discontinuous Galerkin approximations as the penalty parameter tends to infinity. We apply this result to equations of non-negative characteristic form and the non-linear, time dependent system of incompressible miscible displacement. Moreover, we investigate varying the penalty parameter on only a subset of a triangulation and the effects of local super-penalization on the stability of the method, resulting in a partly continuous, partly discontinuous method in the limit. An iterative automatic procedure is also proposed for the determination of the continuous region of the domain without loss of stability of the method.
• ### Implementation of the Continuous-Discontinuous Galerkin Finite Element Method(1201.2878)

Jan. 13, 2012 math.NA, cs.NA
For the stationary advection-diffusion problem the standard continuous Galerkin method is unstable without some additional control on the mesh or method. The interior penalty discontinuous Galerkin method is stable but at the expense of an increased number of degrees of freedom. The hybrid method proposed in [5] combines the computational complexity of the continuous method with the stability of the discontinuous method without a significant increase in degrees of freedom. We discuss the implementation of this method using the finite element library deal.ii and present some numerical experiments.