Featured Researches

Computational Physics

Benchmarking of Numerical Models for Wave Overtopping at Dikes with Shallow Mildly Sloping Foreshores: Accuracy versus Speed

To accurately predict the consequences of nearshore waves, coastal engineers often employ numerical models. A variety of these models, broadly classified as either phase-resolving or phase-averaged, exist; each with strengths and limitations owing to the physical schematization of processes within them. Models which resolve the vertical flow structure or the full wave spectrum (i.e. sea-swell (SS) and infragravity (IG) waves) are considered more accurate, but also more computationally demanding than those with approximations. Here, we assess the speed-accuracy trade-off of six well-known wave models for overtopping (q), under shallow foreshore conditions. The results demonstrate that: i) q is underestimated by an order of magnitude when IG waves are neglected; ii) using more computationally-demanding models does not guarantee more accurate results; and iii) with empirical corrections to account for IG waves, phase-averaged models like SWAN can perform on par, if not better than, phase-resolving models but with far less computational effort.

Read more
Computational Physics

Beyond moments: relativistic Lattice-Boltzmann methods for radiative transport in computational astrophysics

We present a new method for the numerical solution of the radiative-transfer equation (RTE) in multidimensional scenarios commonly encountered in computational astrophysics. The method is based on the direct solution of the Boltzmann equation via an extension of the Lattice Boltzmann (LB) equation and allows to model the evolution of the radiation field as it interacts with a background fluid, via absorption, emission, and scattering. As a first application of this method, we restrict our attention to a frequency independent ("grey") formulation within a special-relativistic framework, which can be employed also for classical computational astrophysics. For a number of standard tests that consider the performance of the method in optically thin, optically thick and intermediate regimes with a static fluid, we show the ability of the LB method to produce accurate and convergent results matching the analytic solutions. We also contrast the LB method with commonly employed moment-based schemes for the solution of the RTE, such as the M1 scheme. In this way, we are able to highlight that the LB method provides the correct solution for both non-trivial free-streaming scenarios and the intermediate optical-depth regime, for which the M1 method either fails or provides inaccurate solutions. When coupling to a dynamical fluid, on the other hand, we present the first self-consistent solution of the RTE with LB methods within a relativistic-hydrodynamic scenario. Finally, we show that besides providing more accurate results in all regimes, the LB method features smaller or comparable computational costs compared to the M1 scheme.

Read more
Computational Physics

Black-box inhomogeneous preconditioning for self-consistent field iterations in density functional theory

We propose a new preconditioner based on the local density of states for computing the self-consistent problem in Kohn-Sham density functional theory. This preconditioner is inexpensive and able to cure the long-range charge sloshing known to hamper convergence in large, inhomogeneous systems such as clusters and surfaces. It is based on a parameter-free and physically motivated approximation to the independent-particle susceptibility operator, appropriate for both metals and insulators. It can be extended to semiconductors by using the macroscopic electronic dielectric constant as a parameter in the model. We test our preconditioner successfully on inhomogeneous systems containing metals, insulators, semiconductors and vacuum.

Read more
Computational Physics

CCDCGAN: Inverse design of crystal structures

Autonomous materials discovery with desired properties is one of the ultimate goals for modern materials science. Applying the deep learning techniques, we have developed a generative model which can predict distinct stable crystal structures by optimizing the formation energy in the latent space. It is demonstrated that the optimization of physical properties can be integrated into the generative model as on-top screening or backwards propagator, both with their own advantages. Applying the generative models on the binary Bi-Se system reveals that distinct crystal structures can be obtained covering the whole composition range, and the phases on the convex hull can be reproduced after the generated structures are fully relaxed to the equilibrium. The method can be extended to multicomponent systems for multi-objective optimization, which paves the way to achieve the inverse design of materials with optimal properties.

Read more
Computational Physics

CIF2WAN: A Tool to Generate Input Files for Electronic Structure Calculations with Wannier90

The generation of input files for density functional theory (DFT) programs must often be manually done by researchers. If one wishes to produce a maximally localized wannier functions (MLWFs) the calculation consists of several separate files that must be formatted correctly in order for the program to work properly. Many of the inputs are repeated throughout the files and can be easily automated. In this work, a program is presented to generate all of the input files needed to produce wannier functions with Wannier90 starting from open source DFT programs such as Quantum Espresso, Abinit, and Siesta. In addition, the input files for WannierTools are also included for those that wish to produce surface green's functions for the generation of surface state bands. The program presented allows for users new to DFT to use the programs with minimal understanding of parameters needed to produce good results, in addition, this program allows for researchers who are advanced DFT users to utilize this program for high throughput wannier calculations.

Read more
Computational Physics

CMInject: Python framework for the numerical simulation of nanoparticle injection pipelines

CMInject simulates nanoparticle injection experiments of particles with diameters in the micrometer to nanometer-regime, e.g., for single-particle-imaging experiments. Particle-particle interactions and particle-induced changes in the surrounding fields are disregarded, due to low nanoparticle concentration in these experiments. CMInject's focus lies on the correct modeling of different forces on such particles, such as fluid-dynamics or light-induced interactions, to allow for simulations that further the scientific development of nanoparticle injection pipelines. To provide a usable basis for this framework and allow for a variety of experiments to be simulated, we implemented first specific force models: fluid drag forces, Brownian motion, and photophoretic forces. For verification, we benchmarked a drag-force-based simulation against a nanoparticle focusing experiment. We envision its use and further development by experimentalists, theorists, and software developers.

Read more
Computational Physics

Calibration of projection-based reduced-order models for unsteady compressible flows

An analysis of calibration for reduced-order models (ROMs) is presented in this work. The Galerkin and least-squares Petrov-Galerkin (LSPG) methods are tested on compressible flows involving a disparity of temporal scales. A novel calibration strategy is proposed for the LSPG method and two test cases are analyzed. The first consists of a subsonic airfoil flow where boundary layer instabilities are responsible for trailing-edge noise generation and the second comprises a supersonic airfoil flow with a transient period where a detached shock wave propagates upstream at the same time that shock-vortex interaction occurs at the trailing edge. Results show that calibration produces stable and long-time accurate for both cases. In order to reduce the computational costs of the LSPG models, an accelerated greedy missing point estimation (MPE) algorithm is employed for hyper-reduction. For the first case investigated, LSPG solutions obtained with hyper-reduction show good comparison with those obtained by the full order model. However, for the supersonic case the transient features of the flow need to be properly captured by the sampled points. Otherwise, the dynamics of the moving shock wave are not fully recovered. The impact of different time-marching schemes is also assessed and, differently than reported in literature, Galerkin models are shown to be more accurate than those computed by LSPG when the non-conservative form of the Navier-Stokes equations are solved. For the supersonic case, the Galerkin and LSPG models (without hyper-reduction) capture the overall dynamics of the detached and oblique shock waves along the airfoil. However, when shock-vortex interaction occurs at the trailing-edge, the Galerkin ROM is able to capture the high-frequency fluctuations from vortex shedding while the LSPG presents a more dissipative solution, not being able to recover the flow dynamics.

Read more
Computational Physics

Carbon nanotubes collapse phase diagram with arbitrary number of walls. Collapse modes and macroscopic analog

Carbon nanotubes tend to collapse when their diameters exceed a certain threshold, or when a sufficiently large external pressure is applied on their walls. Their radial stability of tubes has been studied in each of these cases, however a general theory able to predict collapse is still lacking. Here, we propose a simple model predicting stability limits as a function of the tube diameter, the number of walls and the pressure. The model is supported by atomistic simulations, experiments, and is used to plot collapse phase diagrams. We have identified the most stable carbon nanotube, which can support a maximum pressure of 18 GPa before collapsing. The latter was identified as a multiwall tube with an internal tube diameter of 12nm and 30 walls. This maximum pressure is lowered depending on the internal tube diameter and the number of walls. We then identify a tube diameter domain in which the radial mechanical stability can be treated as equivalent to macroscopic tubes, known to be described by the canonical Lévy-Carrier law. This multiscale behavior is shown to be in good agreement with experiments based on O-ring gaskets collapse, proposed as a simple macroscopic parallel to nanotubes in this domain.

Read more
Computational Physics

Cation interstitial diffusion in lead telluride and cadmium telluride studied by means of neural network potential based molecular dynamics simulations

Using a recently developed approach to represent ab initio based force fields by a neural network potential, we perform molecular dynamics simulations of lead telluride (PbTe) and cadmium telluride (CdTe) crystals. In particular, we study the diffusion of a single cation interstitial in these two systems. Our simulations indicate that the interstitials migrate via two distinct mechanisms: through hops between interstitial sites and through exchanges with lattice atoms. We extract activation energies for both of these mechanisms and show how the temperature dependence of the total diffusion coefficient deviates from Arrhenius behaviour. The accuracy of the neural network approach is estimated by comparing the results for three different independently trained potentials.

Read more
Computational Physics

Central moments multiple relaxation time LBM for hemodynamic simulations in intracranial aneurysms: An in-vitro validation study using PIV and PC-MRI

The lattice Boltzmann method (LBM) has recently emerged as an efficient alternative to classical Navier-Stokes solvers. This is particularly true for hemodynamics in complex geometries. However, in its most basic formulation, {i.e.} with the so-called single relaxation time (SRT) collision operator, it has been observed to have a limited stability domain in the Courant/Fourier space, strongly constraining the minimum time-step and grid size. The development of improved collision models such as the multiple relaxation time (MRT) operator in central moments space has tremendously widened the stability domain, while allowing to overcome a number of other well-documented artifacts, therefore opening the door for simulations over a wider range of grid and time-step sizes. The present work focuses on implementing and validating a specific collision operator, the central Hermite moments multiple relaxation time model with the full expansion of the equilibrium distribution function, to simulate blood flows in intracranial aneurysms. The study further proceeds with a validation of the numerical model through different test-cases and against experimental measurements obtained via stereoscopic particle image velocimetry (PIV) and phase-contrast magnetic resonance imaging (PC-MRI). For a patient-specific aneurysm both PIV and PC-MRI agree fairly well with the simulation. Finally, low-resolution simulations were shown to be able to capture blood flow information with sufficient accuracy, as demonstrated through both qualitative and quantitative analysis of the flow field {while leading to strongly reduced computation times. For instance in the case of the patient-specific configuration, increasing the grid-size by a factor of two led to a reduction of computation time by a factor of 14 with very good similarity indices still ranging from 0.83 to 0.88.}

Read more

Ready to get started?

Join us today