skip to main content

High-order finite element methods for a pressure Poisson equation reformulation of the Navier-Stokes equations with electric boundary conditions.
Pressure Poisson equation (PPE) reformulations of the incompressible Navier–Stokes equations (NSE) replace the incompressibility constraint by a Poisson equation for the pressure and a suitable choice of boundary conditions. This yields a time-evolution equation for the velocity field only, with the pressure gradient acting as a nonlocal operator. Thus, numerical methods based on PPE reformulations are representatives of a class of methods that have no principal limitations in achieving high order. In this paper, it is studied to what extent high-order methods for the NSE can be obtained from a specific PPE reformulation with electric boundary conditions (EBC). To that end, implicit–explicit (IMEX) time-stepping is used to decouple the pressure solve from the velocity update, while avoiding a parabolic time-step restriction; and mixed finite elements are used in space, to capture the structure imposed by the EBC. Via numerical examples, it is demonstrated that the methodology can yield at least third order accuracy in space and time.
Authors:
; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10205749
Journal Name:
Computer methods in applied mechanics and engineering
Volume:
373
Page Range or eLocation-ID:
113451
ISSN:
1879-2138
Sponsoring Org:
National Science Foundation
##### More Like this
1. We introduce a mathematical modeling framework for the conformational dynamics of charged molecules (i.e., solutes) in an aqueous solvent (i.e., water or salted water). The solvent is treated as an incompressible fluid, and its fluctuating motion is described by the Stokes equation with the Landau–Lifschitz stochastic stress. The motion of the solute-solvent interface (i.e., the dielectric boundary) is determined by the fluid velocity together with the balance of the viscous force,hydrostatic pressure, surface tension, solute-solvent van der Waals interaction force, and electrostatic force. The electrostatic interactions are described by the dielectric Poisson–Boltzmann theory.Within such a framework, we derive a generalizedmore »
2. Abstract We consider a collection of Markov chains that model the evolution of multitype biological populations. The state space of the chains is the positive orthant, and the boundary of the orthant is the absorbing state for the Markov chain and represents the extinction states of different population types. We are interested in the long-term behavior of the Markov chain away from extinction, under a small noise scaling. Under this scaling, the trajectory of the Markov process over any compact interval converges in distribution to the solution of an ordinary differential equation (ODE) evolving in the positive orthant. We studymore »
3. Direct numerical simulations of turbulent boundary layers with a nominal free-stream Mach number of $6$ and a Reynolds number of $Re_{\unicode[STIX]{x1D70F}}\approx 450$ are conducted at a wall-to-recovery temperature ratio of $T_{w}/T_{r}=0.25$ and compared with a previous database for $T_{w}/T_{r}=0.76$ in order to investigate pressure fluctuations and their dependence on wall temperature. The wall-temperature dependence of widely used velocity and temperature scaling laws for high-speed turbulent boundary layers is consistent with previous studies. The near-wall pressure-fluctuation intensities are dramatically modified by wall-temperature conditions. At different wall temperatures, the variation of pressure-fluctuation intensities as a function of wall-normal distance is dramatically modifiedmore »
4. Abstract

Many structures are subjected to varying forces, moving boundaries, and other dynamic conditions. Whether part of a vehicle, building, or active energy mitigation device, data on such changes can represent useful knowledge, but also presents challenges in its collection and analysis. In systems where changes occur rapidly, assessment of the system’s state within a useful time span is required to enable an appropriate response before the system’s state changes further. Rapid state estimation is especially important but poses unique difficulties.

In determining the state of a structural system subjected to high-rate dynamic changes, measuring the frequency response is onemore »

In order to facilitate development, the Air Force Research Laboratory created the DROPBEAR, a testbed with an oscillating beam subjected to a continuously variable boundary condition. One end of the beam is held by a stationary fixed support, while a pinned support is able to move along the beam’s length. The free end of the beam structure is instrumented with acceleration, velocity, and position sensors measuring the beam’s vertical axis. Direct position measurement of the pin location is also taken to provide a reference for comparison with numerical models.

This work presents a numerical investigation into methods for extracting the frequency response of a structure in real-time. An FFT based method with a rolling window is used to track the frequency of a data set generated to represent the range of the DROPBEAR, and is run with multiple window lengths. The frequency precision and latency of the FFT method is analyzed in each configuration. A specialized frequency extraction technique, Delayed Comparison Error Minimization, is implemented with parameters optimized for the frequency range of interest. The performance metrics of latency and precision are analyzed and compared to the baseline rolling FFT method results, and applicability is discussed.

5. Abstract The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori , but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surfacemore »