Engineering topics

Morales, J. A.

Publications and source records attributed to Morales, J. A..

4 records

DOE OSTI2024

Limiting Rotation Rate of Neutron Stars from Crust Breaking and Gravitational Waves

Abstract Neutron stars are not observed to spin faster than about half their breakup rate. This limiting rotational frequency may be related to the strength of their crusts. As a star spins up from accretion, centrifugal forces stress the crust. We perform finite-element simulations of rotating neutron stars and find that the crust fails at rotation rates about half the breakup rate. Given uncertainties in microphysics, we have not determined the crust configuration after this failure. Instead, we argue that the crust may fail in an asymmetric way and could produce a configuration with a significant ellipticity (fractional difference in moments of inertia). If the ellipticity is large, a rotating star will radiate gravitational waves that may limit further spin-up. These stars may be promising sources for LIGO/VIRGO and next-generation gravitational-wave detectors.

Morales, J. A.↗

DOE OSTI2024

Anisotropic neutron star crust, solar system mountains, and gravitational waves

"Mountains" or non-axisymmetric deformations of rotating neutron stars (NS) efficiently radiate gravitational waves (GW). We consider analogies between NS mountains and surface features of solar system bodies. Both NS and moons such as Europa or Enceladus have thin crusts over deep oceans while Mercury has a thin crust over a large metallic core. Thin sheets may wrinkle in universal ways. Europa has linear features, Enceladus has "Tiger" stripes, and Mercury has lobate scarps. NS may have analogous features. The innermost inner core of the Earth is anisotropic with a shear modulus that depends on direction. If NS crust material is also anisotropic this will produce an ellipticity, when the crust is stressed, that grows with spin frequency. This yields a braking index (log derivative of spin down rate assuming only GW spin down) very different from n=5 and could explain the maximum spin observed for neutron stars and a possible minimum ellipticity of millisecond pulsars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Neutron star crust can support a large ellipticity

ABSTRACT Non-axisymmetrical deformations of the crust on rapidly rotating neutron stars are one of the main targets of searches for continuous gravitational waves. The maximum ellipticity, or fractional difference in moments of inertia, which can be supported by deformations of the crust (known as ‘mountains’), provides an important upper limit on the strength of these continuous gravitational wave sources. We use the formalism developed by Gittins and Andersson, along with a deforming force that acts mainly in the transverse direction, to obtain a maximum ellipticity of 7.4 × 10−6. This is larger than the original results that Gittins and Andersson obtained but consistent with earlier calculations by Ushomirsky, Cutler, and Bildsten. This suggests that rotating neutron stars could be strong sources of continuous gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗