Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum 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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mineral reactivity determines root effects on soil organic carbon

Modern conceptual models of soil organic carbon (SOC) cycling focus heavily on the microbe-mineral interactions that regulate C stabilization. However, the formation of ‘stable’ (i.e. slowly cycling) soil organic matter, which consists mainly of microbial residues associated with mineral surfaces, is inextricably linked to C loss through microbial respiration. Therefore, what is the net impact of microbial metabolism on the total quantity of C held in the soil? To address this question, we constructed artificial root-soil systems to identify controls on C cycling across the plant-microbe-mineral continuum, simultaneously quantifying the formation of mineral-associated C and SOC losses to respiration. Here we show that root exudates and minerals interacted to regulate these processes: while roots stimulated respiratory C losses and depleted mineral-associated C pools in low-activity clays, root exudates triggered formation of stable C in high-activity clays. Moreover, we observed a positive correlation between the formation of mineral-associated C and respiration. This suggests that the growth of slow-cycling C pools comes at the expense of C loss from the system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

OzDES Reverberation Mapping Programme: Mg ii lags and R − L relation

ABSTRACT The correlation between the broad line region radius and continuum luminosity (R–L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1–2, where AGN activity peaks, the R–L relation is constrained by the reverberation mapping (RM) lags of the Mg ii line. We present 25 Mg ii lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg ii lags and extends the R–L relation to higher redshifts and luminosities. The relative iron line strength $\mathcal {R}_{\rm Fe}$ has little impact on the R–L relation. The best-fitting Mg iiR–L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter $\sigma _{\rm rl} = 0.15^{+0.03}_{-0.02}$ . The slope is consistent with previous measurements and shallower than the H β R–L relation. The intrinsic scatter of the new R–L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R–L relation. Our new R–L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials

Ab initio calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. In this work, we generalize previous work on the ab initio Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for three-dimensional and 2D systems when the interlayer separation tends to zero or infinity, respectively. This expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct ab initio density-functional perturbation theory calculations for monolayer BN and MoS 2 , and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A MeerKAT-meets-LOFAR study of Abell 1413: a moderately disturbed non-cool-core cluster hosting a ~500 kpc ‘mini’-halo

Many relaxed cool-core clusters host diffuse radio emission on scales of hundreds of kiloparsecs: mini-haloes. However, the mechanism responsible for generating them, as well as their connection with central active galactic nuclei, is elusive and many questions related to their physical properties and origins remain unanswered. This paper presents new radio observations of the galaxy cluster Abell 1413 performed with MeerKAT (L-band; 872 to 1712 MHz) and LOFAR HBA (120 to 168 MHz) as part of a statistical and homogeneous census of mini-haloes. Abell 1413 is unique among mini-halo clusters as it is a moderately disturbed non-cool-core cluster. Our study reveals an asymmetric mini-halo up to 584 kpc in size at 1283 MHz, twice as large as first reported at similar frequencies. The spectral index is flatter than previously reported, with an integrated value of α = −1.01 ± 0.06, shows significant spatial variation, and a tentative radial steepening. We studied the point-to-point X-ray/radio surface brightness correlation to investigate the thermal/non-thermal connection: our results show a strong connection between these components, with a superlinear slope of $b = 1.63^{+0.10}_{-0.10}$ at 1283 MHz and $b = 1.20^{+0.13}_{-0.11}$ at 145 MHz. We also explore the X-ray surface brightness/radio spectral index correlation, finding a slope of $b = 0.59^{+0.11}_{-0.11}$. Both investigations support the evidence of spectral steepening. Finally, in the context of understanding the particle acceleration mechanism, we present a simple theoretical model which demonstrates that hybrid scenarios – secondary electrons (re-)accelerated by turbulence – reproduce a superlinear correlation slope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Small Reservoirs Offer a New Perspective on Flood Reduction in Large Basins

The flood reduction potential of individual reservoirs within a large river network continuum remains poorly understood due to the complex interplay between reservoir characteristics and network properties. Here we investigate whether a collection of relatively small reservoirs can play a significant role in mediating downstream floods and assess how that role may be influenced by reservoir network properties compared to traditionally known reservoir characteristics. Our unique contribution was the simulation of downstream flood inundation maps alongside peak flows for each of the 81 major reservoirs (6 × 104 to 8 × 109 m3) across 15,000 river reach segments, integrated into a process-based hydrologic model covering a 415,000 km2 region in the Texas Gulf Coast, United States. The three key takeaways from our study are as follows. (a) Smaller reservoirs can substantially reduce downstream flooding, suggesting that inclusion of large reservoirs—the traditional approach in flood risk management studies—may present only a partial picture. (b) Flood reduction by smaller reservoirs is more effective upstream, although this phenomenon may be linked with aridity and overall water availability. (c) While reservoir size matters, it is not the primary factor determining its downstream flood reduction potential; the influence of network properties, such as catchment area, the count of upstream reservoirs, the cumulative maximum storage capacity of upstream reservoirs, and Stream Order (i.e., location), is equally and often more important. The broader impact of our findings goes beyond just floods, providing foundational insights for addressing emerging challenges such as aging dams and river connectivity.

Patel, Krutikkumar [University of Texas at Arlingt↗
DOE OSTI2023

Path integral suppression of badly behaved causal sets

Abstract Causal set theory is a discrete model of spacetime that retains a notion of causal structure. We understand how to construct causal sets that approximate a given spacetime, but most causal sets are not at all manifold-like, and must be dynamically excluded if something like our Universe is to emerge from the theory. Here we show that the most common of these ‘bad’ causal sets, the Kleitman–Rothschild orders, are strongly suppressed in the gravitational path integral, and we provide evidence that a large class of other ‘bad’ causal sets are similarly suppressed. It thus becomes plausible that continuum behavior could emerge naturally from causal set quantum theory.

Astronomy & Astrophysics↗
DOE OSTI2021

Detecting forest response to droughts with global observations of vegetation water content

Droughts in a warming climate have become more common and more extreme, making understanding forest responses to water stress increasingly pressing. Analysis of water stress in trees has long focused on water potential in xylem and leaves, which influences stomatal closure and water flow through the soil-plant-atmosphere continuum. At the same time, changes of vegetation water content (VWC) are linked to a range of tree responses, including fluxes of water and carbon, mortality, flammability, and more. Unlike water potential, which requires demanding in situ measurements, VWC can be retrieved from remote sensing measurements, particularly at microwave frequencies using radar and radiometry. Here, we highlight key frontiers through which VWC has the potential to significantly increase our understanding of forest responses to water stress. To validate remote sensing observations of VWC at landscape scale and to better relate them to data assimilation model parameters, we introduce an ecosystem-scale analogue of the pressure-volume curve, the non-linear relationship between 44 average leaf or branch water potential and water content commonly used in plant hydraulics. The 45 sources of variability in these ecosystem-scale pressure-volume curves and their relationship to 46 forest response to water stress are discussed. We further show to what extent diel, seasonal, and 47 decadal dynamics of VWC reflect variations in different processes relating the tree response to 48 water stress. VWC can also be used for inferring belowground conditions – which are difficult to 49 impossible to observe directly. Lastly, we discuss how a dedicated geostationary spaceborne 50 observational system for VWC, when combined with existing datasets, can capture diel and 51 seasonal water dynamics to advance the science and applications of global forest vulnerability to 52 future droughts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Shock wave formation in the thermosphere by an earthgrazing fireball: Empirical evidence for volatile-enhanced hydrodynamic shielding

Hydrodynamic shielding is a theoretically well-established but observationally elusive and experimentally difficult-to-replicate phenomenon with implications that extend far beyond meteor physics. Rare earthgrazing meteoroids with infrasound signatures that penetrate to the ground can be used to probe hydrodynamic shielding that leads to strong shock formation at high altitude. Here, we report the first coordinated optical and multi-station infrasound observations of a centimeter-scale earthgrazing fireball that generated sustained cylindrical line shock at thermospheric altitudes near 92 km. The event was recorded by numerous optical stations and three infrasound arrays, allowing trajectory reconstruction, ablation behavior, acoustic source localization, and shock characteristics. Optical observations indicate early mechanical erosion and ablation/evaporation at exceptionally low dynamic pressure, consistent with a cometary or a porous, volatile-bearing CM chondritic object. Independent infrasound detections localize shock generation to multiple points along a 164 km trajectory segment near perigee. Weak-shock modeling yields a consistent blast radius of ∼30 m, implying an acoustic-equivalent source size far exceeding the physical dimensions of the ∼45 g nucleus. We demonstrate that classical gas dynamics and ablation-driven hydrodynamic shielding alone cannot account for these observations under ambient thermospheric conditions. We show that volatile release provides the additional flow-field density enhancement required to amplify hydrodynamic shielding, reduce the effective local Knudsen number, and sustain a shock envelope capable of radiating detectable infrasound. Furthermore, these results demonstrate that small, volatile-rich meteoroids can transiently establish continuum-like flow in rarefied environments.

Astrometry↗
DOE OSTI2023

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗
DOE OSTI2022

Accretion disc sizes from continuum reverberation mapping of AGN selected from the ZTF survey

ABSTRACT We present the accretion disc-size estimates for a sample of 19 active galactic nuclei (AGNs) using the optical g-, r-, and i-band light curves obtained from the Zwicky Transient Facility survey. All the AGNs have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multiband light curves are cross-correlated, and the reverberation lag is estimated using the Interpolated Cross-Correlation Function method and the Bayesian method using the javelin code. As expected from the disc-reprocessing arguments, the g − r band lags are shorter than the g − i band lags for this sample. The interband lags for all, but five sources, are larger than the sizes predicted from the standard Shakura Sunyaev (SS) analytical model. We fit the light curves directly using a thin disc model implemented through the javelin code to get the accretion disc sizes. The disc sizes obtained using this model are on an average 3.9 times larger than the prediction based on the SS disc model. We find a weak correlation between the disc sizes and the known physical parameters, namely the luminosity and the SMBH mass. In the near future, a large sample of AGNs covering broader ranges of luminosity and SMBH mass from large photometric surveys would be helpful in a better understanding of the structure and physics of the accretion disc.

Jha, Vivek Kumar (ORCID:0000000232776335)↗
DOE OSTI2023

Quantum disordered ground state in the triangular-lattice magnet NaRuO 2

It has long been hoped that spin liquid states might be observed in materials that realize the triangular lattice Hubbard model. However, weak spin-orbit coupling and other small perturbations often induce conventional spin freezing or magnetic ordering. Sufficiently strong spin-orbit coupling, however, can renormalize the electronic wave function and induce anisotropic exchange interactions that promote magnetic frustration. Here we show that the cooperative interplay of spin-orbit coupling and correlation effects in the triangular lattice magnet NaRuO 2 produces an inherently fluctuating magnetic ground state. Despite the presence of a charge gap, we find that low-temperature spin excitations generate a metal-like term in the specific heat and a continuum of excitations in neutron scattering, reminiscent of spin liquid states previously found in triangular lattice organic magnets. Further cooling produces a crossover into a different, highly disordered, spin state whose dynamic spin autocorrelation function reflects persistent fluctuations. Furthermore, these findings establish NaRuO 2 as a cousin to organic, Heisenberg spin liquid compounds with a low-temperature crossover in quantum disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effects of emitting surfaces and trapped ions on the sheath physics and current flow in multidimensional plasma systems

Recent one-dimensional simulations of planar sheaths with strong electron emission have shown that trapping of charge-exchange ions causes transitions from space-charge limited (SCL) to inverse sheaths. However, multidimensional emitting sheath phenomena with collisions remained unexplored, due in part to high computational cost. As such, we developed a novel continuum kinetic code to study the sheath physics, current flow and potential distributions in two-dimensional unmagnetized configurations with emitting surfaces. For small negatively biased thermionic cathodes in a plasma, the cathode sheath can exist in an equilibrium SCL state. The SCL sheath carries an immense density of trapped ions, neutralized by thermoelectrons, within the potential well of the virtual cathode. For further increases of emitted flux, the trapped ion cloud expands in space. The trapped ion space charge causes an increase of thermionic current far beyond the saturation limit predicted by conventional collisionless SCL sheath models without ion trapping. For sufficiently strong emission, the trapped ion cloud consumes the entire 2D plasma domain, forming a mode with globally confined ions and an inverse sheath at the cathode. In situations where the emitted flux is fixed and the bias is swept (e.g. emissive probe), the trapped ions cause a large thermionic current to escape for all biases below the plasma potential. Strong suppression of the thermionic emission, required for the probe to float, only occurs when the probe is above the plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Disk-integrated Thermal Properties of Ceres Measured at Millimeter Wavelengths

We observed Ceres at three epochs in 2015 November and 2017 September and October with Atacama Large Millimeter/submillimeter Array (ALMA) 12 m array and in 2017 October with the ALMA Compact Array (ACA), all at ∼265 GHz continuum (wavelengths of ∼1.1 mm) to map the temperatures of Ceres over a full rotation at each epoch. We also used 2017 October ACA observations to search for HCN. The disk-averaged brightness temperature of Ceres is measured to be between 170 and 180 K during our 2017 observations. The rotational light curve of Ceres shows a double-peaked shape with an amplitude of about 4%. Our HCN search returns a negative result with an upper limit production rate of ∼2 × 10{sup 24} molecules s{sup −1}, assuming globally uniform production and a Haser model. A thermophysical model suggests that Ceres’s top layer has higher dielectric absorption than lunar-like materials at a wavelength of 1 mm. However, previous observations showed that the dielectric absorption of Ceres decreases toward longer wavelengths. Such distinct dielectric properties might be related to the hydrated phyllosilicate composition of Ceres and possibly abundant micrometer-sized grains on its surface. The thermal inertia of Ceres is constrained by our modeling as likely being between 40 and 160 thermal inertia units, much higher than previous measurements at infrared wavelengths. Modeling also suggests that Ceres’s light curve is likely dominated by spatial variations in its physical or compositional properties that cause changes in Ceres’s observed thermal properties and dielectric absorption as it rotates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The impact of demographic, clinical, genetic, and imaging variables on tau PET status

Purpose: A substantial proportion of amyloid-β (Aβ)+ patients with clinically diagnosed Alzheimer’s disease (AD) dementia and mild cognitive impairment (MCI) are tau PET–negative, while some clinically diagnosed non-AD neurodegenerative disorder (non-AD) patients or cognitively unimpaired (CU) subjects are tau PET–positive. We investigated which demographic, clinical, genetic, and imaging variables contributed to tau PET status. Methods: We included 2338 participants (430 Aβ+ AD dementia, 381 Aβ+ MCI, 370 non-AD, and 1157 CU) who underwent [ 18 F]flortaucipir ( n = 1944) or [ 18 F]RO948 ( n = 719) PET. Tau PET positivity was determined in the entorhinal cortex, temporal meta-ROI, and Braak V-VI regions using previously established cutoffs. We performed bivariate binary logistic regression models with tau PET status (positive/negative) as dependent variable and age, sex, APOE ε4, Aβ status (only in CU and non-AD analyses), MMSE, global white matter hyperintensities (WMH), and AD-signature cortical thickness as predictors. Additionally, we performed multivariable binary logistic regression models to account for all other predictors in the same model. Results: Tau PET positivity in the temporal meta-ROI was 88.6% for AD dementia, 46.5% for MCI, 9.5% for non-AD, and 6.1% for CU. Among Aβ+ participants with AD dementia and MCI, lower age, MMSE score, and AD-signature cortical thickness showed the strongest associations with tau PET positivity. In non-AD and CU participants, presence of Aβ was the strongest predictor of a positive tau PET scan. Conclusion: We identified several demographic, clinical, and neurobiological factors that are important to explain the variance in tau PET retention observed across the AD pathological continuum, non-AD neurodegenerative disorders, and cognitively unimpaired persons.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Synthesizing land use and demographic change in Southeast Asia’s smaller urbanized areas from 2000–2015

The majority of the human population now reside in urban areas today. The United Nations estimates that nearly half of all urban dwellers currently live in cities smaller than 500 000 persons and the majority of future urban growth will take place in Asia and Africa, likely in these smaller urban areas, not mega cities. Thus, understanding the factors that influence urban demographic trajectories in small urban areas is critical to address sustainable and equitable policy initiatives related to food security, changing climate hazard exposure, and economic opportunities. Here we focus on Southeast Asia—a region historically characterized by lower urban population proportions, yet with a rapidly shifting dynamic demographic—to examine correlates of demographic change among smaller cities. We combine two open-source satellite-informed datasets: GHS urban center database (2015) and age-sex gridded data from WorldPop to calculate socio-demographic characteristics to model drivers of change in annualized urban population growth from 2000–2015 for 505 urbanized places. We find a general pattern of decreasing dependency ratios as city-size increases for most urban areas in Southeast Asia. Higher rates of growth and more variation is observed for smaller cities—those with fewer than 300 000 persons, the lowest population limit for UN data on urbanization. When examining covariates of urban population growth, we find significant statistical associations of population change in smaller urbanized areas with climatic, economic, and land cover/land use variables, but with country-specific variations. Characterizing a continuum of urban population development in the context of changing environmental, economic and climate conditions has been an important sustainable development and equity issue for decades, but newer analysis of city-level drivers allows for systematic inquiry thus moving beyond total population counts for policy-relevant insight.

Southeast Asia synthesis↗
DOE OSTI2023

Probing the rapid formation of black holes and their Galaxy hosts in QSOs

Using the modelling code x-cigale, we reproduced the spectral energy distributions (SEDs) of 1359 SDSS quasi-stellar objects (QSOs) within the redshift range 0 < z < 4, for which we have NIR/MIR fluxes with the highest quality and spectral data characterizing their supermassive black hole (SMBHs). Consistent with a rapid formation of the host galaxies, the star formation histories (SFHs) have small e-folding, at most 750 Myr using an SFH function for spiral or 1000 Myr using one for elliptical. Above z ∼ 1.6, the two solutions are degenerate, the SEDs being dominated by the active galactic nucleus (AGN) continuum and high star formation rates (SFRs), typical of starburst galaxies, while at lower redshifts the starburst nature of the host, independent from its morphology, is better reproduced by a spiral SFH. In general, the SFR increases with the redshift, the mass of the bulge, the AGN luminosity, and Eddington ratio, suggesting there is no evidence of AGN quenching of star formation. Comparing the specific BH accretion rate (BHAR) with specific SFR, all the QSOs at any redshift trace a linear sequence below the Eddington luminosity, in parallel and above the one-to-one relation, implying that QSOs are in a special phase of evolution during which the growth in mass of their SMBH is more rapid than the growth in mass of their galaxy hosts. This particular phase is consistent with a scenario where the galaxy hosts of QSOs in the past grew in mass more rapidly than their SMBHs, suggesting that a high star formation efficiency during their formation was responsible in limiting their masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Prediction of dislocation - grain boundary interactions in FCC aluminum bicrystals using a modified continuum criterion and machine learning methods

Mechanical properties of metals such as strength and toughness are strongly correlated to complex interactions between various defects in the crystalline structure. While elementary interactions between these defects have been investigated using recent micro- and nano-characterization techniques, understanding of the detailed interaction mechanisms has hardly been obtained. To understand defect-driven plasticity at various time and length scales, it is necessary to formulate a general guideline to predict both the interaction type (transmission or reflection) and the dislocation's subsequent slip system after the interaction. Many criteria based on the geometric alignment of the defects have been developed to predict this phenomenon, but these have yet to be found to be accurate when applied to general data sets of grain boundaries (GBs). With this motivation, we conduct a systematic study using molecular dynamics (MD) models of bicrystals to analyze defect interaction process between a prismatic dislocation loop and eleven different grain boundaries of the following character: three tilt, three twist, and five mixed. Based on the MD observations, two new prediction methods are developed: the first is a new data-driven parametric score function based on the classical geometric criteria, and the second is by applying Gaussian process machine learning methods to find the probability distribution of a hidden function. In conclusion, the proposed data-driven prediction methods could pave a new way to predict the unit interaction of dislocations with various GBs, which could show much higher accuracy compared to pre-existing geometric criteria.

36 MATERIALS SCIENCE↗