No Arabic abstract
FFT-based solvers are increasingly used by many researcher groups interested in modelling the mechanical behavior associated to a heterogeneous microstructure. A development is reported here that concerns the viscoelastic behavior of composite structures generally studied experimentally through Dynamic Mechanical Analysis (DMA). A parallelized computation code developed under complex-valued quantities provides virtual DMA experiments directly in the frequency domain on a heterogenous system described by a voxel grid of mechanical properties. The achieved precision and computation times are very good. An effort has been made to show the application of such virtual DMA tool starting from two examples found in the literature: the modelling of glassy/amorphous systems at a small scale and the modelling of experimental data obtained in temperature sweeping mode by DMA on a particulate composite made of glass beads and a polystyrene matrix, at a larger scale. Both examples show how virtual DMA can contribute to question, analyze, understand relaxation phenomena either on the theoretical or experimental point of view.
The geodesic has a fundamental role in physics and in mathematics: roughly speaking, it represents the curve that minimizes the arc length between two points on a manifold. We analyze a basic but misinterpreted difference between the Lagrangian that gives the arc length of a curve and the one that describes the motion of a free particle in curved space. Although they provide the same formal equations of motion, they are not equivalent. We explore this difference from a geometrical point of view, where we observe that the non-equivalence is nothing more than a matter of symmetry. As applications, some distinct models are studied. In particular, we explore the standard free relativistic particle, a couple of spinning particle models and also the forceless mechanics formulated by Hertz.
We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four formed beams along East-West via exact phasing, tiling a sky area of ~250 square degrees. A zero-padding approximation is employed to improve chromatic beam alignment across the wide bandwidth of 400 to 800 MHz. We up-channelize the data in order to achieve fine spectral resolution of $Delta u$=24 kHz and time cadence of 0.983 ms, desirable for detecting transient and dispersed signals such as those from FRBs.
In this short note, we present a new technique to accelerate the convergence of a FFT-based solver for numerical homogenization of complex periodic media proposed by Moulinec and Suquet in 1994. The approach proceeds from discretization of the governing integral equation by the trigonometric collocation method due to Vainikko (2000), to give a linear system which can be efficiently solved by conjugate gradient methods. Computational experiments confirm robustness of the algorithm with respect to its internal parameters and demonstrate significant increase of the convergence rate for problems with high-contrast coefficients at a low overhead per iteration.
The problem of optimization of the rolling dynamics model is considered. That providing safe movement at high frequency when interacting with the railway. Moreover, allowing to evaluate the dynamic parameters when designing new and modernizing existing locomotives. The object of this research is a rail transport dynamic system model. The articles purpose is to increase the efficiency of the digital hardware in the rolling stock loop model by optimizing the organization of the computing process. The mathematical model analysis of the object made it possible to attribute it to the class of hard real-time systems. The computation of the model phase variables with different frequencies is necessary to optimize the simulation time of the train movements and is performed by splitting the original algorithm into parallel threads. The developed planning algorithm and the cyclic schedule implementation for the model of a dynamic real-time object consider microarchitecture solutions of symmetric multiprocessor systems with shared memory and methods for optimizing software tools. The experiments confirmed the operability of the optimized model. Also, allow us to recommend it for practical use in studying objects and determine the dynamic force of trolley structural elements during operation. Analysis of the optimized model simulation results, using cyclic schedules shows the correspondence of the obtained simulation results to the standard. The main advantage of the model is the increase in productivity when performing data processing by reducing the processor time. The optimized cyclic schedule algorithm of the semi-natural modeling platform is used for the subsequent development of the control system in real and accelerated time scales.
Effective medium theory aims to describe a complex inhomogeneous material in terms of a few important macroscopic parameters. To characterise wave propagation through an inhomogeneous material, the most crucial parameter is the effective wavenumber. For this reason, there are many published studies on how to calculate a single effective wavenumber. Here we present a proof that there does not exist a unique effective wavenumber; instead, there are an infinite number of such (complex) wavenumbers. We show that in most parameter regimes only a small number of these effective wavenumbers make a significant contribution to the wave field. However, to accurately calculate the reflection and transmission coefficients, a large number of the (highly attenuating) effective waves is required. For clarity, we present results for scalar (acoustic) waves for a two-dimensional material filled (over a half space) with randomly distributed circular cylindrical inclusions. We calculate the effective medium by ensemble averaging over all possible inhomogeneities. The proof is based on the application of the Wiener-Hopf technique and makes no assumption on the wavelength, particle boundary conditions/size, or volume fraction. This technique provides a simple formula for the reflection coefficient, which can be explicitly evaluated for monopole scatterers. We compare results with an alternative numerical matching method.