Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 127 records · Page 7

Black hole-neutron star binaries near neutron star disruption limit in the mass regime of event GW230529

In May 2023, the LIGO-Virgo-KAGRA (LVK) Collaboration reported the likely black hole-neutron star (BHNS) merger GW230529_181500. The signal was observed with high significance in only one detector, limiting constraints on the black hole spin and motivating our study of disruption in this mass regime. That event is expected to be the merger of a 2.5–4.5 M ⊙ primary with a secondary compact object of mass between 1.2–2.0 M ⊙ . This makes it the first BHNS merger with a significant potential for the production of electromagnetic (EM) counterparts, and provides further evidence for compact objects existing within the suspected lower mass gap. To produce post-merger EM transients, the component of the black hole spin aligned with the orbital angular momentum must be sufficiently high, allowing the neutron star to be tidally disrupted. The disrupting BHNS binary may then eject a few percent of a solar mass of matter, leading to an observable kilonova driven by radioactive decays in ejecta, and/or a compact-binary gamma-ray burst (cbGRB) resulting from the formation of an accretion disk and relativistic jet. Determining which mergers lead to disruption of the neutron star is necessary to predict the prevalence of EM signals from BHNS mergers, yet most BHNS simulations so far have been performed far from the minimum spin required for tidal disruption. Here, we use the Spectral Einstein Code to explore the behavior of BHNS mergers in a mass range consistent with GW230529_181500 close to that critical spin, and compare our results against the mass remnant model currently used by the LVK Collaboration to predict the probability of tidal disruption. Furthermore, our numerical results reveal the emergence of non-zero accretion disks even below the predicted NS disruption limit, of low mass but capable of powering cbGRBs. Our results also demonstrate that the remnant mass model underpredicts the disk mass for the DD2 equation of state, while they are within expected modeling errors for SFHo. The disruption limit itself, however, is not found to significantly differ from the predictions of the analytical model, unless remnant masses M rem ≲ 0.001 M ⊙ prove interesting observationally. In all of our simulations, any kilonova signal would be dim and most likely dominated by post-merger disk outflows.

GW230529_181500↗

Investigating the impact of higher-order phase transitions in binary neutron-star mergers

In this paper we investigate quark deconfinement in neutrons stars and their mergers, focusing on the effects of higher orders for the phase transition between hadronic and quark matter. The different descriptions we use to describe matter microscopically contain varying particle degrees of freedom, including nucleons, hyperons, Delta baryons, and light and strange quarks. We use tabulated equations of state from the CompOSE database in which the quark deconfinement phase transition is described as being first order, and then smooth it out by introducing a percolation, replacing the single first-order phase transition with two transitions of second or third order. We then perform binary neutron-star merger simulations using these new equations of state, focusing on groups of binaries with the same single-star mass, radius, and tidal deformability, but different equations of state. Here, we go on to discuss differences in their evolution, and the ramifications for interpreting future gravitational wave observations and the potential to learn about dense matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Isospin-symmetry implications for nuclear two-body distributions and short-range correlations

Here, we study the implications of isospin symmetry in nuclear systems on two-body distributions and derive relations between different densities or momentum distributions of a given nucleus or of different nuclei in the same isospin multiplet. A connection between neutrinoless double beta decay transition of mirror nuclei and two-body densities is also obtained. The relations are numerically verified by means of quantum Monte Carlo calculations. We discuss the relevance of these results to the analysis of nuclear short-range correlations and also show that isospin symmetry allows us to extract information about the spectator nucleons when a short-range correlated pair is formed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accurate relativistic observables from postprocessing light cone catalogs

We introduce and study a new scheme to construct relativistic observables from postprocessing light cone data. This construction is based on a novel approach, LC-metric, which takes general light cone or snapshot output generated by arbitrary N-body simulations or emulations and solves the linearized Einstein equations to determine the spacetime metric on the light cone. We find that this scheme is able to determine the metric to high precision, and subsequently generate accurate mock cosmological observations sensitive to effects such as post-Born lensing and nonlinear integrated Sachs-Wolfe contributions. By comparing to conventional methods in quantifying those general relativistic effects, we show that this scheme is able to accurately construct the lensing convergence signal. We also find that the accuracy of this method in quantifying the integrated Sachs-Wolfe effects in the highly nonlinear regime outperforms conventional methods by an order of magnitude. Furthermore, this scheme opens a new path for exploring and modeling higher-order and nonlinear general relativistic contributions to cosmological observables, including mock observations of gravitational lensing and the moving lens and Rees-Sciama effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Gravitational waves from binary neutron star mergers with a spectral equation of state

In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

Cocoon shock breakout emission from binary neutron star mergers

Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.

79 ASTRONOMY AND ASTROPHYSICS↗

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Distinction of Electron Dispersion in Time-Resolved Photoemission Spectroscopy

While recent experiments provided compelling evidence for an intricate dependence of attosecond photoemission-time delays on the solid’s electronic band structure, the extent to which electronic transport and dispersion in solids can be imaged in time-resolved photoelectron (PE) spectra remains poorly understood. Emphasizing the distinction between photoemission time delays measured with two-photon, two-color interferometric spectroscopy, and transport times, we demonstrate how the effect of energy dispersion in the solid on photoemission delays can, in principle, be observed in interferometric photoemission. In this work we reveal analytically a scaling relation between the PE transport time in the solid and the observable photoemission delay and confirm this relation in numerical simulations for a model system. We trace photoemission delays to the phase difference the PE accumulates inside the solid and, in particular, predict negative photoemission delays. Based on these findings, we suggest a novel time-domain interferometric solid-state energy-momentum-dispersion imaging method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dust stars in the minimal exponential measure model

Here, we report the existence of horizonless compact object solutions supported by dust in the minimal exponential measure (MEMe) model, a theory which modifies the couplings between gravity and matter without introducing dynamical degrees of freedom. For a perfect fluid source, the field equations for the MEMe model can be rewritten as the Einstein field equations sourced by a perfect fluid with a transformed equation of state, which can endow a sufficiently dense cloud of dust with an effective pressure. The resulting dust-supported horizonless compact objects can have masses below ∼10 −11 ⁢𝑀 ⊙ , making them suitable as massive compact halo objects comprising a significant mass fraction for dark matter. A necessary condition for the existence of these compact object solutions is that the single free parameter in the MEMe model is positive valued. Additionally, we find that this positive sign for the parameter can provide a mechanism for suppressing the formation of (primordial) black holes from the gravitational collapse of matter below a certain mass scale.

alternative gravity theories↗

A Survey of Thermal Sensing Application in Additive Manufacturing

Industry 4.0 represents a significant change in manufacturing with the inclusion of such technologies as additive manufacturing and data analytics. In support of these changes, infrared thermal imaging (IR) and very-near infrared (VNIR) have proven invaluable in numerous applications related to additive manufacturing. In this review, these applications will be presented as they relate to polymer and metal-based additive manufacturing processes. Applications including in-situ characterization, post-build analysis, experimental setups, modeling, and machine learning will be covered. With a general overview complete, the review will conclude with a discussion of challenges identified throughout the literature and thoughts on future trends.

Borish, Michael↗

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Radiometric Testing of Germicidal UV Products, Round 1: UV-C Towers and Whole-Room Luminaires (CALiPER Summary Report)

This summary report analyzes the independently tested performance of 13 germicidal ultraviolet (GUV) products purchased between February and July 2022. A companion full report (DOE 2023) provides additional information and discussion of the tested products, test methods, and results. The products were of three different types: • Seven portable, consumer-oriented GUV towers designed to be placed on the floor or a desk of an unoccupied room to disinfect air and surfaces. Five of these products used LED sources and two products had low-pressure mercury (LPM) sources. • One GUV whole-room luminaire designed to be installed on a ceiling to disinfect air when a room is occupied. This product had LED sources. • Five GUV troffer or high-bay style whole-room luminaires designed to be installed in or suspended from a ceiling to disinfect air and surfaces when a room is unoccupied. All five had LPM sources. Product testing covered radiometric and electrical performance for all 13 products as well as photobiological safety evaluation if product documentation included testable claims. Measurement results enable comparison between products and against manufacturer or vendor claims. Testing identified numerous issues related to the accuracy of claimed GUV product performance. Claims were often untestable, contradictory, ambiguous, or used incorrect units and/or terminology. When claims were testable, they often did not match test results. For example, three LED products that claimed to emit UV-C emitted only UV-A. Product claim issues were more common among consumer-oriented tower products, but all product types exhibited problems with accurate performance claims. The UV-C radiant efficiency (calculated as UV-C output power divided by electrical input power) of the products varied widely, even among similar products using the same source technologies. For example, the UV-C radiant efficiency of LPM products varied by greater than a factor of three for the same product type, indicating a large potential energy savings opportunity for products that are better designed for efficiency. LED products had orders-of-magnitude lower UV-C radiant efficiency than LPM products. This study also identified several testing challenges and limitations. Most significant among these is the capability to accurately test and report the performance of larger GUV products. Whereas integrating spheres are used to quickly measure total radiant flux (i.e., output power) and spectral distribution, goniometers are used to measure radiant intensity distribution (from which radiant flux can be calculated). Integrating spheres require a specialized and costly coating to test UV, and the testing laboratory for this round of products had only a 20-inch diameter hemisphere with this capability. The integrating sphere accommodated just 2 of the 10 UV-C emitting products. Goniometer testing had a different size limitation in that mirrors typically used to increase goniometer test distance to the far field reflect little to no UV. As a result, the study evaluated only 6 of 13 products in the far field. Electronic files of UV-C intensity data for the other 7 products, which would typically be imported into design software for designing GUV applications, may not be reliable for predicting irradiance at arbitrary far-field distances (IES 2022a; CIE 2020). Specifiers and buyers of GUV products need accurate performance claims and data to deploy GUV technology safely and effectively. This CALiPER GUV Round 1 report demonstrates the significant education and training manufacturers and vendors still require to accurately test and report the performance of their GUV products. Further industry standards and guidelines may address testing limitations and improve test methods, product performance, and the accuracy of performance claims.

42 ENGINEERING↗

Radiometric Testing of Germicidal UV Products, Round 1: UV-C Towers and Whole-Room Luminaires (CALiPER Full Report)

This report analyzes the independently tested performance of 13 germicidal ultraviolet (GUV) products purchased between February and July 2022. The products were of three different types: • Seven portable, consumer-oriented GUV towers designed to be placed on the floor or a desk of an unoccupied room to disinfect air and surfaces. Five of these products used LED sources and two products had low-pressure mercury (LPM) sources. • One GUV whole-room luminaire designed to be installed on a ceiling to disinfect air when a room is occupied. This product had LED sources. • Five GUV troffer or high-bay style whole-room luminaires designed to be installed in or suspended from a ceiling to disinfect air and surfaces when a room is unoccupied. All five had LPM sources. Product testing covered radiometric and electrical performance for all 13 products as well as photobiological safety evaluation if product documentation included testable claims. Measurement results enable comparison between products and against manufacturer or vendor claims. Testing identified numerous issues related to the accuracy of claimed GUV product performance. Claims were often untestable, contradictory, ambiguous, or used incorrect units and/or terminology. When claims were testable, they often did not match test results. For example, three LED products that claimed to emit UV-C emitted only UV-A. Product claim issues were more common among consumer-oriented tower products, but all product types exhibited problems with accurate performance claims. The UV-C radiant efficiency (calculated as UV-C output power divided by electrical input power) of the products varied widely, even among similar products using the same source technologies. For example, the UV-C radiant efficiency of LPM products varied by greater than a factor of three for the same product type, indicating a large potential energy savings opportunity for products that are better designed for efficiency. LED products had orders-of-magnitude lower UV-C radiant efficiency than LPM products. This study also identified several testing challenges and limitations. Most significant among these is the capability to accurately test and report the performance of larger GUV products. Whereas integrating spheres are used to quickly measure total radiant flux (i.e., output power) and spectral distribution, goniometers are used to measure radiant intensity distribution (from which radiant flux can be calculated). Integrating spheres require a specialized and costly coating to test UV, and the testing laboratory for this round of products had only a 20-inch diameter hemisphere with this capability. The integrating sphere accommodated just 2 of the 10 UV-C emitting products. Goniometer testing had a different size limitation in that mirrors typically used to increase goniometer test distance to the far field reflect little to no UV. As a result, the study evaluated only 6 of 13 products in the far field. Electronic files of UV-C intensity data for the other 7 products, which would typically be imported into design software for designing GUV applications, may not be reliable for predicting irradiance at arbitrary far-field distances (IES 2022a; CIE 2020). Specifiers and buyers of GUV products need accurate performance claims and data to deploy GUV technology safely and effectively. This CALiPER GUV Round 1 report demonstrates the significant education and training manufacturers and vendors still require to accurately test and report the performance of their GUV products. Further industry standards and guidelines may address testing limitations and improve test methods, product performance, and the accuracy of performance claims.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

3′ RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

3’ RNA-seq↗

Interpreting Mass and Radius Measurements of Neutron Stars with Dark Matter Halos

Abstract The high densities of neutron stars (NSs) could provide astrophysical locations for dark matter (DM) to accumulate. Depending on the DM model, these DM admixed NSs (DANSs) could have significantly different properties than pure baryonic NSs, accessible through X-ray observations of rotation-powered pulsars. We adopt the two-fluid formalism in general relativity to numerically simulate stable configurations of DANSs, assuming a fermionic equation of state (EOS) for the DM with repulsive self-interaction. The distribution of DM in the DANS as a halo affects the path of X-rays emitted from hot spots on the visible baryonic surface, causing notable changes in the pulse profile observed by telescopes such as NICER, compared to pure baryonic NSs. We explore how various DM models affect the DM mass distribution, leading to different types of dark halos. We quantify the deviation in observed X-ray flux from stars with each of these halos. We identify the pitfalls in interpreting mass and radius measurements of NSs inferred from electromagnetic radiation and constraining the baryonic matter EOS if these dark halos exist.

Shawqi, Shafayat (ORCID:0000000210956183)↗

Biobased Flame Retardants Towards Sustainable Building Materials with Low Embodied Carbon

Incorporation of fire retardants in different types of products is one of the essential steps in the material production process to minimize fire risk and meet fire safety requirements. A variety of commonly used flame retardants based on halogen, mineral, and other compounds has gained popularity due to their efficient flame-retardant behaviors. However, especially in the case of halides, there are numerous toxicity related issues and environmental pollution effects that have urged the building sector to deviate from their use and focus on the development of non-toxic alternatives. Therefore, to enhance the safety of flame retardants, the synthesis of flame retardants from waste feedstocks using phosphorous chemistry with a dual flame-retardant mechanism has been established. A range of waste feedstocks that include cardanol, vanillin, and gallic acids has been converted into a series of flame retardants using a one-step approach by incorporating phosphorous moieties into their structures. The established pathways allow to develop a range of flame-retardant materials by converting waste feedstocks into phosphorous-based materials with inherent flame retardancy. The introduction of abundant aromatic structures from biobased feedstocks enables high charring capabilities in materials in which these biobased flame retardants have been incorporated, increasing their char yields which significantly enhances the flame suppressing properties of the designed materials. Studies show that inclusion of phosphoric moieties into structures allows the displacement of flame enhancing radicals, by releasing non-flammable and non-toxic gases, therefore inhibiting the spread of fire in the gas phase. The resulting biobased flame retardants have been incorporated in wood substrates, foam insulation, and hemp fibers where the results show that only 1-5% of loading of biobased flame retardant suppressed the flame completely. This study represents a novel approach for the development of flame retardants with high performance while utilizing waste feedstocks as a source for their design.

Demchuk, Zoriana [ORNL] (ORCID:0000000326292235)↗