Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum scale modeling”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗
NASA NTRS1995

High Performance Parallel Computational Nanotechnology

At a recent press conference, NASA Administrator Dan Goldin encouraged NASA Ames Research Center to take a lead role in promoting research and development of advanced, high-performance computer technology, including nanotechnology. Manufacturers of leading-edge microprocessors currently perform large-scale simulations in the design and verification of semiconductor devices and microprocessors. Recently, the need for this intensive simulation and modeling analysis has greatly increased, due in part to the ever-increasing complexity of these devices, as well as the lessons of experiences such as the Pentium fiasco. Simulation, modeling, testing, and validation will be even more important for designing molecular computers because of the complex specification of millions of atoms, thousands of assembly steps, as well as the simulation and modeling needed to ensure reliable, robust and efficient fabrication of the molecular devices. The software for this capacity does not exist today, but it can be extrapolated from the software currently used in molecular modeling for other applications: semi-empirical methods, ab initio methods, self-consistent field methods, Hartree-Fock methods, molecular mechanics; and simulation methods for diamondoid structures. In as much as it seems clear that the application of such methods in nanotechnology will require powerful, highly powerful systems, this talk will discuss techniques and issues for performing these types of computations on parallel systems. We will describe system design issues (memory, I/O, mass storage, operating system requirements, special user interface issues, interconnects, bandwidths, and programming languages) involved in parallel methods for scalable classical, semiclassical, quantum, molecular mechanics, and continuum models; molecular nanotechnology computer-aided designs (NanoCAD) techniques; visualization using virtual reality techniques of structural models and assembly sequences; software required to control mini robotic manipulators for positional control; scalable numerical algorithms for reliability, verifications and testability. There appears no fundamental obstacle to simulating molecular compilers and molecular computers on high performance parallel computers, just as the Boeing 777 was simulated on a computer before manufacturing it.

Saini, Subhash↗
NASA NTRS2001

The Lyman-Alpha Forest in the Lensed Quasar Q0957+561 and the Characteristic Dimension of the Absorbing Clouds

Far-ultraviolet spectra of the gravitational lens components Q0957+561 A and 9 were obtained with the Hubble Space Telescope Faint Object Spectrograph (FOS) at five equally spaced epochs, one every two weeks. We confirm the flux variability of the quasar's Lyman-alpha and 0 VI emission lines reported by Dolan et al. (1995) in IUE spectra. The fluxes in these lines vary on a time scale of weeks in the observer's rest frame, independently of each other and of the surrounding continuum. The individual spectra of each image were co-added to investigate the properties of the Lyman-alpha forest along the two lines of sight to the quasar. Absorption lines having equivalent width W > 0.3 A in the observer's frame not previously identified by Michalitsianos et al. (1997) as interstellar lines, metal lines, or higher order Lyman lines were taken to be Ly-alpha forest lines. The existence of each line in this consistently selected set was then verified by its presence in two archival FOS spectra with -1.5 times higher signal to noise than our co-added spectra. Ly-alpha forest lines with W > 0.3 A appear at 41 distinct wavelengths in the spectra of the two images. one absorption line in the spectrum of image A has no counterpart in the spectrum of image B and one line in image B has no counterpart in image A. Based on the separation of the lines of sight over the redshift range searched for Ly-forest lines, the density of the absorbing clouds in the direction of Q0957+561 must change significantly over a distance R = 160 (+120, -70)/ h(sub 50) kpc in the simplified model where the absorbers are treated as spherical clouds and the characteristic dimension, R, is the radius. (We adopt H(sub 0) = 50 h(sub 50) km/s/ kpc, q(sub 0) = 1/2, and lambda = 0 throughout the paper.) The 95% confidence interval on R extends from (50 - 950)/h(sub 50) kpc We show in the Appendix that the fraction of Ly-alpha forest lines that appear in only one spectrum can be expressed as a rapidly converging power series in 1/r, where r the ratio of the radius of the cloud to the separation of the two lines of sight at the redshift of the cloud. This power series can be rewritten to give r in terms of the fraction of Ly-forest wavelengths that appear in the spectrum of only one image. A simple linear approximation to the solution that everywhere agrees with the power series solution to better than 0.8% for r > 2 is derived in the Appendix.

Dolan, Joseph↗
NASA NTRS2000

Flow in the Proximity of the Pin-Tool in Friction Stir Welding and Its Relation to Weld Homogeneity

In the Friction Stir Welding (FSW) process a rotating pin inserted into a seam literally stirs the metal from each side of the seam together. It is proposed that the flow in the vicinity of the pin-tool comprises a primary rapid shear over a cylindrical envelope covering the pin-tool and a relatively slow secondary flow taking the form of a ring vortex about the tool circumference. This model is consistent with a plastic characterization of metal flow, where discontinuities in shear flow are allowed but not viscous effects. It is consistent with experiments employing several different kinds of tracer: atomic markers, shot, and wire. If a rotating disc with angular velocity w is superposed on a translating continuum with linear velocity omega, the trajectories of tracer points become circular arcs centered upon a point displaced laterally a distance v/omega from the center of rotation of the disc in the direction of the advancing side of the disc. In the present model a stream of metal approaching the tool (taken as the coordinate system of observation) is sheared at the slip surface, rapidly rotated around the tool, sheared again on the opposite side of the tool, and deposited in the wake of the tool. Local shearing rates are high, comparable to metal cutting in this model. The flow patterns in the vicinity of the pin-tool determine the level of homogenization and dispersal of contaminants that occurs in the FSW process. The approaching metal streams enfold one another as they are rotated around the tool. Neglecting mixing they return to the same lateral position in the wake of the tool preserving lateral tracer positions as if the metal had flowed past the tool like an extrusion instead of being rotated around it. (The seam is, however, obliterated.) The metal stream of thickness approximately that of the tool diameter D is wiped past the tool at elevated temperatures drawn out to a thickness of v/2(omega) in the wiping zone. Mixing distances in the wiping zone are multiplied in the unfolded metal. Inhomogeneities on a smaller scale than the mixing length are obliterated, but structure on a larger scale may be transmitted to the wake of a FSW weld.

Nunes, Arthur C., Jr.↗
DOE OSTI2022

A continuum model for heat and mass transfer in moving-bed reactors for thermochemical energy storage

In this work, a continuum heat and mass transfer model coupling transport phenomena and high-temperature thermochemical reactions is developed for stationary packed-bed and counter-flow moving-bed reactors. After presenting the general modeling framework, we focus on the 2D axisymmetric version of the model for which validation is conducted with experimental results for a packed-bed reactor in the literature for manganese-iron oxide reduction/oxidation and an in-house counter-flow moving-bed reactor for magnesium-manganese oxide reduction up to 1450 °C. Transient simulation results including the local distributions of gas/solid temperatures, oxygen concentration and the extent of reaction, as well as the various energy flow components and energy conversion efficiencies are reported. The results based on the 2D axisymmetric model are also compared with those obtained from a previous 1D model. The comparison shows that capturing the radial variation is critical in reactor modeling and the 2D results demonstrate improved agreement with experiments. Specifically, large temperature variations along the radial direction are observed especially in the reaction zone; this non-uniform radial temperature distribution has a significant effect on the chemical reaction extent due to its strong dependence on temperature; and the overall oxygen concentration at the reactor exit and the predicted system efficiency are slightly lower in the 2D model compared to the 1D model. Finally, the present heat and mass transfer model can provide valuable insights into reactor design, scale-up, and operating conditions selection to maximize system energy storage efficiency.

25 ENERGY STORAGE↗
NASA NTRS1999

UV Spectral Variability and the Lyman-Alpha Forest in the Lensed Quasar Q0957+561

Far-ultraviolet spectra of the gravitational lens components Q0957+561 A and B were obtained with the Hubble Space Telescope Faint Object Spectrograph (FOS) at five equally spaced epochs, one every two weeks. We confirm the flux variability of the quasar's Lyman-alpha and O VI lambda 1037 emission lines reported in IUE (International Ultraviolet Explorer) spectra. The fluxes in these lines vary on a time scale of weeks in the local rest frame (LRF), independently of each other and of the surrounding continuum. The individual spectra of each image were co-added to investigate the properties of the Lyman-alpha forest along the two lines of sight to the quasar. Absorption lines having equivalent width W (sub lambda) greater than or equal to 0.3 Angstroms in the LRF not previously identified as interstellar lines, metal lines, or higher order Lyman lines were taken to be LY-alpha forest lines. The existence of each line in this consistently selected set was then verified by its presence in two archival FOS spectra with approx. 1.5 times higher signal to noise than our co-added spectra. Ly-alpha forest lines with W (sub lambda) greater than or equal to 0.3 Angstroms appear at 42 distinct wavelengths in the spectra of the two images. Two absorption lines in the spectrum of image A have no counterpart at that wavelength in the spectrum of image B, and two lines in image B have no counterpart in image A. Based on the separation of the lines of sight at the redshift of the absorption lines appearing in only one spectrum, the density of the absorbing clouds in the direction of Q0957+561 must change significantly over a distance of 79 (+34, -26) h (sub 50) (sup -1) kpc in the simplified model where the absorbers are treated as spherical clouds and the characteristic dimension is the radius. (We adopt H (sub 0)= 50 h (sub 50) km s (sup -1) kpc (sup -1), q (sub 0) = 1/2, and LAMBDA = 0 throughout the paper.) The two limits define the 68% confidence interval on the characteristic dimension, equivalent to the 1 sigma confidence interval for a Gaussian distribution. The 95% confidence interval extends from (32 - 250) h (sub 50) (sup -1) kpc. We show in the Appendix that the fraction of Ly-alpha forest lines that appear in only one spectrum can be expressed as a rapidly converging power series in 1/r, where r the ratio of the radius of the cloud to the separation of the two LOS at the redshift of the cloud. This power series can be rewritten to give r in terms of the fraction of Ly-alpha forest wavelengths that appear in the spectrum of only one image. A simple linear approximation to the solution which everywhere agrees with the power series solution to better than 0.8% for r greater than or equal to 2 is derived in the Appendix.

Dolan, J. F.↗
DOE OSTI2023

Black hole and galaxy co-evolution in radio-loud active galactic nuclei at z ~ 0.3–4

The relation between the mass of the supermassive black hole (SMBH) in the center of galaxies and their bulge mass or central velocity dispersion is well known. This suggests a coevolution between the SMBHs and their galaxy hosts. Our aim is to study this relation, specifically, for radio loud galaxies, and as a function of redshift z. We selected a sample of 42 radio galaxies and active galactic nuclei (AGN) with broad emission lines and spectroscopic redshifts between z = 0.3 - 4 by cross-matching the low radio frequency sources from Very Large Array (VLA) FIRST with spectroscopically confirmed galaxies from wide-field surveys, including Sloan Digital Sky Survey (SDSS) DR14 ugriz and Dark Energy Survey (DES) DR2 grzY in the optical, Wield Infrared Survey Explorer (WISE), and the Galaxy And Mass Assembly (GAMA) spectroscopic survey. We characterized the stellar mass (M $\star$ ), star formation, and black hole properties (mass of the central SMBH, Eddington ratio η, and jet power, Q jet ). The relation between SMBH mass, M $\star$ , η, and z is placed into context by comparing them with scaling relations (M BH –M $\star$ , M BH /M $\star$ –z, M BH –Q jet , and Q jet –η) from the literature. On the basis of a multiwavelength spectral energy distribution modeling, our radio sources are broadly consistent with being on the star-forming main sequence. They have sub-Eddington accretion rates, η ≃ 1% on average, as typically found in type I AGN, while higher accretion rates favor more powerful jets to be launched by the central engine. We find overmassive SMBHs in (17 ± 5)% of our radio sources, similarly to previous studies on nearby early-type galaxies. Altogether, an evolutionary scenario in which radio-mode AGN feedback regulates the accretion onto the SMBHs and the stellar mass assembly of the radio sources is discussed, which may explain the observed phenomenology. This pilot study represents a benchmark for future studies using wide-field surveys such as those with Euclid and the Vera Rubin Observatory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling the Spectral Diversity of Quasars in the Sixteenth Data Release from the Sloan Digital Sky Survey

Abstract We present a new approach to capturing the broad diversity of emission-line and continuum properties in quasar spectra. We identify populations of spectrally similar quasars through pixel-level clustering on 12,968 high signal-to-noise ratio (S/N) spectra from the Sloan Digital Sky Survey (SDSS) in the redshift range of 1.57 < z < 2.4. Our clustering analysis finds 396 quasar spectra that are not assigned to any population, 15 misclassified spectra, and 6 quasars with incorrect redshifts. We compress the quasar populations into a library of 684 high-S/N composite spectra, anchored in redshift space by the Mg ii emission line. Principal component analysis on the library results in an eigenspectrum basis spanning 1067–4007 Å. We model independent samples of SDSS quasar spectra with the eigenbasis, allowing for a free redshift parameter. Our models achieve a median reduced χ 2 on non–broad absorption line quasar spectra that is reduced by 8.5% relative to models using the eigenspectra from the SDSS spectroscopic pipeline. A significant contribution to the relative improvement is from the ability to reconstruct the range of emission-line variation. The redshift estimates from our model are consistent with the Mg ii emission-line redshift with an average offset that displays 51.4% less redshift-dependent variation relative to the SDSS eigenspectra. Our method for developing quasar spectra models can improve automated classification and predict the intrinsic spectrum in regions affected by intervening absorbers such as Ly α , C iv , and Mg ii , thus benefiting studies of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗
DOE OSTI2023

Diffuse sources, clustering, and the excess anisotropy of the radio synchrotron background

ABSTRACT We present the largest low frequency (120 MHz) arcminute resolution image of the radio synchrotron background (RSB) to date, and its corresponding angular power spectrum of anisotropies (APS) with angular scales ranging from 3° to 0.3 arcmin. We show that the RSB around the north celestial pole has a significant excess anisotropy power at all scales over a model of unclustered point sources based on source counts of known source classes. This anisotropy excess, which does not seem attributable to the diffuse Galactic emission, could be linked to the surface brightness excess of the RSB. To better understand the information contained within the measured APS, we model the RSB varying the brightness distribution, size, and angular clustering of potential sources. We show that the observed APS could be produced by a population of faint clustered point sources only if the clustering is extreme and the size of the Gaussian clusters is ≲1 arcmin. We also show that the observed APS could be produced by a population of faint diffuse sources with sizes ≲1 arcmin, and this is supported by features present in our image. Both of these cases would also cause an associated surface brightness excess. These classes of sources are in a parameter space not well probed by even the deepest radio surveys to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

BASS XXXII: Studying the Nuclear Millimeter-wave Continuum Emission of AGNs with ALMA at Scales ≲100–200 pc

To understand the origin of nuclear (≲100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed subarcsecond resolution Band-6 (211–275 GHz) Atacama Large Millimeter/submillimeter Array data of 98 nearby AGNs (z < 0.05) from the 70 month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling (~1–200 pc), and spans broad ranges of 14–150 keV luminosity {$40\lt \mathrm{log}[{L}_{14-150}/(\mathrm{erg}\,{{\rm{s}}}^{-1})]\lt 45$}, black hole mass $[5\lt \mathrm{log}({M}_{\mathrm{BH}}/{M}_{\odot })\lt 10$], and Eddington ratio ($-4\lt \mathrm{log}{\lambda }_{\mathrm{Edd}}\lt 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14–150 keV luminosities. Its scatter is ≈0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to N H ~ 10 26 cm –2 . While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission: AGN-related dust emission, outflow-driven shocks, and a small-scale (<200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2023

Young but fading radio sources: searching for remnants among compact steep-spectrum radio sources

The incidence of young but fading radio sources provides important information on the life cycle of radio emission in radio-loud active galactic nuclei. Despite its importance for constraining the models of radio source evolution, there are no systematic studies of remnants in complete samples of young radio sources. We report results of the study of 18 compact steep-spectrum (CSS) radio sources, selected from the statistically complete B3-VLA (Very Large Array) CSS sample, characterized by a steep optically thin spectrum (α ≥ 1.0) and no core detection in earlier studies. Our deep multifrequency VLA, pc-scale Very Long Baseline Array, and enhanced Multi Element Remotely Linked Interferometer Network (e-MERLIN) observations allowed us to locate the core component in 10 objects. In three CSS sources, there is no clear evidence of present-time active regions, suggesting that they are likely in a remnant phase. Among sources with core detection, we find three objects that have no clear active regions (hotspots) at the edges of the radio structure, suggesting that the radio emission may have just restarted. Our results support a power-law distribution of the source ages, although the poor statistics prevents us from setting solid constraints on the percentage of remnants and restarted sources in subpopulations of radio sources.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Ultraviolet illuminated molecular cloud boundaries: Extended (C II) 158 micrometer emission toward L1630

The construction of a large scale map of the 158 micrometer C+ line toward the L1630/Orion B molecular cloud, covering an approximately 35' by 45' area which includes the NGC 2024 H II region, zeta-Ori, the reflection nebula NGC 2023, and the Horsehead nebula, is reported. Emission in the C II line is very widespread. The line was detected at levels in excess of a few 0.0001 erg/sq cm/s/sr over almost the entire mapped region. Extended emission associated with the NGC 2024 H II region and its envelope accounts for more than half of the C II flux. Over this approximately 1.5 by 2.5 pc region, the amount of gas phase carbon in the form of C+ is comparable to the amount of carbon in CO. This result, together with the C II distribution implies that C II emission arises on the surface of clumps throughout the cloud rather than in a single layer at the H II region boundary. Away from the H II region, most of the C II emission comes from the western edge of the L1630 cloud and probably results from excitation by external OB stars. The overall extent of the C II emission is comparable to that of millimeter molecular lines but the distributions are different in detail. The difference in C II and molecular line distributions, in particular, the larger extent of the C II emission west of NGC 2024 implies large variations in the ratio of the C II and J = 1 towards 0 intensities. Models of photon dominated regions can explain the relation between C II and CO intensities only if the cloud edges and cloud interior are considered separately. A method for using C II and radio continuum emission to characterize the relationship between OB stars and photon dominated regions is proposed.

Jaffe, D. T.↗
NASA NTRS1986

High resolution observations of the L1551 bipolar outflow

The nearby dark cloud Lynds 1551 contains one of the closest examples of a well-collimated bipolar molecular outflow. This source has the largest angular size of any known outflow and was the first bipolar outflow to be detected. The outflow originates from a low-luminosity young stellar object, IRS-5. Optical and radio continuum observations show the presence of a highly collimated, ionized stellar wind orginating from close to IRS-5 and aligned with the molecular outflow. However, we have little information on the actual mechanism that generates the stellar wind and collimates it into opposed jets. The Very Large Array (VLA) observations indicate that the winds originate within 10(15) cm of IRS-5, unfortunately at a size scale difficult to resolve. For these reasons, observations of the structure and dynamics of the hypersonic molecular gas may provide valuable information on the origin and evolution of these outflows. In addition, the study of the impact of the outflowing gas on the surrounding molecular material is essential to understand the consequence these outflows have on the evolution and star formation history of the entire cloud. Moriarty-Schieven et al. (1986) obtained a oversampled map of the CO emission of a portion of both the blueshifted and redshifted outflows in LI551 using Five College Radio Astronomy Observatory 14 m telescope. The oversampled maps have been reconstructed to an effective angular resolution of 20 arcsec using a maximum entropy algorithm. A continuation of the study of Moriarty-Schieven et al. is presented. The entire L1551 outflow has now been mapped at 12 arcsec sampling requiring roughly 4000 spectra. This data has been constructed to 20 arcsec resolution to provide the first high resolution picture of the entire L1551 outflow. This new data has shown that the blueshifted lobe is more extended than previously thought and has expanded downstream sufficiently to break out of the dense molecular cloud, but the redshifted outflow is still confined within the molecular cloud. Details of the structure and kinematics of the high velocity gas are used to test the various models of the origin and evolution of outflows.

Snell, R.↗
NASA NTRS1994

Ultraviolet-illuminated molecular cloud boundaries: Extended (C II) 158 micron emission toward L1630

We have made a large-scale map of the 158 micrometers C(+) line toward the L1630/Orion B molecular cloud. The map covers a approximately 35 min x 45 min area which includes the NGC 2024 H II region, zeta Ori, the reflection nebula NGC 2023, and the Horsehead nebula. Emission in the (C II) line is very widespread. The line was detected at levels in excess of a few 10(exp -4) ergs/sq cm/s/sr over almost the entire mapped region. Extended emission associated with the NGC 2024 H II region and its envelope accounts for more than half of the (C II) flux. Over this approximately 1.5 x 2.5 pc region, the amount of gas-phase carbon in the form of C(+) is comparable to the amount of carbon in CO. This result, together with the (C II) distribution, implies that (C II) emission arises on the surfaces of clumps throughout the cloud rather than in a single layer at the H II region boundary. Away from the H II region, most of the (C II) emission comes from the western edge of the L1630 cloud and probably results from excitation by external OB stars. The overall extent of the (C II) emission is comparable to that of millimeter molecular lines, but the distributions are different in detail. The difference in (C II) and molecular line distributions, in particular, the larger extent of the (C II) emission west of NGC 2024 implies large variations in the radio of the (C II) and CO J = 1 goes to 0 intensities. Models of photon-dominated regions can explain the relation between (C II) and CO intensities only if one considers the cloud edges and cloud interior separately. We propose a method for using (C II) and radio continuum emission to characterize the relationship between OB stars and photon-dominated regions.

Jaffe, D. T.↗