-
Scalable physical source-to-field inference with hypernetworks
Authors:
Berian James,
Stefan Pollok,
Ignacio Peis,
Jes Frellsen,
Rasmus Bjørk
Abstract:
We present a generative model that amortises computation for the field around e.g. gravitational or magnetic sources. Exact numerical calculation has either computational complexity $\mathcal{O}(M\times{}N)$ in the number of sources and field evaluation points, or requires a fixed evaluation grid to exploit fast Fourier transforms. Using an architecture where a hypernetwork produces an implicit re…
▽ More
We present a generative model that amortises computation for the field around e.g. gravitational or magnetic sources. Exact numerical calculation has either computational complexity $\mathcal{O}(M\times{}N)$ in the number of sources and field evaluation points, or requires a fixed evaluation grid to exploit fast Fourier transforms. Using an architecture where a hypernetwork produces an implicit representation of the field around a source collection, our model instead performs as $\mathcal{O}(M + N)$, achieves accuracy of $\sim\!4\%-6\%$, and allows evaluation at arbitrary locations for arbitrary numbers of sources, greatly increasing the speed of e.g. physics simulations. We also examine a model relating to the physical properties of the output field and develop two-dimensional examples to demonstrate its application. The code for these models and experiments is available at https://github.com/cmt-dtu-energy/hypermagnetics.
△ Less
Submitted 7 May, 2024;
originally announced May 2024.
-
Magnetic levitation by rotation
Authors:
Joachim Marco Hermansen,
Frederik Laust Durhuus,
Cathrine Frandsen,
Marco Beleggia,
Christian R. H. Bahl,
Rasmus Bjørk
Abstract:
A permanent magnet can be levitated simply by placing it in the vicinity of another permanent magnet that rotates in the order of 200 Hz. This surprising effect can be easily reproduced in the lab with off-the-shelf components. Here we investigate this novel type of magnetic levitation experimentally and clarify the underlying physics. Using a 19 mm diameter spherical NdFeB magnet as rotor magnet,…
▽ More
A permanent magnet can be levitated simply by placing it in the vicinity of another permanent magnet that rotates in the order of 200 Hz. This surprising effect can be easily reproduced in the lab with off-the-shelf components. Here we investigate this novel type of magnetic levitation experimentally and clarify the underlying physics. Using a 19 mm diameter spherical NdFeB magnet as rotor magnet, we capture the detailed motion of levitating, spherical NdFeB magnets, denoted floater magnets. We find that as levitation occurs, the floater magnet frequency-locks with the rotor magnet, and, noticeably, that the magnetization of the floater is oriented close to the axis of rotation and towards the like pole of the rotor magnet. This is in contrast to what might be expected by the laws of magnetostatics as the floater is observed to align its magnetization essentially perpendicular to the magnetic field of the rotor. Moreover, we find that the size of the floater has a clear influence on the levitation: the smaller the floater, the higher the rotor speed necessary to achieve levitation, and the further away the levitation point shifts. We verify that magnetostatic interactions between the rotating magnets are responsible for creating the equilibrium position of the floater. Hence, this type of magnetic levitation does not rely on gravity as a balancing force to achieve an equilibrium position. Based on theoretical arguments and a numerical model, we show that a constant, vertical field and eddy-current enhanced damping is sufficient to produce levitation from rest. This enables a gyroscopically stabilised counter-intuitive steady-state moment orientation, and the resulting magnetostatically stable, mid-air equilibrium point. The numerical model display the same trends with respect to rotation speed and the floater magnet size as seen in the experiments.
△ Less
Submitted 23 October, 2023; v1 submitted 24 April, 2023;
originally announced May 2023.
-
Extraordinary magnetometry -- a review on extraordinary magnetoresistance
Authors:
Thierry Desire Pomar,
Ricci Erlandsen,
Bowen Zhou,
Leonid Iliushyn,
Rasmus Bjørk,
Dennis Valbjørn Christensen
Abstract:
Extraordinary magnetoresistance (EMR) is a geometric magnetoresistance effect occurring in hybrid devices consisting of a high-mobility material joined by a metal. The change in resistance can exceed 107% at room temperature when a magnetic field of 5 T is applied. Magnetic field sensors based on EMR hold the potential formeasuring weak magnetic fields with an unprecedented sensitivity, yet, to da…
▽ More
Extraordinary magnetoresistance (EMR) is a geometric magnetoresistance effect occurring in hybrid devices consisting of a high-mobility material joined by a metal. The change in resistance can exceed 107% at room temperature when a magnetic field of 5 T is applied. Magnetic field sensors based on EMR hold the potential formeasuring weak magnetic fields with an unprecedented sensitivity, yet, to date this potential is largely unmet. In this work, we provide an extensive review of the current state-of-the-art in EMR sensors with a focus on the hybrid device geometries, the constituent material properties and applications of EMR. We present a direct comparison of the best devices in literature across magnetoresistance, sensitivity and noise equivalent field for different materials and geometric designs. The compilation of studies collected in this review illustrates the extremely rich possibilities for tuning the magnetoresistive behavior varying the device geometry and material properties. In addition, we aim to improve the understanding of the EMR effect and its interplay with geometry and material properties. Finally, we discuss recent trends in the field and future perspectives for EMR.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
Magnetic Field Prediction Using Generative Adversarial Networks
Authors:
Stefan Pollok,
Nataniel Olden-Jørgensen,
Peter Stanley Jørgensen,
Rasmus Bjørk
Abstract:
Plenty of scientific and real-world applications are built on magnetic fields and their characteristics. To retrieve the valuable magnetic field information in high resolution, extensive field measurements are required, which are either time-consuming to conduct or even not feasible due to physical constraints. To alleviate this problem, we predict magnetic field values at a random point in space…
▽ More
Plenty of scientific and real-world applications are built on magnetic fields and their characteristics. To retrieve the valuable magnetic field information in high resolution, extensive field measurements are required, which are either time-consuming to conduct or even not feasible due to physical constraints. To alleviate this problem, we predict magnetic field values at a random point in space from a few point measurements by using a generative adversarial network (GAN) structure. The deep learning (DL) architecture consists of two neural networks: a generator, which predicts missing field values of a given magnetic field, and a critic, which is trained to calculate the statistical distance between real and generated magnetic field distributions. By minimizing this statistical distance, a reconstruction loss as well as physical losses, our trained generator has learned to predict the missing field values with a median reconstruction test error of 5.14%, when a single coherent region of field points is missing, and 5.86%, when only a few point measurements in space are available and the field measurements around are predicted. We verify the results on an experimentally validated field.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
MagTense: a micromagnetic framework using the analytical demagnetization tensor
Authors:
R. Bjørk,
E. B. Poulsen,
K. K. Nielsen,
A. R. Insinga
Abstract:
We present the open source micromagnetic framework, MagTense, which utilizes a novel discretization approach of rectangular cuboid or tetrahedron geometry "tiles" to analytically calculate the demagnetization field. Each tile is assumed to be uniformly magnetized, and from this assumption only, the demagnetization field can be analytically calculated. Using this novel approach we calculate the sol…
▽ More
We present the open source micromagnetic framework, MagTense, which utilizes a novel discretization approach of rectangular cuboid or tetrahedron geometry "tiles" to analytically calculate the demagnetization field. Each tile is assumed to be uniformly magnetized, and from this assumption only, the demagnetization field can be analytically calculated. Using this novel approach we calculate the solution to the mmag standard micromagnetic problems 2, 3 and 4 and find that the MagTense framework accurately predicts the solution to each of these. Finally, we show that simulation time can be significantly improved by performing the dense demagnetization tensor matrix multiplications using NVIDIA CUDA.
△ Less
Submitted 2 July, 2021;
originally announced July 2021.
-
The magnetic field from a homogeneously magnetized cylindrical tile
Authors:
K. K. Nielsen,
R. Bjørk
Abstract:
The magnetic field of a homogeneously magnetized cylindrical tile geometry, i.e. an angular section of a finite hollow cylinder, is found. The field is expressed as the product between a tensor field describing the geometrical part of the problem and a column vector holding the magnetization of the tile. Outside the tile, the tensor is identical to the demagnetization tensor. We find that four com…
▽ More
The magnetic field of a homogeneously magnetized cylindrical tile geometry, i.e. an angular section of a finite hollow cylinder, is found. The field is expressed as the product between a tensor field describing the geometrical part of the problem and a column vector holding the magnetization of the tile. Outside the tile, the tensor is identical to the demagnetization tensor. We find that four components of the tensor, $N_{xy},N_{xz},N_{yz}$ and $N_{zy}$, can be expressed fully analytically, while the five remaining components, $N_{xx},N_{yx},N_{yy},N_{zx}$ and $N_{zz}$, contain integrals that have to be evaluated numerically. When evaluated numerically the tensor is symmetric. A comparison between the found solution, implemented in the open source magnetic framework MagTense, and a finite element calculation of the magnetic flux density of a cylindrical tile shows excellent agreement.
△ Less
Submitted 15 April, 2020;
originally announced April 2020.
-
The stray- and demagnetizing field from a homogeneously magnetized tetrahedron
Authors:
Kaspar K. Nielsen,
Andrea R. Insinga,
Rasmus Bjørk
Abstract:
The stray- and demagnetization tensor field for a homogeneously magnetized tetrahedron is found analytically. The tetrahedron is a special case of four triangular faces with constant magnetization-charge surface density, for which we also determine the tensor field. The tensor field is implemented in the open source micromagnetic and magnetostatic simulation framework MagTense and compared with th…
▽ More
The stray- and demagnetization tensor field for a homogeneously magnetized tetrahedron is found analytically. The tetrahedron is a special case of four triangular faces with constant magnetization-charge surface density, for which we also determine the tensor field. The tensor field is implemented in the open source micromagnetic and magnetostatic simulation framework MagTense and compared with the obtained magnetic field from an FEM solution, showing excellent agreement. This result is important for modeling magnetostatics in general and for micromagnetism in particular as the demagnetizing field of an arbitrary body discretized using conventional meshing techniques is significantly simplified with this approach.
△ Less
Submitted 15 April, 2020;
originally announced April 2020.
-
The demagnetization factor for randomly packed spheroidal particles
Authors:
R. Bjørk,
Z. Zhou
Abstract:
We investigate if the demagnetization factor for a randomly packed powder of magnetic spheroidal particles depend on the shape of the spheroidal particles and what the internal variation in magnetization is within such a powder. A spheroid is an ellipsoid of revolution, i.e. an ellipsoid with two semi-major axis being equal. The demagnetization factor is calculated as function of particle aspect r…
▽ More
We investigate if the demagnetization factor for a randomly packed powder of magnetic spheroidal particles depend on the shape of the spheroidal particles and what the internal variation in magnetization is within such a powder. A spheroid is an ellipsoid of revolution, i.e. an ellipsoid with two semi-major axis being equal. The demagnetization factor is calculated as function of particle aspect ratio using two independent numerical models for several different packings, and assuming a relative permeability of 2. The calculated demagnetization factor is shown to depend on particle aspect ratio, not because of direct magnetic interaction but because the particle packing depend on the aspect ratio of the particles. The relative standard deviation of the magnetization across the powder was 3\%-8\%, increasing as the particle shape deviates from spherical, while the relative standard deviation within each particle was relatively constant around 5\%.
△ Less
Submitted 10 January, 2019;
originally announced January 2019.
-
A topology optimized switchable permanent magnet system
Authors:
R. Bjørk,
A. R. Insinga
Abstract:
The design of a magnetic field source that can switch from a high field to a low field configuration by rotation by $90^\circ$ of a set of iron pieces is investigated using topology optimization. A Halbach cylinder is considered as the magnetic field source and iron inserts are placed in the air gap of the Halbach cylinder. The ideal shape of these iron inserts is determined as function of the fie…
▽ More
The design of a magnetic field source that can switch from a high field to a low field configuration by rotation by $90^\circ$ of a set of iron pieces is investigated using topology optimization. A Halbach cylinder is considered as the magnetic field source and iron inserts are placed in the air gap of the Halbach cylinder. The ideal shape of these iron inserts is determined as function of the field generated by the Halbach cylinder and as function of the size of the iron segments. The topology optimized structures are parabolic shaped pieces and have a difference in flux density between the high and low positions that is on average 1.29 times higher than optimized regular pole pieces. The maximum increase is a factor of 2.08 times higher than the regular pole pieces.
△ Less
Submitted 10 January, 2019;
originally announced January 2019.
-
The maximum theoretical performance of unconcentrated solar photovoltaic and thermoelectric generator systems
Authors:
Rasmus Bjørk,
Kaspar K. Nielsen
Abstract:
The maximum efficiency for photovoltaic (PV) and thermoelectric generator (TEG) systems without concentration is investigated. Both a combined system where the TEG is mounted directly on the back of the PV and a tandem system where the incoming sunlight is split, and the short wavelength radiation is sent to the PV and the long wavelength to the TEG, are considered. An analytical model based on th…
▽ More
The maximum efficiency for photovoltaic (PV) and thermoelectric generator (TEG) systems without concentration is investigated. Both a combined system where the TEG is mounted directly on the back of the PV and a tandem system where the incoming sunlight is split, and the short wavelength radiation is sent to the PV and the long wavelength to the TEG, are considered. An analytical model based on the Shockley-Queisser efficiency limit for PVs and the TEG figure of merit parameter $zT$ is presented. It is shown that for non-concentrated sunlight, even if the TEG operates at the Carnot efficiency and the PV performance is assumed independent of temperature, the maximum increase in efficiency is 4.5 percentage points (pp.) for the combined case and 1.8 pp. for the tandem case compared to a stand alone PV. For a more realistic case with a temperature dependent PV and a realistic TEG, the gain in performance is much lower. For the combined PV and TEG system it is shown that a minimum $zT$ value is needed in order for the system to be more efficient than a stand alone PV system.
△ Less
Submitted 14 November, 2017; v1 submitted 18 September, 2017;
originally announced October 2017.
-
The magnetic properties of the hollow cylindrical ideal remanence magnet
Authors:
R. Bjørk
Abstract:
We consider the magnetic properties of the hollow cylindrical ideal remanence magnet. This magnet is the cylindrical permanent magnet that generates a uniform field in the cylinder bore, using the least amount of magnetic energy to do so. The remanence distribution of this magnet is derived and the generated field is compared to that of a Halbach cylinder of equal dimensions. The ideal remanence m…
▽ More
We consider the magnetic properties of the hollow cylindrical ideal remanence magnet. This magnet is the cylindrical permanent magnet that generates a uniform field in the cylinder bore, using the least amount of magnetic energy to do so. The remanence distribution of this magnet is derived and the generated field is compared to that of a Halbach cylinder of equal dimensions. The ideal remanence magnet is shown in most cases to generate a significantly lower field than the equivalent Halbach cylinder, although the field is generated with higher efficiency. The most efficient Halbach cylinder is shown to generate a field exactly twice as large as the equivalent ideal remanence magnet.
△ Less
Submitted 27 September, 2016;
originally announced September 2016.
-
Topology optimized permanent magnet systems
Authors:
R. Bjørk,
C. R. H. Bahl,
A. R. Insinga
Abstract:
Topology optimization of permanent magnet systems consisting of permanent magnets, high permeability iron and air is presented. An implementation of topology optimization for magnetostatics is discussed and three examples are considered. First, the Halbach cylinder is topology optimized with iron and an increase of 15% in magnetic efficiency is shown, albeit with an increase of 3.8 pp. in field in…
▽ More
Topology optimization of permanent magnet systems consisting of permanent magnets, high permeability iron and air is presented. An implementation of topology optimization for magnetostatics is discussed and three examples are considered. First, the Halbach cylinder is topology optimized with iron and an increase of 15% in magnetic efficiency is shown, albeit with an increase of 3.8 pp. in field inhomogeneity - a value compared to the inhomogeneity in a 16 segmented Halbach cylinder. Following this a topology optimized structure to concentrate a homogeneous field is shown to increase the magnitude of the field by 111% for the chosen dimensions. Finally, a permanent magnet with alternating high and low field regions is considered. Here a $Λ_\mathrm{cool}$ figure of merit of 0.472 is reached, which is an increase of 100% compared to a previous optimized design.
△ Less
Submitted 3 May, 2017; v1 submitted 26 September, 2016;
originally announced September 2016.
-
An analytical model for the influence of contact resistance on thermoelectric efficiency
Authors:
R. Bjørk
Abstract:
An analytical model is presented that can account for both electrical and hot and cold thermal contact resistances when calculating the efficiency of a thermoelectric generator. The model is compared to a numerical model of a thermoelectric leg, for 16 different thermoelectric materials, as well as the analytical models of Ebling et. al. (2010) and Min \& Rowe (1992). The model presented here is s…
▽ More
An analytical model is presented that can account for both electrical and hot and cold thermal contact resistances when calculating the efficiency of a thermoelectric generator. The model is compared to a numerical model of a thermoelectric leg, for 16 different thermoelectric materials, as well as the analytical models of Ebling et. al. (2010) and Min \& Rowe (1992). The model presented here is shown to accurately calculate the efficiency for all systems and all contact resistances considered, with an average difference in efficiency between the numerical model and the analytical model of $-0.07\pm0.35$ pp. This makes the model more accurate than previously published models. The maximum absolute difference in efficiency between the analytical model and the numerical model is 1.14 pp for all materials and all contact resistances considered.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
A thermoelectric power generating heat exchanger: Part I - Experimental realization
Authors:
R. Bjørk,
A. Sarhadi,
N. Pryds,
N. Lindeburg,
P. Viereck
Abstract:
An experimental realization of a heat exchanger with commercial thermoelectric generators (TEGs) is presented. The power producing capabilities as a function of flow rate and temperature span are characterized for two different commercial heat transfer fluids and for three different thermal interface materials. The device is shown to produce 2 W per TEG or 0.22 W cm$^{-2}$ at a fluid temperature d…
▽ More
An experimental realization of a heat exchanger with commercial thermoelectric generators (TEGs) is presented. The power producing capabilities as a function of flow rate and temperature span are characterized for two different commercial heat transfer fluids and for three different thermal interface materials. The device is shown to produce 2 W per TEG or 0.22 W cm$^{-2}$ at a fluid temperature difference of 175 $^\circ$C and a flow rate per fluid channel of 5 L min$^{-1}$. One experimentally realized design produced 200 W in total from 100 TEGs. For the design considered here, the power production is shown to depend more critically on the fluid temperature span than on the fluid flow rate. Finally, the temperature span across the TEG is shown to be 55% to 75% of the temperature span between the hot and cold fluids.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
The lifetime cost of a magnetic refrigerator
Authors:
R. Bjørk,
C. R. H. Bahl,
K. K. Nielsen
Abstract:
The total cost of a 25 W average load magnetic refrigerator using commercial grade Gd is calculated using a numerical model. The price of magnetocaloric material, magnet material and cost of operation are considered, and all influence the total cost. The lowest combined total cost with a device lifetime of 15 years is found to be in the range \$150-\$400 depending on the price of the magnetocalori…
▽ More
The total cost of a 25 W average load magnetic refrigerator using commercial grade Gd is calculated using a numerical model. The price of magnetocaloric material, magnet material and cost of operation are considered, and all influence the total cost. The lowest combined total cost with a device lifetime of 15 years is found to be in the range \$150-\$400 depending on the price of the magnetocaloric and magnet material. The cost of the magnet is largest, followed closely by the cost of operation, while the cost of the magnetocaloric material is almost negligible. For the lowest cost device, the optimal magnetic field is about 1.4 T, the particle size is 0.23 mm, the length of the regenerator is 40-50 mm and the utilization is about 0.2, for all device lifetimes and material and magnet prices, while the operating frequency vary as function of device lifetime. The considered performance characteristics are based on the performance of a conventional A$^{+++}$ refrigeration unit. In a rough life time cost comparison between the magnetic refrigeration device and such a unit we find similar costs, the former being slightly cheaper, assuming the cost of the magnet can be recuperated at end of life.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Modeling the microstructural evolution during constrained sintering
Authors:
R. Bjørk,
H. L. Frandsen,
N. Pryds
Abstract:
A numerical model able to simulate solid-state constrained sintering is presented. The model couples an existing kinetic Monte Carlo (kMC) model for free sintering with a finite element model (FEM) for calculating stresses on a microstructural level. The microstructural response to the local stress as well as the FEM calculation of the stress field from the microstructural evolution is discussed.…
▽ More
A numerical model able to simulate solid-state constrained sintering is presented. The model couples an existing kinetic Monte Carlo (kMC) model for free sintering with a finite element model (FEM) for calculating stresses on a microstructural level. The microstructural response to the local stress as well as the FEM calculation of the stress field from the microstructural evolution is discussed. The sintering behavior of a sample constrained by a rigid substrate is simulated. The constrained sintering results in a larger number of pores near the substrate, as well as anisotropic sintering shrinkage, with significantly enhanced strain in the central upper part of the sample surface, and minimal strain at the edges near the substrate. All these features have also previously been observed experimentally.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
The performance of a combined solar photovoltaic (PV) and thermoelectric generator (TEG) system
Authors:
R. Bjørk,
K. K. Nielsen
Abstract:
The performance of a combined solar photovoltaic (PV) and thermoelectric generator (TEG) system is examined using an analytical model for four different types of commercial PVs and a commercial bismuth telluride TEG. The TEG is applied directly on the back of the PV, so that the two devices have the same temperature. The PVs considered are crystalline Si (c-Si), amorphous Si (a-Si), copper indium…
▽ More
The performance of a combined solar photovoltaic (PV) and thermoelectric generator (TEG) system is examined using an analytical model for four different types of commercial PVs and a commercial bismuth telluride TEG. The TEG is applied directly on the back of the PV, so that the two devices have the same temperature. The PVs considered are crystalline Si (c-Si), amorphous Si (a-Si), copper indium gallium (di)selenide (CIGS) and cadmium telluride (CdTe) cells. The degradation of PV performance with temperature is shown to dominate the increase in power produced by the TEG, due to the low efficiency of the TEG. For c-Si, CIGS and CdTe PV cells the combined system produces a lower power and has a lower efficiency than the PV alone, whereas for an a-Si cell the total system performance may be slightly increased by the TEG.
△ Less
Submitted 6 August, 2015;
originally announced August 2015.
-
The universal influence of contact resistance on the efficiency of a thermoelectric generator
Authors:
R. Bjørk
Abstract:
The influence of electrical and thermal contact resistance on the efficiency of a segmented thermoelectric generator is investigated. We consider 12 different segmented $p$-legs and 12 different segmented $n$-legs, using 8 different $p$-type and 8 different $n$-type thermoelectric materials. For all systems a universal influence of both the electrical and thermal contact resistance is observed on…
▽ More
The influence of electrical and thermal contact resistance on the efficiency of a segmented thermoelectric generator is investigated. We consider 12 different segmented $p$-legs and 12 different segmented $n$-legs, using 8 different $p$-type and 8 different $n$-type thermoelectric materials. For all systems a universal influence of both the electrical and thermal contact resistance is observed on the leg's efficiency, when the systems are analyzed in terms of the contribution of the contact resistance to the total resistance of the leg. The results are compared with the analytical model of Min and Rowe (1992). In order for the efficiency not to decrease more than 20%, the contact electrical resistance should be less than 30% of the total leg resistance for zero thermal contact resistance, while the thermal contact resistance should be less than 20% for zero electrical contact resistance. The universal behavior also allowed the maximum tolerable contact resistance for a segmented system to be found, i.e. the resistance at which a leg of only the high temperature thermoelectric material has the same efficiency as the segmented leg with a contact resistance at the interface. If e.g. segmentation increases the efficiency by 30% then an electrical contact resistance of 30% or a thermal contact resistance of 20% can be tolerated.
△ Less
Submitted 5 August, 2015;
originally announced August 2015.
-
The efficiency and the demagnetization field of a general Halbach cylinder
Authors:
R. Bjørk,
A. Smith,
C. R. H. Bahl
Abstract:
The maximum magnetic efficiency of a general multipole Halbach cylinder of order $p$ is found as function of $p$. The efficiency is shown to decrease for increasing absolute value of $p$. The optimal ratio between the inner and outer radius, i.e. the ratio resulting in the most efficient design, is also found as function of $p$ and is shown to tend towards smaller and smaller magnet sizes. Finally…
▽ More
The maximum magnetic efficiency of a general multipole Halbach cylinder of order $p$ is found as function of $p$. The efficiency is shown to decrease for increasing absolute value of $p$. The optimal ratio between the inner and outer radius, i.e. the ratio resulting in the most efficient design, is also found as function of $p$ and is shown to tend towards smaller and smaller magnet sizes. Finally, the demagnetizing field in a general $p$-Halbach cylinder is calculated, and it is shown that demagnetization is largest either at $\cos 2pφ=1$ or $\cos 2pφ=-1$. For the common case of a $p=1$ Halbach cylinder the maximum values of the demagnetizing field is either at $φ= 0,π$ at the outer radius, where the field is always equal to the remanence, or at $φ= \pm π/2$ at the inner radius, where it is the magnitude of the field in the bore. Thus to avoid demagnetization the coercivity of the magnets must be larger than these values.
△ Less
Submitted 24 February, 2015;
originally announced February 2015.
-
Metamaterial anisotropic flux concentrators and magnetic arrays
Authors:
R. Bjørk,
A. Smith,
C. R. H. Bahl
Abstract:
A metamaterial magnetic flux concentrator is investigated in detail in combination with a Halbach cylinder of infinite length. A general analytical solution to the field is determined and the magnetic figure of merit is determined for a Halbach cylinder with a flux concentrator. It is shown that an ideal flux concentrator will not change the figure of merit of a given magnet design, while the non-…
▽ More
A metamaterial magnetic flux concentrator is investigated in detail in combination with a Halbach cylinder of infinite length. A general analytical solution to the field is determined and the magnetic figure of merit is determined for a Halbach cylinder with a flux concentrator. It is shown that an ideal flux concentrator will not change the figure of merit of a given magnet design, while the non-ideal will always lower it. The geometric parameters producing maximum figure of merit, i.e. the most efficient devices, are determined. The force and torque between two concentric Halbach cylinders with flux concentrators is determined and the maximum torque is found. Finally, the effect of non-ideal flux concentrators and the practical use of flux concentrators, as well as demagnetization issues, is discussed.
△ Less
Submitted 23 October, 2014;
originally announced October 2014.
-
Determining the minimum mass and cost of a magnetic refrigerator
Authors:
R. Bjørk,
A. Smith,
C. R. H. Bahl,
N. Pryds
Abstract:
An expression is determined for the mass of the magnet and magnetocaloric material needed for a magnetic refrigerator and these are determined using numerical modeling for both parallel plate and packed sphere bed regenerators as function of temperature span and cooling power. As magnetocaloric material Gd or a model material with a constant adiabatic temperature change, representing a infinitely…
▽ More
An expression is determined for the mass of the magnet and magnetocaloric material needed for a magnetic refrigerator and these are determined using numerical modeling for both parallel plate and packed sphere bed regenerators as function of temperature span and cooling power. As magnetocaloric material Gd or a model material with a constant adiabatic temperature change, representing a infinitely linearly graded refrigeration device, is used. For the magnet a maximum figure of merit magnet or a Halbach cylinder is used. For a cost of \$40 and \$20 per kg for the magnet and magnetocaloric material, respectively, the cheapest 100 W parallel plate refrigerator with a temperature span of 20 K using Gd and a Halbach magnet has 0.8 kg of magnet, 0.3 kg of Gd and a cost of \$35. Using the constant material reduces this cost to \$25. A packed sphere bed refrigerator with the constant material costs \$7. It is also shown that increasing the operation frequency reduces the cost. Finally, the lowest cost is also found as a function of the cost of the magnet and magnetocaloric material.
△ Less
Submitted 23 October, 2014;
originally announced October 2014.
-
The influence of the magnetic field on the performance of an active magnetic regenerator (AMR)
Authors:
R. Bjørk,
K. Engelbrecht
Abstract:
The influence of the time variation of the magnetic field, termed the magnetic field profile, on the performance of a magnetocaloric refrigeration device using the active magnetic regeneration (AMR) cycle is studied for a number of process parameters for both a parallel plate and packed bed regenerator using a numerical model. The cooling curve of the AMR is shown to be almost linear far from the…
▽ More
The influence of the time variation of the magnetic field, termed the magnetic field profile, on the performance of a magnetocaloric refrigeration device using the active magnetic regeneration (AMR) cycle is studied for a number of process parameters for both a parallel plate and packed bed regenerator using a numerical model. The cooling curve of the AMR is shown to be almost linear far from the Curie temperature of the magnetocaloric material. It is shown that a magnetic field profile that is 10% of the cycle time out of sync with the flow profile leads to a drop in both the maximum temperature span and the maximum cooling capacity of 20-40\% for both parallel plate and packed bed regenerators. The maximum cooling capacity is shown to depend very weakly on the ramp rate of the magnetic field. Reducing the temporal width of the high field portion of the magnetic field profile by 10% leads to a drop in maximum temperature span and maximum cooling capacity of 5-20%. An increase of the magnetic field from 1 T to 1.5 T increases the maximum cooling capacity by 30-50% but the maximum temperature span by only 20-30%. Finally, it was seen that the influence of changing the magnetic field was more or less the same for the different regenerator geometries and operating parameters studied here. This means that the design of the magnet can be done independently of the regenerator geometry.
△ Less
Submitted 23 October, 2014;
originally announced October 2014.
-
Comparison of adjustable permanent magnetic field sources
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
N. Pryds
Abstract:
A permanent magnet assembly in which the flux density can be altered by a mechanical operation is often significantly smaller than comparable electromagnets and also requires no electrical power to operate. In this paper five permanent magnet designs in which the magnetic flux density can be altered are analyzed using numerical simulations, and compared based on the generated magnetic flux density…
▽ More
A permanent magnet assembly in which the flux density can be altered by a mechanical operation is often significantly smaller than comparable electromagnets and also requires no electrical power to operate. In this paper five permanent magnet designs in which the magnetic flux density can be altered are analyzed using numerical simulations, and compared based on the generated magnetic flux density in a sample volume and the amount of magnet material used. The designs are the concentric Halbach cylinder, the two half Halbach cylinders, the two linear Halbach arrays and the four and six rod mangle. The concentric Halbach cylinder design is found to be the best performing design, i.e. the design that provides the most magnetic flux density using the least amount of magnet material. A concentric Halbach cylinder has been constructed and the magnetic flux density, the homogeneity and the direction of the magnetic field are measured and compared with numerical simulation and a good agreement is found.
△ Less
Submitted 10 October, 2014;
originally announced October 2014.
-
Improving magnet designs with high and low field regions
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
N. Pryds
Abstract:
A general scheme for increasing the difference in magnetic flux density between a high and a low magnetic field region by removing unnecessary magnet material is presented. This is important in, e.g., magnetic refrigeration where magnet arrays has to deliver high field regions in close proximity to low field regions. Also, a general way to replace magnet material with a high permeability soft magn…
▽ More
A general scheme for increasing the difference in magnetic flux density between a high and a low magnetic field region by removing unnecessary magnet material is presented. This is important in, e.g., magnetic refrigeration where magnet arrays has to deliver high field regions in close proximity to low field regions. Also, a general way to replace magnet material with a high permeability soft magnetic material where appropriate is discussed. As an example these schemes are applied to a two dimensional concentric Halbach cylinder design resulting in a reduction of the amount of magnet material used by 42% while increasing the difference in flux density between a high and a low field region by 45%.
△ Less
Submitted 10 October, 2014;
originally announced October 2014.
-
Magnetocaloric properties of LaFe$_{13-x-y}$Co$_x$Si$_y$ and commercial grade Gd
Authors:
R. Bjørk,
C. R. H. Bahl,
M. Katter
Abstract:
The magnetocaloric properties of three samples of LaFe$_{13-x-y}$Co$_x$Si$_y$ have been measured and compared to measurements of commercial grade Gd. The samples have (x=0.86, y=1.08), (x=0.94, y=1.01) and (x=0.97, y=1.07) yielding Curie temperatures in the range 276-288 K. The magnetization, specific heat capacity and adiabatic temperature change have been measured over a broad temperature interv…
▽ More
The magnetocaloric properties of three samples of LaFe$_{13-x-y}$Co$_x$Si$_y$ have been measured and compared to measurements of commercial grade Gd. The samples have (x=0.86, y=1.08), (x=0.94, y=1.01) and (x=0.97, y=1.07) yielding Curie temperatures in the range 276-288 K. The magnetization, specific heat capacity and adiabatic temperature change have been measured over a broad temperature interval. Importantly, all measurements were corrected for demagnetization, allowing the data to be directly compared. In an internal field of 1 T the maximum specific entropy changes were 6.2, 5.1 and 5.0 J/kg K, the specific heat capacities were 910, 840 and 835 J/kg K and the adiabatic temperature changes were 2.3, 2.1 and 2.1 K for the three LaFeCoSi samples respectively. For Gd in an internal field of 1 T the maximum specific entropy change was 3.1 J/kg K, the specific heat capacity was 340 J/kg K and the adiabatic temperature change was 3.3 K. The adiabatic temperature change was also calculated from the measured values of the specific heat capacity and specific magnetization and compared to the directly measured values. In general an excellent agreement was seen.
△ Less
Submitted 8 October, 2014;
originally announced October 2014.
-
An optimized magnet for magnetic refrigeration
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
D. V. Christensen,
N. Pryds
Abstract:
A magnet designed for use in a magnetic refrigeration device is presented. The magnet is designed by applying two general schemes for improving a magnet design to a concentric Halbach cylinder magnet design and dimensioning and segmenting this design in an optimum way followed by the construction of the actual magnet. The final design generates a peak value of 1.24 T, an average flux density of 0.…
▽ More
A magnet designed for use in a magnetic refrigeration device is presented. The magnet is designed by applying two general schemes for improving a magnet design to a concentric Halbach cylinder magnet design and dimensioning and segmenting this design in an optimum way followed by the construction of the actual magnet. The final design generates a peak value of 1.24 T, an average flux density of 0.9 T in a volume of 2 L using only 7.3 L of magnet, and has an average low flux density of 0.08 T also in a 2 L volume. The working point of all the permanent magnet blocks in the design is very close to the maximum energy density. The final design is characterized in terms of a performance parameter, and it is shown that it is one of the best performing magnet designs published for magnetic refrigeration.
△ Less
Submitted 8 October, 2014;
originally announced October 2014.
-
The ideal dimensions of a Halbach cylinder of finite length
Authors:
R. Bjørk
Abstract:
In this paper the smallest or optimal dimensions of a Halbach cylinder of a finite length for a given sample volume and desired flux density are determined using numerical modeling and parameter variation. A sample volume that is centered in and shaped as the Halbach cylinder bore but with a possible shorter length is considered. The external radius and the length of the Halbach cylinder with the…
▽ More
In this paper the smallest or optimal dimensions of a Halbach cylinder of a finite length for a given sample volume and desired flux density are determined using numerical modeling and parameter variation. A sample volume that is centered in and shaped as the Halbach cylinder bore but with a possible shorter length is considered. The external radius and the length of the Halbach cylinder with the smallest possible dimensions are found as a function of a desired internal radius, length of the sample volume and mean flux density. It is shown that the optimal ratio between the outer and inner radius of the Halbach cylinder does not depend on the length of the sample volume. Finally, the efficiency of a finite length Halbach cylinder is considered and compared with the case of a cylinder of infinite length. The most efficient dimensions for a Halbach cylinder are found and it is shown that the efficiency increases slowly with the length of the cylinder.
△ Less
Submitted 2 October, 2014;
originally announced October 2014.
-
The sintering behavior of close packed spheres
Authors:
R. Bjørk,
V. Tikare,
H. L. Frandsen,
N. Pryds
Abstract:
The sintering behavior of close packed spheres is investigated using a numerical model. The investigated systems are the body centered cubic (BCC), face centered cubic (FCC) and hexagonal closed packed spheres (HCP). The sintering behavior is found to be ideal with no grain growth until full density is reached for all systems. During sintering the grains change shape from spherical to tetrakaideca…
▽ More
The sintering behavior of close packed spheres is investigated using a numerical model. The investigated systems are the body centered cubic (BCC), face centered cubic (FCC) and hexagonal closed packed spheres (HCP). The sintering behavior is found to be ideal with no grain growth until full density is reached for all systems. During sintering the grains change shape from spherical to tetrakaidecahedron, similar to the geometry analyzed by Coble (R. L. Coble, J. Appl. Phys. 32 (1961) 787).
△ Less
Submitted 2 October, 2014;
originally announced October 2014.
-
Review and comparison of magnet designs for magnetic refrigeration
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
N. Pryds
Abstract:
One of the key issues in magnetic refrigeration is generating the magnetic field that the magnetocaloric material must be subjected to. The magnet constitutes a major part of the expense of a complete magnetic refrigeration system and a large effort should therefore be invested in improving the magnet design. A detailed analysis of the efficiency of different published permanent magnet designs use…
▽ More
One of the key issues in magnetic refrigeration is generating the magnetic field that the magnetocaloric material must be subjected to. The magnet constitutes a major part of the expense of a complete magnetic refrigeration system and a large effort should therefore be invested in improving the magnet design. A detailed analysis of the efficiency of different published permanent magnet designs used in magnetic refrigeration applications is presented in this paper. Each design is analyzed based on the generated magnetic flux density, the volume of the region where this flux is generated and the amount of magnet material used. This is done by characterizing each design by a figure of merit magnet design efficiency parameter, $Λ_\mathrm{cool}$. The designs are then compared and the best design found. Finally recommendations for designing the ideal magnet design are presented based on the analysis of the reviewed designs.
△ Less
Submitted 29 September, 2014;
originally announced September 2014.
-
Strain in the mesoscale kinetic Monte Carlo model for sintering
Authors:
R. Bjørk,
H. L. Frandsen,
V. Tikare,
E. Olevsky,
N. Pryds
Abstract:
Shrinkage strains measured from microstructural simulations using the mesoscale kinetic Monte Carlo (kMC) model for solid state sintering are discussed. This model represents the microstructure using digitized discrete sites that are either grain or pore sites. The algorithm used to simulate densification by vacancy annihilation removes an isolated pore site at a grain boundary and collapses a col…
▽ More
Shrinkage strains measured from microstructural simulations using the mesoscale kinetic Monte Carlo (kMC) model for solid state sintering are discussed. This model represents the microstructure using digitized discrete sites that are either grain or pore sites. The algorithm used to simulate densification by vacancy annihilation removes an isolated pore site at a grain boundary and collapses a column of sites extending from the vacancy to the surface of sintering compact, through the center of mass of the nearest grain. Using this algorithm, the existing published kMC models are shown to produce anisotropic strains for homogeneous powder compacts with aspect ratios different from unity. It is shown that the line direction biases shrinkage strains in proportion the compact dimension aspect ratios. A new algorithm that corrects this bias in strains is proposed; the direction for collapsing the column is determined by choosing a random sample face and subsequently a random point on that face as the end point for an annihilation path with equal probabilities. This algorithm is mathematically and experimentally shown to result in isotropic strains for all samples regardless of their dimensions. Finally, the microstructural evolution is shown to be similar for the new and old annihilation algorithms.
△ Less
Submitted 29 September, 2014;
originally announced September 2014.
-
Analysis of the internal heat losses in a thermoelectric generator
Authors:
R. Bjørk,
D. V. Christensen,
D. Eriksen,
N. Pryds
Abstract:
A 3D thermoelectric numerical model is used to investigate different internal heat loss mechanisms for a thermoelectric generator with bismuth telluride p- and n-legs. The model considers all thermoelectric effects, temperature dependent material parameters and simultaneous convective, conductive and radiative heat losses, including surface to surface radiation. For radiative heat losses it is sho…
▽ More
A 3D thermoelectric numerical model is used to investigate different internal heat loss mechanisms for a thermoelectric generator with bismuth telluride p- and n-legs. The model considers all thermoelectric effects, temperature dependent material parameters and simultaneous convective, conductive and radiative heat losses, including surface to surface radiation. For radiative heat losses it is shown that for the temperatures considered here, surface to ambient radiation is a good approximation of the heat loss. For conductive heat transfer the module efficiency is shown to be comparable to the case of radiative losses. Finally, heat losses due to internal natural convection in the module is shown to be negligible for the millimetre sized modules considered here. The combined case of radiative and conductive heat transfer resulted in the lowest efficiency. The optimized load resistance is found to decrease for increased heat loss. The leg dimensions are varied for all heat losses cases and it is shown that the ideal way to construct a TEG module with minimal heat losses and maximum efficiency is to either use a good insulating material between the legs or evacuate the module completely, and use small and wide legs closely spaced.
△ Less
Submitted 19 September, 2014;
originally announced September 2014.
-
Demagnetization factor for a powder of randomly packed spherical particles
Authors:
R. Bjørk,
C. R. H. Bahl
Abstract:
The demagnetization factors for randomly packed spherical particle powders with different porosities, sample aspect ratios and monodisperse, normal and log-normal particle size distributions have been calculated using a numerical model. For a relative permeability of 2, comparable to room temperature Gd, the calculated demagnetization factor is close to the theoretical value. The normalized standa…
▽ More
The demagnetization factors for randomly packed spherical particle powders with different porosities, sample aspect ratios and monodisperse, normal and log-normal particle size distributions have been calculated using a numerical model. For a relative permeability of 2, comparable to room temperature Gd, the calculated demagnetization factor is close to the theoretical value. The normalized standard deviation of the magnetization in the powder was 6.0%-6.7%. The demagnetization factor decreased significantly, while the standard deviation of the magnetization increased, for increasing relative permeability.
△ Less
Submitted 19 September, 2014;
originally announced September 2014.
-
Optimization and improvement of Halbach cylinder design
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
N. Pryds
Abstract:
In this paper we describe the results of a parameter survey of a 16 segmented Halbach cylinder in three dimensions in which the parameters internal radius, $r_{\mathrm{\scriptsize{in}}}$, external radius, $r_{\mathrm{\scriptsize{ex}}}$, and length, $L$, have been varied. Optimal values of $r_{\mathrm{\scriptsize{ex}}}$ and $L$ were found for a Halbach cylinder with the least possible volume of mag…
▽ More
In this paper we describe the results of a parameter survey of a 16 segmented Halbach cylinder in three dimensions in which the parameters internal radius, $r_{\mathrm{\scriptsize{in}}}$, external radius, $r_{\mathrm{\scriptsize{ex}}}$, and length, $L$, have been varied. Optimal values of $r_{\mathrm{\scriptsize{ex}}}$ and $L$ were found for a Halbach cylinder with the least possible volume of magnets with a given mean flux density in the cylinder bore. The volume of the cylinder bore could also be significantly increase by only slightly increasing the volume of the magnets, for a fixed mean flux density. Placing additional blocks of magnets on the end faces of the Halbach cylinder also improved the mean flux density in the cylinder bore, especially so for short Halbach cylinders with large $r_{\mathrm{\scriptsize{ex}}}$. Moreover magnetic cooling as an application for Halbach cylinders was considered. A magnetic cooling quality parameter, $Λ_{\mathrm{cool}}$, was introduced and results showed that this parameter was optimal for long Halbach cylinders with small $r_{\mathrm{\scriptsize{ex}}}$. Using the previously mentioned additional blocks of magnets can improve the parameter by as much as 15% as well as improve the homogeneity of the field in the cylinder bore.
△ Less
Submitted 8 September, 2014;
originally announced September 2014.
-
In situ characterization of delamination and crack growth of a CGO-LSM multi-layer ceramic sample investigated by X-ray tomographic microscopy
Authors:
R. Bjørk,
V. Esposito,
E. M. Lauridsen,
P. S. Jørgensen,
J. L. Fife,
K. B. Andersen,
S. P. V. Foghmoes,
N. Pryds
Abstract:
The densification, delamination and crack growth behavior in a Ce$_{0.9}$Gd$_{0.1}$O$_{1.95}$ (CGO) and (La$_{0.85}$Sr$_{0.15})_{0.9}$MnO$_{3}$ (LSM) multi-layer ceramic sample was studied using in situ X-ray tomographic microscopy (microtomography), to investigate the critical dynamics of crack propagation and delamination in a multilayered sample. Naturally occurring defects, caused by the sampl…
▽ More
The densification, delamination and crack growth behavior in a Ce$_{0.9}$Gd$_{0.1}$O$_{1.95}$ (CGO) and (La$_{0.85}$Sr$_{0.15})_{0.9}$MnO$_{3}$ (LSM) multi-layer ceramic sample was studied using in situ X-ray tomographic microscopy (microtomography), to investigate the critical dynamics of crack propagation and delamination in a multilayered sample. Naturally occurring defects, caused by the sample preparation process, are shown not to be critical in sample degradation. Instead defects are nucleated during the debinding step. Crack growth is significantly faster along the material layers than perpendicular to them, and crack growth and delamination only accelerates when sintering occurs.
△ Less
Submitted 8 September, 2014;
originally announced September 2014.
-
Analysis of the magnetic field, force, and torque for two-dimensional Halbach cylinders
Authors:
R. Bjørk,
C. R. H. Bahl,
A. Smith,
N. Pryds
Abstract:
The Halbach cylinder is a construction of permanent magnets used in applications such as nuclear magnetic resonance apparatus, accelerator magnets and magnetic cooling devices. In this paper the analytical expression for the magnetic vector potential, magnetic flux density and magnetic field for a two dimensional Halbach cylinder are derived. The remanent flux density of a Halbach magnet is charac…
▽ More
The Halbach cylinder is a construction of permanent magnets used in applications such as nuclear magnetic resonance apparatus, accelerator magnets and magnetic cooling devices. In this paper the analytical expression for the magnetic vector potential, magnetic flux density and magnetic field for a two dimensional Halbach cylinder are derived. The remanent flux density of a Halbach magnet is characterized by the integer $p$. For a number of applications the force and torque between two concentric Halbach cylinders are important. These quantities are calculated and the force is shown to be zero except for the case where $p$ for the inner magnet is one minus $p$ for the outer magnet. Also the force is shown never to be balancing. The torque is shown to be zero unless the inner magnet $p$ is equal to minus the outer magnet $p$. Thus there can never be a force and a torque in the same system.
△ Less
Submitted 5 September, 2014;
originally announced September 2014.
-
Magnetic cooling at Risoe DTU
Authors:
K. K. Nielsen,
R. Bjoerk,
J. B. Jensen,
C. R. H Bahl,
N. Pryds,
A. Smith,
A. Nordentoft,
J. Hattel
Abstract:
Magnetic refrigeration at room temperature is of great interest due to a long-term goal of making refrigeration more energy-efficient, less noisy and free of any environmentally hostile materials. A refrigerator utilizing an active magnetic regenerator (AMR) is based on the magnetocaloric effect, which manifests itself as a temperature change in magnetic materials when subjected to a varying mag…
▽ More
Magnetic refrigeration at room temperature is of great interest due to a long-term goal of making refrigeration more energy-efficient, less noisy and free of any environmentally hostile materials. A refrigerator utilizing an active magnetic regenerator (AMR) is based on the magnetocaloric effect, which manifests itself as a temperature change in magnetic materials when subjected to a varying magnetic field. In this work we present the current state of magnetic refrigeration research at Risoe DTU with emphasis on the numerical modeling of an existing AMR test machine. A 2D numerical heat-transfer and fluid-flow model that represents the experimental setup is presented. Experimental data of both no-heat load and heat load situations are compared to the model. Moreover, results from the numerical modeling of the permanent magnet design used in the system are presented.
△ Less
Submitted 4 February, 2009;
originally announced February 2009.
-
Exploring the Galaxy using space probes
Authors:
Rasmus Bjoerk
Abstract:
This paper investigates the possible use of space probes to explore the Milky Way, as a means both of finding life elsewhere in the Galaxy and as finding an answer to the Fermi paradox. I simulate exploration of the Galaxy by first examining how long time it takes a given number of space probes to explore 40,000 stars in a box from -300 to 300 pc above the Galactic thin disk, as a function of Ga…
▽ More
This paper investigates the possible use of space probes to explore the Milky Way, as a means both of finding life elsewhere in the Galaxy and as finding an answer to the Fermi paradox. I simulate exploration of the Galaxy by first examining how long time it takes a given number of space probes to explore 40,000 stars in a box from -300 to 300 pc above the Galactic thin disk, as a function of Galactic radius. I then model the Galaxy to consist of $\sim{}260,000$ of these 40,000 stellar systems all located in a defined Galactic Habitable Zone and show how long time it takes to explore this zone. The result is that with 8 probes, each with 8 subprobes $\sim{}4%$ of the Galaxy can be explored in $2.92\cdot{}10^{8}$ years. Increasing the number of probes to 200, still with 8 subprobes each, reduces the exploration time to $1.52\cdot{}10^{7}$ years.
△ Less
Submitted 23 April, 2007; v1 submitted 9 January, 2007;
originally announced January 2007.
-
Theoretical Uncertainties in Red Giant Branch Evolution: The Red Giant Branch Bump
Authors:
Stephan R. Bjork,
Brian Chaboyer
Abstract:
A Monte Carlo simulation exploring uncertainties in standard stellar evolution theory on the red giant branch of metal-poor globular clusters has been conducted. Confidence limits are derived on the absolute V-band magnitude of the bump in the red giant branch luminosity function (M_v,b) and the excess number of stars in thebump, R_b. The analysis takes into account uncertainties in the primordi…
▽ More
A Monte Carlo simulation exploring uncertainties in standard stellar evolution theory on the red giant branch of metal-poor globular clusters has been conducted. Confidence limits are derived on the absolute V-band magnitude of the bump in the red giant branch luminosity function (M_v,b) and the excess number of stars in thebump, R_b. The analysis takes into account uncertainties in the primordial helium abundance, abundance of alpha-capture elements, radiative and conductive opacities, nuclear reaction rates, neutrino energy losses, the treatments of diffusion and convection, the surface boundary conditions, and color transformations.
The uncertainty in theoretical values for the red giant bump magnitude varies with metallicity between +0.13/-0.12 mag at [Fe/H] = -2.4 and +0.23/-0.21 mag at [Fe/H] = -1.0$. The dominant sources of uncertainty are the abundance of the alpha-capture elements, the mixing length, and the low-temperature opacities. The theoretical values of M_v,b are in good agreement with observations. The uncertainty in the theoretical value of R_b is +/-0.01 at all metallicities studied. The dominant sources of uncertainty are the abundance of the alpha-capture elements, the mixing length, and the high-temperature opacities. The median value of R_b varies from 0.44 at [Fe/H] = -2.4$ to 0.50 at [Fe/H] = -1.0. These theoretical values for R_b are in agreement with observations.
△ Less
Submitted 20 December, 2005;
originally announced December 2005.