SEARCH · Engineering Papers
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.
A new strategy is presented for computing anharmonic partition functions for the motion of adsorbates relative to a catalytic surface. Importance sampling is compared with conventional Monte Carlo.
The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.
Abstract After large-magnitude earthquakes, a crucial task for impact assessment is to rapidly and accurately estimate the ground shaking in the affected region. To satisfy real-time constraints, intensity measures are traditionally evaluated with empirical Ground Motion Models that can drastically limit the accuracy of the estimated values. As an alternative, here we present Machine Learning strategies trained on physics-based simulations that require similar evaluation times. We trained and validated the proposed Machine Learning-based Estimator for ground shaking maps with one of the largest existing datasets (<100M simulated seismograms) from CyberShake developed by the Southern California Earthquake Center covering the Los Angeles basin. For a well-tailored synthetic database, our predictions outperform empirical Ground Motion Models provided that the events considered are compatible with the training data. Using the proposed strategy we show significant error reductions not only for synthetic, but also for five real historical earthquakes, relative to empirical Ground Motion Models.
Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.
We revisit the Kibble–Robinson theory, first proposed in 2014 by Kibble and Robinson. This theory significantly simplifies the construction and operation of Kibble balances. We conducted a theoretical investigation of the theory’s assumptions, using a corner cube as the optical target in the interferometer for velocity measurement. We find that it is advantageous to build a mechanism whose output has minimal rotation and horizontal motion. For balances with relative uncertainty targets below $1\times10^{-6}$, the mass pan and the optical target should be suspended from a common gimbal so they have the same vertical velocity and no rotation. In this case, the measurement bias due to Abbe offset is minimized, and the bias due to corner loading is repeatable.
Resonant drag instabilities (RDIs) in protoplanetary discs are driven by the aerodynamic back-reaction of dust on gas and occur when the relative dust-gas motion resonate with a wave mode intrinsic to the gas fluid. Axisymmetric models indicate that the RDI generates filamentary perturbations, leading to grain clumping and planetesimal formation. Motivated by these findings, we investigate the dust–gas interaction in a non-axisymmetric inviscid protoplanetary disc with an embedded low-mass planet ($M_{\mathrm{p}}\in [0.3, 3] \, \rm M_{\oplus }$, here $\, \rm M_{\oplus }$ is the Earth mass). We conduct global 3D high-resolution two-fluid simulations, with the dust being parametrized by the Stokes number $\mathrm{St}\in [0.01,0.5]$. We find that planetary waves (PWs; also known as Rossby waves), which propagate along the downstream separatrices of the horseshoe region, resonate with the streaming motion and trigger the RDI. The consequent development of a global-scale filamentary dust distribution does not sensitively depend on the Stokes number, nor does it depend on the fast dust settling that takes place in an inviscid disc. The rapid onset of this instability, which is comparable to the dynamical orbital time-scale, suppresses the formation of asymmetric structures in the dust in the vicinity of the planet (such as dust voids and filaments). Additionally, we find that the dust feedback enables buoyancy resonances in an otherwise non-buoyant (globally isothermal) disc. Therefore, our results provide the first numerical evidence of RDIs generation driven by planetary waves.
We discover a surprising relation between the collective motion of nucleons within atomic nuclei, traditionally understood to be driven by long-range correlations, and short-range nucleon-nucleon interactions. Specifically, we find that quadrupole collectivity in low-lying states of 6 Li and 12 C, calculated with state-of-the-art ab initio techniques, is significantly influenced by two opposing 𝑆-wave contact couplings that subtly alter the surface oscillations of one largely deformed nuclear shape, without changing that shape’s overall contribution within the nucleus. The results offer new insights into the nature of emergent nuclear collectivity and its link to the underlying nucleon-nucleon interaction at short distances.
Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.
The Idaho National Laboratory (INL) Seismic Monitoring Network is located in eastern Idaho and monitors a portion of the intermountain seismic belt. It has been in place for 50 yr and has undergone several major changes, the most recent of which has been the transition to the Antelope real‐time acquisition system and the implementation of automatic phase picking algorithms to aid in analysis. This study discusses the efforts to evaluate the performance of the INL seismic monitoring network (and other surrounding stations) using the new real‐time acquisition system. The method outlined by Wilson et al. (2021) is used to develop an empirical relationship between the observability of local earthquakes as a function of magnitude and distance. This relationship is used to produce detection thresholds for Pwaves for all stations of interest. The INL seismic network has two main goals: monitor tectonic‐and volcanic‐related events and measure ground motions for input into seismic hazard analysis. Because of these two overall objectives, several seismic stations have been installed near critical facilities and, therefore, are not as quiet as stations that are used primarily for earthquake detection. This is reflected in their detection thresholds, which are much smaller for stations away from facilities. This study shows that the INL Seismic Monitoring Network is able to detect earthquakes near INL facilities with M L > 1.2, with redundancies built in to ensure this sensitivity even if data became unavailable from some stations. This study also shows “holes” in the monitoring network where the detection of smaller earthquakes is highly dependent on sparsely placed seismic stations. In conclusion, the results of this study will be used to govern plans for expansion of earthquake monitoring in Idaho and the surrounding region and to fine‐tune the detection thresholds for individual stations.
This dataset provides long-duration large-eddy simulations (LES) of Hurricane Laura (2020), capturing the turbulent boundary-layer wind field during the storm's passage at horizontal grid spacings of 55.55 m and 10 m. Hurricane Laura made landfall as a Category 4 storm in southwestern Louisiana at 06:00 UTC on 27 August 2020, producing widespread damage to coastal infrastructure. The simulations characterize the spatial and temporal variability of the turbulent wind at heights spanning the lowest 300 m of the boundary layer, throughout the approach and passage of the storm's eyewall over a fixed location. To resolve the spatial heterogeneity of wind and turbulence conditions across the storm, seven independent refined LES runs are performed at radial locations spanning the storm relative to its direction of motion. The refinement positions are defined on a normalized radial coordinate r ̂=r/R, where R≈21 km is Laura's radius of maximum wind, and lie between r ̂=-1.5 and r ̂=+1.5 in ∆r ̂=0.5 increments along a line perpendicular to the storm's track at the midpoint of the simulation. Negative r ̂ corresponds to locations south-southwest of the storm center and positive r ̂ to locations north-northeast. The dataset is segmented into seven radial refinement groups.
During testing of the Stress Managed Cosine-Theta dipole mirror magnet SMCTM1, magnet quenches were observed following large voltage spikes in the half-coil voltage taps which preceded normal quench initiation. Following some recent work at CERN measuring power current converter transients in relation to flux jumps, an additional differential current probe was added to the magnet instrumentation. Measurements of this probe, as well as the half coil voltage taps show a fast and substantial current change preceding the quench, hinting at coil motion. A simple model based on relative shift of the inner and outer layer coils was developed which suggests the magnitude of coil shift events is on the estimated relative coil shift from the magnet constraint conditions. The behavior is shown to be reversible and characteristic of stick-slip.
Using the equations of motion of QCD and Lorentz invariance relations, we show a new, more direct way of obtaining the decomposition of the proton's quark transverse angular momentum into its spin and orbital components. The new decomposition can be understood as a twist-three relation, with the quark transverse spin described by the parton distribution g T . We find that the orbital component can be expressed in terms of a moment in the quark transverse momentum, k T , of the generalized transverse momentum-dependent distribution F 12 , or alternatively in terms of the twist-three generalized parton distributions H 2T and $\tilde{E}$ 2T .
The American Bureau of Shipping performed an evaluation of national maritime nuclear testing capabilities and determined there currently is an inability to test maritime nuclear power plants in a ship motion environment. A Ship Motion Test Facility (SMTF) is proposed to fill this gap. To inform design of this facility, this report documents an evaluation performed by a panel of technical experts that identifies key physical phenomena associated with maritime nuclear power plant operation that are important to safety and reliability and have a low state of knowledge. This evaluation concluded that all of the major test needs to understand a nuclear plant's behavior in a ship motion environment are for thermal-hydraulics related phenomena. Based on this finding, the evaluation panel concluded that the SMTF would not be required to use nuclear fuel to provide a heat source for the fluids. This would significantly decrease the cost, timeline, and complexity for this test facility to enable the advancement of maritime nuclear technology.
This TEAMER project had three main objectives related to characterizing the performance of a vortex-induced vibration (VIV) marine energy harvesting device called the Whatever Input To Torsion (WITT). The first was to characterize the motion of the WITT device as it vibrates on top of a pipe in a natural tidal flow. The second was to measure how much electricity the WITT device generates when it is deployed on top of a pipe that is secured to the seabed with a lander. The third objective was to measure how much electricity the WITT device generates in the same configuration but with a longer pipe. This TEAMER project succeeded in satisfying the first objective by measuring the motion of the WITT device with an inertial measurement unit (IMU) mounted in the WITT housing. The motion was described by examining changes in the acceleration, pitch, roll, and yaw, and displacement. The frequency of the motion was characterized using spectral analysis. The TEAMER project succeeded in reaching the second objective by measuring the power produced by the WITT. It is reported here as the average and maximum power in hourly intervals. WITT Energy decided to not have us pursue the third objective of testing a longer pipe because the first test was unexpectedly successful at having the system vibrate at the target frequency. A longer pipe would have resulted in a decreased frequency and therefore not have been a useful test. Based on these findings, the decision was made to keep the same pipe length for both tests and improve on other aspects of the test structure that had failed during the first test. The key finding of this TEAMER project is that the WITT is capable of generating electricity in moderate flow speeds when it is in a housing that is attached to a 2.5 m pipe and secured to the seabed by a lander. This report contains the first power values for a VIV instrument that generates electricity from a pendulum swinging at the top of a long pipe. The experiment demonstrated that the vibration frequency of the device varied with the tidal current speed, which changed throughout a given tidal cycle. This meant that the desired peak frequency was not sustained for more than one hour. The success of this TEAMER test provides evidence that continued research in the area of electricity generation from VIV holds promise as a source of marine energy for offshore applications. Future tests should focus on refining the pipe length equations and improving the power take off system of the WITT power electronics. In addition, experiments aimed at expanded testing and evaluation (e.g. longer than one month) are needed to measure the power production at higher tidal current speeds and test the system's durability.
In this project we use ground-based DOE ARM observations, satellite observations, and Weather Research and Forecasting (WRF) model simulations to arrive at a more complete understanding (or, new theory) of how atmospheric heating profiles in thunderstorm cores drives how quickly air rises in the storms. The rate at which air motion increases or decreases with altitude directly relates to how much cloud mass is ejected laterally into the adjacent atmosphere, which then influences the development and maintenance of “stratiform raining cloud anvils”. Our work has led to new equations connecting the heating in these thunderstorms to these stratiform clouds, with the connection modulated by the stability structure of the atmosphere (i.e., if more stable, then a more rapid change in thunderstorm heating with altitude results in more stratiform cloud, independent of whether over land or ocean, or how organized the thunderstorms are). Improving knowledge of stratiform clouds is important, since they contribute to substantial accumulated rainfall (important to agriculture and local hydrology) as well as radiation, given their large areal extents and lifetime. The radiation impacts likely impact simulated cloud feedback in all Earth System models (ESMs).
This study reports the first use of in-situ synchrotron X-ray diffraction (SXRD) to study the effects of ultrasonic melt processing (USMP) on phase and grain size evolution during solidification in a commercial Al–Si–Mg casting alloy. USMP is a technique that, when applied to aluminum as it solidifies, can be used to refine the local microstructure of large-scale castings. Analysis of the in-situ SXRD data to estimate the average grain size of primary α-Al grains during USMP demonstrates that USMP slows the growth rate of α-Al grains and reduces grain size by 36 pct. Furthermore, there is also evidence that USMP causes the primary α-Al grains to move relative to the X-ray beam; such motion increases the probability of primary α-Al grains colliding and fragmenting. This movement becomes constrained at the onset of the Al–Si binary eutectic, suggesting that USMP ceases to effectively refine the microstructure once the Al–Si binary eutectic begins to form. Complementary laboratory-scale X-ray diffraction (XRD) data were used to correlate the lattice parameters of the α-Al and Si (D-A4) phases with temperature to estimate cooling rate during solidification. Thus, this study can guide the design of novel castings with spatially distributed fine-grained regions produced using local ultrasonic processing.
Here, this paper develops a novel method for reconstructing the full-field response of structural dynamic systems using sparse measurements. The singular value decomposition is applied to a frequency response matrix relating the structural response to physical loads, base motion, or modal loads. The left singular vectors form a non-physical reduced basis that can be used for response reconstruction with far fewer sensors than existing methods. The contributions of the singular vectors to measured response are termed singular-vector loads (SVLs) and are used in a regularized Bayesian framework to generate full-field response estimates and confidence intervals. The reconstruction framework is applicable to the estimation of single data records and power spectral densities from multiple records. Reconstruction is successfully performed in configurations where the number of SVLs to identify is less than, equal to, and greater than the number of sensors used for reconstruction. In a simulation featuring a seismically excited shear structure, SVL reconstruction significantly outperforms modal FRF-based reconstruction and successfully estimates full-field responses with as few as two uniaxial accelerometers. SVL reconstruction is further verified in a simulation featuring an acoustically excited cylinder. Finally, response reconstruction and uncertainty quantification are performed on an experimental structure with three shaker inputs and 27 triaxial accelerometer outputs.
This project focussed on two aspects of computational modeling with a view toward application for photovoltaic design: (i) increased reliability of exchange-correlation functionals in time-dependent density functional theory (TDDFT) (a method of choice of the calculation of electronic spectra and dynamics), especially for time-resolved non-perturbative dynamics, and (ii) the development of a practical but rigorously-based method for coupling electronic and nuclear motion via the exact-factorization (EF) approach (a relatively new framework for developing approximations).