Verification of BOUT++ by the method of manufactured solutions

BOUT++ is a software package designed for solving plasma fluid models. It has been used to simulate a wide range of plasma phenomena ranging from linear stability analysis to 3D plasma turbulence and is capable of simulating a wide range of drift-reduced plasma fluid and gyro-fluid models. A verification exercise has been performed as part of a EUROfusion Enabling Research project, to rigorously test the correctness of the algorithms implemented in BOUT++, by testing order-of-accuracy convergence rates using the Method of Manufactured Solutions (MMS). We present tests of individual components including time-integration and advection schemes, non-orthogonal toroidal field-aligned coordinate systems and the shifted metric procedure which is used to handle highly sheared grids. The flux coordinate independent approach to differencing along magnetic field-lines has been implemented in BOUT++ and is here verified using the MMS in a sheared slab configuration. Finally, we show tests of three complete models: 2-field Hasegawa-Wakatani in 2D slab, 3-field reduced magnetohydrodynamics (MHD) in 3D field-aligned toroidal coordinates, and 5-field reduced MHD in slab geometry

Preliminary Evidence for cue-induced Alcohol Craving Modulated by Serotonin Transporter Gene Polymorphism rs1042173

We previously have shown that cue-induced alcohol craving and propensity for higher drinking are modulated by allelic differences in SLC6A4 associated with serotonin transporter (5-HTT) expression level alterations. In an independent study, we characterized another polymorphism, SNP rs1042173, in 3′-untranslated region (3′-UTR) of the same gene, which also altered 5-HTT expression levels; the T allele of rs1042173 was associated with lower mRNA and protein levels. In subsequent analyses, the TT genotype was found to be associated with higher drinking intensity in alcohol-dependent (AD) individuals of Caucasian descent. Building upon these findings, we hypothesized that the low-expressing TT genotype associated with intense drinking would predict higher craving for alcohol in AD individuals. In this pilot study, we sought to test our hypothesis by examining 34 Hispanic AD volunteers (mean age, 34.8 years) for rs1042173 genotype-based [i.e., TT versus TG/GG (Gx)] differences in subjective response to alcohol. We employed a human laboratory paradigm and analyzed the data using a linear mixed-effects model (SAS® PROC MIXED) to assess treatment, cue procedures, and genotype main effects as well as the two-way interaction effects between them. On subjective “urge to drink” and “crave for a drink,” we found a significant main effect of the cue experiment (p ≤ 0.01) and an interaction effect between genotype and cue effects (p < 0.05). TT genotype was associated with higher urge to drink (p = 0.002) and crave for a drink (p = 0.005) when exposed to alcohol cue. Our results not only support the hypothesis that rs1042173 is a genetic marker for cue-induced alcohol craving among AD males but also are suggestive of a neurobiological mechanism associated with the rs1042173-TT genotype that triggers a disproportionate craving in response to alcohol consumption, which in turn may lead to more intense drinking. Future studies with larger sample sizes are needed to characterize the interactive effects of the serotonin transporter-linked polymorphic region (5′-HTTLPR)-L-allele reported in our previous study and of the rs1042173-TT genotype on cue-induced alcohol craving

The synergistic integration of computational fluid dynamics and experimental fluid dynamics for ground effect aerodynamics studies

This article highlights the ‘synergistic’ use of experimental fluid dynamics (EFD) and computational fluid dynamics (CFD), where the two sets of simulations are performed concurrently and by the same researcher. In particular, examples from the area of ground effect aerodynamics are discussed, where the major facility used was also designed through a combination of CFD and EFD. Three examples are than outlined, to demonstrate the insight that can be obtained from the integration of CFD and EFD studies. The case studies are the study of dimple flow (to enhance aerodynamic performance), the analysis of a Formula-style front wing and wheel, and the study of compressible flow ground effect aerodynamics. In many instances, CFD has been used to not only provide complementary information to an experimental study, but to design the experiments. Laser-based, non-intrusive experimental techniques were used to provide an excellent complement to CFD. The large datasets found from both experimental and numerical simulations have required a new methodology to correlate the information; a new post-processing method has been developed, making use of the kriging and co-kriging estimators, to develop correlations between the often disparate data types

Sensitivity analysis of unsteady RANS flows

In this paper, we presents a continuous sensitivity equation method (SEM) for Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations with an adaptative Finite Element Method. The development is performed for value parameters that do not affect the geometry of the computational domain. We use the standard \u3ba - \u3b5 model of turbulence with wall functions and some of its usual variants such as Kato-Launder modification and RNG. The logarithmic form of transport equations for \u3ba - \u3b5 are used. An adaptative remeshing algorithm is developed for time integration of the URANS equations. Formulation and implementation are first verified on a problem with a closed form solution : the decay grid turbulence. This is followed by applications to transient flow reaching a steady state for a backward facing step and vortex shedding behind a cylinder. Flow and sensitivity analyses are performed. Serveral uses of sensitivity information are demonstrated : identification of key parameters controlling the flow, assessing the influence of closure coefficients, and nearby solutions.Dans le pr\ue9sent article, on pr\ue9sente une m\ue9thode avec \ue9quation de sensibilit\ue9 continue pour des simulations Reynolds-Averaged Navier-Stokes instables (URANS) avec une m\ue9thode adaptative \ue0 \ue9l\ue9ments finis. Le d\ue9veloppement est r\ue9alis\ue9 pour des valeurs de param\ue8tres qui n\u2019affectent pas la g\ue9om\ue9trie du domaine de calcul. On utilise le mod\ue8le standard k 12e\u2013epsilon de turbulence avec fonctions de paroi et certaines de ses variantes habituelles comme une modification Kato-Launder ou RNG. On utilise la forme logarithmique des \ue9quations de transport pour k et e -epsilon. Un algorithme adaptatif de remaillage a \ue9t\ue9 d\ue9velopp\ue9 pour l\u2019int\ue9gration en temps des \ue9quations URANS. La formulation et la mise en \u153uvre sont d\u2019abord v\ue9rifi\ue9es au moyen d\u2019un probl\ue8me ayant une solution \ue0 forme ferm\ue9e : la turbulence de la grille de d\ue9croissance. Ceci est suivi d\u2019applications \ue0 des \ue9coulements transitoires atteignant un \ue9tat d\u2019\ue9quilibre pour une \ue9tape orient\ue9e vers l\u2019arri\ue8re et un d\ue9collement de tourbillon derri\ue8re un cylindre. On a r\ue9alis\ue9 des analyses d\u2019\ue9coulement et de sensibilit\ue9. On a \ue9tabli plusieurs utilisations des renseignements sur la sensibilit\ue9 : identification de param\ue8tres cl\ue9s contr\uf4lant l\u2019\ue9coulement, \ue9valuation de l\u2019influence des coefficients de fermeture et solutions de proximit\ue9.Peer reviewed: YesNRC publication: Ye

Processes and Procedures for Application of CFD to Nuclear Reactor Safety Analysis

Traditionally, nuclear reactor safety analysis has been performed using systems analysis codes such as RELAP5, which was developed at the INL. However, goals established by the Generation IV program, especially the desire to increase efficiency, has lead to an increase in operating temperatures for the reactors. This increase pushes reactor materials to operate towards their upper temperature limits relative to structural integrity. Because there will be some finite variation of the power density in the reactor core, there will be a potential for local hot spots to occur in the reactor vessel. Hence, it has become apparent that detailed analysis will be required to ensure that local ‘hot spots’ do not exceed safety limits. It is generally accepted that computational fluid dynamics (CFD) codes are intrinsically capable of simulating fluid dynamics and heat transport locally because they are based on ‘first principles.’ Indeed, CFD analysis has reached a fairly mature level of development, including the commercial level. However, CFD experts are aware that even though commercial codes are capable of simulating local fluid and thermal physics, great care must be taken in their application to avoid errors caused by such things as inappropriate grid meshing, low-order discretization schemes, lack of iterative convergence and inaccurate time-stepping. Just as important is the choice of a turbulence model for turbulent flow simulation. Turbulence models model the effects of turbulent transport of mass, momentum and energy, but are not necessarily applicable for wide ranges of flow types. Therefore, there is a well-recognized need to establish practices and procedures for the proper application of CFD to simulate flow physics accurately and establish the level of uncertainty of such computations. The present document represents contributions of CFD experts on what the basic practices, procedures and guidelines should be to aid CFD analysts to obtain accurate estimates of the flow and energy transport as applied to nuclear reactor safety. However, it is expected that these practices and procedures will require updating from time to time as research and development affect them or replace them with better procedures. The practices and procedures are categorized into five groups. These are: 1.Code Verification 2.Code and Calculation Documentation 3.Reduction of Numerical Error 4.Quantification of Numerical Uncertainty (Calculation Verification) 5.Calculation Validation. These five categories have been identified from procedures currently required of CFD simulations such as those required for publication of a paper in the ASME Journal of Fluids Engineering and from the literature such as Roache [1998]. Code verification refers to the demonstration that the equations of fluid and energy transport have been correctly coded in the CFD code. Code and calculation documentation simply means that the equations and their discretizations, etc., and boundary and initial conditions used to pose the fluid flow problem are fully described in available documentation. Reduction of numerical error refers to practices and procedures to lower numerical errors to negligible or very low levels as is reasonably possible (such as avoiding use of first-order discretizations). The quantification of numerical uncertainty is also known as calculation verification. This means that estimates are made of numerical error to allow the characterization of the numerica

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

We present a study of $\gamma$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $\gamma$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $\gamma$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{\mu G}$ is deduced.Comment: 33 pages, 9 Figures, 6 Table

Compressive properties of min-mod-type limiters in modelling shockwave-containing flows

The long-ignored compressive properties of Min-mod-type limiter is investigated in this manuscript by demonstrating its potential in numerically modelling shockwave-containing flows, especially in shock wave/boundary layer interaction (SWBLI) problems. Theoretical studies were firstly performed based on Sweby’s total variation diminishing (TVD) limiter region and Spekreijse’s monotonicity-preserving limiter region to indicate Min-mod-type limiters’ compressive properties. The influence of limiters on the solution accuracy was evaluated using a hybrid-order analysis method based on the grid-independent study in three typical shockwave-containing flows. The conclusions are that, Min-mod-type limiter can be utilized as a dissipative and/or compressive limiter, but depending on the reasonable value of the compression parameter. The compressive Min-mod limiter tends to be more attractive in modelling shockwave-containing flows as compared to other commonly preferred limiters because of its stable computational process and its high-resolution predictions. However, the compressive Min-mod limiter may suffer from its slightly poor convergence, as that observed in other commonly accepted smooth limiters in modelling SWBLI problems. © 2020, The Author(s)

Measurement of Cosmic-ray Electrons at TeV Energies by VERITAS

Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering processes while propagating within the Galaxy and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, USA, is primarily utilized for gamma-ray astronomy, but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710 $\pm$ 40$_{stat}$ $\pm$ 140$_{syst}$ GeV.Comment: 17 pages, 2 figures, accepted for publication in PR