Stochastic turbulence modeling in RANS simulations via Multilevel Monte Carlo

A multilevel Monte Carlo (MLMC) method for quantifying model-form uncertainties associated with the Reynolds-Averaged Navier-Stokes (RANS) simulations is presented. Two, high-dimensional, stochastic extensions of the RANS equations are considered to demonstrate the applicability of the MLMC method. The first approach is based on global perturbation of the baseline eddy viscosity field using a lognormal random field. A more general second extension is considered based on the work of [Xiao et al.(2017)], where the entire Reynolds Stress Tensor (RST) is perturbed while maintaining realizability. For two fundamental flows, we show that the MLMC method based on a hierarchy of meshes is asymptotically faster than plain Monte Carlo. Additionally, we demonstrate that for some flows an optimal multilevel estimator can be obtained for which the cost scales with the same order as a single CFD solve on the finest grid level.Comment: 40 page

Coherent Pion Production in Neutrino Nucleus Scattering

In this article, we study the coherent pion production in neutrino-nucleus interaction in the resonance region using the formalism based on partially conserved axial current (PCAC) theorem which relates the neutrino-nucleus cross section to the pion-nucleus elastic cross section. The pion nucleus elastic cross section is calculated using the Glauber model in terms of pion-nucleon cross sections obtained by parameterizing the experimental data. We calculate the differential and integrated cross sections for charged current coherent pion production in neutrino carbon scattering. The results of integrated cross section calculations are compared with the measured data. Predictions for the differential and integrated cross sections for coherent pion productions in neutrino iron scattering using above formalism are also made