Independent Tensor Components Research Articles

Fast normalized cross-correlation for template matching with rotations

AbstractNormalized cross-correlation is the reference approach to carry out template matching on images. When it is computed in Fourier space, it can handle efficiently template translations but it cannot do so with template rotations. Including rotations requires sampling the whole space of rotations, repeating the computation of the correlation each time.This article develops an alternative mathematical theory to handle efficiently, at the same time, rotations and translations. Our proposal has a reduced computational complexity because it does not require to repeatedly sample the space of rotations. To do so, we integrate the information relative to all rotated versions of the template into a unique symmetric tensor template -which is computed only once per template-. Afterward, we demonstrate that the correlation between the image to be processed with the independent tensor components of the tensorial template contains enough information to recover template instance positions and rotations. Our proposed method has the potential to speed up conventional template matching computations by a factor of several magnitude orders for the case of 3D images.

ON THE RECONSTRUCTION OF PRESTRESS FIELDS IN A HOLLOW CYLINDER

The present research is devoted to the development of the theoretical foundations of nondestructive acoustic method for identifying inhomogeneous prestress fields in a hollow cylinder, depending on the probing loading type. A linearized model of steady oscillations of an elastic body in the presence of an inhomogeneous prestress field of arbitrary nature is considered in the standard and weak formulations. On the basis of this model, we formulate a problem for a cantilever-clamped prestressed hollow cylinder that performs steady axisymmetric vibrations under three types of probing loading. A corresponding weak formulation of the problem in the cylindrical coordinate system is presented, in which six independent components of the prestress tensor are taken into account. At that, a case of prestress fields obtained by applying some initial mechanical external static load is considered. In the presence and absence of prestresses of various types, amplitude-frequency dependences are analyzed, and resonant and natural frequencies are found in a wide frequency range. Numerical calculations were carried out using the FEM on a non-uniform grid; mesh refinement is carried out in the vicinity of the boundary points, where the type of boundary conditions changes. Based on the numerical solution of an auxiliary set of direct problems, seven types of prestress fields are constructed, differing in the types of initial loading, most often encountered in practice. To assess the possibility of implementing the procedure for reconstructing prestresses of each of the considered types, a sensitivity analysis was additionally performed, which showed that for some prestress types there are frequencies and types of probing loading for which the presence of prestress is practically not manifested. The sensitivity analysis performed made it possible to implement the optimal method of probing loading when solving the inverse coefficient problem. The statement of the new inverse problem on the restoration of arbitrary inhomogeneous prestress fields in the considered finite hollow cylinder is formulated. When restoring the prestress of a given structure, the inverse problem is reduced to finding a set of parameters from an ill-conditioned algebraic system, which was studied with the help of the A.N. Tikhonov regularization method. Additional data for solving the inverse problem was obtained on the basis of probing both via a single load and via combined probing modes. It has been found that it is most effective to use a combined loading mode and use a sufficiently wide frequency range when selecting sounding frequencies. The results of computational experiments on the reconstruction of six components of the prestress tensor are presented and analyzed, and recommendations are proposed for choosing the optimal modes of acoustic sounding.

Onsager Reciprocal Relation between Anomalous Transverse Coefficients of an Anisotropic Antiferromagnet.

Whenever two irreversible processes occur simultaneously, time-reversal symmetry of microscopic dynamics gives rise, on a macroscopic level, to Onsager's reciprocal relations, which impose constraints on the number of independent components of any transport coefficient tensor. Here, we show that in the antiferromagnetic YbMnBi_{2}, which displays a strong temperature-dependent anisotropy, Onsager's reciprocal relations are strictly satisfied for anomalous electric (σ_{ij}^{A}) and anomalous thermoelectric (α_{ij}^{A}) conductivity tensors. In contradiction with what was recently reported by Pan etal. [Nat. Mater. 21, 203 (2022)NMAACR1476-112210.1038/s41563-021-01149-2], we find that σ_{ij}^{A}(H)=σ_{ji}^{A}(-H) and α_{ij}^{A}(H)=α_{ji}^{A}(-H). This equality holds in the whole temperature window irrespective of the relative weights of the intrinsic or extrinsic mechanisms. The α_{ij}^{A}/σ_{ij}^{A} ratio is close to k_{B}/e at room temperature but peaks to an unprecedented magnitude of 2.9k_{B}/e at ∼150 K, which may involve nondegenerate carriers of small Fermi surface pockets.

Invariance embedded physics-infused deep neural network-based sub-grid scale models for turbulent flows

In this paper, we present two novel physics-infused neural network (NN) architectures that satisfy various invariance conditions for constructing efficient and robust Deep Learning (DL)-based sub-grid scale (SGS) turbulence models for use in Large Eddy Simulation (LES) procedures widely used in various fluid engineering applications. The first architecture, called tensor basis neural networks (TBNN), recently proposed in the context of Reynolds-averaged Navier–Stokes (RANS) modeling, and introduced herein for SGS modeling of wall-bounded turbulence, embeds the analytical expansion of the SGS turbulence stresses into integrity basis tensors composed of the symmetric and anti-symmetric parts of the resolved velocity gradient tensor, thus incorporating the Galilean, rotational and reflectional invariances. In our second approach, a relatively simple yet powerful architecture, called Galilean Invariance embedded Neural Network (GINN) incorporates the Galilean invariance, and takes in as inputs the independent components of the integrity basis tensors in addition to the invariant inputs in a single input layer. Explicit filtering of data from direct simulations of the canonical channel flow at two friction Reynolds numbers Reτ≈395 and 590 provided accurate data for training and testing these models. Both sets of models are used to predict the SGS stresses for feature datasets generated with different filter widths, and at different Reynolds numbers. The GINN model yields less prediction error on test datasets (mean squared error ∼0.4) compared to the TBNN model (mean squared error ∼0.5). Upon comparison, it is revealed that GINN has better feature extraction capacity owing to its ability to establish relations between and extract information from cross-components of the integrity basis tensors and the SGS stresses. Based on their predictive performance, both models proposed herein have shown great promise for use in a posteriori actual LESs. The present work illustrates the significance of physics infusion as well as invariance embedding into NN architecture for constructing efficient and robust DL-based SGS turbulent models in achieving superior predictive efficacy.

ELECTROMAGNETIC COUPLING OF PIEZOELECTRIC/FERRITE COMPOSITE WITH INITIAL STRESS STATE

A mathematical model of electromagnetic thermoelasticity for an initial-stressed transversal-isotropic composite with piezoelectric magnetostrictive phases has been developed. To solve the related boundary value problem of electromagnetic thermoe-lasticity, the Green function method was used as part of a generalized singular approximation of the statistical mechanics of composites, taking into account the initial stress state of the representative domain of the composite at micro- and macro-levels. The solution of the problem of "effective module" for tensors of effective elastic, piezomechanical, magnetostrictive properties, dielectric permittivity and magnetic permit-tivity, temperature, pyroelectric, pyromagnetic coefficients and (which appeared only at the macro level) electromagnetic and magnetoelectric couplings of a quasi-periodic composite with an initial electromagnetic elastic state was obtained. The solutions for the desired tensors of effective properties of the quasi-periodic composite are presented in the form of analytical formulas of simple linear decompositions by solutions for tensors of effective properties of the periodic structure and statistical mixture, decomposition coefficients are the coefficient of "periodicity" (correlation of quasi-periodic and periodic structures) p and "disordering" 1-p, respectively. Results of calculation of all independent components of tensors of effective coefficients of electromagnetic and magnetoelectric couplings of different structures (periodic, quasi-periodic and statistical mixture) are presented in unidirectional direction of fibrous composite "PZT-4/ferrite" with axisymmet-ric tensor of initial macrostrain of composite. For a quasi-periodic composite (with initial macrostrain), the significantly non-monotonic nature of the dependencies of relative (to values in the absence of an initial stress state) values of effective coefficients of electromagnetic and magnetoelectric couplings from the volume fraction of ferrite fibers was revealed. It was revealed that we have an increase in the absolute values of the electromagnetic and magnetoelectric coupling coefficients of composite at negative values of its initial axisymmetric axial macrostrains. The most significant effect is the initial comprehensive macrostrain in the transversal plane on the coefficients of the transversal magnetoelectric coupling of the composite.

From moments of the distribution function to hydrodynamics: The nonconformal case

We study the one-dimensional boost-invariant Boltzmann equation in the relaxation-time approximation using special moments of the distribution function for a system with a finite particle mass. The infinite hierarchy of moments can be truncated by keeping only the three lowest moments that correspond to the three independent components of the energy-momentum tensor. We show that such a three-moment truncation reproduces accurately the exact solution of the kinetic equation after a simple renormalization that takes into account the effects of the neglected higher moments. We derive second-order Israel-Stewart hydrodynamic equations from the three-moment equations, and show that, for most physically relevant initial conditions, these equations yield results comparable to those of the three-moment truncation, albeit less accurate. We attribute this feature to the fact that the structure of Israel-Stewart equations is similar to that of the three-moment truncation. In particular, the presence of the relaxation term in the Israel-Stewart equations, yields an early-time regime that mimics approximately the collisionless regime. A detailed comparison of the three-moment truncation with second-order non-conformal hydrodynamics reveals ambiguities in the definition of second-order transport coefficients. These ambiguities affect the ability of Israel-Stewart hydrodynamics to reproduce results of kinetic theory.

Picosecond acoustics: a new way to access elastic properties of materials at pressure and temperature conditions of planetary interiors

Picosecond acoustics is an optical pump-probe technique allowing to access thermoelastic properties and sound velocities of a large variety of materials under extreme conditions. Coupled with diamond anvil cells and laser heating, picosecond acoustics measurements offer the possibility to probe materials over a pressure and temperature range directly pertinent for the deep planetary interiors. In this paper we highlight the capabilities and versatility of this technique by presenting some recent applications on materials of geophysical interest. All the independent components of the elastic tensor of MgO are simultaneously determined by measurements on a single crystal at ambient conditions. Compressional sound velocity is measured at high pressure on an iron-carbon alloy and on polycrystalline argon. First laser heating test measurements performed on molybdenum at high pressure are also presented. These examples demonstrate that picosecond acoustics is a valuable alternative to already existing techniques for determining the physical properties of samples under extreme pressure and temperature conditions.

DETERMINATION OF KINEMATIC CHARACTERISTICS OF METAL PLASTIC FLOW UNDER CONDITIONS OF AXIAL SYMMETRY OF THE DEFORMATION PROCESS

With axisymmetric deformation, the deformable body and the load have a common axis of symmetry. Such deformation occurs in numerous technological operations. About 70% of parts obtained by cold extrusion are deformed under conditions of axial symmetry. The plane problem of the theory of plasticity is reduced to a solution in the framework of a two-dimensional statement, when the motion of points in the section of a workpiece is analyzed. Each point can move only in the section plane, and its velocity can be decomposed into two mutually perpendicular directions along the coordinate axes. For plane deformation, it is assumed that the velocity in the direction of the third coordinate axis is equal to zero. The presence of axial symmetry allows us to confine ourselves to studying the behavior of points located on the plane of the meridional section of the workpiece. In this case, each point can move only in the section plane and its velocity can be decomposed into two orthogonal directions: along the axis and along the radius. The component of the velocity vector in the circumferential direction is equal to zero, so only four of the six independent components of the strain rate tensor remain. In this regard, axisymmetric problems of the theory of plasticity are of considerable interest from the point of view of solving applied problems. The article discusses the possibility of using mixed Euler and Lagrange coordinates to determine the components of the strain rate tensor in plastic deformation processes characterized by the axisymmetric nature of metal plastic flow. The vector field of displacements at each point in space reflects the transformation of the initial (non-deformed) configuration into the current one, and therefore determines the configuration of the deformed body in space in a certain frame of reference. The deformation process is first considered in a single Cartesian coordinate system, and then in a cylindrical coordinate system, the use of which is more appropriate for axisymmetric deformation. It is assumed that the functions of the Euler coordinates from the Lagrange coordinates are obtained by approximating the experimental data.

Second-harmonic generation in the surface layer of a dielectric spheroidal particle: II. Analysis of the solution

The vectors used in the solution of the problem of second-harmonic generation in the surface layer of a dielectric spheroidal particle are explicitly expressed in terms of the basis vectors of the spherical, cylindrical, and Cartesian coordinate systems. Three-dimensional directivity patterns characterizing the spatial distribution of the generated radiation in the far-field region and its polarization are constructed. It is found that for a small size of a spheroidal particle, the directivity pattern due to each of the non-chiral components of the nonlinear dielectric susceptibility tensor has its own individual shape. A proportional increase in the linear dimensions of the particle leads to separation of several directions of predominant radiation with a high directivity in the directivity pattern. If the exciting radiation has a linear polarization, then the generated radiation due to one (any) of the independent components of the tensor is also linearly polarized. Mathematical properties characterizing the spatial distribution of the generated radiation in the far-field region and properties associated with the change of problem parameters are found for the functions used in the solution. A relationship between the symmetries of the directivity patterns of doubled-frequency radiation and the indicated properties is revealed. The conditions under which the generation of radiation does not occur and the conditions under which the generated radiation has a linear polarization are found. The above conditions are related to the features of the spatial distribution of the generated radiation and its polarization, illustrated in the directivity patterns. Methods for estimating the independent components of the nonlinear dielectric susceptibility tensor using these conditions are proposed. Keywords: second-harmonic generation, spheroidal dielectric particle, symmetry of the spatial distribution of radiation, conditions for the absence of generation, conditions for the generation of linearly polarized radiation.

Генерация второй гармоники в поверхностном слое диэлектрической сфероидальной частицы: II. Анализ решения

The vectors used in the solution of the problem of second-harmonic generation in the surface layer of a dielectric spheroidal particle are explicitly expressed in terms of the basis vectors of the spherical, cylindrical, and Cartesian coordinate systems. Three-dimensional directivity patterns characterizing the spatial distribution of the generated radiation in the far-field region and its polarization are constructed. It is found that for a small size of a spheroidal particle, the directivity pattern due to each of the non-chiral components of the nonlinear dielectric susceptibility tensor has its own individual shape. A proportional increase in the linear dimensions of the particle leads to separation of several directions of predominant radiation with a high directivity in the directivity pattern. If the exciting radiation has a linear polarization, then the generated radiation due to one (any) of the independent components of the tensor is also linearly polarized. Mathematical properties characterizing the spatial distribution of the generated radiation in the far-field region and properties associated with the change of problem parameters are found for the functions used in the solution. A relationship between the symmetries of the directivity patterns of doubled-frequency radiation and the indicated properties is revealed. The conditions under which the generation of radiation does not occur and the conditions under which the generated radiation has a linear polarization are found. The above conditions are related to the features of the spatial distribution of the generated radiation and its polarization, illustrated in the directivity patterns. Methods for estimating the independent components of the nonlinear dielectric susceptibility tensor using these conditions are proposed.

Error Calibration for Full Tensor Magnetic Gradiometer Probe Based on Coordinate Transformation Method

In order to calibrate the zero offset error and triaxial non-orthogonal error of full tensor magnetic gradiometer (FTMG) probe composed of superconducting quantum interference device (SQUID) magnetometers, we propose a new calibration method in this article. With the help of the triaxial Helmholtz coil, we can generate a magnetic field with known parameters. By establishing the objective functions containing the error parameters for the outputs of the SQUIDs and solving them with the least-squares fitting algorithm, the included angle between the SQUID sensor and the Helmholtz coil can be obtained, then all the required error parameters can be calculated by coordinate transformation method. Compared with the traditional method, this method does not need to rotate the SQUID magnetometers-based probe to collect the data set for calibration, thus avoiding the influence of motion induced noise on the calibration results. In order to verify the effectiveness of the proposed method, we adopted the widely used ellipsoid fitting method as a comparison method and conducted field comparison experiments. The results show that the proposed calibration method is better than the traditional ellipsoid fitting calibration method under the same experimental environment. The root mean square errors of five independent magnetic gradient tensor components are reduced to 0.39, 0.35, 0.41, 0.31 and 0.57 nT/m from 1630.55, 756.59, 3097.56, 1051.09 and 1245.64 nT/m by the proposed method.

A systematic approach to reduce the independent tensor components by symmetry transformations: a commented translation of \u201cTensors and Crystal Symmetry\u201d by Carl Hermann

We here present a faithful translation of Carl Hermann’s important 1934 paper “Tensoren und Kristallsymmetrie”. This work, originally published in German language, is transferred into English while the preceding foreword summarizes Hermann’s achievements.

Addressing biases in spectral databases for increasing accuracy and computational efficiency of crystal plasticity computations

Spectral databases of crystal plasticity solutions stored as discrete Fourier transforms (DFTs) have now been used for over a decade to speed-up the crystal plasticity simulations. These approaches have demonstrated impressive computational advantages over the conventional numerical approaches that explicitly solve the governing constitutive equations. However, their implementation in finite element models (FEM) referred to as spectral crystal plasticity finite element method (SCPFEM) framework has encountered significant hurdles mainly related to instabilities in jacobian calculations leading to convergence issues in solving the governing equilibrium equations. In this work, these main hurdles encountered are tracked to certain biases introduced in the formulation of the spectral databases. These biases are related to the choices made in the selection of the independent diagonal components of the deviatoric stress tensor and the ordering of the eigenvalues of the plastic stretching tensor. Specific strategies have been designed and implemented to address these biases such as, generating separate DFT database for all three diagonal components of deviatoric stress tensor, and standardizing the assignment of eigenvalues of the plastic stretching tensor among the other improvements. These improvements to the spectral databases have been shown to improve markedly the accuracy of the database solutions for both the stress-strain responses as well as the jacobians needed in finite element computations. The benefits of the improved databases are demonstrated with selected case studies.

Bayesian optimization for inverse problems in time-dependent quantum dynamics.

We demonstrate an efficient algorithm for inverse problems in time-dependent quantum dynamics based on feedback loops between Hamiltonian parameters and the solutions of the Schrödinger equation. Our approach formulates the inverse problem as a target vector estimation problem and uses Bayesian surrogate models of the Schrödinger equation solutions to direct the optimization of feedback loops. For the surrogate models, we use Gaussian processes with vector outputs and composite kernels built by an iterative algorithm with the Bayesian information criterion (BIC) as a kernel selection metric. The outputs of the Gaussian processes are designed to model an observable simultaneously at different time instances. We show that the use of Gaussian processes with vector outputs and the BIC-directed kernel construction reduces the number of iterations in the feedback loops by, at least, a factor of 3. We also demonstrate an application of Bayesian optimization for inverse problems with noisy data. To demonstrate the algorithm, we consider the orientation and alignment of polyatomic molecules SO2 and propylene oxide (PPO) induced by strong laser pulses. We use simulated time evolutions of the orientation or alignment signals to determine the relevant components of the molecular polarizability tensors. We show that, for the five independent components of the polarizability tensor of PPO, this can be achieved with as few as 30 quantum dynamics calculations.

Global Navigation Process Simulation Based on Different Types of Gravity Data.

A theoretical study on the feasibility of global navigation based on three different types of gravity data was performed. A computer simulation of gravity-aided navigation was performed for three models of sections of the Earth’s surface with gravity anomalies distributed as specified. For navigation, three types of data sources were used, e.g., the gravity vector magnitude, three orthogonal projections of the gravity vector, and five independent components of the full gravity tensor. For each data source, when searching a specified route, the dependencies of the number of the identified true and false points were determined in accordance with the measurement error specified. The problem of determining the true route on the set of the identified points is briefly reviewed. General conclusions are presented regarding the practical applicability of the reviewed data sources to the problem of global navigation.

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

The 3D stress state from geomechanical\u2013numerical modelling and its uncertainties: a case study in the Bavarian Molasse Basin

Knowledge of the undisturbed stress state is a key parameter for borehole stability, productivity and induced seismicity hazard assessment. Due to the sparse and incomplete availability of data records, 3D geomechanical–numerical modelling is applied to estimate the stress state in a volume. To assess the quality of the model results, the model uncertainties have to be quantified. We present an approach that provides the uncertainties of the six independent stress tensor components at each location in the model volume by an average value and its standard deviation. We explore the uncertainties introduced to the model results with respect to the stress magnitude data used for model calibration. We test our approach on a model of the Bavarian Molasse Basin which includes the area around Munich with many geothermal projects. In the test area, we find large uncertainties in the modelled magnitude of the maximum horizontal stress (text{S}_text {Hmax}) in the order of 15–30%. The uncertainties in the magnitude of the minimum horizontal stress (text{S}_text {hmin}) are smaller between 5 and 20%. In connection with an adequate failure criterion, we compare our model results to the seismological observations at two neighbouring geothermal projects. Whilst Aschheim/Feldkirchen/Kirchheim remained seismically quiet, induced seismicity was recorded in Poing. Our modelled undisturbed stress state was confirmed by these observations by showing in average a stable stress state in Aschheim/Feldkirchen/Kirchheim and a critical stress state in Poing. If the uncertainties in terms of two standard deviations are considered, this does not change. This demonstrates the increase in model significance when uncertainties are available.

Higher order elastic tensors of crystal structure under non-linear deformation

The higher-order elastic tensors can be used to characterize the linear and non-linear mechanical properties of crystals at ultra-high pressures. Besides the widely studied second-order elastic constants, the third- and fourth-order elastic constants are sixth and eighth tensors, respectively. The independent tensor components of them are completely determined by the symmetry of the crystal. Using the relations between elastic constants and sound velocity in solid, the independent elastic constants can be measured experimentally. The anisotropy in elasticity of crystal structures is directly determined by the independent elastic constants.

Single-crystal elasticity of iron-bearing phase E and seismic detection of water in Earth's upper mantle

AbstractThe elastic properties of Mg2.12(2)Fe0.21(2)Ni0.01Si1.15(1)O6H2.67(8) phase E single crystals with Fe3+/ΣFe = 0.25(3) have been determined by Brillouin spectroscopy at ambient conditions. We find that that the elasticity of iron-bearing phase E is described by the six independent stiffness tensor components (all in units of GPa): C11 = 192.2(6), C12 = 56.4(8), C13 = 43.5(8), C14 = –4.3(3), C33 = 192.1(7), C44 = 46.4(3). The Voigt-Reuss-Hill averages of bulk and shear moduli are 95.9(4) and 59.6(2) GPa, respectively. The aggregate velocities of iron-bearing phase E are νP = 7.60(2) and νS = 4.43(1) km/s, markedly lower than those of major mantle minerals at ambient conditions. Modeling based on our results suggests that the presence of iron-bearing phase E may reduce the sound wave velocities in upper mantle and transition zone rocks, making it a possible target for future seismological investigations aiming to map hydration in subducting slabs.

A method of using geomagnetic anomaly to recognize objects based on HOG and 2D-AVMD

In order to identify the shape of underground small magnetic anomaly objects, we use Support Vector Machines (SVM) to identify the underground magnetic anomaly targets. Firstly, as the SVM needs a lot of training data, and we also need to make full use of the magnetic field signal, nine component signals including total magnetic intensity (TMI) and five independent components of tensor are calculated from the original detected magnetic signal. Secondly, the nine component signals are subjected respectively to two-dimensional adaptively variational mode decomposition (2D-AVMD), which is advanced based on the two indicators, namely Mutual information (MI) and empirical entropy (EE), and we can get the nine primary signals from the decomposition results of nine component signals called the Intrinsic Mode Function (IMF). Then, the Histogram of Oriented Gradients (HOG) of the nine primary signals is extracted, and the feature data would be constructed into feature vectors. In the end, Support Vector Machines (SVM) are adopted to process these feature vectors. The output of the SVM can indicate the result of small objects’ shape recognition under the ground. Experiments prove that the shape recognition accuracy of underground small magnetic anomaly object recognition reaches 90%.