Isotropic Function Research Articles

Homogenization of quasi-periodic conformal architectured materials and applications to chiral lattices

In this study, we propose to extend asymptotic periodic homogenization for non-periodic continuous microstructured media, assuming that the non-periodic geometry (called quasi-periodic) can be designed by a conformal planar transformation of a periodic parent domain architectured media with periodically disposed unit cells. Conformal transformations are shown to play a privileged role in the design of circular macroscopic heterogeneous domains tessellated with non-periodic unit cells, obtained from a periodic parent domain architectured with these unit cells. The conditions for conformal invariance are established, leading to the general form of conformal transformation in their dependencies upon the periodic coordinates. It is shown that any conformal map can be decomposed into the product of an isotropic dilatation function of the first periodic spatial position of decreasing exponential type and a rotation characterized by an angular function linear in the second periodic position. A general theory of quasi-periodic homogenization in the framework of conformal transformations is established for the first time, leading to an expression of the tensor of quasi-periodic moduli which is fully evaluated from the solution of the elasticity boundary value problem posed over the periodic unit cell. The influence of microcurvature distortion of individual unit cells on their effective properties is evaluated. Closed-form solutions are confronted to numerical examples issued from the implementation of circular periodicity in a finite element solver, showing overall a good agreement with the identified homogenized moduli.

Earing Prediction of AA5754-H111 (Al-Alloy) with Linear Transformation-Based Anisotropic Drucker Yield Function under Non-Associated Flow Rule (Non-AFR).

The yield behavior of aluminum alloy 5754-H111 under different stress conditions for three kinds of plastic work is studied using an anisotropic Drucker model. It is found that when the plastic work is 30 MPa, the anisotropic Drucker model has the most accurate prediction. Comparing the Hill48 and Yld91 models with the Drucker model, the results show that both the anisotropic Drucker and Yld91 models can accurately predict the yield behavior of the alloy. Cylinder drawing finite element analysis is performed under the AFR, but it is not possible to accurately predict the position and height of earing appearance. The anisotropic Drucker model is used to predict the earing behavior under the non-AFR, which can accurately predict the earing phenomenon. Numerical simulation is conducted using three different combinations of yield functions: the anisotropic yield function and the anisotropic plastic potential function (AYAPP), the anisotropic yield function and the isotropic plastic potential function (AYIPP), and the isotropic yield function and the anisotropic plastic potential function (IYAPP). It is concluded that the influence of the plastic potential function on predicting earing behavior is more critical than that of the yield function.

A model for capturing the rate-dependent mechanical behaviour of liquid crystal elastomers

A modelling framework is adopted and adapted here to model the rate-dependent deformation behaviour of liquid crystal elastomers (LCEs). On using an isotropic hyperelastic strain energy function WI1,I2 of binomial form, previously employed for capturing the (quasi-static) deformation of polydomain LCEs, an extension is presented for application to the rate-dependent behaviour of LCEs by incorporating the deformation rate into the model parameters. That is, the model parameters of the core hyperelastic strain energy function W are considered to evolve with, i.e., be a function of, the deformation rate. A specific measure of deformation rate ɛ̇ is utilised, and a particular functional dependency of the hyperelastic model parameters on ɛ̇, of a skewed exponential form, is devised and presented. The ensuing rate-dependent model Wr is then applied to a range of extant experimental datasets, encompassing various LCE specimens and across different deformation rates, under uniaxial loading and including both tension and compression deformations. The model is shown to capture and predict the data favourably. Given the simplic of devising and implementing the presented modelling framework, the flexibility it provides for incorporating different choices of strain energy W function as well as the functional form(s) of the dependency of the model parameters on the deformation rate, and the accuracy of the modelling results, the application of this modelling approach to capturing the rate-dependent deformation behaviour of LCEs is hereby presented.

Anisotropy of Magnetohydrodynamic and Kinetic Scale Fluctuations through Correlation Tensor in Solar Wind at 0.8 au

Space plasma turbulence is inherently characterized by anisotropic fluctuations. The generalized k-th order correlation tensor of magnetic field increments allow us to separate the mixed isotropic and anisotropic structure functions from the purely anisotropic ones. In this work, we quantified the relative importance of anisotropic fluctuations in solar wind turbulence using two Alfvénic data samples gathered by the Solar Orbiter at 0.8 astronomical units. The results based on the joined statistics suggest that the anisotropic fluctuations are ubiquitous in solar wind turbulence and persist at kinetic scales. Using the RTN coordinate system, we show that their presence depends on the anisotropic sector under consideration, e.g., the RN and RT sectors exhibit enhanced anisotropy toward kinetic scales, in contrast with the TN. We then study magnetic field fluctuations parallel and perpendicular to the local mean magnetic field separately. We find that perpendicular fluctuations are representative of the global statistics, resembling the typical picture of magnetohydrodynamic turbulence, whereas parallel fluctuations exhibit a scaling law with slope ∼1 for all the joined isotropic and anisotropic components. These results are in agreement with predictions based on the critical balance phenomenology. This topic is potentially of interest for future space missions measuring kinetic and MHD scales simultaneously in a multi-spacecraft configuration.

Self-consistent dynamical models with a finite extent – IV. Wendland models based on compactly supported radial basis functions

ABSTRACT We present a new step in our systematic effort to develop self-consistent dynamical models with a finite radial extent. The focus is on models with simple analytical density profiles allowing for analytical calculations of many dynamical properties. In this paper, we introduce a family of models, termed Wendland models, based on compactly supported radial basis functions. The family of models is characterized by a parameter k that controls the smoothness of the transition at the truncation radius. In the limit $k\rightarrow \infty$, the Wendland model reduces to a non-truncated model with a Gaussian density profile. For each Wendland model, the density, mass and gravitational potential are simple truncated polynomial functions of radius. Via the sphecow tool we demonstrate that all Wendland models can be supported by isotropic distribution functions. Surprisingly, the isotropic distribution function exhibits varied behaviour across different Wendland models. Additionally, each model can be supported by a continuum of Osipkov–Merritt orbital structures, ranging from radially anisotropic to completely tangential at the truncation radius. To the best of our knowledge, the Wendland models presented here are the first family of models accommodating both radial and tangential Osipkov–Merritt distribution functions. Using linear superposition, these models can easily be combined to generate Wendland models with even more diverse orbital structures. While the Wendland models are not fully representative of real dynamical systems due to their Gaussian-like density profile, this study lays important groundwork for constructing more realistic models with truncated density profiles that can be supported by a range of orbital structures.

Anisotropic remixing of a phase separated binary colloidal system with particles of different sizes in an external modulation.

We explore the phase behavior of a binary colloidal system under external spatially periodic modulation. We perform Monte Carlo simulations on a binary mixture of big and small repulsive Lennard-Jones particles with a diameter ratio of 2:1. We characterize structure by isotropic and anisotropic pair correlation functions, cluster size distribution, bond angle distribution, order parameter, and specific heat. We observe the demixing of the species in the absence of external modulation. However, the mixing of the species gets enhanced with increasing potential strength along with the alignment of the particles transverse to the modulation. The de-mixing order parameter shows discontinuity with increasing modulation strength, characterizing a first order phase transition. The peak in specific heat increases linearly with the size of the system. We also look into the dynamical behavior of the system via computing Mean Square Displacement (MSD) along both parallel and perpendicular directions to the modulation. We observe a decrease in the diffusion coefficient for both types of particles as we increase the strength of the modulation.

Statistical Modeling of Right-Censored Spatial Data Using Gaussian Random Fields

Consider a fixed number of clustered areas identified by their geographical coordinates that are monitored for the occurrences of an event such as a pandemic, epidemic, or migration. Data collected on units at all areas include covariates and environmental factors. We apply a probit transformation to the time to event and embed an isotropic spatial correlation function into our models for better modeling as compared to existing methodologies that use frailty or copula. Composite likelihood technique is employed for the construction of a multivariate Gaussian random field that preserves the spatial correlation function. The data are analyzed using counting process and geostatistical formulation that led to a class of weighted pairwise semiparametric estimating functions. The estimators of model parameters are shown to be consistent and asymptotically normally distributed under infill-type asymptotic spatial statistics. Detailed small sample numerical studies that are in agreement with theoretical results are provided. The foregoing procedures are applied to the leukemia survival data in Northeast England. A comparison to existing methodologies provides improvement.

Superconductivity properties of Pd-Heusler rare earth from ab initio and isotropic Eliashberg function

Ab initio was employed to investigate the structural, electronic, and superconducting properties of the Heusler alloys Pd2TmSn, Pd2YbSn, and Pd2LuSn. The results reveal that these alloys exhibit both structural and dynamic stability in a cubic face-centered lattice. The isotropic Eliashberg function was used to calculate the superconductor properties. The phonon–electron couplings fall within the range characteristic of medium-coupling superconductors. Analysis of the Fermi surfaces identifies the presence of necks connecting the first Brillouin zones, which facilitate electron pairing. The calculated superconducting critical temperatures using the Allen–Dynes equation are in agreement with previous experimental and theoretical findings.

Cox processes driven by transformed Gaussian processes on linear networks—A review and new contributions

Cox processes driven by transformed Gaussian processes on linear networks—A review and new contributions

On Itskov (2024) counterexample to the functional basis of isotropic vector and tensor functions by Shariff (2023)

On Itskov (2024) counterexample to the functional basis of isotropic vector and tensor functions by Shariff (2023)

On the functional basis of isotropic vector and tensor functions by Shariff (2023)

In the paper Shariff (Q. J. Mech. Appl. Math. 76:143–161, 2023) a functional basis of a system of vectors and symmetric tensors is proposed. The functional basis is expressed in terms of eigenvalues and eigenvectors of the first tensor and includes a smaller number of terms in comparison to the classical irreducible representation (see, e.g., Boehler, J. Appl. Math. Mech. 57:323–327, 1977; Pennisi and Trovato, Int. J. Eng. Sci. 25:1059–1065, 1987). In the present contribution, we show that elements of the functional basis by Shariff (Q. J. Mech. Appl. Math. 76:143–161, 2023) do not represent isotropic invariants of the vector and tensor arguments and cannot thus be referred to as the functional basis. To this end, a counterexample with two symmetric tensors is considered. Under an arbitrary orthogonal transformation the functional basis (Shariff, Q. J. Mech. Appl. Math. 76:143–161, 2023) of these two tensors should remain constant but it does change in contrast to the classical representation.

Monte Carlo simulation of spatial resolution of lens-coupled LYSO scintillator for intense pulsed gamma-ray imaging system with large field of view

In this paper, we use a Monte Carlo (MC) simulation based on Geant4 to investigate the influence of four parameters on the spatial resolution of the lens-coupled lutetium yttrium orthosilicate (LYSO) scintillator, including the thickness of the LYSO scintillator, the F-number and minification factor of the lens, and the incident position of the gamma-rays. Simulation results show that when the gamma-rays are incident along the lens axis, the smaller the thickness, the larger the F-number, the larger the minification factor, the higher the spatial resolution, with an isotropic point spread function (PSF). As the incident position of the gamma-rays deviates from the lens axis, the spatial resolution decreases, and the PSF becomes anisotropic. In addition, by analyzing the whole physical process of the lens-coupled LYSO scintillator from gamma-rays to secondary electrons to fluorescence photons, we aim to provide a detailed analysis of the influence of each parameter on the spatial resolution. The results show that the PSF of the secondary electrons energy deposition is almost constant in the simulation, which determines the upper limit of the spatial resolution. Meanwhile, the dispersion process of the fluorescence photons can explain the reason why each parameter affects the spatial resolution.

An anisotropic finite strain elastoplastic model considering different plastic spin effects on the intermediate configuration

Large deformation anisotropic elastoplastic models have important application background and research value in engineering and materials fields. In this paper, a macroscopic phenomenological model of elastoplastic anisotropic deformation considering plastic spin is proposed. The model is based on the multiplication decomposition of deformation gradient. The free energy function is expressed as the isotropic function of the strain and the structural tensor on the intermediate configuration, which is a push-forward of initial configuration by using plastic deformation gradient. Advantageously, the non-equilibrium free energy function remains invariant under the superimposed rigid body rotation on the intermediate configuration, which is due to the non-uniqueness of the multiplicative decomposition of deformation gradient. However, the rate of superimposed rigid body rotation has effects on the model. The effects are discussed by considering three different spin assumptions. Numerical simulations show that the plastic spin has influence on the calculation results. Therefore, the plastic spin assumption should be carefully selected in the practical application of elastoplastic anisotropic model.

Expansion and Collapse of Spherically Symmetric Isotropic Elastic Bodies Surrounded by Vacuum

A class of isotropic and scale-invariant strain energy functions is given for which the corresponding spherically symmetric elastic motion includes bodies whose diameter becomes infinite with time or collapses to zero in finite time, depending on the sign of the residual pressure. The bodies are surrounded by vacuum so that the boundary surface forces vanish, while the density remains strictly positive. The body is subject only to internal elastic stress.

A novel and fully automated coordinate system transformation scheme for near optimal surrogate construction

This work develops a novel coordinate system transformation scheme to improve the performance of common radial basis function surrogate models. This coordinate system transformation scheme is based on the fact that commonly used basis functions are isotropic.Three main empirical findings are established in this study. Firstly, in general isotropic functions are inadequate to describe anisotropic data manifolds due to a mismatch between the functional form and the form of the data manifold resulting in poor generative performance. Counter-intuitively, utilising additional gradients during surrogate training often worsens the generative capability.Secondly, component-wise scaling of isotropic model forms during cross-validation is inadequate to enhance the functional form of the data manifold form as anisotropic coupling in the data manifold remains coupled. Improving the match between the functional form and the data manifold form requires both rotation and scaling.Thirdly, the coordinate system transformation scheme should predominantly be based on a collection of local curvature estimations and not on global curvature approximations. Gradients are critical to estimating the local curvature for identifying a near-optimal reference frame for surrogate construction, which then translates to additional benefits of gradients in gradient-enhanced surrogates.Based on the above observations, this paper proposes an isotropic transformation for the data coordinate system that performs near-optimal transformations on lower dimensional data without requiring any cross-validation. The method is compared against commonly applied component-wise cross-validation data coordinate system scaling as well as the more modern Active Subspace Method on a carefully crafted decomposable test problem, which has a known optimal coordinate system, that varies between 2 and 16 dimensions.The paper concludes after demonstrating that the developed transformation scheme, as well as the other common methods, will offer little benefit on non-decompose problems and offers some suggestions on future work to create a more general isotropic transformation.

Matrix-valued isotropic covariance functions with local extrema

Multivariate random fields are commonly used in spatial statistics and natural science to model coregionalized variables. In this context, the matrix-valued covariance function plays a central role in capturing their spatial continuity and interdependence. This study aims to contribute to the literature on covariance modeling by proposing new parametric families of isotropic matrix-valued functions exhibiting non-monotonic behaviors, namely hole effects and cross-dimples. The benefit of the proposed models is shown on a bivariate data set consisting of concentrations of airborne particulate matter.

Histology-informed multiscale modeling of human brain white matter

In this study, we propose a novel micromechanical model for the brain white matter, which is described as a heterogeneous material with a complex network of axon fibers embedded in a soft ground matrix. We developed this model in the framework of RVE-based multiscale theories in combination with the finite element method and the embedded element technique for embedding the fibers. Microstructural features such as axon diameter, orientation and tortuosity are incorporated into the model through distributions derived from histological data. The constitutive law of both the fibers and the matrix is described by isotropic one-term Ogden functions. The hyperelastic response of the tissue is derived by homogenizing the microscopic stress fields with multiscale boundary conditions to ensure kinematic compatibility. The macroscale homogenized stress is employed in an inverse parameter identification procedure to determine the hyperelastic constants of axons and ground matrix, based on experiments on human corpus callosum. Our results demonstrate the fundamental effect of axon tortuosity on the mechanical behavior of the brain’s white matter. By combining histological information with the multiscale theory, the proposed framework can substantially contribute to the understanding of mechanotransduction phenomena, shed light on the biomechanics of a healthy brain, and potentially provide insights into neurodegenerative processes.

Spherical harmonic coefficients of isotropic polynomial functions with applications to gravity field modeling

Various aspects of gravity field modeling rely upon analytical mathematical functions for calculating spherical harmonic coefficients. Such functions allow quick and efficient evaluation of cumbersome convolution integrals defined on the sphere. In this work, we present a new analytical method for determining spherical harmonic coefficients of isotropic polynomial functions. This method in computationally flexible and efficient, since it makes use of recurrence relations. Also, its use is universal and could be extended to piecewise polynomials and polynomials with compact support. Our numerical investigation of the proposed method shows that certain recurrence relations lose accuracy as the order of implemented polynomials increases because of accumulation of numerical errors. Propagation of these errors could be mitigated by hybrid methods or using extended precision arithmetic. We demonstrate the relevance of our method in gravity field modeling and discuss two areas of application. The first one is the design of B-spline windows and filter kernels for the low-pass filtering of gravity field functionals (e.g., GRACE Follow-On monthly geopotential solutions). The second one is the calculation of spherical harmonic coefficients of isotropic polynomial covariance functions.

Polarized Anisotropic Synchrotron Emission and Absorption and Its Application to Black Hole Imaging

Low-collisionality plasma in a magnetic field generically develops anisotropy in its distribution function with respect to the magnetic field direction. Motivated by the application to radiation from accretion flows and jets, we explore the effect of temperature anisotropy on synchrotron emission. We derive analytically and provide numerical fits for the polarized synchrotron emission and absorption coefficients for a relativistic bi-Maxwellian plasma (we do not consider Faraday conversion/rotation). Temperature anisotropy can significantly change how the synchrotron emission and absorption coefficients depend on observing angle with respect to the magnetic field. The emitted linear polarization fraction does not depend strongly on anisotropy, while the emitted circular polarization does. We apply our results to black hole imaging of Sgr A* and M87* by ray tracing a GRMHD simulation and assuming that the plasma temperature anisotropy is set by the thresholds of kinetic-scale anisotropy-driven instabilities. We find that the azimuthal asymmetry of the 230 GHz images can change by up to a factor of 3, accentuating (T ⊥ > T ∥) or counteracting (T ⊥ < T ∥) the image asymmetry produced by Doppler beaming. This can change the physical inferences from observations relative to models with an isotropic distribution function, e.g., by allowing for larger inclination between the line of sight and spin direction in Sgr A*. The observed image diameter and the size of the black hole shadow can also vary significantly due to plasma temperature anisotropy. We describe how the anisotropy of the plasma can affect future multifrequency and photon ring observations. We also calculate kinetic anisotropy-driven instabilities (mirror, whistler, and firehose) for relativistically hot plasmas.

Controllable Deformations of Unconstrained Ideal Nematic Elastomers

We establish that, for ideal unconstrained uniaxial nematic elastomers described by a homogeneous isotropic strain-energy density function, the only smooth deformations that can be controlled by the application of surface tractions only and are universal in the sense that they are independent of the strain-energy density are those for which the deformation gradient is constant and the liquid crystal director is either aligned uniformly or oriented randomly in Cartesian coordinates. This result generalizes the classical Ericksen’s theorem for nonlinear homogeneous isotropic hyperelastic materials. While Ericksen’s theorem is directly applicable to liquid crystal elastomers in an isotropic phase where the director is oriented randomly, in a nematic phase, the constitutive strain-energy density must account also for the liquid crystal orientation which leads to significant differences in the analysis compared to the purely elastic counterpart.