Related papers: A Singularity-Free Coordinate System for X-point G…
A nonhydrostatic dynamical core has been developed by using the multi-moment finite volume method that ensures the rigorous numerical conservation. To represent the spherical geometry free of polar problems, the cubed-sphere grid is…
We propose a new discrete FFT-based method for computational homogenization of micromechanics on a regular grid that is simple, fast and robust. The discretization scheme is based on a tetrahedral stencil that displays three crucial…
The article proposes a Caputo fractional conjugate gradient (CFCG) method for unconstrained optimization problems which is applicable to smooth as well as non-smooth problmes. The proposed method uses a non-adaptive version of the Caputo…
This paper proposes a nonlinear guidance algorithm for fuel-optimal impulsive trajectories for rendezvous operations close to a reference orbit. The approach involves overparameterized monomial coordinates and a high-order approximation of…
Full waveform inversion (FWI) is an iterative identification process that serves to minimize the misfit of model-based simulated and experimentally measured wave field data, with the goal of identifying a field of parameters for a given…
Stability and stabilization analysis of fractional-order linear time-invariant (FO-LTI) systems with different derivative orders is studied in this paper. First, by using an appropriate linear matrix function, a single-order equivalent…
In this paper, we present a monolithic multigrid method for the efficient solution of flow problems in fractured porous media. Specifically, we consider a mixed-dimensional model which couples Darcy flow in the porous matrix with…
To identify under what conditions guiding-centre or full-orbit tracing should be used, an estimation of the spatial variation of the magnetic field is proposed, not only taking into account gradient and curvature terms but also parallel…
The magnetic phase diagram at zero external field of an ensemble of dipoles with uniaxial anisotropy on a FCC lattice is investigated from tempered Monte Carlo simulations. The uniaxial anisotropy is characterized by a random distribution…
We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry,…
This paper presents a direct data-driven approach for computing robust control invariant (RCI) sets and their associated state-feedback control laws for linear time-invariant systems affected by bounded disturbances. The proposed method…
We present a simple yet powerful framework for solving inverse problems by leveraging automatic differentiation. Our method is broadly applicable whenever a smooth cost function can be defined near the true solution, and a numerical…
Hamilton-Jacobi theory for general relativity provides an elegant covariant formulation of the gravitational field. A general `coordinate-free' method of integrating the functional Hamilton-Jacobi equation for gravity and matter is…
We use the Bianchi-I spacetime to study the local dynamics of a magnetized self-gravitating Fermi gas. The set of Einstein-Maxwell field equations for this gas becomes a dynamical system in a 4-dimensional phase space. We consider a…
Codes were written to simulate the propagation of monochromatic light through a bare optical resonator, using a computational Fourier method to solve the Huygens-Fresnel integral. This was used, in the Fox-Li method, to find the lowest-loss…
The classical Helmholtz problem is applied for modelling the axisymmetric inviscid cusp-ended separated flow around a sphere. Two coordinate systems are employed: polar for initial calculations and parabolic the latter being more suitable…
A semi-implicit finite difference time domain (FDTD) numerical Maxwell solver is developed for full electromagnetic Particle-in-Cell (PIC) codes for the simulations of plasma-based acceleration. The solver projects the volumetric Yee…
Harmonic coordinate conditions in stationary asymptotically flat spacetimes with matter sources have more than one solution. The solutions depend on the degree of smoothness of the metric and its first derivatives, which we wish to impose…
Efficient performance of a number of engineering systems is achieved through different modes of operation - yielding systems described as "hybrid", containing both real-valued and discrete decision variables. Prominent examples of such…
We introduce a general purpose algorithm for rapidly computing certain types of oscillatory integrals which frequently arise in problems connected to wave propagation and general hyperbolic equations. The problem is to evaluate numerically…