BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

John A. Evans

John A. Evans contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.flu-dynComputational Geometryphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
5topics
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

9 published item(s)

preprint2022arXiv

A Divergence-Conforming Hybridized Discontinuous Galerkin Method for the Incompressible Magnetohydrodynamics Equations

We introduce a new hybridized discontinuous Galerkin method for the incompressible magnetohydrodynamics equations. If particular velocity, pressure, magnetic field, and magnetic pressure spaces are employed for both element and trace solution fields, we arrive at an energy stable method which returns pointwise divergence-free velocity fields and magnetic fields and properly balances linear momentum. We discretize in time using a second-order-in-time generalized-$α$ method, and we present a block iterative method for solving the resulting nonlinear system of equations at each time step. We numerically examine the effectiveness of our method using a manufactured solution and observe our method yields optimal convergence rates in the $L_2$ norm for the velocity field, pressure field, magnetic field, and magnetic pressure field. We further find our method is pressure robust. We then apply our method to a selection of benchmark problems and numerically confirm our method is energy stable.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A note on the conservation properties of the generalized-$α$ method

We show that the second-order accurate generalized-$α$ method on a uniform temporal mesh may be viewed as an implicit midpoint method on a shifted temporal mesh. With this insight, we demonstrate generalized-$α$ time integration of a finite element spatial discretization of a conservation law system results in a fully-discrete method admitting discrete balance laws when (i) the time integration is second-order accurate, (ii) a uniform temporal mesh is employed, (iii) the spatial discretization is conservative, and (iv) conservation variables are discretized.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constructing Nitsche's method for variational problems

Nitsche's method is a well-established approach for weak enforcement of boundary conditions for partial differential equations (PDEs). It has many desirable properties, including the preservation of variational consistency and the fact that it yields symmetric, positive-definite discrete linear systems that are not overly ill-conditioned. In recent years, the method has gained in popularity in a number of areas, including isogeometric analysis, immersed methods, and contact mechanics. However, arriving at a formulation based on Nitsche's method can be a mathematically arduous process, especially for high-order PDEs. Fortunately, the derivation is conceptually straightforward in the context of variational problems. To facilitate the process, we devised an abstract framework for constructing Nitsche's method for these types of problems in [J. Benzaken, J. A. Evans, S. McCormick, and R. Tamstorf, Nitsche's method for linear Kirchhoff-Love shells: Formulation, error analysis, and verification, Comput. Methods Appl. Mech. Eng., 374 (2021), p. 113544]. The goal of this paper is to elucidate the process through a sequence of didactic examples. First, we show the derivation of Nitsche's method for Poisson's equation to gain an intuition for the various steps. Next, we present the abstract framework and then revisit the derivation for Poisson's equation to use the framework and add mathematical rigor. In the process, we extend our derivation to cover the vector-valued setting. Armed with a basic recipe, we then show how to handle a higher-order problem by considering the vector-valued biharmonic equation and the linearized Kirchhoff-Love plate. In the end, the hope is that the reader will be able to apply Nitsche's method to any problem that arises from variational principles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-accuracy mesh-free quadrature for trimmed parametric surfaces and volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. The scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Invariant Data-Driven Subgrid Stress Modeling in the Strain-Rate Eigenframe for Large Eddy Simulation

We present a new approach for constructing data-driven subgrid stress models for large eddy simulation of turbulent flows. The key to our approach is representation of model input and output tensors in the filtered strain rate eigenframe. Provided inputs and outputs are selected and non-dimensionalized in a suitable manner, this yields a model form that is symmetric, Galilean invariant, rotationally invariant, reflectionally invariant, and unit invariant. We use this model form to train a simple and efficient neural network model using only one time step of filtered direct numerical simulation data from a forced homogeneous isotropic turbulence simulation. We demonstrate the accuracy of this model as well as the model's ability to generalize to previously unseen filter widths, Reynolds numbers, and flow physics using a priori and a posteriori tests.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Skeleton-stabilized divergence-conforming B-spline discretizations for highly advective incompressible flow problems

We consider a stabilization method for divergence-conforming B-spline discretizations of the incompressible Navier--Stokes problem wherein jumps in high-order normal derivatives of the velocity field are penalized across interior mesh facets. We prove that this method is pressure robust, consistent, and energy stable, and we show how to select the stabilization parameter appearing in the method so that excessive numerical dissipation is avoided in both the cross-wind direction and in the diffusion-dominated regime. We examine the efficacy of the method using a suite of numerical experiments, and we find the method yields optimal $\textbf{L}^2$ and $\textbf{H}^1$ convergence rates for the velocity field, eliminates spurious small-scale structures that pollute Galerkin approximations, and is effective as an Implicit Large Eddy Simulation (ILES) methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Direct Numerical Simulation of a Turbulent Boundary Layer on a Flat Plate Using Synthetic Turbulence Generation

The turbulent boundary layer over a flat plate is computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations as a test bed for a synthetic turbulence generator (STG) inflow boundary condition. The inlet momentum thickness Reynolds number is approximately 1,000. The study provides validation of the ability of the STG to develop accurate turbulence in 5 to 7 boundary layer thicknesses downstream of the boundary condition. Also tested was the effect of changes in the stabilization scheme on the development of the boundary layer. Moreover, the grid resolution required for both the development region and the downstream flow is investigated when using a stabilized finite element method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

S-Frame Discrepancy Correction Models for Data-Informed Reynolds Stress Closure

Despite their well-known limitations, RANS models remain the most commonly employed tool for modeling turbulent flows in engineering practice. RANS models are predicated on the solution of the RANS equations, but these equations involve an unclosed term, the Reynolds stress tensor, which must be modeled. The Reynolds stress tensor is often modeled as an algebraic function of mean flow field variables and turbulence variables. This introduces a discrepancy between the Reynolds stress tensor predicted by the model and the exact Reynolds stress tensor. This discrepancy can result in inaccurate mean flow field predictions. In this paper, we introduce a data-informed approach for arriving at Reynolds stress models with improved predictive performance. Our approach relies on learning the components of the Reynolds stress discrepancy tensor associated with a given Reynolds stress model in the mean strain-rate tensor eigenframe. These components are typically smooth and hence simple to learn using state-of-the-art machine learning strategies and regression techniques. Our approach automatically yields Reynolds stress models that are symmetric, and it yields Reynolds stress models that are both Galilean and frame invariant provided the inputs are themselves Galilean and frame invariant. To arrive at computable models of the discrepancy tensor, we employ feed-forward neural networks and an input space spanning the integrity basis of the mean strain-rate tensor, the mean rotation-rate tensor, the mean pressure gradient, and the turbulent kinetic energy gradient, and we introduce a framework for dimensional reduction of the input space to further reduce computational cost. Numerical results illustrate the effectiveness of the proposed approach for data-informed Reynolds stress closure for a suite of turbulent flow problems of increasing complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spectral mesh-free quadrature for planar regions bounded by rational parametric curves

This work presents spectral, mesh-free, Green's theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region's boundary curves corresponding to Green's theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green's theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree k with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time.

0 citations0 reviewsNo code linkedOpen access