Computational Physics
Convolutional neural networks (CNN) are utilized to encode the relation between initial configurations of obstacles and three fundamental quantities in porous media: porosity ($\varphi$), permeability $k$, and tortuosity ($T$). The…
Large High Energy Physics (HEP) experiments adopted a distributed computing model more than a decade ago. WLCG, the global computing infrastructure for LHC, in partnership with the US Open Science Grid, has achieved data management at the…
This work outlines a new multi-physics-compatible immersed rigid body method for Eulerian finite-volume simulations. To achieve this, rigid bodies are represented as a diffuse scalar field and an interface seeding method is employed to…
There is a growing consensus that solutions to complex science and engineering problems require novel methodologies that are able to integrate traditional physics-based modeling approaches with state-of-the-art machine learning (ML)…
In this paper we propose a computational framework for the investigation of the correlated motion between positive and negative ions exposed to the attraction of a bubble surface that mimics the (oscillating) cell membrane. The correlated…
In this paper, three-dimensional numerical simulations of ballooning in spiders using multiple silk threads are performed using the discrete elastic rods method. The ballooning of spiders is hypothesised to be caused by the presence of the…
In this article we present a new formulation and an associated algorithm for the simultaneous numerical simulation of multiple condensed phase explosives in direct contact with each other, which may also be confined by (or interacting with…
It is well known that collective electronic excitations in fullerene C$_{60}$ are manifested as Mie plasmons, and in graphene (the limiting case of an infinitely large fullerene) collective excitations are of plasmon-polariton type. How the…
We study the nodal properties of many-body eigenstates of stationary Schr\"odinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a…
In this paper, a unified gas-kinetic wave-particle scheme (UGKWP) for the disperse dilute gas-particle multiphase flow is proposed. The gas phase is always in the hydrodynamic regime. However, the particle phase covers different flow…
The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gr\'ehan [Part. Part. Syst. Charact. 17…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…
The capillary time-step constraint is the dominant limitation on the applicable time-step in many simulations of interfacial flows with surface tension and, consequently, governs the execution time of these simulations. We propose a…
For decades, there is no unified model developed for Silicon carriers from 4K to above room temperature. In this paper, a unified undoped Silicon low field and high field mobility model for both electron in the <100> and <111> directions…
This article proposes a new statistical numerical method to address gas kinetics problems obeying the Boltzmann equation. This method is inspired from some Monte-Carlo algorithms used in linear transport physics, where virtual particles are…
We develop a computationally and numerically efficient method to calculate binding energies and corresponding wave functions of quantum mechanical three-body problems in low dimensions. Our approach exploits the tensor structure of the…
We develop the resonant mode coupling approximation to calculate the optical spectra of a stack of two photonic crystal slabs. The method is based on a derivation of the input and output resonant vectors in each slab in terms of the Fourier…
In our recent work, we introduced the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) and demonstrated enhanced solver accuracy, scalability, and performance compared to conventional approaches. We…
A high fidelity flow simulation for complex geometries for high Reynolds number ($Re$) flow is still very challenging, which requires more powerful computational capability of HPC system. However, the development of HPC with traditional CPU…
Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR…