skip to main content

Title: 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.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Computer methods in applied mechanics and engineering
Page Range or eLocation-ID:
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 »Rayleigh–Plesset equation, a nonlinear stochastic ordinary differential equation (SODE), for the radius of a spherical charged molecule, such as anion. The spherical average of the stochastic stress leads to a multiplicative noise. We design and test numerical methods for solving the SODE and use the equation, together with explicit solvent molecular dynamics simulations, to study the effective radius of a single ion. Potentially, our general modeling framework can be used to efficiently determine the solute-solvent interfacial structures and predict the free energies of more complex molecular systems.« less
  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 »the asymptotic behavior of the quasi-stationary distributions (QSD) in this scaling regime. Our main result shows that, under conditions, the limit points of the QSD are supported on the union of interior attractors of the flow determined by the ODE. We also give lower bounds on expected extinction times which scale exponentially with the system size. Results of this type when the deterministic dynamical system obtained under the scaling limit is given by a discrete-time evolution equation and the dynamics are essentially in a compact space (namely, the one-step map is a bounded function) have been studied by Faure and Schreiber (2014). Our results extend these to a setting of an unbounded state space and continuous-time dynamics. The proofs rely on uniform large deviation results for small noise stochastic dynamical systems and methods from the theory of continuous-time dynamical systems. In general, QSD for Markov chains with absorbing states and unbounded state spaces may not exist. We study one basic family of binomial-Poisson models in the positive orthant where one can use Lyapunov function methods to establish existence of QSD and also to argue the tightness of the QSD of the scaled sequence of Markov chains. The results from the first part are then used to characterize the support of limit points of this sequence of QSD.« less
  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 »in the near-wall region but remains almost intact away from the wall. Wall cooling also has a strong effect on the frequency spectrum of wall-pressure fluctuations, resulting in a higher dominant frequency and a sharper spectrum peak with a faster roll-off at both the high- and low-frequency ends. The effect of wall cooling on the free-stream noise spectrum can be largely accounted for by the associated changes in boundary-layer velocity and length scales. The pressure structures within the boundary layer and in the free stream evolve less rapidly as the wall temperature decreases, resulting in an increase in the decorrelation length of coherent pressure structures for the colder-wall case. The pressure structures propagate with similar speeds for both wall temperatures. Due to wall cooling, the generated pressure disturbances undergo less refraction before they are radiated to the free stream, resulting in a slightly steeper radiation wave front in the free stream. Acoustic sources are largely concentrated in the near-wall region; wall cooling most significantly influences the nonlinear (slow) component of the acoustic source term by enhancing dilatational fluctuations in the viscous sublayer while damping vortical fluctuations in the buffer and log layers.« less
  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 »method that can be used to draw inferences, provided the system is adequately understood and defined. The work presented here is the result of an investigation into methods to determine the frequency response, and thus state, of a structure subjected to high-rate boundary changes in real-time.

    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.

    « less
  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 »fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. Our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.« less