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.


  • NSF Public Access
  • Search Results
  • An efficient tree-topological local mesh refinement on Cartesian grids for multiple moving objects in incompressible flow
Title: An efficient tree-topological local mesh refinement on Cartesian grids for multiple moving objects in incompressible flow
This paper develops a tree-topological local mesh refinement (TLMR) method on Cartesian grids for the simulation of bio-inspired flow with multiple moving objects. The TLMR nests refinement mesh blocks of structured grids to the target regions and arrange the blocks in a tree topology. The method solves the time-dependent incompressible flow using a fractional-step method and discretizes the Navier-Stokes equation using a finite-difference formulation with an immersed boundary method to resolve the complex boundaries. When iteratively solving the discretized equations across the coarse and fine TLMR blocks, for better accuracy and faster convergence, the momentum equation is solved on all blocks simultaneously, while the Poisson equation is solved recursively from the coarsest block to the finest ones. When the refined blocks of the same block are connected, the parallel Schwarz method is used to iteratively solve both the momentum and Poisson equations. Convergence studies show that the algorithm is second-order accurate in space for both velocity and pressure, and the developed mesh refinement technique is benchmarked and demonstrated by several canonical flow problems. The TLMR enables a fast solution to an incompressible flow problem with complex boundaries or multiple moving objects. Various bio-inspired flows of multiple moving objects show that the solver can save over 80% computational time, proportional to the grid reduction when refinement is applied.  more » « less
Award ID(s):
1931929
PAR ID:
10473217
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
ELSEVIER
Date Published:
Journal Name:
Journal of Computational Physics
Volume:
479
Issue:
C
ISSN:
0021-9991
Page Range / eLocation ID:
111983
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, American Institute of Aeronautics and Astronautics)
    This paper presents a robust mesh moving solver developed to address moving boundary problems. Crucially, the resulting deformed mesh retains the same topology as the original mesh without being overly distorted. The mesh is treated as an elastic material, and the deformation of the computational domain resulting from moving boundaries is determined by solving the equilibrium linear elasticity equations. The linear elasticity equations are discretized by the classic Galerkin finite element method and solved by the block conjugate gradient iterative method. To maintain the quality of the mesh after motion, the Young's modulus of each element is weighted by the reciprocal of the distance between the element center and the moving boundaries. The effectiveness of this approach is demonstrated through a set of 2D and 3D test cases featuring prescribed translational and/or rotational motion of the embedded object. The method is now ready for integration into our existing in-house CFD solvers. 
    more » « less
  2. ; (, International Journal of Computer Mathematics)
    Nonlinear Kawarada equations have been used to model solid fuel combustion processes in the oil industry. An effective way to approximate solutions of such equations is to take advantage of the finite difference configurations. Traditionally, the nonlinear term of the equation is linearized while the numerical stability of a difference scheme is investigated. This leaves certain ambiguity and uncertainty in the analysis. Based on nonuniform grids generated through a quenching-seeking moving mesh method in space and adaptation in time, this paper introduces a completely new stability analysis of the approximation without freezing the nonlinearity involved. Pointwise orders of convergence are investigated numerically. Simulation experiments are carried out to accompany the mathematical analysis to strengthen our conclusions. 
    more » « less
  3. ; ; ; ; (, Numerical Linear Algebra with Applications)
    Abstract It is known that the solutions to space‐fractional diffusion equations exhibit singularities near the boundary. Therefore, numerical methods discretized on the composite mesh, in which the mesh size is refined near the boundary, provide more precise approximations to the solutions. However, the coefficient matrices of the corresponding linear systems usually lose the diagonal dominance and are ill‐conditioned, which in turn affect the convergence behavior of the iteration methods.In this work we study a finite volume method for two‐sided fractional diffusion equations, in which a locally refined composite mesh is applied to capture the boundary singularities of the solutions. The diagonal blocks of the resulting three‐by‐three block linear system are proved to be positive‐definite, based on which we propose an efficient block Gauss–Seidel method by decomposing the whole system into three subsystems with those diagonal blocks as the coefficient matrices. To further accelerate the convergence speed of the iteration, we use T. Chan's circulant preconditioner31as the corresponding preconditioners and analyze the preconditioned matrices' spectra. Numerical experiments are presented to demonstrate the effectiveness and the efficiency of the proposed method and its strong potential in dealing with ill‐conditioned problems. While we have not proved the convergence of the method in theory, the numerical experiments show that the proposed method is convergent. 
    more » « less
  4. ; ; ; (, arXiv but under review in Monthly Weather Review)
    This paper examines the application of adaptive mesh refinement (AMR) in the field of numerical weather prediction (NWP). We implement and assess two distinct AMR approaches and evaluate their performance through standard NWP benchmarks. In both cases, we solve the fully compressible Euler equations, fundamental to many non-hydrostatic weather models. The first approach utilizes oct-tree cell-based mesh refinement coupled with a high-order discontinuous Galerkin method for spatial discretization. In the second approach, we employ level-based AMR with the finite difference method. Our study provides insights into the accuracy and benefits of employing these AMR methodologies for the multi-scale problem of NWP. Additionally, we explore essential properties including their impact on mass and energy conservation. Moreover, we present and evaluate an AMR solution transfer strategy for the tree-based AMR approach that is simple to implement, memory-efficient, and ensures conservation for both flow in the box and sphere. Furthermore, we discuss scalability, performance portability, and the practical utility of the AMR methodology within an NWP framework -- crucial considerations in selecting an AMR approach. The current de facto standard for mesh refinement in NWP employs a relatively simplistic approach of static nested grids, either within a general circulation model or a separately operated regional model with loose one-way synchronization. It is our hope that this study will stimulate further interest in the adoption of AMR frameworks like AMReX in NWP. These frameworks offer a triple advantage: a robust dynamic AMR for tracking localized and consequential features such as tropical cyclones, extreme scalability, and performance portability. 
    more » « less
  5. (, Transactions of the American Mathematical Society)
    Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email