Local government authority attitudes to road traffic CO₂ emissions modelling: a British case study

Local government authorities (LGAs) play a key role in facilitating mitigation of road traffic CO2 emissions and must engage in emissions modelling to quantify the impact of transport interventions. Existing Emissions Model (EM) methodologies range from aggregate to disaggregate approaches, with more detail normally entailing more resources. However, it is not clear which approaches LGAs actually utilise. This article reports results of a survey designed to discover the level of detail considered practical by British LGAs (n = 34). Results show that resource scarcity is important, with particular importance attached to EM reusability and convenient input data sources. Most LGA EMs use traffic variable inputs (predominantly traffic flow and traffic average speed), with this approach being the best-fit for LGA resources. Link-by-link sources of data rated highly for convenience are road traffic models and urban traffic control systems

In-Use Emissions Testing of Diesel-Driven Buses in Southampton:Is Selective Catalytic Reduction as Effective as Fleet Operators Think?

Despite the continuously tightening emissions legislation, urban concentrations of nitrogen oxides (NOx) remain at harmful levels. Road transport is responsible for a large fraction, wherein diesel engines are the principal culprits. Turbocharged diesel engines have long been preferred in heavy duty applications, due to their torque delivery and low fuel consumption. Fleet operators are under pressure to understand and control the emissions of their vehicles, yet the performance of emissions abatement technology in real-world driving is largely unquantified. The most popular NOx abatement technology for heavy duty diesel vehicles is selective catalytic reduction. In this work, we empirically determine the efficiency of a factory-fitted SCR system in realworld driving by instrumenting passenger buses with both a portable emissions measurement system (PEMS) and a custom built telematics unit to record key parameters from the vehicle diagnostics systems. Wefindthateveninrelativelyfavourableconditions, while there is some improvement due to the use of SCR, the vehicles operate far from the design emissions targets. The archival value of this paper is in quantification of real world emissions versus design levels and the factors responsible for the discrepancy, as well as in examination of technologies to reduce this difference

Stabilized leapfrog based local time-stepping method for the wave equation

Local time-stepping methods permit to overcome the severe stability constraint on explicit methods caused by local mesh refinement without sacrificing explicitness. In \cite{DiazGrote09}, a leapfrog based explicit local time-stepping (LF-LTS) method was proposed for the time integration of second-order wave equations. Recently, optimal convergence rates were proved for a conforming FEM discretization, albeit under a CFL stability condition where the global time-step, $\Delta t$, depends on the smallest elements in the mesh \cite{grote_sauter_1}. In general one cannot improve upon that stability constraint, as the LF-LTS method may become unstable at certain discrete values of $\Delta t$. To remove those critical values of $\Delta t$, we apply a slight modification (as in recent work on LF-Chebyshev methods \cite{CarHocStu19}) to the original LF-LTS method which nonetheless preserves its desirable properties: it is fully explicit, second-order accurate, satisfies a three-term (leapfrog like) recurrence relation, and conserves the energy. The new stabilized LF-LTS method also yields optimal convergence rates for a standard conforming FE discretization, yet under a CFL condition where $\Delta t$ no longer depends on the mesh size inside the locally refined region

Uncertainty Quantification by MLMC and Local Time-stepping For Wave Propagation

Because of their robustness, efficiency and non-intrusiveness, Monte Carlo methods are probably the most popular approach in uncertainty quantification to computing expected values of quantities of interest (QoIs). Multilevel Monte Carlo (MLMC) methods significantly reduce the computational cost by distributing the sampling across a hierarchy of discretizations and allocating most samples to the coarser grids. For time dependent problems, spatial coarsening typically entails an increased time-step. Geometric constraints, however, may impede uniform coarsening thereby forcing some elements to remain small across all levels. If explicit time-stepping is used, the time-step will then be dictated by the smallest element on each level for numerical stability. Hence, the increasingly stringent CFL condition on the time-step on coarser levels significantly reduces the advantages of the multilevel approach. By adapting the time-step to the locally refined elements on each level, local time-stepping (LTS) methods permit to restore the efficiency of MLMC methods even in the presence of complex geometry without sacrificing the explicitness and inherent parallelism

Uncertainty quantification by multilevel Monte Carlo and local time-stepping

Because of their robustness, efficiency, and non intrusiveness, Monte Carlo methods are probablythe most popular approach in uncertainty quantification for computing expected values of quantitiesof interest. Multilevel Monte Carlo (MLMC) methods significantly reduce the computational costby distributing the sampling across a hierarchy of discretizations and allocating most samples tothe coarser grids. For time dependent problems, spatial coarsening typically entails an increasedtime step. Geometric constraints, however, may impede uniform coarsening thereby forcing someelements to remain small across all levels. If explicit time-stepping is used, the time step will thenbe dictated by the smallest element on each level for numerical stability. Hence, the increasinglystringent CFL condition on the time step on coarser levels significantly reduces the advantages of themultilevel approach. To overcome that bottleneck we propose to combine the multilevel approach ofMLMC with local time-stepping. By adapting the time step to the locally refined elements on eachlevel, the efficiency of MLMC methods is restored even in the presence of complex geometry withoutsacrificing the explicitness and inherent parallelism. In a careful cost comparison, we quantify thereduction in computational cost for local refinement either inside a small fixed region or towards areentrant corner

Formation of Jets and Equatorial Superrotation on Jupiter

The zonal flow in Jupiter's upper troposphere is organized into alternating retrograde and prograde jets, with a prograde (superrotating) jet at the equator. Existing models posit as the driver of the flow either differential radiative heating of the atmosphere or intrinsic heat fluxes emanating from the deep interior; however, they do not reproduce all large-scale features of Jupiter's jets and thermal structure. Here it is shown that the difficulties in accounting for Jupiter's jets and thermal structure resolve if the effects of differential radiative heating and intrinsic heat fluxes are considered together, and if upper-tropospheric dynamics are linked to a magnetohydrodynamic (MHD) drag that acts deep in the atmosphere. Baroclinic eddies generated by differential radiative heating can account for the off-equatorial jets; meridionally propagating equatorial Rossby waves generated by intrinsic convective heat fluxes can account for the equatorial superrotation. The zonal flow extends deeply into the atmosphere, with its speed changing with depth, up to depths at which the MHD drag acts. The theory is supported by simulations with an energetically consistent general circulation model of Jupiter's outer atmosphere. A simulation that incorporates differential radiative heating and intrinsic heat fluxes reproduces Jupiter's observed jets and thermal structure and makes testable predictions about as-yet unobserved aspects thereof. A control simulation that incorporates only differential radiative heating but not intrinsic heat fluxes produces off-equatorial jets but no equatorial superrotation; another control simulation that incorporates only intrinsic heat fluxes but not differential radiative heating produces equatorial superrotation but no off-equatorial jets.Comment: 23 pages, 10 figure