Eigensolution analysis of spectral/hp continuous Galerkin approximations to advection-diffusion problems: insights into spectral vanishing viscosity

AbstractThis study addresses linear dispersion–diffusion analysis for the spectral/hp continuous Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion curves are obtained for the advection–diffusion problem and the role of multiple eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves' behaviour, we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion characteristics for under-resolved simulations of problems strongly dominated by either convection or diffusion. Subsequently, the linear advection equation augmented with spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of CG with SVV-based stabilization are verified to display similar non-smooth features in flow regions where convection is much stronger than dissipation or vice-versa, owing to a dependency of the standard SVV operator on a local Péclet number. First a modification is proposed to the traditional SVV scaling that enforces a globally constant Péclet number so as to avoid the previous issues. In addition, a new SVV kernel function is suggested and shown to provide a more regular behaviour for the eigencurves along with a consistent increase in resolution power for higher-order discretizations, as measured by the extent of the wavenumber range where numerical errors are negligible. The dissipation characteristics of CG with the SVV modifications suggested are then verified to be broadly equivalent to those obtained through upwinding in the discontinuous Galerkin (DG) scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is found to have a slightly higher resolution power for the same dissipation levels. These results show that improved CG-SVV characteristics can be pursued via different kernel functions with the aid of optimization algorithms

Tearing of free-standing graphene

We examine the fracture mechanics of tearing graphene. We present a molecular dynamics simulation of the propagation of cracks in clamped, free-standing graphene as a function of the out-of-plane force. The geometry is motivated by experimental configurations that expose graphene sheets to out-of-plane forces, such as back-gate voltage. We establish the geometry and basic energetics of failure and obtain approximate analytical expressions for critical crack lengths and forces. We also propose a method to obtain graphene's toughness. We observe that the cracks' path and the edge structure produced are dependent on the initial crack length. This work may help avoid the tearing of graphene sheets and aid the production of samples with specific edge structures.CAPESNational Science Foundation DMR 1002428Physic

Linear dispersion-diffusion analysis and its application to under-resolved turbulence simulations using discontinuous Galerkin spectral/hp methods

We investigate the potential of linear dispersion–diffusion analysis in providing direct guidelines for turbulence simulations through the under-resolved DNS (sometimes called implicit LES) approach via spectral/hp methods. The discontinuous Galerkin (DG) formulation is assessed in particular as a representative of these methods. We revisit the eigensolutions technique as applied to linear advection and suggest a new perspective to the role of multiple numerical modes, peculiar to spectral/hp methods. From this new perspective, “secondary” eigenmodes are seen to replicate the propagation behaviour of a “primary” mode, so that DG's propagation characteristics can be obtained directly from the dispersion–diffusion curves of the primary mode. Numerical dissipation is then appraised from these primary eigencurves and its effect over poorly-resolved scales is quantified. Within this scenario, a simple criterion is proposed to estimate DG's effective resolution in terms of the largest wavenumber it can accurately resolve in a given hp approximation space, also allowing us to present points per wavelength estimates typically used in spectral and finite difference methods. Although strictly valid for linear advection, the devised criterion is tested against (1D) Burgers turbulence and found to predict with good accuracy the beginning of the dissipation range on the energy spectra of under-resolved simulations. The analysis of these test cases through the proposed methodology clarifies why and how the DG formulation can be used for under-resolved turbulence simulations without explicit subgrid-scale modelling. In particular, when dealing with communication limited hardware which forces one to consider the performance for a fixed number of degrees of freedom, the use of higher polynomial orders along with moderately coarser meshes is shown to be the best way to translate available degrees of freedom into resolution power

Charging ahead on the transition to electric vehicles with standard 120 v wall outlets

Electrification of transportation is needed soon and at significant scale to meet climate goals, but electric vehicle adoption has been slow and there has been little systematic analysis to show that today's electric vehicles meet the needs of drivers. We apply detailed physics-based models of electric vehicles with data on how drivers use their cars on a daily basis. We show that the energy storage limits of today's electric vehicles are outweighed by their high efficiency and the fact that driving in the United States seldom exceeds 100 km of daily travel. When accounting for these factors, we show that the normal daily travel of 85-89% of drivers in the United States can be satisfied with electric vehicles charging with standard 120 V wall outlets at home only. Further, we show that 77-79% of drivers on their normal daily driving will have over 60 km of buffer range for unexpected trips. We quantify the sensitivities to terrain, high ancillary power draw, and battery degradation and show that an extreme case with all trips on a 3% uphill grade still shows the daily travel of 70% of drivers being satisfied with electric vehicles. These findings show that today's electric vehicles can satisfy the daily driving needs of a significant majority of drivers using only 120 V wall outlets that are already the standard across the United States

Service scheduling in garden maintenance

Neoturf is a Portuguese company working in the area of project, building and garden’s maintenance. Neoturf would like to have a procedure for scheduling and routing efficiently the clients from garden maintenance services. The company has two teams available during the whole year and an additional team during summer to handle all the maintenance jobs. Each team consists of two or three employees with a vehicle fully equipped with the tools that allow to carry out every kind of maintenance service. In the beginning of each year, the number and frequency of maintenance interventions to conduct during the year, on each client, are accorded. Each client is assigned to the same team and, usually, time windows are established so that visits to the client should occur only within these periods. As the Neoturf costumers’ are geographically spread over a wide region, the total distance on visiting clients is a factor that has a heavy weight on the costs of the company. Neoturf is concerned with reducing these costs, while satisfying the agreements with the clients

Four dimensional R^4 superinvariants through gauge completion

We fully compute the N=1 supersymmetrization of the fourth power of the Weyl tensor in d=4 x-space with the auxiliary fields. In a previous paper, we showed that their elimination requires an infinite number of terms; we explicitely compute those terms to order \kappa^4 (three loop). We also write, in superspace notation, all the possible N=1 actions, in four dimensions, that contain pure R^4 terms (with coupling constants). We explicitely write these actions in terms of the \theta components of the chiral density \epsilon and the supergravity superfields R, G_m, W_{ABC}. Using the method of gauge completion, we compute the necessary \theta components which allow us to write these actions in x-space. We discuss under which circumstances can these extra R^4 correction terms be reabsorbed in the pure supergravity action, and their relevance to the quantum supergravity/string theory effective actions.Comment: 20 pages, no figures. Sec. 3 clarified; typos correcte

Energy in an Expanding Universe in the Teleparallel Geometry

The main purpose of this paper is to explicitly verify the consistency of the energy-momentum and angular momentum tensor of the gravitational field established in the Hamiltonian structure of the Teleparallel Equivalent of General Relativity (TEGR). In order to reach these objectives, we obtained the total energy and angular momentum (matter plus gravitational field) of the closed universe of the Friedmann-Lemaitre-Robertson-Walker (FLRW). The result is compared with those obtained from the pseudotensors of Einstein and Landau-Lifshitz. We also applied the field equations (TEGR) in an expanding FLRW universe. Considering the stress energy-momentum tensor for a perfect fluid, we found a teleparallel equivalent of Friedmann equations of General Relativity (GR).Comment: 19 pages, no figures. Revised in view of Referee's comments. Version to appear in the Brazilian Journal of Physic

Effects of industrial heat treatments on the kinetics of inactivation of antimicrobial bovine milk xanthine oxidase.

Milk is a source of antimicrobial systems such as xanthine oxidoreductase, which has been proposed to modulate the oral and gut microbiota of infants. Heat treatments are applied to milk to ensure its microbial safety, however, the effects of heat on this antimicrobial enzyme are not known. The effects of batch pasteurization (BP), high-temperature short time (HTST), and ultra high temperature (UHT) on kinetics of inactivation of xanthine oxidase and its antimicrobial properties were determined. Xanthine oxidase activity was preserved by HTST (100%). Partial (8%) and nearly complete (95%) enzyme inactivation were observed for BP and UHT milks, respectively. K m values of 100 μM and V max values of 6.85, 5.12, 6.31, and 0.40 μmol/min/mg were determined for xanthine oxidase in raw, BP, HTST, and UHT milks, respectively. These results demonstrate that xanthine oxidase maintains apparent affinity and activity for its substrate when milk is treated by BP and HTST and yet the enzyme is inactivated with UHT. To investigate heat treatment-induced alterations in the biological activity of xanthine oxidase, heat treated milks were compared to raw milk for their ability to inhibit the growth of S. aureus. Raw, BP, and HTST milk xanthine oxidase efficiently inhibited S. aureus growth. However, these antibacterial properties were lost when milk was subjected to UHT. These results demonstrate that HTST and BP preserves bovine milk xanthine oxidase activity compared with UHT and that, the judicious selection of thermal treatments could be exploited to preserve the antimicrobial properties of bovine milk