Dynamics of a driven magneto-martensitic ribbon

We investigate the dynamics of a sinusoidally driven ferromagnetic martensitic ribbon by adopting a recently introduced model that involves strain and magnetization as order parameters. Retaining only the dominant mode of excitation we reduce the coupled set of partial differential equations for strain and magnetization to a set of coupled ordinary nonlinear equations for the strain and magnetization amplitudes. The equation for the strain amplitude takes the form of parametrically driven oscillator. Finite strain amplitude can only be induced beyond a critical value of the strength of the magnetic field. Chaotic response is seen for a range of values of all the physically interesting parameters. The nature of the bifurcations depends on the choice of temperature relative to the ordering of the Curie and the martensite transformation temperatures. We have studied the nature of response as a function of the strength and frequency of the magnetic field, and magneto-elastic coupling. In general, the bifurcation diagrams with respect to these parameters do not follow any standard route. The rich dynamics exhibited by the model is further illustrated by the presence of mixed mode oscillations seen for low frequencies. The geometric structure of the mixed mode oscillations in the phase space has an unusual deep crater structure with an outer and inner cone on which the orbits circulate. We suggest that these features should be seen in experiments on driven magneto-martensitic ribbons.     

On the dynamics of deformable ferromagnets

The paper begins with a fairly detailed presentation of a general mathematical model for the dynamics of ferromagnetic bodies undergoing arbitrarily large deformations. Next, a mathematical study is presented of the evolution problem for the magnetization field in a «soft» ferromagnetic body which is «mechanically at rest». No matter how special and simple this problem within the frame-work of the full theory, the governing equation is interesting: it is identical to the dynamic version of the harmonic-map equation usually referred to in the mathematical literature as the Gilbert form of the Landau-Lifshitz equation. Motivated by recent nonuniqueness results for the dynamic harmonicmap equation, we give a new proof of global existence of weak solutions to the Gilbert-Landau-Lifshitz equation.     

The two-time Green’s function and the diagram technique

Based on constructing the equations of motion for the two-time Green’s functions, we discuss calculating the dynamical spin susceptibility and correlation functions in the Heisenberg model. Using a Mori-type projection, we derive an exact Dyson equation with the self-energy operator in the form of a multiparticle Green’s function. Calculating the self-energy operator in the mode-coupling approximation in the ferromagnetic phase, we reproduce the results of the temperature diagram technique, including the correct formula for low-temperature magnetization. We also consider calculating the spin fluctuation spectrum in the paramagnetic phase in the framework of the method of equations of motion for the relaxation function.     

Determination of precession and dissipation parameters in micromagnetism

The precession β and the dissipation parameter α of a ferromagnetic material can be considered microscopically space dependent. Their space distribution is difficult to obtain by direct measurements. In this article we consider an inverse problem, where we aim at recovering α and β from space measurements of the magnetization. The evolution of the magnetization in micromagnetism is governed by the Landau-Lifshitz (LL) equation. We first study the sensitivity of the LL equation. We derive the existence, uniqueness and stability results for the LL equation and the corresponding sensitivity equations. On the basis of the results we analyze the inverse problem. We employ the energy method and we minimize the underlying cost functional by means of the steepest descent method. We derive a convergence result for the proposed algorithm. The presented numerical examples support the theoretical results.     

A Magneto-Mechanically Coupled Phase-Field Model at Large Deformation

The aggregate magneto-mechanical behavior of magneto-rheological elastomers (MREs) stems from the magnetic properties of the ferromagnetic inclusion and the mechanical properties of the matrix material. We propose a large deformation micro-magnetic theory, to predict the behavior and interaction of ferromagnetic particles inside an elastomeric matrix. A rate-type variational principle, with the magnetization as the order parameter is proposed. A large deformation Landau-Lifshitz-Gilbert equation for the time evolution of the magnetization, is obtained directly from the proposed variational principle. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Band Model and Ferromagnetism

The subject is reviewed from an elementary point of view. Topicstreated include: band structure of ferromagnetic metals, theband model at 0°K, magnetization at low temperatures andat higher temperatures, the Curie point, the susceptibilityabove T_c, interactions and spin waves in metals. The literaturecited goes up to and includes 1966.     

Heat Conduction in a Phase Field Model for Martensitic Transformation

We study the martensitic transformation with a phase field model, where we consider the Bain transformation path in a small strain setting. For the order parameter, interpolating between an austenitic parent phase and martensitic phases, we use a Ginzburg-Landau evolution equation, assuming a constant mobility. In [1], a temperature dependent separation potential is introduced. We use this potential to extend the model in [2], by considering a transient temperature field, where the temperature is introduced as an additional degree of freedom. This leads to a coupling of both the evolution equation of the order parameter and the mechanical field equations (in terms of thermal expansion) with the heat equation. The model is implemented in FEAP as a 4-node element with bi-linear shape functions. Numerical examples are given to illustrate the influence of the temperature on the evolution of the martensitic phase. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Existence and uniqueness for a mathematical model in superfluidity

In this paper we propose a model to study superfluidity by considering as state variables the order parameter, describing the concentration of the superfluid phase, the velocity of the superfluid and the absolute temperature. We assume that the order parameter satisfies a Ginzburg–Landau equation and that the velocity is decomposed as the sum of a normal and a superfluid component. The heat equation provides the evolution equation for the temperature. We prove that this model is consistent with the principles of thermodynamics. Well‐posedness of the resulting initial and boundary value problem is shown. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical methods for the modeling of the magnetization vector in multiferroic heterostructures

Multiferroic heterostructures are commonly used to obtain electro-magnetic coupling effects. Thereby, the ferroelectric layer is used to control the magnetization in the ferromagnetic layer. The coupling between the layers is obtained by the mechanical coupling between the layers, which have well-defined interfaces. Within this contribution we use phase field models to define the polarization and magnetization in the ferroelectric and ferromagnetic layers, respectively. A coupling between polarization/magnetization and strains in each layer in combination with coherent deformations at the interface yields an electromagnetic coupling within the entire heterostructure. Numerical formulations for the interpolation of the polarization vector are well-defined in the literature. However, the establishment of a consistent numerical formulation for the ferromagnetic layer, where the length of the magnetization vector has to be constant, remains a difficult task. We propose a new numerical approach for the consistent treatment of the ferromagnetic layer and provide numerical simulations which illustrate the electromagnetic coupling effect. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of dust particles on a layer of micropolar ferromagnetic fluid heated from below saturating a porous medium

The problem of the effect of dust particles on the thermal convection in micropolar ferromagnetic fluid saturating a porous medium subject to a transverse uniform magnetic field has been investigated theoretically. Linear stability analysis and normal mode analysis methods are used to find an exact solution for a flat micropolar ferromagnetic fluid layer contained between two free boundaries. In case of stationary convection, the effect of various parameters like medium permeability, dust particles, non-buoyancy magnetization, coupling parameter, spin-diffusion parameter and micropolar heat conduction parameter are analyzed. For sufficiently large values of magnetic parameter M₁, the critical magnetic thermal Rayleigh number for the onset of instability is determined numerically and results are depicted graphically. It is also observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the micropolar ferromagnetic fluid saturating a porous medium heated from below in the absence of micropolar viscous effect, microinertia and dust particles.     

A dual approach to regularity in thin film micromagnetics

We prove a regularity result in the two-dimensional theory of soft ferromagnetic films. The associated Euler–Lagrange equation is given by a nonlocal degenerate variational inequality involving fractional derivatives. A difference quotient type argument based on a dual formulation in terms of magnetostatic potentials yields a Hölder estimate for the uniquely determined gradient projection of the magnetization field.     

Effect of magnetic field dependent viscosity on ferroconvection in the presence of dust particles

This paper deals with the theoretical investigation of the effect of magnetic field dependent (MFD) viscosity on the thermal convection in a ferromagnetic fluid in the presence of dust particles. For a flat ferromagnetic fluid layer contained between two free boundaries, the exact solution is obtained using a linear stability analysis and a normal mode analysis method. For the case of stationary convection, dust particles always have a destabilizing effect, whereas the MFD viscosity has a stabilizing effect on the onset of convection. In the absence of MFD viscosity, the destabilizing effect of magnetization is depicted but in the presence of MFD viscosity, non-buoyancy magnetization may have a destabilizing or a stabilizing effect on the onset of convection. The critical wave number and critical magnetic thermal Rayleigh number for the onset of stationary convection are also determined numerically for sufficiently large values of buoyancy magnetization parameter M ₁. Graphs have been plotted by giving numerical values to the parameters to depict the stability characteristics. It is observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the ferromagnetic fluid heated from below in the absence of dust particles. The oscillatory modes are introduced due to the presence of the dust particles, which were non-existent in their absence. A sufficient condition for the non-existence of overstability is also obtained.     

Self-consistent approximation in the Ising model of pure and dilute magnets using a pair correlation

Theoretical and Mathematical Physics - We construct a self-consistent approximation for calculating the magnetization and Curie temperature in the Ising model on an arbitrary lattice based on using...     

Martingale weak solutions of the stochastic Landau–Lifshitz–Bloch equation

The present paper studies the stochastic Landau–Lifshitz–Bloch equation which is recommended as the only valid model at temperature around the Curie temperature and is especially important for the simulation of heat-assisted magnetic recording. We study the stochastic Landau–Lifshitz–Bloch equation in the case that the temperature is raised higher than the Curie temperature. The global existence of martingale weak solutions is proved by using a new argument and regularity properties of the weak solutions are discussed.     

A Second-Order Semi-Implicit Method for the Inertial Landau-Lifshitz-Gilbert Equation

Electron spins in magnetic materials have preferred orientations collectively and generate the macroscopic magnetization. Its dynamics spans over a wide range of timescales from femtosecond to picosecond, and then to nanosecond. The Landau-Lifshitz-Gilbert (LLG) equation has been widely used in micromagnetics simulations over decades. Recent theoretical and experimental advances have shown that the inertia of magnetization emerges at sub-picosecond timescales and contributes significantly to the ultrafast magnetization dynamics, which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial LLG (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves as a hyperbolic system at sub-picosecond timescales, while behaves as a parabolic system at larger timescales spanning from picosecond to nanosecond. Such hybrid behaviors impose additional difficulties on designing efficient numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second-order temporal derivative of magnetization is approximated by the standard centered difference scheme, and the first-order temporal derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with second-order accuracy. At each time step, the unconditionallyunique solvability of the unsymmetric linear system is proved with detailed discussions on the condition number. Numerically, the second-order accuracy of the proposed method in both time and space is verified. At sub-picosecond timescales, the inertial effect of ferromagnetics is observed in micromagnetics simulations, in consistency with the hyperbolic property of the iLLG model; at nanosecond timescales, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model.     

A stereographic projection based phase field model for ferromagnetics

By using the stereographic projection to satisfy the constant magnetization magnitude, a constraint-free phase field model is developed for ferromagnetic materials. Implemented by finite element method, the model is shown to be capable of simulating magnetic domain, vortex, and hysteresis curves. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A thermodynamically consistent Ginzburg–Landau model for superfluid transition in liquid helium

In this paper, we propose a thermodynamically consistent model for superfluid-normal phase transition in liquid helium, accounting for variations of temperature and density. The phase transition is described by means of an order parameter, according to the Ginzburg–Landau theory, emphasizing the analogies between superfluidity and superconductivity. The normal component of the velocity is assumed to be compressible, and the usual phase diagram of liquid helium is recovered. Moreover, the continuity equation leads to a dependence between density and temperature in agreement with the experimental data.     

Computational Homogenization in Micro-Magneto-Elasticity

The overall macroscopic response of magneto-mechanically coupled materials stems from complex magnetization evolution and corresponding domain wall motion occurring on a lower length scale. In order to account for such effects we propose a computational homogenization approach that incorporates a ferromagnetic phase-field formulation into a macroscopic Boltzmann continuum. This scale-bridging is obtained by rigorous definition of rate-type and incremental variational principles. An extended version of the classical Hill-Mandel macro-homogeneity condition is obtained as a consequence. In order to satisfy the unity constraint of the magnetization on the micro-scale, an efficient operator-split method is proposed. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Large-Deformation Phase-Field Approach for the Modeling of Magneto-Sensitive Elastomers

Magneto-sensitive materials show magneto-mechanical coupled response and are thus of increasing interest in the recent age of smart functional materials. Ferromagnetic particles suspended in an elastomeric matrix show realignment under the influence of an external applied field, in turn causing large deformations of the substrate material. The magneto-mechanical coupling in this case is governed by the magnetic properties of the inclusion and the mechancial properties of the matrix. The magnetic phenomenon in ferromagnetic materials is governed by the formation and evolution of domains on the micro scale. A better understanding of the behavior of these particles under the influence of an external applied field is required to accurately predict the behavior of such materials. In this context it is of particular importance to model the macro scopic magneto-mechanically coupled behavior based on the micro-magnetic domain evolution. The key aspect of this work is to develop a large-deformation micro-magnetic model that can accurately capture the microscopic response of such materials. Rigorous exploitation of appropriate rate-type variational principles and consequent incremental variational principles directly give us field equations including the time evolution equation of the magnetization, which acts as the order parameter in our formulation. The theory presented here is the continuation of the work presented in [1, 7] for small deformations. A summary of magneto-mechanical theories spanning over multiple scales has been presented in [4]. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

WAVES IN FERROMAGNETIC MEDIA

ABSTRACT

It is shown that small perturbations of equilibrium states in ferromagnetic media give rise to standing and traveling waves that are stable for long times. The evolution of the wave profiles is governed by semilinear heat equations. The mathematical model underlying these results consists of the Landau–Lifshitz equation for the magnetization vector and Maxwell's equations for the electromagnetic field variables. The model belongs to a general class of hyperbolic equations for vector-valued functions, whose asymptotic properties are analyzed rigorously. The results are illustrated with numerical examples.