Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dense”

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 19 records

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shear Stress Dependence of Force Networks in 3D Dense Suspensions

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicate a uniformly spreading network that rapidly forms across empty domains, however the growth patterns have structure that are significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Edens, Lance E.↗

Detection of random noise and anatomy of continuous seismic waveforms in dense array data near Anza California

We develop a methodology to separate continuous seismic waveforms into random noise (RN), not random noise (NRN) produced by earthquakes, wind, traffic and other sources of ground motions, and an undetermined mixture of signals. The analysis is applied to continuous records from a dense seismic array on the San Jacinto fault zone. To detect RN signals, we cut hourly waveforms into non-overlapping 1 s time windows and apply cross-correlations to separate RN candidates from outliers. The cross-correlation coefficients between different RN candidates fall into a tight range (i.e. 0.09–0.35), while cross-correlation coefficients of RN candidates with NRN signals (e.g. seismic or air-traffic events) are lower. The amplitude spectra of RN candidates have a well-defined level, while the amplitude spectra of other signals deviate from that level. Using these properties, we examine the amplitude spectra of moving time windows and cross-correlation coefficients with RN templates in each hour. The hourly RN is quasi-stationary and the results cluster tightly in the parameter space of cross-correlation coefficients and L2 norm deviations from the mean spectra of RN candidates. Time windows with parameters in this tight cluster are identified as RN, windows that deviate significantly from the RN cluster are identified as NRN and windows with values in between are identified as mixed signals. Several iterations on each hourly data are used to update and stabilize the selection of RN templates and mean noise spectra. For the days examined, the relative fractions of RN, NRN and mixed signals in local day (night) times are about 26 (42 percent), 40 (33 percent) and 34 percent (25 percent), respectively.

Geochemistry & Geophysics↗

Fast truncated SVD of sparse and dense matrices on graphics processors

We investigate the solution of low-rank matrix approximation problems using the truncated singular value decomposition (SVD). For this purpose, we develop and optimize graphics processing unit (GPU) implementations for the randomized SVD and a blocked variant of the Lanczos approach. Our work takes advantage of the fact that the two methods are composed of very similar linear algebra building blocks, which can be assembled using numerical kernels from existing high-performance linear algebra libraries. Furthermore, the experiments with several sparse matrices arising in representative real-world applications and synthetic dense test matrices reveal a performance advantage of the block Lanczos algorithm when targeting the same approximation accuracy.

Computer Science↗

Phase transformation of hydrous ringwoodite to the lower-mantle phases and the formation of dense hydrous silica

In order to understand the effects of H 2 O on the mineral phases forming under the pressure-temperature conditions of the lower mantle, we have conducted laser-heated diamond-anvil cell experiments on hydrous ringwoodite (Mg 2 SiO 4 with 1.1wt% H 2 O) at pressures between 29 and 59 GPa and temperatures between 1200 and 2400 K. Our experimental results show that hydrous ringwoodite (hRw) converts to crystalline dense hydrous silica, stishovite (Stv) or CaCl 2 -type SiO 2 (mStv), containing 1 wt% H 2 O together with Brd and MgO at the pressure-temperature conditions expected for shallow lower-mantle depths between approximately 660 to 1600 km. Considering lack of sign for melting in our experiments, our preferred interpretation of the observation is that Brd partially breaks down to dense hydrous silica and periclase (Pc), forming Brd + Pc + Stv mineralogy. Our experiments may provide an explanation for the enigmatic coexistence of Stv and Fp inclusions in lower-mantle diamonds

58 GEOSCIENCES↗

Enhanced Visibility of Subduction Slabs by the Formation of Dense Hydrous Phase A

Abstract Phase A (Mg 7 Si 2 O 8 (OH) 6 ) is one of the important dense hydrous magnesium silicates in subducting slab, because it forms after the breakdown of antigorite serpentine and could be the dominant hydrous phase in the upper‐mantle deep slab for the water transportation into deep Earth. In this study, the compressional ( P ) and shear ( S ) wave velocities of phase A were measured at simultaneous pressure and temperature conditions up to 11 GPa and 1073 K. Combined with elastic properties of olivine and pyroxenes, we calculate the hydration effect on the velocities of harzburgite lithology in cold subduction zones throughout the depth ranges where phase A is thermodynamically stable. Our calculations suggest that the hydration increases both P ‐ and S ‐wave velocities of harzburgite; at ∼5 wt% hydration, its seismic detectability is enhanced by 1%–1.5% in velocity contrasts relative to its anhydrous counterpart.

Cai, Nao↗

4D Proxy Imaging of Fracture Dilation and Stress Shadowing Using Electrical Resistivity Tomography During High Pressure Injections into a Dense Rock Formation

Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.

electrical resisitivity tomography, stress, 4D Tim↗

Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas

A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.

Solmaz, Evrim (ORCID:0000000211830881)↗

Kiloparsec-Scale Variations in the Star Formation Efficiency of Dense Gas: The Antennae Galaxies (NGC 4038/39)

We study the relationship between dense gas and star formation in the Antennae galaxies by comparing Atacama large millimeter/submillimeter array (ALMA) observations of dense-gas tracers (HCN, HCO{sup +}, and HNC J = 1–0) with the total infrared luminosity (L {sub TIR}) calculated using data from the Herschel Space Observatory and the Spitzer Space Telescope. We compare the luminosities of our star formation rate (SFR) and gas tracers using aperture photometry and employing two methods for defining apertures. We taper the ALMA data set to match the resolution of our L {sub TIR} maps and present new detections of dense-gas emission from complexes in the overlap and western arm regions. Using Owens Valley Radio Observatory CO J = 1–0 data, we compare with the total molecular gas content, M(H{sub 2}){sub tot}, and calculate star formation efficiencies and dense-gas mass fractions for these different regions. We derive HCN, HCO{sup +}, and HNC upper limits for apertures where emission was not significantly detected, because we expect that emission from dense gas should be present in most star-forming regions. The Antennae extends the linear L{sub TIR}--L{sub HCN} relationship found in previous studies. The L{sub TIR}--L{sub HCN} ratio varies by up to a factor of ∼10 across different regions of the Antennae, implying variations in the star formation efficiency of dense gas, with the nuclei, NGC 4038 and NGC 4039, showing the lowest SFE{sub dense} (0.44 and 0.70 × 10{sup −8} yr{sup −1}). The nuclei also exhibit the highest dense-gas fractions (∼9.1% and ∼7.9%).

79 ASTRONOMY AND ASTROPHYSICS↗