Engineering topics

Cook, William

Publications and source records attributed to Cook, William.

7 records

DOE OSTI2023

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Long-term general relativistic magnetohydrodynamics simulations of magnetic field in isolated neutron stars

ABSTRACT Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

Space Qualification of the Optical Filter Assemblies for the ICESat-2/ATLAS Instrument

The Advanced Topographic Laser Altimeter System (ATLAS) will be the only instrument on the Ice, Cloud, and Land Elevation Satellite -2 (ICESat-2). ICESat-2 is the 2nd-generation of the orbiting laser altimeter ICESat, which will continue polar ice topography measurements with improved precision laser-ranging techniques. In contrast to the original ICESat design, ICESat-2 will use a micro-pulse, multi-beam approach that provides dense cross-track sampling to help scientists determine a surface's slope with each pass of the satellite. The ATLAS laser will emit visible, green laser pulses at a wavelength of 532 nm and a rate of 10 kHz and will be split into 6 beams. A set of six identical, thermally-tuned etalon filter assemblies will be used to remove background solar radiation from the collected signal while transmitting the laser light to the detectors. A seventh etalon assembly will be used to monitor the laser center wavelength during the mission. In this paper, we present the design and optical performance measurements of the ATLAS optical filter assemblies (OFA) in air and in vacuum before integration on the ATLAS instrument.

Etalons↗

NASA NTRS2012

First Assessments of Predicted ICESat-2 Performance Using Aircraft Data

The Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) is a next-generation laser altimeter designed to continue key observations of ice sheet elevation change, sea ice freeboard, vegetation canopy height, earth surface elevation, and sea surface height. Scheduled for launch in mid-2016, ICESat-2 will use a high repetition rate (10 kHz), small footprint (10 m nominal ground diameter) laser, and a single-photon-sensitive detection strategy (photon counting) to measure precise range to the earth's surface. Using green light (532 nm), the six beams of ICESat-2 will provide improved spatial coverage compared with the single beam of ICESat, while the differences in transmit energy among the beams provide a large dynamic range. The six beams are arranged into three pairs of beams which allow slopes to measured on an orbit-by-orbit basis. In order to evaluate models of predicted ICESat-2 performance and provide ICESat-2-like data for algorithm development, an airborne ICESat-2 simulator was developed and first flown in 2010. This simulator, the Multiple Altimeter Beam Experimental Lidar (MABEL) was most recently deployed to Iceland in April 2012 and collected approx 85 hours of science data over land ice, sea ice, and calibration targets. MABEL uses a similar photon-counting measurement strategy to what will be used on ICESat-2. MABEL collects data in 16 green channels and an additional 8 channels in the infrared aligned across the direction of flight. By using NASA's ER-2 aircraft flying at 20km altitude, MABEL flies as close to space as is practical, and collects data through approx 95% of the atmosphere. We present background on the MABEL instrument, and data from the April 2012 deployment to Iceland. Among the 13 MABEL flights, we collected data over the Greenland ice sheet interior and outlet glaciers in the southwest and western Greenland, sea ice data over the Nares Strait and Greenland Sea, and a number of small glaciers and ice caps in Iceland and Svalbard. Several of the flights were coincident in time and space with NASA's Operation IceBridge, which provides an independent data set for validation. MABEL also collected data along CryoSat track 10482 in north central Greenland approximately one month after CryoSat passed overhead.

Neumann, Thomas↗

NASA NTRS1999

Development of the Double Etalon Fabry-Perot Interferometer for Determining Total and Tropospheric Ozone Concentrations

Measuring and understanding the distribution of ozone through the lower levels of Earth's atmosphere are high priorities in global change and climate research. Of particular interest now is the global distribution of ozone in the upper troposphere and lower stratosphere. Global coverage of the stratospheric ozone is feasible only via remote sensing instruments on a space-based platform. And though extensive monitoring tropospheric ozone is possible using instruments flown aboard conventional aircraft, a space-based system would be significantly less costly and provide information over a much broader area and produce more uniform coverage. Here we describe the prototype of an instrument being developed to monitor, from an orbiting spacecraft, the ozone found in Earth's upper troposphere and lower stratosphere. Our new spectrometer is an infrared Fabry-Perot interferometer which uses two synchrounously tuned etalons: a high resolution narrow band device and a lower resolution broader band filtering etalon. The prototype is a scanning device making use of nearly collimated input radiation and a single element detector. As presently configured, it is capable of providing a resolution better than 0.07/cm with a spectral band width approximately 5/cm wide and centered at 1054.7/cm. For the future space-based emission device a modification of the the prototype was to be made to employ innovative circle-to-line detector optics, those developed or in development at UM/SPRL, and a focal plane array detector. These enhancements would enable a simultaneous recording of the entire spectral range of interest, but with simple detection electronics and a significant gain in signal-to-noise over that of the scanning version.

Cook, William↗