Source author record

Robert D. Moser

Robert D. Moser appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

A pseudospectral method for direct numerical simulation of low-Mach, variable-density, turbulent flows

A novel algorithm for the direct numerical simulation of the variable-density, low-Mach Navier-Stokes equations extending the method of Kim, Moin, and Moser (1987) for incompressible flow is presented here. A Fourier representation is employed in the two homogeneous spatial directions and a number of discretizations can be used in the inhomogeneous direction. The momentum is decomposed into divergence- and curl-free portions which allows the momentum equations to be rewritten, removing the need to solve for the pressure. The temporal discretization is based on an explicit, segregated Runge-Kutta method and the scalar equations are reformulated to directly address the redundancy of the equation of state and the mass conservation equation. An efficient, matrix-free, iterative solution of the resulting equations allows for second-order accuracy in time and numerical stability for large density ratios, which is demonstrated for ratios up to $\sim 25.7$.

preprint2021arXiv

Active model split hybrid RANS/LES

Reliably predictive simulation of complex flows requires a level of model sophistication and robustness exceeding the capabilities of current Reynolds-averaged Navier-Stokes (RANS) models. The necessary capability can often be provided by well-resolved large eddy simulation (LES), but, for many flows of interest, such simulations are too computationally intensive to be performed routinely. In principle, hybrid RANS/LES (HRL) models capable of transitioning through arbitrary levels of modeled and resolved turbulence would ameliorate both RANS deficiencies and LES expense. However, these HRL approaches have led to a host of unique complications, in addition to those already present in RANS and LES. This work proposes a modeling approach aimed at overcoming such challenges. The approach presented here relies on splitting the turbulence model into three distinct components: two responsible for the standard subgrid model roles of either providing the unresolved stress or dissipation and a third which reduces the model length scale by creating resolved turbulence. This formulation renders blending functions unnecessary in HRL. Further, the split-model approach both reduces the physics-approximation burden on simple eddy-viscosity-based models and provides convenient flexibility in model selection. In regions where the resolution is adequate to support additional turbulence, fluctuations are generated at the smallest locally resolved scales of motion. This active forcing drives the system towards a balance between RANS and grid-resolved LES for any combination of resolution and flow while the split-model formulation prevents local disruption to the total stress. The model is demonstrated on fully-developed, incompressible channel flow [1] and the periodic hill [2], in which it is shown to produce accurate results and avoid common HRL shortcomings, such as model stress depletion.

preprint2021arXiv

Bayesian inference of an uncertain generalized diffusion operator

This paper defines a novel Bayesian inverse problem to infer an infinite-dimensional uncertain operator appearing in a differential equation, whose action on an observable state variable affects its dynamics. Inference is made tractable by parametrizing the operator using its eigendecomposition. The plausibility of operator inference in the sparse data regime is explored in terms of an uncertain, generalized diffusion operator appearing in an evolution equation for a contaminant's transport through a heterogeneous porous medium. Sparse data are augmented with prior information through the imposition of deterministic constraints on the eigendecomposition and the use of qualitative information about the system in the definition of the prior distribution. Limited observations of the state variable's evolution are used as data for inference, and the dependence on the solution of the inverse problem is studied as a function of the frequency of observations, as well as on whether or not the data is collected as a spatial or time series.

preprint2021arXiv

Near wall patch representation of wall bounded turbulence

Recent experimental and computational studies indicate that near wall turbulent flows can be characterized by universal small scale autonomous dynamics that are modulated by large scale structures. We formulate numerical simulations of near wall turbulence in a small domain localized to the boundary, whose size scales in viscous units. To mimic the environment in which the near wall turbulence evolves, our formulation accounts for the flux of mean momentum through the upper boundary of the domain. Comparisons of the model's two dimensional energy spectra and low order single-point statistics with the corresponding quantities computed from zero and mild favorable pressure gradient direct numerical simulations indicate it successfully captures the dynamics of the small scale near wall turbulence.

preprint2015arXiv

The Parallel C++ Statistical Library for Bayesian Inference: QUESO

The Parallel C++ Statistical Library for the Quantification of Uncertainty for Estimation, Simulation and Optimization, Queso, is a collection of statistical algorithms and programming constructs supporting research into the quantification of uncertainty of models and their predictions. Queso is primarily focused on solving statistical inverse problems using Bayes's theorem, which expresses a distribution of possible values for a set of uncertain parameters (the posterior distribution) in terms of the existing knowledge of the system (the prior) and noisy observations of a physical process, represented by a likelihood distribution. The posterior distribution is not often known analytically, and so requires computational methods. It is typical to compute probabilities and moments from the posterior distribution, but this is often a high-dimensional object and standard Reimann-type methods for quadrature become prohibitively expensive. The approach Queso takes in this regard is to rely on Markov chain Monte Carlo (MCMC) methods which are well suited to evaluating quantities such as probabilities and moments of high-dimensional probability distributions. Queso's intended use is as tool to assist and facilitate coupling uncertainty quantification to a specific application called a forward problem. While many libraries presently exist that solve Bayesian inference problems, Queso is a specialized piece of software primarily designed to solve such problems by utilizing parallel environments demanded by large-scale forward problems. Queso is written in C++, uses MPI, and utilizes libraries already available to the scientific community.

preprint2014arXiv

Validating Predictions of Unobserved Quantities

The ultimate purpose of most computational models is to make predictions, commonly in support of some decision-making process (e.g., for design or operation of some system). The quantities that need to be predicted (the quantities of interest or QoIs) are generally not experimentally observable before the prediction, since otherwise no prediction would be needed. Assessing the validity of such extrapolative predictions, which is critical to informed decision-making, is challenging. In classical approaches to validation, model outputs for observed quantities are compared to observations to determine if they are consistent. By itself, this consistency only ensures that the model can predict the observed quantities under the conditions of the observations. This limitation dramatically reduces the utility of the validation effort for decision making because it implies nothing about predictions of unobserved QoIs or for scenarios outside of the range of observations. However, there is no agreement in the scientific community today regarding best practices for validation of extrapolative predictions made using computational models. The purpose of this paper is to propose and explore a validation and predictive assessment process that supports extrapolative predictions for models with known sources of error. The process includes stochastic modeling, calibration, validation, and predictive assessment phases where representations of known sources of uncertainty and error are built, informed, and tested. The proposed methodology is applied to an illustrative extrapolation problem involving a misspecified nonlinear oscillator.

preprint2013arXiv

Estimating Uncertainties in Statistics Computed from DNS

Rigorous assessment of uncertainty is crucial to the utility of DNS results. Uncertainties in the computed statistics arise from two sources: finite statistical sampling and the discretization of the Navier-Stokes equations. Due to the presence of non-trivial sampling error, standard techniques for estimating discretization error (such as Richardson extrapolation) fail or are unreliable. This work provides a systematic and unified approach for estimating these errors. First, a sampling error estimator that accounts for correlation in the input data is developed. Then, this sampling error estimate is used as part of a Bayesian extension of Richardson extrapolation in order to characterize the discretization error. These methods are tested using the Lorenz equations and are shown to perform well. These techniques are then used to investigate the sampling and discretization errors in the DNS of a wall-bounded turbulent flow. For both cases, it is found that while the sampling uncertainty is large enough to make the order of accuracy difficult to determine, the estimated discretization errors are quite small. This indicates that the commonly used heuristics provide ad- equate resolution for this class of problems. However, it is also found that, for some quantities, the discretization error is not small relative to sampling error, indicating that the conventional wisdom that sampling error dominates discretization error for this class of simulations needs to be reevaluated.