Researcher profile

Oana Marin

Oana Marin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

6 published item(s)

preprint2022arXiv

Fitting Matérn Smoothness Parameters Using Automatic Differentiation

The Matérn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $ν$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Matérn covariance function with respect to $ν$ require derivatives of the modified second-kind Bessel function $\mathcal{K}_ν$ with respect to $ν$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $ν$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $ν$ use finite difference estimates for $\partial_ν\mathcal{K}_ν$. In this work, we introduce a new implementation of $\mathcal{K}_ν$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. We provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

preprint2022arXiv

On automatic differentiation for the Matérn covariance

To target challenges in differentiable optimization we analyze and propose strategies for derivatives of the Matérn kernel with respect to the smoothness parameter. This problem is of high interest in Gaussian processes modelling due to the lack of robust derivatives of the modified Bessel function of second kind with respect to order. In the current work we focus on newly identified series expansions for the modified Bessel function of second kind valid for complex orders. Using these expansions we obtain highly accurate results using the complex step method. Furthermore, we show that the evaluations using the recommended expansions are also more efficient than finite differences.

preprint2022arXiv

The PETSc Community Is the Infrastructure

The communities who develop and support open source scientific software packages are crucial to the utility and success of such packages. Moreover, these communities form an important part of the human infrastructure that enables scientific progress. This paper discusses aspects of the PETSc (Portable Extensible Toolkit for Scientific Computation) community, its organization, and technical approaches that enable community members to help each other efficiently.

preprint2020arXiv

A scalable matrix-free spectral element approach for unsteady PDE constrained optimization using PETSc/TAO

We provide a new approach for the efficient matrix-free application of the transpose of the Jacobian for the spectral element method for the adjoint based solution of partial differential equation (PDE) constrained optimization. This results in optimizations of nonlinear PDEs using explicit integrators where the integration of the adjoint problem is not more expensive than the forward simulation. Solving PDE constrained optimization problems entails combining expertise from multiple areas, including simulation, computation of derivatives, and optimization. The Portable, Extensible Toolkit for Scientific computation (PETSc) together with its companion package, the Toolkit for Advanced Optimization (TAO), is an integrated numerical software library that contains an algorithmic/software stack for solving linear systems, nonlinear systems, ordinary differential equations, differential algebraic equations, and large-scale optimization problems and, as such, is an ideal tool for performing PDE-constrained optimization. This paper describes an efficient approach in which the software stack provided by PETSc/TAO can be used for large-scale nonlinear time-dependent problems. While time integration can involve a range of high-order methods, both implicit and explicit. The PDE-constrained optimization algorithm used is gradient-based and seamlessly integrated with the simulation of the physical problem.

preprint2020arXiv

Direct Numerical Simulation of Fluid Flow in a 5x5 Square Rod Bundle Using Nek5000

Rod bundle flows are commonplace in nuclear engineering, and are present in light water reactors (LWRs) as well as other more advanced concepts. Inhomogeneities in the bundle cross section can lead to complex flow phenomena, including varying local conditions of turbulence. Despite the decades of numerical and experimental investigations regarding this topic, and the importance of elucidating the physics of the flow field, to date there are few publicly available direct numerical simulations (DNS) of the flow in multiple-pin rod bundles. Thus a multiple-pin DNS study can provide significant value toward reaching a deeper understanding of the flow physics, as well as a reference simulation for development of various reduced-resolution analysis techniques. To this end, DNS of the flow in a square 5x5 rod bundle at Reynolds number of 19,000 has been performed using the highly-parallel spectral element code Nek5000. The geometrical dimensions were representative of typical LWR fuel designs. The DNS was designed using microscales estimated from an advanced Reynolds-Averaged Navier-Stokes (RANS) model. Characteristics of the velocity field, Reynolds stresses, and anisotropy are presented in detail for various regions of the bundle. The turbulent kinetic energy budget is also presented and discussed

preprint2020arXiv

Towards Direct Numerical Simulation of a 5x5 Rod Bundle

Rod bundle flows are prevalent in nuclear engineering for both light water reactors (LWR) and advanced reactor concepts. Unlike canonical channel flow, the flow in rod bundles presents some unique characteristics, notably due to the inhomogeneous cross section which can present different local conditions of turbulence as well as localized effects characteristic of external flows. Despite the ubiquity of rod bundle flows and the decades of experimental and numerical knowledge acquired in this field, there are no publicly available direct numerical simulations (DNS) of the flow in multiple pin rod bundles with heat transfer. A multiple pin DNS study is of great value as it would allow for assessment of the reliability of various turbulence models in the presence of heat transfer, as well as allow for a deeper understanding of the flow physics. We present work towards DNS of the flow in a square 5x5 rod bundle representative of LWR fuel. We consider standard configurations as well as configurations where the central pin is replaced with a guide thimble. We perform simulations in STAR-CCM+ to design the numerical DNS, which is to be conducted using the open source spectral element code Nek5000. Large Eddy Simulations are also performed in Nek5000 to confirm that the resolution requirements are adequate. We compare results from STAR-CCM+ and Nek5000, which show very good agreement in the wide gaps with larger discrepancies in the narrow gaps. In particular, evidence of a gap vortex street is seen in the edge subchannels in LES but is not predicted by STAR-CCM+.