Search Results

We present a new hybrid discrete exterior calculus (DEC) and finite difference (FD) method to simulate fully three-dimensional Boussinesq convection in spherical shells subject to internal heating and basal heating, relevant in the ...

• Boussinesq convection in spherical shells.
• Hybrid discrete exterior calculus ...

We develop a set of highly efficient and effective computational algorithms and simulation tools for fluid simulations on a network. The mathematical models are a set of hyperbolic conservation laws on edges of a network, as well as coupling conditions ...

The aim of this work is to introduce a unified description of topology optimization (TO) methods, which modularizes and generalizes all TO methods, both density based and boundary based. This unified description allows for the implementation of a ...

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one ...

• A preconditioning strategy was developed for low-temperature streamer discharges.

Gaussian processes are a widely used alternative to neural networks for non-linear system identification. The method requires computing the inversion of a large covariance matrix. In this work, we introduce our new task-based asynchronous ...

Co-simulation is a widely used solution to enable global simulation of a modular system via the composition of black-boxed simulators. Among co-simulation methods, the IFOSMONDI implicit iterative algorithm, previously introduced by the authors, ...

In this paper, a chain of scientific tools is applied to provide a flexible, efficient, and extensible toolchain that is capable of addressing limitations in corrosion simulation. A previously developed computational model of galvanic ...

• Simulating complex corrosion models yields many challenges in the implementation process.

Particle-in-cell methods employ particle representations of unknown fields but also employ continuum fields for other parts of the problem. Thus projection between particle and continuum bases is required. Moreover, we often need to enforce conservation ...

We present a new software system, PETSc TSAdjoint, for first-order and second-order adjoint sensitivity analysis of time-dependent nonlinear differential equations. The derivative calculation in PETSc TSAdjoint is essentially a high-level algorithmic ...

This paper presents a Python wrapper and extended functionality of the parallel topology optimization framework introduced by Aage et al. (Topology optimization using PETSc: an easy-to-use, fully parallel, open source topology optimization ...$\mathtt{}\mathtt{}\mathtt{}\mathtt{}\mathtt{}\mathtt{}$

In large, collaborative open-source projects, developers must follow good coding standards to ensure the quality and sustainability of the resulting software. This is especially a challenge in high-performance computing projects, which admit a diverse ...

Battery performance is strongly correlated with electrode microstructure. Electrode materials for lithium-ion batteries have complex microstructure geometries that require millions of degrees of freedom to solve the electrochemical system at the ...

Large-scale PDE simulations using high-order finite-element methods on unstructured meshes are an indispensable tool in science and engineering. The widely used open-source PETSc library offers an efficient representation of generic unstructured meshes ...

A computational bottleneck during the solution to multiphase formulations of the incompressible Navier–Stokes equations is often during the implicit solution of the pressure-correction equation that results from operator-splitting ...

• CFD code MFiX was interfaced with the PETSc linear solver library.
• Block Jacobi ...

The computational complexity of solving nonlinear support vector machine (SVM) is prohibitive on large-scale data. In particular, this issue becomes very sensitive when the data represents additional difficulties such as highly imbalanced class ...

We herein describe the performance evaluation of a modular implementation of the MGRIT (MultiGrid-In-Time) algorithm within the context of the PETSc (the Portable, Extensible Toolkit for Scientific computing) library. Our aim is to give the PETSc ...

Emerging many-core CPU architectures with high degrees of single-instruction, multiple data (SIMD) parallelism promise to enable increasingly ambitious simulations based on partial differential equations (PDEs) via extreme-scale computing. However, such ...

The paper deals with handling multiple subdomains per computational core in the PERMON toolbox, namely in the PermonFLLOP module, to fully exploit the potential of the Total Finite Element Tearing and Interconnecting (TFETI) domain decomposition method (...

OpenFOAM is a widely used open source framework for simulation in several areas of computational fluid dynamics and engineering. As a partial differential equation (PDE)-based framework, OpenFOAM suffers from a performance bottleneck in solving large-...

Multilevel balancing domain decomposition by constraints (BDDC) deluxe algorithms are developed for the saddle point problems arising from mixed formulations of Darcy flow in porous media. In addition to the standard no-net-flux constraints on each face, ...