187 research outputs found

    A survey of high level frameworks in block-structured adaptive mesh refinement packages

    Get PDF
    pre-printOver the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a representative set of SAMR packages and SAMR-based codes that have been in existence for half a decade or more, have a reasonably sized and active user base outside of their home institutions, and are publicly available. The set consists of a mix of SAMR packages and application codes that cover a broad range of scientific domains. We look at their high-level frameworks, their design trade-offs and their approach to dealing with the advent of radical changes in hardware architecture. The codes included in this survey are BoxLib, Cactus, Chombo, Enzo, FLASH, and Uintah

    AMRA: An Adaptive Mesh Refinement Hydrodynamic Code for Astrophysics

    Get PDF
    Implementation details and test cases of a newly developed hydrodynamic code, AMRA, are presented. The numerical scheme exploits the adaptive mesh refinement technique coupled to modern high-resolution schemes which are suitable for relativistic and non-relativistic flows. Various physical processes are incorporated using the operator splitting approach, and include self-gravity, nuclear burning, physical viscosity, implicit and explicit schemes for conductive transport, simplified photoionization, and radiative losses from an optically thin plasma. Several aspects related to the accuracy and stability of the scheme are discussed in the context of hydrodynamic and astrophysical flows.Comment: 41 pages, 21 figures (9 low-resolution), LaTeX, requires elsart.cls, submitted to Comp. Phys. Comm.; additional documentation and high-resolution figures available from http://www.camk.edu.pl/~tomek/AMRA/index.htm

    Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media

    Get PDF
    This paper introduces a new discrete fracture model accounting for non-isothermal compositional multiphase Darcy flows and complex networks of fractures with intersecting, immersed and non immersed fractures. The so called hybrid-dimensional model using a 2D model in the fractures coupled with a 3D model in the matrix is first derived rigorously starting from the equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully implicit time integration combined with the Vertex Approximate Gradient (VAG) finite volume scheme which is adapted to polyhedral meshes and anisotropic heterogeneous media. The fully coupled systems are assembled and solved in parallel using the Single Program Multiple Data (SPMD) paradigm with one layer of ghost cells. This strategy allows for a local assembly of the discrete systems. An efficient preconditioner is implemented to solve the linear systems at each time step and each Newton type iteration of the simulation. The numerical efficiency of our approach is assessed on different meshes, fracture networks, and physical settings in terms of parallel scalability, nonlinear convergence and linear convergence

    ISCR annual report FY 1998

    Full text link

    Multiphysics simulations: challenges and opportunities.

    Full text link

    Software for Exascale Computing - SPPEXA 2016-2019

    Get PDF
    This open access book summarizes the research done and results obtained in the second funding phase of the Priority Program 1648 "Software for Exascale Computing" (SPPEXA) of the German Research Foundation (DFG) presented at the SPPEXA Symposium in Dresden during October 21-23, 2019. In that respect, it both represents a continuation of Vol. 113 in Springer’s series Lecture Notes in Computational Science and Engineering, the corresponding report of SPPEXA’s first funding phase, and provides an overview of SPPEXA’s contributions towards exascale computing in today's sumpercomputer technology. The individual chapters address one or more of the research directions (1) computational algorithms, (2) system software, (3) application software, (4) data management and exploration, (5) programming, and (6) software tools. The book has an interdisciplinary appeal: scholars from computational sub-fields in computer science, mathematics, physics, or engineering will find it of particular interest

    Techniques of High Performance Reservoir Simulation for Unconventional Challenges

    Get PDF
    The quest to improve the performance of reservoir simulators has been evolving with the newly encountered challenges of modeling more complex recovery mechanisms and related phenomena. Reservoir subsidence, fracturing and fault reactivation etc. require coupled flow and poroelastic simulation. These features, in turn, bring a heavy burden on linear solvers. The booming unconventional plays such as shale/tight oil in North America demand reservoir simulation techniques to handle more physics (or more hypotheses). This dissertation deals with three aspects in improving the performance of reservoir simulation toward these unconventional challenges. Compositional simulation is often required for many reservoir studies with complex recovery mechanisms such as gas inject. But, it is time consuming and its parallelization often suffers sever load imbalance problems. In the first section, a novel approach based on domain over-decomposition is investigated and implemented to improve the parallel performance of compositional simulation. For a realistic reservoir case, it is shown the speedup is improved from 29.27 to 62.38 on 64 processors using this technique. Another critical part that determines the performance of a reservoir simulator is the linear solver. In the second section, a new type of linear solver based the combinatorial multilevel method (CML) is introduced and investigated for several reservoir simulation applications. The results show CML has better scalability and performance empirically and is well-suited for coupled poroelastic problems. These results also suggest that CML might be a promising way of precondition for flow simulation with and without coupled poroelastic calculations. In order to handle unconventional petroleum fluid properties for tight oil, the third section incorporates a simulator with extended vapor-liquid equilibrium calculations to consider the capillarity effect caused by the dynamic nanopore properties. The enhanced simulator can correctly capture the pressure dependent impact of the nanopore on rock and fluid properties. It is shown inclusion of these enhanced physics in simulation will lead to significant improvements in field operation decision-making and greatly enhance the reliability of recovery predictions

    Time integration for complex fluid dynamics

    Get PDF
    2021 Fall.Includes bibliographical references.Efficient and accurate simulation of turbulent combusting flows in complex geometry remains a challenging and computationally expensive proposition. A significant source of computational expense is in the integration of the temporal domain, where small time steps are required for the accurate resolution of chemical reactions and long solution times are needed for many practical applications. To address the small step sizes, a fourth-order implicit-explicit additive Runge-Kutta (ARK4) method is developed to integrate the stiff chemical reactions implicitly while advancing the convective and diffusive physics explicitly in time. Applications involving complex geometry, stiff reaction mechanisms, and high-order spatial discretizations are challenged by stability issues in the numerical solution of the nonlinear problem that arises from the implicit treatment of the stiff term. Techniques for maintaining a physical thermodynamic state during the numerical solution of the nonlinear problem, such as placing constraints on the nonlinear solver and the use of a nonlinear optimizer to find valid thermodynamic states, are proposed and tested. Verification and validation are performed for the new adaptive ARK4 method using lean premixed flames burning hydrogen, showing preservation of 4th-order error convergence and recovery of literature results. ARK4 is then applied to solve lean, premixed C3H8-air combustion in a bluff-body combustor geometry. In the two-dimensional case, ARK4 provides a 70× speedup over the standard explicit four-stage Runge-Kutta method and, for the three-dimensional case, three-orders-of-magnitude-larger time step sizes are achieved. To further increase the computational scaling of the algorithms, parallel-in-time (PinT) techniques are explored. PinT has the dual benefit of providing parallelization to long temporal domains as well as taking advantage of hardware trends towards more concurrency in modern high-performance computing platforms. Specifically, the multigrid reduction-in-time (MGRIT) method is adapted and enhanced by adding adaptive mesh refinement (AMR) in time. This creates a space-time algorithm with efficient solution-adaptive grids. The new MGRIT+AMR algorithm is first verified and validated using problems dominated by diffusion or characterized by time periodicity, such as Couette flow and Stokes second problem. The adaptive space-time parallel algorithm demonstrates up to a 13.7× speedup over a time-sequential algorithm for the same solution accuracy. However, MGRIT has difficulties when applied to solve practical fluid flows, such as turbulence, governed by strong hyperbolic partial differential equations. To overcome this challenge, the multigrid operations are modified and applied in a novel way by exploiting the space-time localization of fine turbulence scales. With these new operators, the coarse-scale errors are advected out of the temporal domain while the fine-scale dynamics iterate to equilibrium. This leads to rapid convergence of the bulk flow, which is important for computing macroscopic properties useful for engineering purposes. The novel multigrid operations are applied to the compressible inviscid Taylor-Green vortex flow and the convergence of the low-frequency modes is achieved within a few iterations. Future work will be focused on a performance study for practical highly turbulent flows

    A summary of my twenty years of research according to Google Scholars

    Get PDF
    I am David Pardo, a researcher from Spain working mainly on numerical analysis applied to geophysics. I am 40 years old, and over a decade ago, I realized that my performance as a researcher was mainly evaluated based on a number called \h-index". This single number contains simultaneously information about the number of publications and received citations. However, dif- ferent h-indices associated to my name appeared in di erent webpages. A quick search allowed me to nd the most convenient (largest) h-index in my case. It corresponded to Google Scholars. In this work, I naively analyze a few curious facts I found about my Google Scholars and, at the same time, this manuscript serves as an experiment to see if it may serve to increase my Google Scholars h-index
    corecore