Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MADNESS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Automated selection of the optimal cardiac phase for single-beat coronary CT angiography reconstruction

Purpose: Reconstructing a low-motion cardiac phase is expected to improve coronary artery visualization in coronary computed tomography angiography (CCTA) exams. This study developed an automated algorithm for selecting the optimal cardiac phase for CCTA reconstruction. The algorithm uses prospectively gated, single-beat, multiphase data made possible by wide cone-beam imaging. The proposed algorithm differs from previous approaches because the optimal phase is identified based on vessel image quality (IQ) directly, compared to previous approaches that included motion estimation and interphase processing. Because there is no processing of interphase information, the algorithm can be applied to any sampling of image phases, making it suited for prospectively gated studies where only a subset of phases are available. Methods: An automated algorithm was developed to select the optimal phase based on quantitative IQ metrics. For each reconstructed slice at each reconstructed phase, an image quality metric was calculated based on measures of circularity and edge strength of through-plane vessels. The image quality metric was aggregated across slices, while a metric of vessel-location consistency was used to ignore slices that did not contain through-plane vessels. The algorithm performance was evaluated using two observer studies. Fourteen single-beat cardiac CT exams (Revolution CT, GE Healthcare, Chalfont St. Giles, UK) reconstructed at 2% intervals were evaluated for best systolic (1), diastolic (6), or systolic and diastolic phases (7) by three readers and the algorithm. Pairwise inter-reader and reader-algorithm agreement was evaluated using the mean absolute difference (MAD) and concordance correlation coefficient (CCC) between the reader and algorithm-selected phases. A reader-consensus best phase was determined and compared to the algorithm selected phase. In cases where the algorithm and consensus best phases differed by more than 2%, IQ was scored by three readers using a five point Likert scale. Results: There was no statistically significant difference between inter-reader and reader-algorithm agreement for either MAD or CCC metrics (p > 0.1). The algorithm phase was within 2% of the consensus phase in 15/21 of cases. The average absolute difference between consensus and algorithm best phases was 2.29% ± 2.47%, with a maximum difference of 8%. Average image quality scores for the algorithm chosen best phase were 4.01 ± 0.65 overall, 3.33 ± 1.27 for right coronary artery (RCA), 4.50 ± 0.35 for left anterior descending (LAD) artery, and 4.50 ± 0.35 for left circumflex artery (LCX). Average image quality scores for the consensus best phase were 4.11 ± 0.54 overall, 3.44 ± 1.03 for RCA, 4.39 ± 0.39 for LAD, and 4.50 ± 0.18 for LCX. There was no statistically significant difference (p > 0.1) between the image quality scores of the algorithm phase and the consensus phase. Conclusions: The proposed algorithm was statistically equivalent to a reader in selecting an optimal cardiac phase for CCTA exams. When reader and algorithm phases differed by >2%, image quality as rated by blinded readers was statistically equivalent. By detecting the optimal phase for CCTA reconstruction, the proposed algorithm is expected to improve coronary artery visualization in CCTA exams.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

SU-F-J-108: TMR Correction Factor Based Online Adaptive Radiotherapy for Stereotactic Radiosurgery (SRS) of L-Spine Tumors Using Cone Beam CT

Purpose: To investigate the feasibility of using TMR ratio correction factors for a fast online adaptive plan to compensate for anatomical changes in stereotactic radiosurgery (SRS) of L-spine tumors. Methods: Three coplanar treatment plans were made for 11 patients: Uniform (9 IMRT beams equally distributed around the patient); Posterior (IMRT with 9 posterior beams every 20 degree) and VMAT (2 360° arcs). For each patient, the external body and bowel gas were contoured on the planning CT and pre-treatment CBCT. After registering CBCT and the planning CT by aligning to the tumor, the CBCT contours were transferred to the planning CT. To estimate the actual delivered dose while considering patient’s anatomy of the treatment day, a hybrid CT was created by overriding densities in planning CT using the differences between CT and CBCT external and bowel gas contours. Correction factors (CF) were calculated using the effective depth information obtained from the planning system using the hybrid CT: CF = TMR (delivery)/TMR (planning). The adaptive plan was generated by multiplying the planned Monitor Units with the CFs. Results: The mean absolute difference (MAD) in V16Gy of the target between planned and estimated delivery with and without TMR correction was 0.8 ± 0.7% vs. 2.4 ± 1.3% for Uniform and 1.0 ± 0.9% vs. 2.6 ± 1.3% for VMAT plans(p<0.05), respectively. For V12Gy of cauda-equina with and without TMR correction, MAD was 0.24 ± 0.19% vs. 1.2 ± 1.02% for Uniform and 0.23 ± 0.20% vs. 0.78 ± 0.79% for VMAT plans(p<0.05), respectively. The differences between adaptive and original plans were not significant for posterior plans. Conclusion: The online adaptive strategy using TMR ratios and pre-treatment CBCT information was feasible strategy to compensate for anatomical changes for the patients treated for L-spine tumors, particularly for equally spaced IMRT and VMAT plans.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

A Public-Private-Acadmic Partnership to Advance Solar Power Forecasting

The National Center for Atmospheric Research (NCAR) is pleased to have led a partnership to advance the state-of-the-science of solar power forecasting by designing, developing, building, deploying, testing, and assessing the SunCast™ Solar Power Forecasting System. The project has included cutting edge research, testing in several geographically- and climatologically-diverse high penetration solar utilities and Independent System Operators (ISOs), and wide dissemination of the research results to raise the bar on solar power forecasting technology. The partners include three other national laboratories, six universities, and industry partners. This public-private-academic team has worked in concert to perform use-inspired research to advance solar power forecasting through cutting-edge research to advance both the necessary forecasting technologies and the metrics for evaluating them. The project has culminated in a year-long, full-scale demonstration of provide irradiance and power forecasts to utilities and ISOs to use in their operations. The project focused on providing elements of a value chain, beginning with the weather that causes a deviation from clear sky irradiance and progresses through monitoring and observations, modeling, forecasting, dissemination and communication of the forecasts, interpretation of the forecast, and through decision-making, which produces outcomes that have an economic value. The system has been evaluated using metrics developed specifically for this project, which has provided rich information on model and system performance. Research was accomplished on the very short range (0-6 hours) Nowcasting system as well as on the longer term (6-72 hour) forecasting system. The shortest range forecasts are based on observations in the field. The shortest range system, built by Brookhaven National Laboratory (BNL) and based on Total Sky Imagers (TSIs) is TSICast, which operates on the shortest time scale with a latency of only a few minutes and forecasts that currently go out to about 15 min. This project has facilitated research in improving the hardware and software so that the new high definition cameras deployed at multiple nearby locations allow discernment of the clouds at varying levels and advection according to the winds observed at those levels. Improvements over “smart persistence” are about 29% for even these very short forecasts. StatCast is based on pyranometer data measured at the site as well as concurrent meteorological observations and forecasts. StatCast is based on regime-dependent artificial intelligence forecasting techniques and has been shown to improve on “smart persistence” forecasts by 15-50%. A second category of short-range forecasting systems employ satellite imagery and use that information to discern clouds and their motion, allowing them to project the clouds, and the resulting blockage of irradiance, in time. CIRACast (the system produced by the Cooperative Institute for Atmospheric Research [CIRA] at Colorado State University) was already one of the more advanced cloud motion systems, which is the reason that team was brought to this project. During the project timeframe, the CIRA team was able to advance cloud shadowing, parallax removal, and implementation of better advecting winds at different altitudes. CIRACast shows generally a 25-40% improvement over Smart Persistence between sunrise and approximately 1600 UTC (Coordinated Universal Time) . A second satellite-based system, MADCast (Multi-sensor Advective Diffusive foreCast system), assimilates data from multiple satellite imagers and profilers to assimilate a fully three-dimensional picture of the cloud into the dynamic core of WRF. During 2015, MADCast (provided at least 70% improvement over Smart Persistence, with most of that skill being derived during partly cloudy conditions. That allows advection of the clouds via the Weather Research and Forecasting (WRF) model dynamics directly. After WRF-Solar™ showed initial success, it was also deployed in nowcasting mode with coarser runs out to 6 hours made hourly. It provided improvements on the order of 50-60% over Smart Persistence for forecasts up to 1600 UTC. The advantages of WRF-Solar-Nowcasting and MADCast were then blended to develop the new MAD-WRF model that incorporates the most important features of each of those models, both assimilating satellite cloud fields and using WRF-So far physics to develop and dissipate clouds. MAE improvements for MAD-WRF for forecasts from 3-6 hours are improved over WRF-Solar-Now by 20%. While all the Nowcasting system components by themselves provide improvement over Smart Persistence, the largest benefit is derived when they are smartly blended together by the Nowcasting Integrator to produce an integrated forecast. The development of WRF-Solar™ under this project has provided the first numerical weather prediction (NWP) model specifically designed to meet the needs of irradiance forecasting. The first augmentation improved the solar tracking algorithm to account for deviations associated with the eccentricity of the Earth’s orbit and the obliquity of the Earth. Second, WRF-Solar™ added the direct normal irradiance (DNI) and diffuse (DIF) components from the radiation parameterization to the model output. Third, efficient parameterizations were implemented to either interpolate the irradiance in between calls to the expensive radiative transfer parameterization, or to use a fast radiative transfer code that avoids computing three-dimensional heating rates but provides the surface irradiance. Fourth, a new parameterization was developed to improve the representation of absorption and scattering of radiation by aerosols (aerosol direct effect). A fifth advance is that the aerosols now interact with the cloud microphysics, altering the cloud evolution and radiative properties, an effect that has been traditionally only implemented in atmospheric computationally costly chemistry models. A sixth development accounts for the feedbacks that sub-grid scale clouds produce in shortwave irradiance as implemented in a shallow cumulus parameterization Finally, WRF-Solar™ also allows assimilation of infrared irradiances from satellites to determine the three dimensional cloud field, allowing for an improved initialization of the cloud field that increases the performance of short-range forecasts. We find that WRF-Solar™ can improve clear sky irradiance prediction by 15-80% over a standard version of WRF, depending on location and cloud conditions. In a formal comparison to the NAM baseline, WRF-Solar™ showed improvements in the Day-Ahead forecast of 22-42%. The SunCast™ system requires substantial software engineering to blend all of the new model components as well as existing publically available NWP model runs. To do this we use an expert system for the Nowcasting blender and the Dynamic Integrated foreCast (DICast®) system for the NWP models. These two systems are then blended, we use an empirical power conversion method to convert the irradiance predictions to power, then apply an analog ensemble (AnEn) approach to further tune the forecast as well as to estimate its uncertainty. The AnEn module decreased RMSE (root mean squared error) by 17% over the blended SunCast™ power forecasts and provided skill in the probabilistic forecast with a Brier Skill Score of 0.55. In addition, we have also developed a Gridded Atmospheric Forecast System (GRAFS) in parallel, leveraging cost share funds. An economic evaluation based on Production Cost Modeling in the Public Service Company of Colorado showed that the observed 50% improvement in forecast accuracy will save their customers $819,200 with the projected MW deployment for 2024. Using econometrics, NCAR has scaled this savings to a national level and shown that an annual expected savings for this 50% forecast error reduction ranges from $11M in 2015 to $43M expected in 2040 with increased solar deployment. This amounts to a $455M discounted savings over the 26 year period of analysis.

14 SOLAR ENERGY↗
DOE OSTI2017

The MARS Code System User’s Guide Version 15(2016)

This paper is a user’s guide to the current version of the MARS15 Monte Carlo code. MARS allows inclusive and, in many cases, exclusive simulations of three-dimensional hadronic and electromagnetic cascades and modeling of heavy ion, muon and low energy neutron-photon transport in accelerator, detector, spacecraft and shielding components. Particles can have energies from a fraction of an electronvolt up to about 100 TeV. The code has undergone substantial improvements since the last documented version MARS13(95) and since previous releases of MARS15; all these as well as other specific features of the MARS code system are explained in detail. Descriptions of general input and output with commentary and recommendations are given. Examples are given for running the program with distributed sources, complex compounds, arbitrary geometries, and magnetic fields. A built-in graphical-user interface allows the visualization and detailed analysis of the geometry description. Enhanced options are described for the MARS code being coupled with: the DPMJET hadron event generator, the MCNP code for low-energy neutron transport, the ANSYS system for thermal and stress analysis, physics analysis and graphics packages, the MAD accelerator lattice description via a universal MAD-MARS beam line builder, the STRUCT program for tracking particles in accelerator lattices with beam loss recording, the DETRA code for nuclide decay and transmutation calculations, the EGS5 code for precise modeling of low-energy electromagnetic showers, and the ROOT system for description and visualization of arbitrary geometries. Use of the code in a multistage mode, coupled with event generators (DPMJET), with the STRUCT program for tracking particles in accelerator lattices with beam loss recording, and with physics analysis and graphics packages is demonstrated with typical input and output examples.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Beam Optics Studies for Larger Beampipe Diameters for the LAMP-PSR

The proposed beampipe diameter increase of the PSR beamlines for the LAMP PSR upgrade is expected decrease beam loss due to tuning errors and misalignment. However, a larger beampipe diameter would require increasing the pole-to-pole gap height within all magnets in the PSR (dipole gap, quadrupole bore, etc.). This would in turn affect the extent and profile of the fringe fields at the ends of the magnets, altering the effective length of the magnets and thus the magnetic optics. In this work, we used the simulation code MAD-X to model the PSR lattice for different beampipe diameters. The results of these simulations will be later used to benchmark particle simulations using the code pyORBIT. The MAD-X simulations are described, and the simulation results are discussed in this technote.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Modelling Beam Loss Within the LANSCE Proton Storage Ring

Several upgrades are being considered for the proton storage ring (PSR) at the Los Alamos Neutron Science Center (LANSCE) to reduce beam loss and thereby reduce the cooldown period of the PSR following a beam run. First, we have considered an increased beam pipe diameter would reduce beam loss due to beam scraping caused by misalignments and tuning errors. However, this would require increased pole-to-pole gap height within the dipoles and quadrupoles, which would change their effective length and alter their fringe fields. The effect of different magnet gaps on the beam optical parameters and on beam loss was studied using the simulation codes MAD-X and PyORBIT. Second, we are developing a detailed particle tracking model within the framework of the simulation package General Particle Tracer (GPT) for the PSR H - stripping system. This model will be used to study the effect of the stripper foil parameters (position, composition, areal density, depth, etc.) on first-turn losses, where most of the observed beam loss and emittance growth occurs due to foil scattering, foil stripping, and Lorentz stripping. The model implements C++ custom elements created to model foil scattering, foil stripping, and Lorentz stripping, as these are not built-in features of GPT. The preliminary results of the MAD-X and PyORBIT simulations, as well as the GPT simulation model, is described in this tech note.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

The Sloan Digital Sky Survey Reverberation Mapping Project: An investigation of biases in C iv emission line properties

We investigate the dependence on data quality of quasar properties measured from the C iv emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch (SE) and coadded spectra, focusing on the C iv λ1549 emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ∼2–4), we find increasing systematic offsets with decreasing signal-to-noise ratio (S/N). The systematic difference (measurement uncertainty) in our lowest-S/N (<5) subsample between the SE and coadded spectrum (i) C iv equivalent width is 17 Å (31 Å), (ii) centroid wavelength is <1 Å (2 Å), and fractional velocity widths, ΔV/V, characterized by (iii) the line dispersion, σ {sub l}, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1σ measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the C iv FWHM velocity width and the C iv profile characterized by FWHM/σ {sub l} are only smaller than the statistical uncertainties when S/N > 10, although offsets in lower-S/N spectra exceed the statistical uncertainties by only a factor of ∼1.5 and may depend on the type of functional fit to the line. Characterizing the C iv line profile by the kurtosis is the least robust property investigated, as the median systematic coadded–SE measurement differences are larger than the statistical uncertainties for all S/N subsamples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Extending the Calibration of C iv-based Single-epoch Black Hole Mass Estimators for Active Galactic Nuclei

We provide an updated calibration of C iv λ1549 broad emission line–based single-epoch (SE) black hole (BH) mass estimators for active galactic nuclei (AGNs) using new data for six reverberation-mapped AGNs at redshift z=0.005--0.028 with BH masses (bolometric luminosities) in the range 10{sup 6.5}--10{sup 7.5} M{sub ⊙} (10{sup 41.7}--10{sup 43.8} erg s{sup −1}). New rest-frame UV-to-optical spectra covering 1150–5700 Å for the six AGNs were obtained with the Hubble Space Telescope (HST). Multicomponent spectral decompositions of the HST spectra were used to measure SE emission-line widths for the C iv, Mg ii, and Hβ lines, as well as continuum luminosities in the spectral region around each line. We combine the new data with similar measurements for a previous archival sample of 25 AGNs to derive the most consistent and accurate calibrations of the C iv-based SE BH mass estimators against the Hβ reverberation-based masses, using three different measures of broad-line width: full width at half maximum (FWHM), line dispersion (σ{sub line}), and mean absolute deviation (MAD). The newly expanded sample at redshift z=0.005--0.234 covers a dynamic range in BH mass (bolometric luminosity) of logM{sub BH}/M{sub ⊙}=6.5--9.1 (logL{sub bol}/ erg s{sup −1} = 41.7--46.9), and we derive the new C iv-based mass estimators using a Bayesian linear regression analysis over this range. We generally recommend the use of σ{sub line} or MAD rather than FWHM to obtain a less biased velocity measurement of the C iv emission line, because its narrow-line component contribution is difficult to decompose from the broad-line profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Blazar Variability from Turbulence in Jets Launched by Magnetically Arrested Accretion Flows

Blazars show variability on timescales ranging from minutes to years, the former being comparable to and in some cases even shorter than the light-crossing time of the central black hole. The observed γ -ray light curves can be described by a power-law power density spectrum (PDS), with a similar index for both BL Lacs and flat-spectrum radio quasars. We show that this variability can be produced by turbulence in relativistic jets launched by magnetically arrested accretion flows (MADs). We perform radiative transport calculations on the turbulent, highly magnetized jet launching region of a MAD with a rapidly rotating supermassive black hole. The resulting synchrotron and synchrotron self-Compton emission, originating from close to the black hole horizon, is highly variable. This variability is characterized by PDS, which is remarkably similar to the observed power-law spectrum at frequencies less than a few per day. Furthermore, turbulence in the jet launching region naturally produces fluctuations in the plasma on scales much smaller than the horizon radius. We speculate that similar turbulent processes, operating in the jet at large radii (and therefore a high bulk Lorentz factor), are responsible for blazar variability over many decades in frequency, including on minute timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Winds and Disk Turbulence Exert Equal Torques on Thick Magnetically Arrested Disks

Abstract The conventional accretion disk lore is that magnetized turbulence is the principal angular momentum transport process that drives accretion. However, when dynamically important large-scale magnetic fields thread an accretion disk, they can produce mass and angular momentum outflows, known as winds,that also drive accretion. Yet, the relative importance of turbulent and wind-driven angular momentum transport is still poorly understood. To probe this question, we analyze a long-duration (1.2 × 10 5 r g /c) simulation of a rapidly rotating (a= 0.9) black hole feeding from a thick (H/r∼ 0.3), adiabatic, magnetically arrested disk (MAD), whose dynamically important magnetic field regulates mass inflow and drives both uncollimated and collimated outflows (i.e., winds and jets, respectively). By carefully disentangling the various angular momentum transport processes within the system, we demonstrate the novel result that disk winds and disk turbulence both extract roughly equal amounts of angular momentum from the disk. We find cumulative angular momentum and mass accretion outflow rates of L ̇ ∝ r 0.9 and M ̇ ∝ r 0.4 , respectively. This result suggests that understanding both turbulent and laminar stresses is key to understanding the evolution of systems where geometrically thick MADs can occur, such as the hard state of X-ray binaries, low-luminosity active galactic nuclei, some tidal disruption events, and possibly gamma-ray bursts.

Astronomy & Astrophysics↗
DOE OSTI2022

First Sagittarius A* Event Horizon Telescope Results. V. Testing Astrophysical Models of the Galactic Center Black Hole

In this paper we provide a first physical interpretation for the Event Horizon Telescope's (EHT) 2017 observations of Sgr A*. Our main approach is to compare resolved EHT data at 230 GHz and unresolved non-EHT observations from radio to X-ray wavelengths to predictions from a library of models based on time-dependent general relativistic magnetohydrodynamics simulations, including aligned, tilted, and stellar-wind-fed simulations; radiative transfer is performed assuming both thermal and nonthermal electron distribution functions. We test the models against 11 constraints drawn from EHT 230 GHz data and observations at 86 GHz, 2.2μm, and in the X-ray. All models fail at least one constraint. Light-curve variability provides a particularly severe constraint, failing nearly all strongly magnetized (magnetically arrested disk (MAD)) models and a large fraction of weakly magnetized models. A number of models fail only the variability constraints. We identify a promising cluster of these models, which are MAD and have inclination i ≤ 30°. They have accretion rate (5.2–9.5) × 10$^{−9}$ M $_{⊙}$ yr$^{−1}$, bolometric luminosity (6.8–9.2) × 10$^{35}$ erg s$^{−1}$, and outflow power (1.3–4.8) × 10$^{38}$ erg s$^{−1}$. We also find that all models with i ≥ 70° fail at least two constraints, as do all models with equal ion and electron temperature; exploratory, nonthermal model sets tend to have higher 2.2μm flux density; and the population of cold electrons is limited by X-ray constraints due to the risk of bremsstrahlung overproduction. Finally, we discuss physical and numerical limitations of the models, highlighting the possible importance of kinetic effects and duration of the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Unified Picture of Short and Long Gamma-Ray Bursts from Compact Binary Mergers

Abstract The recent detections of the ∼10 s longγ-ray bursts (GRBs) 211211A and 230307A followed by softer temporally extended emission (EE) and kilonovae point to a new GRB class. Using state-of-the-art first-principles simulations, we introduce a unifying theoretical framework that connects binary neutron star (BNS) and black hole–NS (BH–NS) merger populations with the fundamental physics governing compact binary GRBs (cbGRBs). For binaries with large total masses,M tot ≳ 2.8M ⊙ , the compact remnant created by the merger promptly collapses into a BH surrounded by an accretion disk. The duration of the pre-magnetically arrested disk (MAD) phase sets the duration of the roughly constant power cbGRB and could be influenced by the disk mass,M d . We show that massive disks (M d ≳ 0.1M ⊙ ), which form for large binary mass ratiosq≳ 1.2 in BNS orq≲ 3 in BH–NS mergers, inevitably produce 211211A-like long cbGRBs. Once the disk becomes MAD, the jet power drops with the mass accretion rate as M ̇ ∼ t − 2 , establishing the EE decay. Two scenarios are plausible for short cbGRBs. They can be powered by BHs with less massive disks, which form for otherqvalues. Alternatively, for binaries withM tot ≲ 2.8M ⊙ , mergers should go through a hypermassive NS (HMNS) phase, as inferred for GW170817. Magnetized outflows from such HMNSs, which typically live for ≲1 s, offer an alternative progenitor for short cbGRBs. The first scenario is challenged by the bimodal GRB duration distribution and the fact that the Galactic BNS population peaks at sufficiently low masses that most mergers should go through an HMNS phase.

Astronomy & Astrophysics↗
DOE OSTI2023

Synchrotron Polarization Signatures of Surface Waves in Supermassive Black Hole Jets

Abstract Supermassive black holes in active galactic nuclei are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of nonthermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched from the vicinity of the black hole and propagate along the interface between the jet and the wind. At this interface, a steep gradient in velocity is present between the mildly relativistic wind and the highly relativistic jet. The interface is, therefore, a shear layer, and due to the shear, the waves generate roll-ups that alter the magnetic field configuration and the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The telltale polarization signatures of the waves could be observable by future very long baseline interferometric observations, e.g., the next-generation Event Horizon Telescope.

Astronomy & Astrophysics↗
DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗
NASA NTRS1988

Computer modeling of polymers

A Polymer Molecular Analysis Display System (p-MADS) was developed for computer modeling of polymers. This method of modeling allows for the theoretical calculation of molecular properties such as equilibrium geometries, conformational energies, heats of formations, crystal packing arrangements, and other properties. Furthermore, p-MADS has the following capabilities: constructing molecules from internal coordinates (bonds length, angles, and dihedral angles), Cartesian coordinates (such as X-ray structures), or from stick drawings; manipulating molecules using graphics and making hard copy representation of the molecules on a graphics printer; and performing geometry optimization calculations on molecules using the methods of molecular mechanics or molecular orbital theory.

Green, Terry J.↗
NASA NTRS1989

Mesoscale acid deposition modeling studies

The work performed in support of the EPA/DOE MADS (Mesoscale Acid Deposition) Project included the development of meteorological data bases for the initialization of chemistry models, the testing and implementation of new planetary boundary layer parameterization schemes in the MASS model, the simulation of transport and precipitation for MADS case studies employing the MASS model, and the use of the TASS model in the simulation of cloud statistics and the complex transport of conservative tracers within simulated cumuloform clouds. The work performed in support of the NASA/FAA Wind Shear Program included the use of the TASS model in the simulation of the dynamical processes within convective cloud systems, the analyses of the sensitivity of microburst intensity and general characteristics as a function of the atmospheric environment within which they are formed, comparisons of TASS model microburst simulation results to observed data sets, and the generation of simulated wind shear data bases for use by the aviation meteorological community in the evaluation of flight hazards caused by microbursts.

Kaplan, Michael L.↗
NASA NTRS2011

Thrust Stand Measurements of the Microwave Assisted Discharge Inductive Plasma Accelerator

Pulsed inductive plasma thrusters [1-3] are spacecraft propulsion devices in which electrical energy is capacitively stored and then discharged through an inductive coil. This type of pulsed thruster is electrodeless, with a time-varying current in the coil interacting with a plasma covering the face of the coil to induce a plasma current. Propellant is accelerated and expelled at a high exhaust velocity (O(10-100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, pulsed inductive plasma thrusters require high pulse energies to inductively ionize propellant. The Microwave Assisted Dis- charge Inductive Plasma Accelerator (MAD-IPA), shown in Fig. 1, is a pulsed inductive plasma thruster that addressees this issue by partially ionizing propellant inside a conical inductive coil before the main current pulse via an electron cyclotron resonance (ECR) discharge. The ECR plasma is produced using microwaves and a static magnetic field from a set of permanent magnets arranged to create a thin resonance region along the inner surface of the coil, restricting plasma formation, and in turn current sheet formation, to a region where the magnetic coupling between the plasma and the theta-pinch coil is high. The use of a conical theta-pinch coil also serves to provide neutral propellant containment and plasma plume focusing that is improved relative to the more common planar geometry of the Pulsed Inductive Thruster (PIT) [1, 2]. In this paper, we describe thrust stand measurements performed to characterize the performance (specific impulse, thrust efficiency) of the MAD-IPA thruster. Impulse data are obtained at various pulse energies, mass flow rates and inductive coil geometries. Dependencies on these experimental parameters are discussed in the context of the current sheet formation and electromagnetic plasma acceleration processes.

Hallock, Ashley K.↗
NASA NTRS2011

Thrust Stand Measurements Using Alternative Propellants in the Microwave Assisted Discharge Inductive Plasma Accelerator

Storable propellants (for example water, ammonia, and hydrazine) are attractive for deep space propulsion due to their naturally high density at ambient interplanetary conditions, which obviates the need for a cryogenic/venting system. Water in particular is attractive due to its ease of handling and availability both terrestrially and extra-terrestrially. While many storable propellants are reactive and corrosive, a propulsion scheme where the propellant is insulated from vulnerable (e.g. metallic) sections of the assembly would be well-suited to process these otherwise incompatible propellants. Pulsed inductive plasma thrusters meet this criterion because they can be operated without direct propellant-electrode interaction. During operation of these devices, electrical energy is capacitively stored and then discharged through an inductive coil creating a time-varying current in the coil that interacts with a plasma covering the face of the coil to induce a plasma current. Propellant is accelerated and expelled at a high exhaust velocity (O(10-100 km/s)) by the Lorentz body force arising from the interaction of the magnetic field and the induced plasma current. While this class of thruster mitigates the life-limiting issues associated with electrode erosion, many pulsed inductive plasma thrusters require high pulse energies to inductively ionize propellant. The Microwave Assisted Discharge Inductive Plasma Accelerator (MAD-IPA) is a pulsed inductive plasma thruster that addressees this issue by partially ionizing propellant inside a conical inductive coil before the main current pulse via an electron cyclotron resonance (ECR) discharge. The ECR plasma is produced using microwaves and a static magnetic field from a set of permanent magnets arranged to create a thin resonance region along the inner surface of the coil, restricting plasma formation, and in turn current sheet formation, to a region where the magnetic coupling between the plasma and the theta-pinch coil is high. The use of a conical theta-pinch coil also serves to provide neutral propellant containment and plasma plume focusing that is improved relative to the more common planar geometry of the Pulsed Inductive Thruster (PIT). In this paper, we describe thrust stand measurements performed to evaluate the specific impulse and thrust efficiency of the MAD-IPA for a variety of propellants. Propellants tested include both widely-used, non-reactive noble gases like argon, and rarely-used propellants such as water, hydrazine and ammonia. Dependencies of impulse data on propellant species are discussed in the context of the current sheet formation and electromagnetic plasma acceleration processes.

Hallock, Ashley K.↗