skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Stable nodal projection method on octree grids
We propose a novel collocated projection method for solving the incompressible Navier-Stokes equations with arbitrary boundaries. Our approach employs non-graded octree grids, where all variables are stored at the nodes. To discretize the viscosity and projection steps, we utilize supra-convergent finite difference approximations with sharp boundary treatments. We demonstrate the stability of our projection on uniform grids, identify a sufficient stability condition on adaptive grids, and validate these findings numerically. We further demonstrate the accuracy and capabilities of our solver with several canonical two- and three-dimensional simulations of incompressible fluid flows. Overall, our method is second-order accurate, allows for dynamic grid adaptivity with arbitrary geometries, and reduces the overhead in code development through data collocation.  more » « less
Award ID(s):
1840265
PAR ID:
10505313
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Elsevier
Date Published:
Journal Name:
Journal of Computational Physics
Volume:
499
Issue:
C
ISSN:
0021-9991
Page Range / eLocation ID:
112695
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, SPIE)
    Zelnio, Edmund; Garber, Frederick D (Ed.)
    Ground Penetrating Radar (GPR) is essential for subsurface exploration. Conventional GPR 3D imaging demands dense spatial sampling along regular grids, which is both time-consuming and impractical in complex environments. In this work, we propose a novel method that combines sparse recovery techniques with a placement matrix to merge arbitrarily and sparsely sampled measurements into a regular grid framework. By exploiting the inherent sparsity of subsurface targets and using the Dantzig Selector with cross-validation, our method reconstructs the target reflectivity vector from random spatial sampling. The recovered data is then processed via the Back-Projection Algorithm (BPA) to generate high-resolution 3D images. Simulations demonstrate that our approach not only improves imaging quality under reduced sampling conditions but also efficiently handles arbitrary scanning paths by mapping irregular measurements onto the desired grid. 
    more » « less
  2. ; (, Journal of Computational Physics)
    Hypo-elastoplasticity is a framework suitable for modeling the mechanics of many hard materials that have small elastic deformation and large plastic deformation. In laboratory tests for these materials the Cauchy stress is often in quasi-static equilibrium. Rycroft et al. discovered a mathematical correspondence between this physical system and the incompressible Navier–Stokes equations, and developed a projection method similar to Chorin's projection method (1968) for incompressible Newtonian fluids. Here, we improve the original projection method to simulate quasi-static hypo-elastoplasticity, by making three improvements. First, drawing inspiration from the second-order projection method for incompressible Newtonian fluids, we formulate a second-order in time numerical scheme for quasi-static hypo-elastoplasticity. Second, we implement a finite element method for solving the elliptic equations in the projection step, which provides both numerical benefits and flexibility. Third, we develop an adaptive global time-stepping scheme, which can compute accurate solutions in fewer timesteps. Our numerical tests use an example physical model of a bulk metallic glass based on the shear transformation zone theory, but the numerical methods can be applied to any elastoplastic material. 
    more » « less
  3. ; ; ; ; ; ; (, ACM Transactions on Graphics)
    We propose Leapfrog Flow Maps (LFM) to simulate incompressible fluids with rich vortical flows in real time. Our key idea is to use a hybrid velocityimpulse scheme enhanced with leapfrog method to reduce the computational workload of impulse-based flow map methods, while possessing strong ability to preserve vortical structures and fluid details. In order to accelerate the impulse-to-velocity projection, we develop a fast matrix-free Algebraic Multigrid Preconditioned Conjugate Gradient (AMGPCG) solver with customized GPU optimization, which makes projection comparable with impulse evolution in terms of time cost. We demonstrate the performance of our method and its efficacy in a wide range of examples and experiments, such as real-time simulated burning fire ball and delta wingtip vortices. 
    more » « less
  4. ; ; (, SIAM journal on scientific computing)
    In this paper, we consider the numerical approximations for a hydrodynamical model of smectic-A liquid crystals. The model, derived from the variational approach of the modified Oseen– Frank energy, is a highly nonlinear system that couples the incompressible Navier–Stokes equations and a constitutive equation for the layer variable. We develop two linear, second order time marching schemes based on the Invariant Energy Quadratization method for nonlinear terms in the constitutive equation, the projection method for the Navier–Stokes equations, and some subtle implicit-explicit treatments for the convective and stress terms. Moreover, we prove the well-posedness of the linear system and their unconditionally energy stabilities rigorously. Various numerical experiments are presented to demonstrate the stability and the accuracy of the numerical schemes in simulating the dynamics under shear flow and the magnetic field. 
    more » « less
  5. ; ; (, International Journal for Numerical Methods in Engineering)
    Summary We present a spatially varying Robin interface condition for solving fluid‐structure interaction problems involving incompressible fluid flows and nonuniform flexible structures. Recent studies have shown that for uniform structures with constant material and geometric properties, a constant one‐parameter Robin interface condition can improve the stability and accuracy of partitioned numerical solution procedures. In this work, we generalize the parameter to a spatially varying function that depends on the structure's local material and geometric properties, without varying the exact solution of the coupled fluid‐structure system. We present an algorithm to implement the Robin interface condition in an embedded boundary method for coupling a projection‐based incompressible viscous flow solver with a nonlinear finite element structural solver. We demonstrate the numerical effects of the spatially varying Robin interface condition using two example problems: a simplified model problem featuring a nonuniform Euler‐Bernoulli beam interacting with an inviscid flow and a generalized Turek‐Hron problem featuring a nonuniform, highly flexible beam interacting with a viscous laminar flow. Both cases show that a spatially varying Robin interface condition can clearly improve numerical accuracy (by up to two orders of magnitude in one instance) for the same computational cost. Using the second example problem, we also demonstrate and compare two models for determining the local value of the combination function in the Robin interface condition. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email