Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “exponential integrators”

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 343 records · Page 19

A consumer-centric approach to quantify efficiency of receiving goods purchased via online

Virtual participation in shopping activities has increased exponentially in the past four years compared to the last couple of decades. E-tailing or online shopping offers the convenience of goods reaching a consumer instead of a consumer traveling to a store, but it has downsides like geographical service variability and negative social externalities such as increased energy consumption and emissions. This study proposes a novel approach to quantify e-tailing efficiency from the consumers’ viewpoint. The methodology is innovative in its integration of accessibility theory with energy and cost impedance factors and consideration of delivering and picking up goods purchased via online. The methodology is implemented for the San Francisco Bay Area and is subject to scenarios that see enhancements to various facets of online shopping delivery. Results demonstrate that increasing the frequency of e-commerce deliveries helps improve e-tailing efficiency in rural locations, while improvements in energy and cost aspects of delivery modes are seen to improve e-tailing efficiencies in the central parts of the region. The approach proposed can provide valuable insights on where people have limited benefits from online shopping and how emerging delivery mechanisms can change the quality of the e-commerce experience within a city. This research offers a replicable framework for assessing e-commerce systems in diverse geographic contexts, contributing to the development of equitable and environmentally sustainable urban freight systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Microscopic Origin of the Entropy of Black Holes in General Relativity

We construct an infinite family of microstates with geometric interiors for eternal black holes in general relativity with a negative cosmological constant in any dimension. Wormholes in the Euclidean path integral for gravity cause these states to have small, but nonzero, quantum mechanical overlaps that have a universal form. The overlaps have a dramatic consequence: The microstates span a Hilbert space of log dimension equal to the Bekenstein-Hawking entropy. The semiclassical microstates we construct contain Einstein-Rosen bridges of arbitrary size behind their horizons. Our results imply that all these bridges can be interpreted as quantum superpositions of wormholes of size at most exponential in the entropy.

Physics↗

Some simple properties of stellar pulsation modes

Except for the lowest harmonics, small-amplitude stellar pulsation modes possess many simple properties whose evaluation does not require numerical integration of the fourth-order equations of motion. All antinodes tend to have the same total kinetic energy except for those lying near physical or geometric boundaries. However, when kinetic energy per unit volume is considered, order-of-magnitude enhancements are seen in antinodes lying near the center of the star, and factor-of-2 enhancements occur near the polar axis. The nodes are distributed very regularly along the radius. They follow an exponential law in g-regions, and their separation is proportional to the sound travel time in p-regions. A simple graphical procedure is described for surveying the oscillation frequencies of a new stellar model. A precise condition is derived giving the division of energy between radial and angular motion. Another condition gives the fractional contribution to the velocity field of its two sources, the divergence and the curl. Certain simplifying results of weak coupling among the linear modes are briefly described.

Wolff, C. L.↗

The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z = 4 to 8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Hubble Frontier Field deep parallel observations near the Abell 2744 and MACS J0416.1- 2403 clusters. The combination of these surveys provides an effective volume of 0.6-1.2 ×10(exp 6) Mpc(exp 3) over this epoch, allowing us to perform a robust search for bright (M(sub UV) less than −21) and faint (M(sub UV) = −18) galaxies. We select galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5 less than z less than 8.5, with more than 1000 galaxies at z of approximately 6 - 8. We measure both a stepwise luminosity function for galaxies in our redshift samples, as well as a Schechter function, using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our UV luminosity functions agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z of greater than or equal to 6. Our bestfit value of the characteristic magnitude M* is consistent with −21 at z of greater than or equal to 5, different than that inferred based on previous trends at lower redshift. At z = 8, a single power-law provides an equally good fit to the UV luminosity function, while at z = 6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare our luminosity functions to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift, though a decreasing impact of feedback may also be possible. We measure the evolution of the cosmic star-formation rate (SFR) density by integrating our observed luminosity functions to M(sub UV) = −17, correcting for dust attenuation, and find that the SFR density declines proportionally to (1 + z)((exp −4.3)(+/-)(0.5)) at z greater than 4, consistent with observations at z greater than or equal to 9. Our observed luminosity functions are consistent with a reionization history that starts at redshift of approximately greater than 10, completes at z greater than 6, and reaches a midpoint (x(sub HII) = 0.5) at 6.7 less than z less than 9.4. Finally, using a constant cumulative number density selection and an empirically derived rising star-formation history, our observations predict that the abundance of bright z = 9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z similar to 10 galaxies.

Galaxy↗

Improved Constraints for the XUV Luminosity Evolution of Trappist-1

We re-examine the XUV luminosity evolution of TRAPPIST-1 utilizing new observational constraints of the current stellar parameters (XUV and bolometric luminosity) from multi-epoch X-ray/UV photometry. Following the formalism presented on Fleming et al. (2020), we infer that TRAPPIST-1maintained a saturated XUV luminosity, relative to the bolometric luminosity, of log10(LXUV/Lbol)=−3.03+0.25−0.23at early times for a period of tsat= 3.14+2.22−1.46Gyr. After the saturation phase, we find LXUV decayed over time by an exponential rate ofβXUV=−1.17+0.27−0.28. Compared to our inferred age of the system, age = 7.96+1.78−1.87Gyr, our result for tsat suggests that there is only a∼4% chance thatTRAPPIST-1 still remains in the saturated phase today, which is significantly lower than the previous estimate of 40%. Despite this reduction in tsat, our results remain consistent in the conclusion that the TRAPPIST-1 planets likely received an extreme amount XUV energy—an estimated integrated XUV energy of∼1030−1032erg over the star’s lifetime—that is∼15% lower than the original result.

Stellar evolution↗

Modeling Cyber Supply Chain Incidents with Multilayered Graph Motifs

As noted within the literature, supply chain includes people and organizations---manufacturers, integrators, and third-party vendors---that are involved in one or more stages of a product lifecycle. Since supply chains, by definition, include organizations and people, supply chain risk management activities must consider dependencies between an organization's business processes and third-party resources. Just as adversarial tactics can be implemented via techniques implemented via networked computer systems, so can such tactics be expressed via legal business relationships. A cyber incident may have an exponential impact downstream, for example, by leveraging a product's distribution channel (e.g. malicious updates in SolarWinds, buggy updates in CrowdStrike). Similarly, legitimate and legal business relationships also affect the attack surface exposure of systems, enabling long-term persistence and/or unknown impacts to product quality that are hard to detect. This paper catalogs several recent digital supply chain incidents and applies a multilayered network formalism to develop structural indicators (graph motifs) that reflect potentially-adversarial behavior. Finally, we compare and contrast the characteristics of adversarial tactics (e.g. Loss of Availability, Data Collection) that leverage cyber-physical dependencies to those that leverage legal organizational relationships.

97 - MATHEMATICS AND COMPUTING↗

Sensitivity analysis of the effect of wind and wake characteristics on wind turbine loads in a small wind farm

Abstract. Wind turbines are designed using a set of simulations to determine the fatigue and ultimate loads, which are typically focused solely on unwaked wind turbine operation. These structural loads can be significantly influenced by the wind inflow conditions. Turbines experience altered inflow conditions when placed in the wake of upstream turbines, which can additionally influence the fatigue and ultimate loads. It is important to understand the impact of uncertainty on the resulting loads of both unwaked and waked turbines. The goal of this work is to assess which wind-inflow-related and wake-related parameters have the greatest influence on fatigue and ultimate loads during normal operation for turbines in a three-turbine wind farm. Twenty-eight wind inflow and wake parameters are screened using an elementary effects sensitivity analysis approach to identify the parameters that lead to the largest variation in the fatigue and ultimate loads of each turbine. This study uses the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine, simulated with OpenFAST and synthetically generated inflow based on the International Electrotechnical Commission (IEC) Kaimal turbulence spectrum with the IEC exponential coherence model using the NREL tool TurbSim. The focus is on sensitivity to individual parameters, though interactions between parameters are considered, and how sensitivity differs between waked and unwaked turbines. The results of this work show that for both waked and unwaked turbines, ambient turbulence in the primary wind direction and shear are the most sensitive parameters for turbine fatigue and ultimate loads. Secondary parameters of importance for all turbines are identified as yaw misalignment, streamwise integral length, and the exponent and streamwise components of the IEC coherence model. The tertiary parameters of importance differ between waked and unwaked turbines. Tertiary effects account for up to 9.0 % of the significant events for waked turbine ultimate loads and include veer, non-streamwise components of the IEC coherence model, Reynolds stresses, wind direction, air density, and several wake calibration parameters. For fatigue loads, tertiary effects account for up to 5.4 % of the significant events and include vertical turbulence standard deviation, lateral and vertical wind integral lengths, non-streamwise components of the IEC coherence model, Reynolds stresses, wind direction, and all wake calibration parameters. This information shows the increased importance of non-streamwise wind components and wake parameters in the fatigue and ultimate load sensitivity of downstream turbines.

17 WIND ENERGY↗

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils from Combined Active-Passive Microwave SMAP Observations

An active-passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height (s) and horizontal correlation length (l)] is presented for bare soils. The algorithm is based on active-passive microwave covariation, including the improved Integral Equation Method (I²EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for εₛ > 10 [-]. Furthermore, the physical model basis of this approach (I²EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

correlation length↗

Implicit–explicit multirate infinitesimal stage-restart methods

Implicit–Explicit (IMEX) methods are flexible numerical time integration methods which solve an initial-value problem (IVP) that is split into stiff and nonstiff processes with the goal of lower computational costs than a purely implicit or explicit approach. A complementary form of flexible IVP solvers are multirate infinitesimal methods for problems split into fast- and slow-changing dynamics, that solve a multirate IVP by evolving a sequence of “fast” IVPs using any suitably accurate algorithm. This article introduces a new class of high-order implicit–explicit multirate methods that are designed for multirate IVPs in which the slow-changing dynamics are further split in an IMEX fashion. This new class, which we call implicit–explicit multirate infinitesimal stage-restart (IMEX-MRI-SR), both improves upon the previous implicit–explicit multirate infinitesimal generalized-structure additive Runge Kutta (IMEX-MRI-GARK) methods by allowing for far easier creation of new embedded methods, and extends multirate exponential Runge Kutta (MERK) methods by allowing the fast-changing dynamics to be nonlinear and the methods to be implicit. We leverage GARK theory to derive conditions for orders of accuracy up to four, and we provide second- and third-order accurate example methods, which are the first known embedded MRI methods with IMEX structure. We then perform numerical simulations demonstrating convergence rates and computational performance in both fixed-step and adaptive-step settings.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Blockchain-Enabled Secure Device-to-Device Communication in Software-Defined Networking

The Internet of Things (IoT) continues to increase the demand for seamless communication among IoT devices. The rapid growth of IoT devices has led to an exponential increase in device-to-device (D2D) communication within the Software-Defined Networking (SDN), though it enables a flexible archi-tecture for managing network resources. However, traditional security models face challenges (e.g., Security, privacy, and trust) in addressing the dynamic and decentralized nature of these communications. Despite of these challenges, this paper proposes a novel approach that leverages blockchain technology to enhance the security, privacy, and trustworthiness of D2D communication within an SDN environment. The proposed approach integrates blockchain nodes in sDN components to establish a decentralized ledger for transparent and verifiable records. Smart contracts enforce authentication rules to ensure that only authenticated devices can access the network and engage in transactions securely. It also automates the security policies to ensure temper resistance execution using the cryptographic mechanism for data integrity and authentic communication. The Implementation of the proposed algorithms validates the resilience of the proposed approach against cyberattacks. Overall, the proposed approach enables efficient and secure D2D communication for resilient SDN infrastructure in IoT ecosystems.

Das, Debashis↗

Analyses Made to Order: Using Transformation to Rapidly Configure a Multidisciplinary Environment

Aerospace problems are highly multidisciplinary. Four or more major disciplines are involved in analyzing any particular vehicle. Moreover, the choice of implementation technology of various subsystems can lead to a change of leading domain or reformation of the driving equations. An excellent example is the change of expertise required to consider aircraft built from composite or metallic structures, or those propelled by chemical or electrical thrusters. Another example is in the major reconfiguration of handling and stability equations with different control surface configuration (e.g., canards, t-tail v four-post tail). Combinatorial problems are also commonplace anytime that a major system is to be designed. If there are only 5 attributes of a design to consider with 4 different options, this is already 1024 options. Adding just 5 more dimensions to the study explodes the space to over one million. Even generous assumptions like the idea that only 10% of the combinations are physically feasible can only contain the problem for so long. To make matters worse, the simple number of combinations is only the beginning. Combining the issue of trade space size with the need to reformulate the design problem for many of the possibilities makes life exponentially more difficult. Advances in software modeling approaches have led to the development of model-driven architecture. This approach uses the transformation of models into inferred models (e.g. inferred execution traces from state machines) or the skeletons for code generation. When the emphasis on transformation is applied to aerospace, it becomes possible to exploit redundancy in the information specified in multiple domain models into a unified system model. F1urther, it becomes possible to overcome the combinatorial nature of specifying integrated system behavior by manually combining the equations governing a given component technology. Transformations from a system specification combined with a system-analysis mapping specification enable one-click combination of domain analyses. This is a flexibility that has been missing from many engineering codes, which often entangle design specification and physical examination much more than is required to conduct the analysis. This capability has been investigated and cultivated within the DARPA F6 program by a team of JPL and Phoenix Integration engineers building the Adapatable Systems Design and Analysis (ASDA) framework. By embracing system modeling with SysML and the Query-View-Transformation (QVT) language, the ASDA team has been able to build a flexible, easily reconfigurable framework for building up and solving large tradespaces. Examples of application and lessons learned in building the framework will be described in this paper. In addition, the motivation will be laid for various tool vendors to develop open model description standards while being able to maintain competitive advantage through proprietary algorithms and approaches. These standards will also be compared to the underpinnings of model-driven architecture and the OMG standards of the Meta-Object Facility (MOF), SysML, and QVT.

Cole, Bjorn↗

Convective instability in differentially rotating disks

A normal mode analysis for nonaxisymmetric perturbations in a thin, differentially rotating disk with a vertical structure that is isothermal and convectively unstable is performed. The vertical gravity is assumed to be external and constant. The perturbation scale is assumed to be much shorter than the radius of the disk but comparable to or less than the thickness. The initial value problem is formulated in shearing coordinates. Dispersion relations are obtained for the three limiting cases of zero shear, axisymmetric perturbations, and small radial wavelengths. The full effects of shear are studied by integrating numerically the initial value problem. Nonaxisymmetric local Fourier modes are found to have a radial wavenumber that increases linearly with time in proportion to the shear times the azimuthal wavenumber. While Coriolis forces exert stabilizing effects on the convective modes, reducing their growth rate and the range of unstable wavelengths, shear has destablizing effects inasmuch as it reduces the epicyclic frequency at a given angular velocity. In a Keplerian disk, perturbations with azimuthal wavelengths about 2 times smaller than vertical wavelengths grow exponentially.

Ryu, Dongsu↗

Reconfigurable quantum phononic circuits via piezo-acoustomechanical interactions

Abstract We show that piezoelectric strain actuation of acoustomechanical interactions can produce large phase velocity changes in an existing quantum phononic platform: aluminum nitride on suspended silicon. Using finite element analysis, we demonstrate a piezo-acoustomechanical phase shifter waveguide capable of producing ± π phase shifts for GHz frequency phonons in 10s of μm with 10s of volts applied. Then, using the phase shifter as a building block, we demonstrate several phononic integrated circuit elements useful for quantum information processing. In particular, we show how to construct programmable multi-mode interferometers for linear phononic processing and a dynamically reconfigurable phononic memory that can switch between an ultra-long-lifetime state and a state strongly coupled to its bus waveguide. From the master equation for the full open quantum system of the reconfigurable phononic memory, we show that it is possible to perform read and write operations with over 90% quantum state transfer fidelity for an exponentially decaying pulse.

Taylor, Jeffrey C. (ORCID:0000000215982436)↗

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↗

Reduced-order modeling of near-field THMC coupled processes for nuclear waste repositories in shale

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke’s law models) are implemented to describe mechanical characteristics of the system. Overall, our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Integration of EDWARD readout architecture in full-field fluorescence imaging detector

Data bandwidth, timing resolution and resource utilization in readouts of radiation detectors are a constant challenge. Event driven solutions are pushing against well-trenched framed solutions. The idea for an asynchronous readout architecture called EDWARD ( E vent- D riven W ith A ccess and R eset D ecoder) was presented at the TWEPP 2021 conference. Here we show the progress of our work which resulted in two chip prototypes. The first one, named 3FI65P1, is a full device with the analog pixel circuitry suited for full-field fluorescence imaging. It is already manufactured, and preliminary results are presented. Finally, the second chip, named EDWARD65P1, contains digital pulse generators with Poisson-exponential distribution in each pixel for extraction of the performance matrix of the EDWARD architecture alone.

47 OTHER INSTRUMENTATION↗

Nanoparticle-Enhanced Plasma Discharge Using Nanosecond High-Voltage Pulses

By discharging nanosecond high-voltage (5 kV) pulses across an insulating substrate containing Au, Pt, or Cu nanoparticles, a 3 order of magnitude (1000×) enhancement in the generation of plasma can be achieved through local field enhancement on the surface of the nanoparticles. The low-temperature nature of this transient plasma is crucial to maintaining the structural integrity of these delicate nanoparticles. These nanoparticles provide up to a 1000-fold enhancement in the generation of the plasma, which is localized to the surface of the nanoparticles where it is potentially useful (e.g., for catalysis). Additionally, we performed both time-domain and frequency-domain calculations of the electromagnetic response of the nanoparticles based on high-resolution transmission electron microscope (HRTEM) images, which show local field enhancement of the nanosecond high-voltage pulse on the order of 3×. Since the plasma initiation depends exponentially on the peak electric field strength, this 3-fold increase in the local electric field can result in a several orders of magnitude increases in the generation of plasma at a given applied external field strength. In order to rule out plasmon-resonance enhancement, which is often associated with small metal nanoparticles, we performed finite difference time domain (FDTD) simulations in the optical frequency domain, which show that the effect of plasmon resonance is negligible for Pt nanoparticles. We therefore attribute the nanoparticle-based enhancement to the generation of plasma (an electrostatic effect) rather than enhanced coupling of light from the near field to the far field via the plasmon resonance phenomenon (an optical effect).

36 MATERIALS SCIENCE↗

Use of Diffusion Bonded Cu Strap and Integrated MLI for Thermal Control of 100 K IR Detector on L’Ralph Instrument

Passive cooling of cryogenic instruments is one of the most challenging aspects of spaceflight thermal control systems. They are highly sensitive to parasitic heat leaks from their warmer environment, especially the spacecraft components. The L’Ralph instrument on the Lucy mission is an example of such a system, with its LEISA detector requiring passive cooling to temperatures below 112 K. MLI performance can be measured in terms of estar, and analysts apply bias between a high and low estar value. This approach works well for conventional MLI at room temperature, however, it can be dangerous when a mission goes through a wide range of temperatures during its lifetime, with estar values increasing exponentially as temperatures get colder. Various estar correlations are presented as well as how L’Ralph is approaching the MLI problem with the use of IMLI. IMLI uses discrete spacers to isolate each layer, and doesn’t require intermediate layers such as dacron netting. This results in estar performance that can be more precisely estimated, and yields valuesaround 0.004 at the 112K operating temperature. L’Ralph is using IMLI on the backside of the radiator, which permits a more effective use of increased radiator area, by minimizing the area dependent heat leak impact from the instrument backloading. The detector is thermally coupled to the radiator via a diffusion bonded Cu strap. The diffusion bonding process chemically bonds the Cu foil ends in order to maximize heat transfer effectiveness across the multiple foils used. Employing a pressure compensation design via the use of Ti blocks and a clamp at the strap ends ensures that as temperatures go colder, pre-load is maintained in order to minimize resistance at the strap ends. The use of both of these technologies, and how they work together are crucial to the success of the L’Ralph thermal control system.

Thermal↗