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.
This project explored the feasibility of using fiber-optic sensing technologies for condition monitoring of offshore synthetic mooring lines used in marine renewable energy systems. The work focused on evaluating optical fiber materials, interrogation methods, and rope-embedment approaches capable of enabling distributed strain or elongation measurements along mooring lines.
Two-dimensional (2D) materials inherently exhibit instabilities. Structurally, this may lead to modulations along the third dimension, e.g., wrinkles. Electronically, 2D instabilities can manifest themselves as charge density waves (CDWs). Although wrinkles can alter anisotropic electronic structures susceptible to forming CDWs, less is known about their impact on broken-symmetry ground states. Here, in this study, using scanning tunneling microscopy and spectroscopy, we investigate the CDW states on the wrinkled surface of DyTe 3 . We identify elongated, parallel nanoscale wrinkles stabilized by ribbon-shaped defects. Interestingly, the CDW order persists across the nanowrinkles with a gradual phase shift but is locally suppressed near the defects, where phase windings occur. In addition, these defects induce quantum confinement effects along the nanowrinkles, indicating the presence of one-dimensional metallic states with hole- like dispersion, while angle-resolved photoemission spectroscopy identifies a gap along the wrinkle direction. We ascribe this discrepancy to strain-induced changes in the Fermi surface, which lead to the closure of the gap at the sites of the nanowrinkles. Taken together, our results underscore the complex interplay between structural features and Fermi surface topology, allowing for the deliberate manipulation of quantum states in strongly correlated systems via local crystal deformations.
Shock wave (Hugoniot), shock temperature, and release data are presented for several geophysically important, refractory materials. A sensitive multichannel optical pyrometer was developed to measure shock temperatures (2500 to 5600 K at pressures from 48 to 117 GPa) in anorthite (CaAl2Si2O8) glass. Shock temperatures of 3750 to 6000 K at pressures from 140 to 182 GPa were measured in calcium oxide (CaO). Temperature data were used to constrain the energetics of the B1-B2 phase transition at 70 GPa in CaO, and to construct a finite strain equation of state for CaO consistent with previous Hugoniot data. The CaO equation of state was used with equation of state parameters of other oxides to construct a theoretical mixed oxide Hugoniot of anorthite, which is in agreement with new Hugoniot data above about 50 GPa, determined using experimental techniques developed. The mixed oxide model, however, overestimates the shock temperatures, and does not accurately predict measured release paths.
Abstract Earthquake triggering by seismic waves has been recognized as a phenomenon for nearly 30 years. However, our ability to study dynamic triggering has been limited by our ability to capture the triggering stresses accurately and record the resultant earthquakes. Here we use full waveforms from a dense seismic network and a modern, high‐resolution seismic catalog to measure triggering in Southern California from 2008 to 2017 based on interevent time ratios. We find that the fractional seismicity rate change, which we term triggering intensity or triggerability, as a function of peak strain change for the period of ∼20 s due to distant earthquakes is monotonically increasing and compatible with earlier measurements made with a disjoint data set from 1984 to 2008. A triggering strain of 1 microstrain is equivalent to the local productivity generated by an M 1.8 earthquakes. This result implies that a prediction of seismicity rate changes can be made based on recorded ground shaking using the same formalism as currently used for aftershock prediction. For a teleseismic event, this small level of triggering occurs throughout the region and thus aggregates to a regional effect. We find that the triggering rate decays after the triggerer follows an Omori‐Utsu law, but at a much slower rate than a typical aftershock sequence. The slow decay rate suggests that an ancillary process such as creep or fluid flow must be part of dynamic triggering. The prevalence of triggering in areas of creep or fluid involvement reinforces this inference. A triggering cascade of secondary earthquakes is insufficient to explain the data.
"Prediction and subsequent discovery of topological insulators is considered to be one of the main results in condensed matter physics in the last decades. Not surprisingly, it has received major attention of both researchers and funding agencies. This attention is well-deserved; yet, one cannot but note that nearly all this research revolves about essentially the same concept: electronic excitations with linear dispersion, covering, of course, such diverse and intriguing phenomena as topological insulator, bulk Dirac states (or Weyl, if not spin degenerate), Mayorana fermions. In this project, we will address, mainly, other topological phenomena, such as topologically nontrivial magnetic patterns (as, for instance, topological Hall and related phenomena). Specifically, we propose three interrelated trusts: (1) Time-reversal symmetry breaking nonrelativistic antiferromagnets, called altermagnets. These are materials that break Kramers degeneracy of electronic bands, despite having zero net magnetization by symmetry and being fully collinear, and not necessarily non-centrosymmetric. The corresponding band structure is very similar to the band structure in non-centrosymmetric spin-orbital materials, but materials that we propose to study are distinctly different, first and foremost in the sense that despite sharing many aspects of their electronic properties with the latter, they break the time-reversal symmetry without either spin-orbit coupling or lack of inversion symmetry. (2) Topologically nontrivial magnetic spirals. The PI has been engaged with the experimental group of Dr. Ghimire at GMU investigating Dirac materials with helical magnetism, based on stacked magnetic Kagome layers, with a generic formula of RMn6Sn6. In particular, Y Mn6Sn6 demonstrates a component of the Hall effect that is naturally interpreted in terms of a topological spin texture, as well as linear magnetoresistance. Our calculation identify Dirac states that are robust with respect to the spiral formation, and let us derive an advance mean-field model explaining the observed phase diagram. This model predicted four distinct phases, with very distinct properties, which have now been seen in neutron experiments. The same compound is known to demonstrate topological Hall effect in a particular magnetic phase, and only at elevated temperature. Based on our understanding of the phase diagram, we have worked out a phenomenological theory of a chiral (skyrmionic) response to an external magnetic field, similar to the nematic response to external strain in Fe-based superconductors, which is possible in a centrosymmetric lattice and without interplanar Dzyaloshinskii-Moriya interaction. This phenomenological theory agrees quantitatively with the experiment. It is in our plans to research other similar materials for this effect. (3) Search for 3D analogues of Fe-based superconductors. We want to investigate materials that can be viewed as 3D analogues of FeSe. Specifically, we want materials that are good metals and host antiferromagnetism, which can be suppressed by pressure and generate an s-wave superconductivity, as in Fe-based superconductors. We have in mind some candidates already. This work will proceed in close collaboration with the experimental group of Prof. Nirmal Ghimire in the same department, whose expertise lies in sample making, magnetometry and transport measurements of materials with complex magnetic structures."
A hybrid atomic-scale and continuum modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline (nc) Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. Here, the simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.
Here, we discuss the effect of a non-constant dislocation drag coefficient on the very high strain rate regime within an analytic model describing mobile-immobile dislocation intersections applicable to fcc polycrystals. Based on previous work on dislocation drag, we estimate its temperature and pressure dependence and its effects on stress–strain rate relations. In the high temperature regime, we show that drag can remain the dominating effect even down to intermediate strain rates. We also discuss the consequences of having a limiting dislocation velocity, a feature which is typically predicted by analytic models of dislocation drag, but which is somewhat under debate because a number of MD simulations predict supersonic dislocations.
Future NASA missions for earth observation and planetary science require large photovoltaic detector arrays with high performance in the long wavelength region to 18 microns and at operating temperatures above 65 K where single-cycle long-life cryocoolers are being developed. Since these detector array requirements exceed the state of current HgCdTe technology, alternative detector materials are being investigated as a possible option for future missions. Advanced growth techniques (e.g., MBE and MOCVD) of column III-V semiconductors have opened opportunities for engineering new detector materials and device structures. The technical approaches under investigation at JPL (with university and industry participation) include: quantum well infrared photodetectors, heterojunction internal photoemission (HIP) photodetectors, type-II strained layer superlattices, and nipi doping superlattices. Each of these options are briefly described with some of their pros and cons. A more detailed description is given for the HIP approach being pioneered at JPL.
Compressional (VP) and shear (VS) wave velocities of polycrystalline vanadium were measured simultaneously up to 11.5 GPa at room temperature using ultrasonic interferometry in a multi-anvil press. Complex softening behavior in VS and resulting shear moduli are discovered, possibly revealing a precursor to the reported phase transition within 30–60 GPa. The current data enables a comprehensive assessment of the elastic and mechanical properties of vanadium at high pressures, including bulk and shear moduli, Young’s modulus, Poisson’s ratio, and Pugh’s ratio. Through fitting to the 3rd-order finite strain equations, the elastic moduli and their pressure derivatives were determined to be K S0 = 151 (2) GPa, G 0 = 46.9 (8) GPa, K$^{’}_{S0}$ = 3.47 (5), and G$^{’}_{0}$ = 0.62 (1). These experimental results allow us to compare with and benchmark the existing Steinberg–Guinan models for extrapolations to extreme pressure and temperature conditions.
Abstract High-density phase change memory (PCM) storage is proposed for materials with multiple intermediate resistance states, which have been observed in 1 T -TaS 2 due to charge density wave (CDW) phase transitions. However, the metastability responsible for this behavior makes the presence of multistate switching unpredictable in TaS 2 devices. Here, we demonstrate the fabrication of nanothick verti-lateral H -TaS 2 /1 T -TaS 2 heterostructures in which the number of endotaxial metallic H -TaS 2 monolayers dictates the number of resistance transitions in 1 T -TaS 2 lamellae near room temperature. Further, we also observe optically active heterochirality in the CDW superlattice structure, which is modulated in concert with the resistivity steps, and we show how strain engineering can be used to nucleate these polytype conversions. This work positions the principle of endotaxial heterostructures as a promising conceptual framework for reliable, non-volatile, and multi-level switching of structure, chirality, and resistance.
Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.
It is well known that macroscopically normalizable zero-energy wave functions of spin-$\frac{1}{2}$ particles in a two-dimensional inhomogeneous magnetic field are spin-polarized and exactly calculable with degeneracy equaling the number of flux quanta linking the whole system. Here, extending this argument to massless Dirac fermions subjected to magnetic fields that have zero net flux but are doubly periodic in real space, we show that there exist only two Bloch-normalizable zero-energy eigenstates, one for each spin flavor. This result is immediately relevant to graphene multilayer systems subjected to doubly periodic strain fields, which at low energies enter the Hamiltonian as periodic pseudogauge vector potentials. Furthermore, we explore various related settings including nonlinearly dispersing band structure models and systems with singly periodic magnetic fields.
Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.
Homogeneous barotropic beta-plane turbulence is investigated, taking into account the effects of spatial inhomogeneity in the form of a zonal shear flows. Attention is given to the case of zonal flows that are barotropically stable and of larger scale than the resulting transient eddy field. Numerical simulations reveal that large-scale zonal flows alter the picture of classical beta-plane turbulence. It is found that the disturbance field penetrates to the largest scales of motion, that the larger disturbance scales show a tendency to meridional rather than zonal anisotropy, and that the initial spectral transfer rate away from an isotropic intermediate-scale source is enhanced by the shear-induced transfer associated with straining by the zonal flow.
Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).
Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.
When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.
In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.