Fast computation of magnetostatic fields by Non-uniform Fast Fourier Transforms

The bottleneck of micromagnetic simulations is the computation of the long-ranged magnetostatic fields. This can be tackled on regular N-node grids with Fast Fourier Transforms in time N logN, whereas the geometrically more versatile finite element methods (FEM) are bounded to N^4/3 in the best case. We report the implementation of a Non-uniform Fast Fourier Transform algorithm which brings a N logN convergence to FEM, with no loss of accuracy in the results

A Quasi-Hamiltonian Discretization of the Thermal Shallow Water Equations

International audienceThe rotating shallow water (RSW) equations are the usual testbed for the development of numerical methods for three-dimensional atmospheric and oceanic models. However, an arguably more useful set of equations are the thermal shallow water equations (TSW), which introduce an additional thermodynamic scalar but retain the single layer, two-dimensional structure of the RSW. As a stepping stone towards a three-dimensional atmospheric dynamical core, this work presents a quasi-Hamiltonian discretization of the thermal shallow water equations using compatible Galerkin methods, building on previous work done for the shallow water equations. Structure-preserving or quasi-Hamiltonian discretizations methods, that discretize the Hamiltonian structure of the equations of motion rather than the equations of motion themselves, have proven to be a powerful tool for the development of models with discrete conservation properties. By combining these ideas with an energy-conserving Poisson time integrator and a careful choice of Galerkin spaces, a large set of desirable properties can be achieved. In particular, for the first time total mass, buoyancy and energy are conserved to machine precision in the fully discrete model

Dimensionality cross-over in magnetism: from domain walls (2D) to vortices (1D)

Dimensionality cross-over is a classical topic in physics. Surprisingly it has not been searched in micromagnetism, which deals with objects such as domain walls (2D) and vortices (1D). We predict by simulation a second-order transition between these two objects, with the wall length as the Landau parameter. This was conrmed experimentally based on micron-sized ux-closure dots

Innovative Weak Formulation for The Landau-Lifshitz-Gilbert Equations

A non-conventional finite element formalism is proposed to solve the dynamic Landau-Lifshitz-Gilbert micromagnetic equations. Two bidimensional test problems are treated to estimate the validity and the accuracy of this finite element approachComment: 4 pages, proceedings for Intermag Madrid, May 2008 (oral contribution

Modelling of magnetisation dynamics

On présente ici un ensemble de méthodes numériques performantes pour lasimulation micromagnétique 3D reposant sur l’équation de Landau-Lifchitz-Gilbert, constituantun code nommé feeLLGood. On a choisi l’approche éléments finis pour sa flexibilitégéométrique. La formulation adoptée respecte la contrainte d’orthogonalité entre l’aimantationet sa dérivée temporelle, contrairement à la formulation classique sur-dissipative.On met au point un schéma de point milieu pour l’équation Landau-Lifchitz-Gilbert quiest stable et d’ordre deux en temps. Cela permet de prendre, à précision égale, des pas detemps beaucoup plus grands (typiquement un ordre de grandeur) que les schémas classiques.Un véritable enjeu numérique est le calcul du champ démagnétisant, non local. Oncompare plusieurs techniques de calcul rapide pour retenir celles, inédites dans le domaine,des multipôles rapides (FMM) et des transformées de Fourier hors-réseau (NFFT). Aprèsavoir validé le code sur des cas-tests et établi son efficacité, on présente les applications àla simulation des nanostructures : sélection de chiralité et résonance ferromagnétique d’unplot monovortex de cobalt, hystérésis des chapeaux de Néel dans un plot allongé de fer.Enfin, l’étude d’un oscillateur spintronique prouve l’évolutivité du code.Here is presented a set of efficient numerical methods for 3D micromagneticsimulation based on the Landau-Lifchitz-Gilbert equation, making up a code named feeLLGood.The finite element approach was chosen for its geometrical flexibility. The adoptedformulation meets the orthogonality constraint between the magnetization and its time derivative,unlike the over-dissipative classical formulation. A midoint rule was developed forthe Landau-Lifchitz-Gilbert equation which is stable and second order in time. This allowsfor much bigger time steps (typically an order of magnitude) than classical schemes at thesame precision. Computing the nonlocal demagnetizing interaction is a real numerical challenge.Several fast computation techniques are compared. Those selected are novel to thefield : the Fast Multipole Method (FMM) and Non-uniform Fast Fourier Transforms (NFFT).After the code is validated on test cases and its efficiency established, applications to the simulationof nanostructures are presented : chirality selection and ferromagnetic resonanceof a cobalt monovortex dot, Neel caps hysteresis in an iron dot. Finally, the study of a spintronicoscillator proves the code’s upgradability

Modélisation de la dynamique de l’aimantation par éléments finis

Here is presented a set of efficient numerical methods for 3D micromagneticsimulation based on the Landau-Lifchitz-Gilbert equation, making up a code named feeLLGood.The finite element approach was chosen for its geometrical flexibility. The adoptedformulation meets the orthogonality constraint between the magnetization and its time derivative,unlike the over-dissipative classical formulation. A midoint rule was developed forthe Landau-Lifchitz-Gilbert equation which is stable and second order in time. This allowsfor much bigger time steps (typically an order of magnitude) than classical schemes at thesame precision. Computing the nonlocal demagnetizing interaction is a real numerical challenge.Several fast computation techniques are compared. Those selected are novel to thefield : the Fast Multipole Method (FMM) and Non-uniform Fast Fourier Transforms (NFFT).After the code is validated on test cases and its efficiency established, applications to the simulationof nanostructures are presented : chirality selection and ferromagnetic resonanceof a cobalt monovortex dot, Neel caps hysteresis in an iron dot. Finally, the study of a spintronicoscillator proves the code’s upgradability.On présente ici un ensemble de méthodes numériques performantes pour lasimulation micromagnétique 3D reposant sur l’équation de Landau-Lifchitz-Gilbert, constituantun code nommé feeLLGood. On a choisi l’approche éléments finis pour sa flexibilitégéométrique. La formulation adoptée respecte la contrainte d’orthogonalité entre l’aimantationet sa dérivée temporelle, contrairement à la formulation classique sur-dissipative.On met au point un schéma de point milieu pour l’équation Landau-Lifchitz-Gilbert quiest stable et d’ordre deux en temps. Cela permet de prendre, à précision égale, des pas detemps beaucoup plus grands (typiquement un ordre de grandeur) que les schémas classiques.Un véritable enjeu numérique est le calcul du champ démagnétisant, non local. Oncompare plusieurs techniques de calcul rapide pour retenir celles, inédites dans le domaine,des multipôles rapides (FMM) et des transformées de Fourier hors-réseau (NFFT). Aprèsavoir validé le code sur des cas-tests et établi son efficacité, on présente les applications àla simulation des nanostructures : sélection de chiralité et résonance ferromagnétique d’unplot monovortex de cobalt, hystérésis des chapeaux de Néel dans un plot allongé de fer.Enfin, l’étude d’un oscillateur spintronique prouve l’évolutivité du code

Modélisation de la dynamique de l’aimantation par éléments finis

Here is presented a set of efficient numerical methods for 3D micromagneticsimulation based on the Landau-Lifchitz-Gilbert equation, making up a code named feeLLGood.The finite element approach was chosen for its geometrical flexibility. The adoptedformulation meets the orthogonality constraint between the magnetization and its time derivative,unlike the over-dissipative classical formulation. A midoint rule was developed forthe Landau-Lifchitz-Gilbert equation which is stable and second order in time. This allowsfor much bigger time steps (typically an order of magnitude) than classical schemes at thesame precision. Computing the nonlocal demagnetizing interaction is a real numerical challenge.Several fast computation techniques are compared. Those selected are novel to thefield : the Fast Multipole Method (FMM) and Non-uniform Fast Fourier Transforms (NFFT).After the code is validated on test cases and its efficiency established, applications to the simulationof nanostructures are presented : chirality selection and ferromagnetic resonanceof a cobalt monovortex dot, Neel caps hysteresis in an iron dot. Finally, the study of a spintronicoscillator proves the code’s upgradability.On présente ici un ensemble de méthodes numériques performantes pour lasimulation micromagnétique 3D reposant sur l’équation de Landau-Lifchitz-Gilbert, constituantun code nommé feeLLGood. On a choisi l’approche éléments finis pour sa flexibilitégéométrique. La formulation adoptée respecte la contrainte d’orthogonalité entre l’aimantationet sa dérivée temporelle, contrairement à la formulation classique sur-dissipative.On met au point un schéma de point milieu pour l’équation Landau-Lifchitz-Gilbert quiest stable et d’ordre deux en temps. Cela permet de prendre, à précision égale, des pas detemps beaucoup plus grands (typiquement un ordre de grandeur) que les schémas classiques.Un véritable enjeu numérique est le calcul du champ démagnétisant, non local. Oncompare plusieurs techniques de calcul rapide pour retenir celles, inédites dans le domaine,des multipôles rapides (FMM) et des transformées de Fourier hors-réseau (NFFT). Aprèsavoir validé le code sur des cas-tests et établi son efficacité, on présente les applications àla simulation des nanostructures : sélection de chiralité et résonance ferromagnétique d’unplot monovortex de cobalt, hystérésis des chapeaux de Néel dans un plot allongé de fer.Enfin, l’étude d’un oscillateur spintronique prouve l’évolutivité du code

A convergent finite element approximation for Landau-Lifschitz-Gilbert equation

International audienc

Conservative interpolation between general spherical meshes

International audienceAn efficient, local, explicit, second-order, conservative interpolation algorithm between spherical meshes is presented. The cells composing the source and target meshes may be either spherical polygons or latitude–longitude quadrilaterals. Second-order accuracy is obtained by piece-wise linear finite-volume reconstruction over the source mesh. Global conservation is achieved through the introduction of a supermesh, whose cells are all possible intersections of source and target cells. Areas and intersections are computed exactly to yield a geometrically exact method. The main efficiency bottleneck caused by the construction of the supermesh is overcome by adopting tree-based data structures and algorithms, from which the mesh connectivity can also be deduced efficiently.The theoretical second-order accuracy is verified using a smooth test function and pairs of meshes commonly used for atmospheric modelling. Experiments confirm that the most expensive operations, especially the supermesh construction, have O(NlogN) computational cost. The method presented is meant to be incorporated in pre- or post-processing atmospheric modelling pipelines, or directly into models for flexible input/output. It could also serve as a basis for conservative coupling between model components, e.g., atmosphere and ocean

Modélisation de la dynamique de l'aimantation par éléments finis

On présente ici un ensemble de méthodes numériques performantes pour lasimulation micromagnétique 3D reposant sur l équation de Landau-Lifchitz-Gilbert, constituantun code nommé feeLLGood. On a choisi l approche éléments finis pour sa flexibilitégéométrique. La formulation adoptée respecte la contrainte d orthogonalité entre l aimantationet sa dérivée temporelle, contrairement à la formulation classique sur-dissipative.On met au point un schéma de point milieu pour l équation Landau-Lifchitz-Gilbert quiest stable et d ordre deux en temps. Cela permet de prendre, à précision égale, des pas detemps beaucoup plus grands (typiquement un ordre de grandeur) que les schémas classiques.Un véritable enjeu numérique est le calcul du champ démagnétisant, non local. Oncompare plusieurs techniques de calcul rapide pour retenir celles, inédites dans le domaine,des multipôles rapides (FMM) et des transformées de Fourier hors-réseau (NFFT). Aprèsavoir validé le code sur des cas-tests et établi son efficacité, on présente les applications àla simulation des nanostructures : sélection de chiralité et résonance ferromagnétique d unplot monovortex de cobalt, hystérésis des chapeaux de Néel dans un plot allongé de fer.Enfin, l étude d un oscillateur spintronique prouve l évolutivité du code.Here is presented a set of efficient numerical methods for 3D micromagneticsimulation based on the Landau-Lifchitz-Gilbert equation, making up a code named feeLLGood.The finite element approach was chosen for its geometrical flexibility. The adoptedformulation meets the orthogonality constraint between the magnetization and its time derivative,unlike the over-dissipative classical formulation. A midoint rule was developed forthe Landau-Lifchitz-Gilbert equation which is stable and second order in time. This allowsfor much bigger time steps (typically an order of magnitude) than classical schemes at thesame precision. Computing the nonlocal demagnetizing interaction is a real numerical challenge.Several fast computation techniques are compared. Those selected are novel to thefield : the Fast Multipole Method (FMM) and Non-uniform Fast Fourier Transforms (NFFT).After the code is validated on test cases and its efficiency established, applications to the simulationof nanostructures are presented : chirality selection and ferromagnetic resonanceof a cobalt monovortex dot, Neel caps hysteresis in an iron dot. Finally, the study of a spintronicoscillator proves the code s upgradability.SAVOIE-SCD - Bib.électronique (730659901) / SudocGRENOBLE1/INP-Bib.électronique (384210012) / SudocGRENOBLE2/3-Bib.électronique (384219901) / SudocSudocFranceF