Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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,117 records · Page 62

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗

Modeling Nanoscale Ohmics in Carbon Supports of Fuel Cell Cathodes

Here, reducing platinum (Pt) loading in polymer electrolyte fuel cells (PEFCs) while meeting performace requirements is critical to their widespread deployment. However, significant polarization losses manifest at higher current densities in cathodes with lower Pt content. The morphology of the carbon supports in PEFC cathodes affects the location of Pt deposition into the micropits or onto the surface of the carbon support, translating into different kinetic and transport resistances. In this work, we derive an agglomerate scale model that differentiates the sink terms for Pt on the surface and in the pits of carbon supports. We develop an approach to assess the impact of nanoscale ohmic resistance to Pt in the micropits arising from weakly ionic solution in the carbon support on PEFC performance. Effectiveness factors relating the actual reaction rate to the maximum reaction rate (had all the Pt been exposed) are derived and embedded into a one-dimensional catalyst layer model. Parameters in the catalyst layer model are tuned based on experimental local oxygen transport resistances. Subsequently, we estimate bounds for the micropore resistances based on geometric and physical arguments. Lastly, polarization curves are simulated to assess the effect of the micropore resistance in fully-humidified and oxygen-rich environments.

25 ENERGY STORAGE↗

Observed Temperature and Moisture Advection Characteristics in the Southern Great Plains

Until recently, advection profiles in the planetary boundary layer have been difficult to quantify from observations. As a result, not much is known about its basic characteristics, like magnitude and diurnal cycles, or how these differ as a function of large-scale environments. Here, to provide insight into advective characteristics, a Green’s Theorem–based method for calculating profiles of advection in the lowest 3 km from an array of ground-based thermodynamic and kinematic profiling instruments has been applied to 2 years of observations from the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observing site. Since advection largely depends on synoptic-scale forcing, a self-organizing map (SOM) was used to classify days based on mean sea level pressure analyses; from that, the diurnal evolution of vertical profiles of advection for different synoptic environments was quantified. The overall mean magnitude of the potential temperature advection is approximately ±3 K h −1 , while the moisture advection is approximately ±1.5 g kg −1 h −1 , with substantial variability in the sign and magnitude of the advective tendencies at different heights and throughout the diurnal cycle. Advection magnitude is strongly connected to the strength of synoptic-scale forcing, which varies depending on the time of year. Variability in advection is larger for potential temperature than it is for moisture, and strongly forced environments, like midlatitude cyclones, feature greater variability in the advection magnitude than weakly forced or quiescent environments.

advection↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Multi-device Studies of Pedestal Physics and Confinement in the I-mode Regime

This paper describes joint ITPA studies of the I-mode regime, which features an edge thermal barrier but with L-mode-like particle and impurity transport, and without ELMs. The regime has been demonstrated on Alcator C-Mod, ASDEX Upgrade and DIII-D tokamaks, over a wide range of device parameters and pedestal conditions, including low collisionality. Access to I-mode is independent of heating method (neutral beam injection, ICRH and/or ECRH). Normalized energy confinement H98y,2 ≥ 1 has been achieved for a range of 3≤q95≤4.9 and scales favourably with power. Changes in turbulence in the pedestal region accompany the transition from L-mode to I-mode. The L-I threshold increases with plasma density and current, and with device size, but weakly with BT. The upper limit of power for I-modes increases with toroidal field, and is largest on Alcator C-Mod at B > 5 T. Issues for extrapolation to ITER and other future fusion devices are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of a weak shock on interstellar gas toward the Rho Ophiuchi cloud

The effect of a low-velocity shock on depletion in interstellar clouds is studied. High-resolution Copernicus observations of interstellar absorption lines toward four stars in the Rho Ophiuchi cloud complex were used to measure differential depletion in interstellar clouds separated by velocities of 10-15 km/s. Optical observations of CH and CH(+) were used as indicators of shock strength and direction of propagation. Observations indicate the presence of a shock with velocity about 10 km/s expanding away from the sun into the Rho Oph cloud. The existence of this shock is supported by the other observations quoted in the literature. Two distinct regions in the lines of sight have been found. A low-density, less depleted, predominantly atomic region is associated with preshock gas, while postshock gas accounts for a predominantly molecular, more highly depleted region. Apparently a weak shock, has the effect of enhancing grain formation or grain growth as a result of increased density in the postshock gas. It is possible for grain growth by accretion to occur in the shocked gas on time scales short compared with the cloud-crossing time of the shock.

Meyers, K. A.↗

Radial diffusion in magnetodiscs

The orbits of charged particles in magnetodiscs are considered. The bounce motion is assumed adiabatic except for transits of a small equatorial region of weak magnetic field strength and high field curvature. Previous theory and modeling have shown that particles scatter randomly in pitch angle with each passage through the equator. A peaked distribution thus diffuses in pitch angle on the time scale of many bounces. It is argued in this paper that spatial diffusion is a further consequence when the magnetodisc has a longitudinal asymmetry. A general expression for DLL, the diffusion of equatorial crossing radii, is derived. DLL is evaluated explicitly for ions in Jupiter's 20-35 radii magnetodisc, assumed to be represented by Connerney et al.'s (1982) Voyager model plus a small image dipole asymmetry. Rates are energy, species, and space dependent but can average as much as a few tenths of a planetary radius per bounce period.

Birmingham, T. J.↗

A comparative study of rapidly and slowly rotating dynamical regimes in a terrestrial general circulation model

As a preliminary step in the development of a general circulation model for general planetary use, a simplified vesion of thef GISS Model I GCM has been run at various rotation periods to investigate differences between the dynamical regimes of rapidly and slowly rotating planets. To isolate the dynamical processes, the hydrologic cycle is suppressed and the atmosphere is forced with perpetual annual mean solar heating. All other parameters except the rotation period remain fixed at their terrestrial values. Experiments were conducted for rotation periods of 2/3, 1, 2, 4, 8, 16, 64 and 256 days. The results are in qualitative agreement with similar experiments carried out previously with other GCMs and with certain aspects of one Venus GCM simulation. As rotation rate decreases, the energetics shifts from baroclinc to quasi-barotropic when the Rossby radius of deformation reaches planetary scale. The Hadley cell expands poleward and replaces eddies as the primary mode of large-scale heat transport. Associated with this is a poleward shift of the baroclinic zone and jet stream and a reduction of the equator-pole temperature contrast. Midlatitude jet strength peaks at 8 days period, as does the weak positive equatorial zonal wind which occurs at upper levels at all rotation periods. Eddy momentum transport switches from poleward to equatorward at the same period. Tropospheric mean static stability generally increases in the tropics and decreases in midlatitudes as rotation rate decreases, but the global mean static stability is independent of rotation rate. The peak in the eddy kinetic energy spectrum shifts toward lower wavenumbers, reaching wavenumber 1 at a period of 8 days. Implications of these results for the dynamics of Venus and Titan are discussed. Specifically, it is suggested that the extent of low-level convection determines whether the Gierasch mechanism contributes significantly to equatorial superrotation on these planets.

Del Genio, Anthony D.↗

Technologies and Designs for Electronic Nanocomputers

Diverse space-related applications have been proposed for microscopic and sub-microscopic structures, mechanisms, and 'organisms'. To govern their functions, many of these tiny systems will require even smaller, nanometer-scale programmable computers, i.e. 'nanocomputers' on-board. This paper provides an overview of the results of a nearly two-year study of the technologies and designs that presently are in development for electronic nanocomputers. Strengths and weaknesses of the various technologies and designs are discussed, as well as promising directions for remedying some of the present research issues in this area. The presentation is a synopsis of a longer MITRE review article on the same subject.

Montemerlo, Michael S.↗

Fracture in Phenolic Impregnated Carbon Ablator

The thermal protection materials used for spacecraft heat shields are subjected to various thermal-mechanical loads during an atmospheric entry which can threaten the structural integrity of the system. This paper discusses the development of a novel technique to understand the failure mechanisms inside thermal protection materials. The focus of research is Phenolic Impregnated Carbon Ablator (PICA). It has successfully flown on the Stardust spacecraft and is the TPS material chosen for the Mars Science Laboratory (MSL) and Dragon spacecraft. Although PICA has good thermal properties, structurally, it is a weak material. In order to thoroughly understand failure in PICA, fracture tests were performed on FiberForm* (precursor of PICA), virgin and charred PICA materials. Several samples of these materials were tested to investigate failure mechanisms at a microstructural scale. Stress-strain data were obtained simultaneously to estimate the fracture toughness. It was found that cracks initiated and grew in the FiberForm when a critical stress limit was reached such that the carbon fibers separated from the binder. However, both for virgin and charred PICA, crack initiation and growth occurred in the matrix (phenolic) phase. Both virgin and charred PICA showed greater strength values compared to FiberForm coupons, confirming that the presence of the porous matrix helps in absorbing the fracture energy.

Agrawal, Parul↗

Lower Urban Humidity Moderates Outdoor Heat Stress

Abstract Surface temperature is often used to examine heat exposure in multi‐city studies and for informing urban heat mitigation efforts due to scarcity of urban air temperature measurements. Cities also have lower relative humidity, traditionally not accounted for in large‐scale observational urban heat risk assessments. Here, using crowdsourced measurements from over 40,000 weather stations in ≈600 urban clusters in Europe, we show the moderating effect of this urbanization‐induced humidity reduction on outdoor heat stress during the 2019 heatwave. We demonstrate that daytime differences in heat index between urban clusters and their surroundings are weak, and associations of this urban‐rural difference with background climate, generally examined from the surface temperature perspective, are diminished due to moisture feedbacks. We also examine the spatial variability of surface temperature, air temperature, and heat index within these clusters—relevant for detecting hotspots and potential disparities in heat exposure—and find that surface temperature is a poor proxy for the intra‐urban distribution of heat index during daytime. Finally, urban vegetation shows much weaker (∼1/6th as strong) associations with heat index than with surface temperature, which has broad implications for optimizing urban heat stress mitigation strategies. These findings are valid for operational metrics of heat stress for shaded conditions (apparent temperature and humidex), thermodynamic proxies (wet‐bulb temperature), and empirical heat indices. Based on this large‐scale empirical evidence, surface temperature, used due to the lack of better alternatives, may not be suitable for accurately informing heat mitigation strategies within and across cities, necessitating more urban‐scale observations and better urban‐resolving models.

54 ENVIRONMENTAL SCIENCES↗

Space–time structure of weak magnetohydrodynamic turbulence

The two-time energy spectrum of weak magnetohydrodynamic turbulence is found by applying a wave-turbulence closure to the cumulant hierarchy constructed from the dynamical equations. Solutions are facilitated via asymptotic expansions in terms of the small parameter $\varepsilon$ , describing the ratio of time scales corresponding to Alfvénic propagation and nonlinear interactions between counter-propagating Alfvén waves. The strength of nonlinearity at a given spatial scale is further quantified by an integration over all possible delta-correlated modes compliant in a given set of three-wave interactions that are associated with energy flux through the said scale. The wave-turbulence closure for the two-time spectrum uncovers a secularity occurring on a time scale of order $\varepsilon ^{-2}$ , and the asymptotic expansion for the spectrum is reordered in a manner comparable to the one-time case. It is shown that for the regime of stationary turbulence, the two-time energy spectrum exponentially decays on a lagged time scale $(\varepsilon ^2 \gamma _k^s)^{-1}$ in proportion to the strength of the associated three-wave interactions, characterized by nonlinear decorrelation frequency $\gamma _k^s$ . The scaling of the form $k_{\perp } v_0 \chi _0$ exhibited by this frequency is reminiscent of random sweeping by the outer scale with characteristic fluctuation velocity $v_0$ that is modified due to competition with Alfvénic propagation (characterized by $\chi _0$ ) at the said scale. A brief calculation of frequency broadening of the power spectrum due to nonlinear interactions is also presented.

Physics↗

Towards a Unified Setup to Simulate Mid‐Latitude and Tropical Mesoscale Convective Systems at Kilometer‐Scales

Abstract Mesoscale convective systems (MCSs) are the main source of precipitation in the tropics and parts of the mid‐latitudes and are responsible for high‐impact weather worldwide. Studies showed that deficiencies in simulating mid‐latitude MCSs in state‐of‐the‐art climate models can be alleviated by kilometer‐scale models. However, whether these models can also improve tropical MCSs and whether we can find model settings that perform well in both regions is understudied. We take advantage of high‐quality MCS observations collected over the Atmospheric Radiation Measurement (ARM) facilities in the US Southern Great Plains (SGP) and the Amazon basin near Manaus (MAO) to evaluate a perturbed physics ensemble of simulated MCSs with 4 km horizontal grid spacing. A new model evaluation method is developed that enables to distinguish biases stemming from spatiotemporal displacements of MCSs from biases in their reflectivity and cloud shield. Amazon MCSs are similarly well simulated across these evaluation metrics than SGP MCSs despite the challenges anticipated from weaker large‐scale forcing in the tropics. Generally, SGP MCSs are more sensitive to the choice of model microphysics, while Amazon cases are more sensitive to the planetary boundary layer (PBL) scheme. Although our tested model physics combinations had strengths and weaknesses, combinations that performed well for SGP simulations result in worse results in the Amazon basin and vice versa. However, we identified model settings that perform well at both locations, which include the Thompson and Morrison microphysics coupled with the Yonsei University (YSU) PBL scheme and the Thompson scheme coupled with the Mellor‐Yamada‐Janjic PBL scheme.

54 ENVIRONMENTAL SCIENCES↗

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The H I content of red geyser galaxies

Red geysers are a specific type of quiescent galaxy, denoted by twin jets emerging from their galactic centres. These bisymmetric jets possibly inject energy and heat into the surrounding material, effectively suppressing star formation by stabilizing cool gas. In order to confirm the presence and evolutionary consequences of these jets, this paper discusses the scaling, stacking, and conversion of 21-cm H I flux data sourced from the H I -MaNGA survey into H I gas-to-stellar mass (G/S) spectra. Our samples were dominated by non-detections, or galaxies with weak H I signals, and consequently by H I upper limits. The stacking technique discussed successfully resolved emission features in both the red geyser G/S spectrum and the control sample G/S spectrum. From these stacked spectra, we find that on average, red geyser galaxies have G/S of 0.086 ± 0.011 (random) + 0.029 (systematic), while non-red geyser galaxies of similar stellar mass have a G/S ratio of 0.039 ± 0.018 (random) + 0.013 (systematic). Therefore, we find no statistically significant evidence that the H I content of red geysers is different from the general quiescent population.

79 ASTRONOMY AND ASTROPHYSICS↗

QCD thermalization: Ab initio approaches and interdisciplinary connections

Heavy-ion collisions at BNL's Relativistic Heavy Ion Collider and CERN's Large Hadron Collider provide strong evidence for the formation of a quark-gluon plasma, with temperatures extracted from relativistic viscous hydrodynamic simulations shown to be well above the transition temperature from hadron matter. Outstanding problems in QCD include how the strongly correlated quark gluon matter forms in a heavy-ion collision, its properties off-equilibrium, and the thermalization process in the plasma. In this work, we review the theoretical progress in this field in weak coupling QCD effective field theories and in strong coupling holographic approaches based on gauge-gravity duality. We outline the interdisciplinary connections of different stages of the thermalization process to non equilibrium dynamics in other systems across energy scales ranging from inflationary cosmology, to strong field QED, to ultracold atomic gases, with emphasis on the universal dynamics of non-thermal and hydrodynamic attractors. We survey measurements in heavy-ion collisions that are sensitive to the early non-equilibrium stages of the collision and discuss the potential for future measurements. We summarize the current state-of-the art in thermalization studies and identify promising avenues for further progress.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DEReliction: A Cybersecurity Vulnerability Assessment Methodology for Distributed Energy Resources

With the increasing integration of Distributed Energy Resources (DER) into the electric grid, maintaining grid reliability and resilience requires that these devices remain secure. This paper discusses a cybersecurity vulnerability assessment methodology that incorporates best practices from Sandia National Laboratories, SANS Institute, OWASP Foundation, and other web and Internet of Things (IoT) penetration testing (“pen testing”) programs, courses, and frameworks for assessing the security posture of devices. The methodology involves five sequential steps: (1) Collect Public Information, (2) Extract Hardware Details, (3) Inventory Software Components, (4) Identify Vulnerabilities, and (5) Test Vulnerabilities. Each step uncovers potential weaknesses in both hardware and software components of DER devices, considering adversary tactics, techniques, and procedures (TTPs), and potential attack vectors along the way. The results from the execution of this method on multiple residential- and small commercial-scale photovoltaic (PV) inverters reveled hardware and software vulnerabilities, which highlight the benefit of taking a methodical approach to discover vulnerabilities. While the specific vulnerability details are not shared here, a generalized overview of findings underscore the importance of robust security assessments for DER devices. Adoption of an assessment framework of this kind will identify and mitigate cybersecurity threats and bolster the resilience of DER-integrated electric grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗

HectoMAP: The Complete Redshift Survey (Data Release 2)

HectoMAP is a dense redshift survey of 95,403 galaxies based primarily on MMT spectroscopy with a median redshift z = 0.345. The survey covers 54.64 deg 2 in a 1$_{.}^{°}$5 wide strip across the northern sky centered at a decl. of 43$_{.}^{°}$25. We report the redshift, the spectral indicator D n 4000, and the stellar mass. The red-selected survey is 81% complete for 55,962 galaxies with (g - r) > 1 and r < 20.5; it is 72% complete for 32,908 galaxies with (g - r > 1, (r - i) > 0.5 and 20.5 < r < 21.3. Comparison of the survey basis SDSS photometry with the HSC-SSP photometry demonstrates that HectoMAP provides complete magnitude-limited surveys based on either photometric system. We update the comparison between the HSC-SSP photometric redshifts with HectoMAP spectroscopic redshifts; the comparison demonstrates that the HSC-SSP photometric redshifts have improved between the second and third data releases. HectoMAP is a foundation for examining the quiescent galaxy population (63% of the survey), clusters of galaxies, and the cosmic web. HectoMAP is completely covered by the HSC-SSP survey, thus enabling a variety of strong- and weak-lensing investigations.

79 ASTRONOMY AND ASTROPHYSICS↗