Search NASA⌕ Search

SEARCH · Search NASA

Results for “relative phase angle”

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 217 records · Page 12

Electronic structure and charge-density wave transition in monolayer $\mathrm{VS}_2$

Vanadium disulfide (VS 2 ) attracts elevated interests for its charge-density wave (CDW) phase transition, ferromagnetism, and catalytic reactivity, but the electronic structure of monolayer has not been well understood yet. Here we report synthesis of epitaxial 1T VS 2 monolayer on bilayer graphene grown by molecular-beam epitaxy (MBE). Angle-resolved photoemission spectroscopy (ARPES) measurements reveal that Fermi surface with six elliptical pockets centered at the M points shows gap opening at low temperature. Temperature-dependence of the gap size suggests existence of CDW phase transition above room temperature. Here, our observations provide important evidence to understand the strongly correlated electron physics and the related surface catalytic properties in two-dimensional transition-metal dichalcogenides (TMDCs).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Io's sodium cloud - Explanation of the east-west asymmetries. II

A three- dimensional model for Io's sodium cloud incorporating the full effects of solar radiation acceleration associated with solar resonance scattering of atoms in the D-1 and D-2 lines is developed and applied to interpret two different but related east-west asymmetries exhibited in the cloud, using data obtained from ground-based telescopes. The physical changes produced in the cloud by the solar radiation acceleration as a function of satellite phase angle are documented, and how these changes explain the asymmetries is clarified. The model results suggest an east-west phase lag asymmetry of 10-15 degrees; results for the observed east-west intensity asymmetry of Bergstrahl (1977) are also presented. The east-west shape and intensity asymmetries are found to follow directly from the lack of circular symmetry about Jupiter of the forces acting on the cloud atoms.

Smyth, W. H.↗

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully non-linear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel-dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel-dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗

Lunar phase function effects on spectral ratios used for resource assessment

Groundbased telescopic CCD images of 36 selected locations on the moon were obtained in five 'standard' bandpasses at 12 phase angles ranging from -78 deg to +75 deg to measure phase function effects on the ratio values used to quantify the abundance of TiO2 and qualitatively indicate soil maturity. Consistent with previous studies, we find that the moon is 'bluer' at small phase angles, but that the effect on the ratio values for TiO2 abundance for the phase angles of our data is on the order of the measurement uncertainties throughout the range of abundances found in the mare. The effect is more significant as seen from orbiting spacecraft over a range of selenographic latitude. Spectral ratio images (400/560 and 400/730 nm) were used to map the abundance of TiO2 using the empirical relation found by Charlette et al from analysis of returned lunar soils. Additionally, the 950/560 and 950/730 nm image ratios were used to define the regions of mature mare soil in which the relation is valid. Although the phase function dependence on wavelength was investigated and quantified for small areas and the integrated disc, the effect specifically on TiO2 mapping was not rigorously determined. For consistency and convenience in observing the whole lunar front side, our mapping utilized images taken -15 deg less than alpha less than 15 deg when the moon was fully illuminated from earth; however, this includes the strong opposition peak.

Larson, S. M.↗

Phase Function Determination in Support of Orbital Debris Size Estimation

To recover the size of a space debris object from photometric measurements, it is necessary to determine its albedo and basic shape: if the albedo is known, the reflective area can be calculated; and if the shape is known, the shape and area taken together can be used to estimate a characteristic dimension. Albedo is typically determined by inferring the object s material type from filter photometry or spectroscopy and is not the subject of the present study. Object shape, on the other hand, can be revealed from a time-history of the object s brightness response. The most data-rich presentation is a continuous light-curve that records the object s brightness for an entire sensor pass, which could last for tens of minutes to several hours: from this one can see both short-term periodic behavior as well as brightness variations with phase angle. Light-curve interpretation, however, is more art than science and does not lend itself easily to automation; and the collection method, which requires single-object telescope dedication for long periods of time, is not well suited to debris survey conditions. So one is led to investigate how easily an object s brightness phase function, which can be constructed from the more survey-friendly point photometry, can be used to recover object shape. Such a recovery is usually attempted by comparing a phase-function curve constructed from an object s empirical brightness measurements to analytically-derived curves for basic shapes or shape combinations. There are two ways to accomplish this: a simple averaged brightness-versus phase curve assembled from the empirical data, or a more elaborate approach in which one is essentially calculating a brightness PDF for each phase angle bin (a technique explored in unpublished AFRL/RV research and in Ojakangas 2011); in each case the empirical curve is compared to analytical results for shapes of interest. The latter technique promises more discrimination power but requires more data; the former can be assembled in its essentials from fewer measurements but will be less definitive in its assignments. The goal of the present study is to evaluate both techniques under debris survey conditions to determine their relative performance and, additionally, to learn precisely how a survey should be conducted in order to maximize their performance. Because the distendedness of objects has more of an effect than their precise shape in calculating a characteristic dimension, one is interested in the techniques discrimination ability to distinguish between an elongated rectangular prism and a short rectangular prism or cube, or an elongated cylinder from a squat cylinder or sphere. Sensitivity studies using simulated data will be conducted to determine discrimination power for both techniques as a function of amount of data collected and range (and specific region) of phase angles sampled. Empirical GEODSS photometry data for distended objects (dead payloads with solar panels, rocket bodies) and compact objects (cubesats, calibration spheres, squat payloads) will also be used to test this discrimination ability. The result will be a recommended technique and data collection paradigm for debris surveys in order to maximize this type of discrimination.

Hejduk, M. D.↗

Trajectory Design of the Lunar Impactor Mission Concept

The National Aeronautics and Space Administration (NASA) solicited proposals in 2006 for an opportunity to include a small secondary payload with the launch of the Lunar Reconnaissance Orbiter (LRO) scheduled for October 2008. The cost cap of the proposal was between $50 and $80M, and the mass cap was 1,000 kilograms. JPL proposed a Lunar Impactor (LI) concept for this solicitation. The mission objective of LI was to impact the permanently shadowed region of a South polar crater ultimately to detect the presence of water. The detection of water ice would prove to be an important factor on future lunar exploration. NASA Ames Research Center also proposed a similar concept, the Lunar Crater observation and Sensing Satellite (LCROSS), which was selected by NASA for the mission. However, in this paper, the trajectory design of the LI proposed by JPL is considered. Since the LI spacecraft was to be launched on the LRO launch vehicle as a secondary payload, its initial trajectory must be diverted at some later time from the LRO trans-lunar trajectory for the subsequent impact. Several such trajectories have been considered, where each trajectory option fields some specific values for the mission parameters. The mission parameters include the availability of LRO instruments at the time of impact for the observation by LRO, the mission duration, the impact velocity, the impact angle, etc. It is possible for the LI to be deflected with a relatively low delta-V to impact a South polar crater at a reasonable impact velocity and impact angle directly with no delay. However, the instruments on-board LRO may not be ready for observation. Thus, several delayed trajectory options have been considered further. The lunar phase at the time of impact may also play an important factor for observation, especially from Earth. Several lunar flyby trajectory maneuvers have been identified to arrive at the Moon for impact at the desired lunar phase. By using a combination of these successive lunar flyby maneuvers, the impact lunar phase may be adjusted to the desired location. A few such trajectories have been suggested. Also, some attempts have been made to maximize the impact velocity by converting the impact trajectory into a retrograde orbit with respect to Earth. Since these types of trajectories take advantage of the Sun-Earth three-body region to minimize the delta-V, the mission duration is relatively long. A few such trajectories are suggested. Also, an attempt has been made to adjust the lunar impact within a desired time period for the optimum Earth observation for the above trajectories. The mission parameters resulting from each trajectory option above are considered and weighed against the cost and robustness of the mission in a brief summary.

trajectory design↗

A New Machine Learning Based Analysis for Improving Satellite Retrieved Atmospheric Composition Data: OMI SO2 as an Example

Despite recent progress, satellite retrievals of anthropogenic SO2 still suffer from relatively low signal-tonoise ratios. In this study, we demonstrate a new machine learning data analysis method to improve the quality of satellite SO2 products. In the absence of large ground-truth datasets for SO2, we start from SO2 slant column densities (SCDs) retrieved from the Ozone Monitoring Instrument (OMI) using a data-driven, physically based algorithm and calculate the ratio between the SCD and the root mean square (rms) of the fitting residuals for each pixel. To build the training data, we select presumably clean pixels with small SCD / rms ratios (SRRs) and set their target SCDs to zero. For polluted pixels with relatively large SRRs, we set the target to the original retrieved SCDs. We then train neural networks (NNs) to reproduce the target SCDs using predictors including SRRs for individual pixels, solar zenith, viewing zenith and phase angles, scene reflectivity, and O3 column amounts, as well as the monthly mean SRRs. For data analysis, we employ two NNs: (1) one trained daily to produce analyzed SO2 SCDs for polluted pixels each day and (2) the other trained once every month to produce analyzed SCDs for less polluted pixels for the entire month. Test results for 2005 show that our method can significantly reduce noise and artifacts over background regions. Over polluted areas, the monthly mean NN-analyzed and original SCDs generally agree to within ±15 %, indicating that our method can retain SO2 signals in the original retrievals except for large volcanic eruptions. This is further confirmed by running both the NN-analyzed and original SCDs through a topdown emission algorithm to estimate the annual SO2 emissions for ∼ 500 anthropogenic sources, with the two datasets yielding similar results. We also explore two alternative approaches to the NN-based analysis method. In one, we employ a simple linear interpolation model to analyze the original SCD retrievals. In the other, we develop a PCA–NN algorithm that uses OMI measured radiances, transformed and dimension-reduced with a principal component analysis (PCA) technique, as inputs to NNs for SO2 SCD retrievals. While the linear model and the PCA–NN algorithm can reduce retrieval noise, they both underestimate SO2 over polluted areas. Overall, the results presented here demonstrate that our new data analysis method can significantly improve the quality of existing OMI SO2 retrievals. The method can potentially be adapted for other sensors and/or species and enhance the value of satellite data in air quality research and applications.

Can Li↗

Acoustic Phased Array Quantification of Quiet Technology Demonstrator 3 Advanced Inlet Liner Noise Component

Acoustic phased array flyover noise measurements were acquired as part of the Boeing 737 MAX-7 NASA Advanced Inlet Liner segment of the Quiet Technology Demonstrator 3 (QTD3) flight test program. This paper reports on the processes used for separating and quantifying the engine inlet, exhaust and airframe noise source components and provides sample phased array-based comparisons of the component noise source levels associated with the inlet liner treatment configurations. Full scale flyover noise testing of NASA advanced inlet liners was conducted as part of the Quiet Technology Demonstrator 3 flight test program in July and August of 2018. Details on the inlet designs and testing are provided in the companion paper of Reference 1. The present paper provides supplemental details relating to the acoustic phased array portion of the analyses provided in Ref. 1. In brief, the test article was a Boeing 737MAX-7 aircraft with a modified right hand (starboard side) engine inlet, which consisted of either a production inlet liner, a NASA designed inlet liner or a simulated hard wall configuration (accomplished by applying speed tape over the inlet acoustic treatment areas). In all three configurations, the engine forward fan case acoustic panel was replaced with a unperforated (hardwall) panel. No other modifications to any other acoustic treatment areas were made. The left hand (port side) engine was a production engine and was flown at idle thrust for all measurements in order to isolate the effects of the inlet liners to the right hand engine. As described in Ref. 1, the NASA inlet treatment consists of laterally cut slots (cut perpendicular to the flow direction) which are designed to reduce excrescence drag while maintaining or exceeding the liner acoustic noise reduction capabilities. The NASA inlet liner consists of a Multi-Degree of Freedom (MDOF) design with two breathable septum layers inserted into each honeycomb cell [1]. The aircraft noise measurements were acquired for both takeoff (flaps 1 setting, gear up) and approach (flaps 30 gear up and gear down) configurations. The inlet and flight test configurations are summarized in Table 1. Table 1: Inlet Treatment and Flight Configurations Inlet Forward Fan Case Aircraft Production Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down NASA Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down Hardwall Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down III.Test Description and Hardware The flight testing was conducted at the Grant County airport in Moses Lake, WA, between 27 July and 6 August 2018. The noise measurement instrumentation included 8 flush dish microphones arranged in a noise certification configuration as well as an 840 microphone phased array. The flush dish microphones were used to quantify the levels and differences in levels between the various inlet treatments. The phased array was used to separate and quantify the narrowband (tonal) and broadband noise component levels from the engine inlet/exhaust and from the airframe. Phased array extraction of the broadband component was critical to this study because it allowed for the separation of the inlet component from the total airplane level noise even when it was significantly below the total level. Figure 1 provides an overview of the phased array microphone layout as well as a detailed image of an individual phased array microphone mounted in a plate holder (the microphone sensor is the dot in the center of the plate). The ground plane ensemble array microphones (referred to as “ensemble array” in this paper) were mounted in plates with “flower petal” edges designed to minimize edge scattering effects. Fig. 1 Flyover test microphone layout. The phased array configuration was the result of a progressive development of concepts originally implemented in Ref. 2 and refined over the following years, consisting namely of multiple multi-arm logarithmic spiral subarrays designed to cover overlapping frequency ranges and optimized for various aircraft emission angles. For the present case, the signals from all 840 microphones were acquired on a single system. The 840 microphones were parsed into 11 primary subarray sets spanning from smallest to largest aperture size and labeled accordingly as “a, b, …, k”, where “a” corresponds to the smallest fielded subarray and “k” corresponds to the largest aperture subarray. The apertures ranged from approximately 10 ft to 427 ft in size (in the flight direction) with the subarrays consisting of between 215 and 312 microphones. Figure 2 shows three such subarrays, k, h and a. As done in Ref. 2, microphones were shared between subarrays in order to reduce total channel count. Fig. 2 Sample subarray sizes (20° from overhead – refer to Figure 3a discussion). In addition to the above, each of the 11 primary subarray sets consisted of four subarrays optimized to provide near equivalent array spatial resolution in both the flight and lateral directions within 30 degrees of overhead (i.e., airplane directly above the center of the array), namely, at angles  of 0, ±10, ±20 and ±30 degrees relative to overhead where angle  is defined as shown in Figure 3a. This allowed for optimized aircraft noise measurements from 60 to 120 degree emission angle.6 An example of this pletharray design is shown in Figure 3b for the k subarray. When the aircraft is at overhead, the microphones indicated by the blue markers are used for beamforming. When the aircraft is at angles ±10 degrees from overhead, both the blue and red colored microphones are used, and so on for the ±20 and ±30 degree aircraft locations. See Ref. 3 for extensive details on pletharray design for aeroacoustic phased array testing. 6 In the discussions that follow, emission angle values are used. These are the angles at the time sound is emitted relative to the engine axis and are calculated based on flight path angle, body aircraft body angle with respect to the relative wind direction, and engine axis angle relative to aircraft body angle.

Brusniak, Leon↗

Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors

The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.

36 MATERIALS SCIENCE↗

Discrete ordinate theory of radiative transfer. 2: Scattering from maritime haze

Discrete ordinate theory was used to calculate the reflected and transmitted radiance of photons which have interacted with plane parallel maritime haze layers. The results are presented for three solar zenith angles, three values of the surface albedo, and a range of optical thicknesses from very thin to very thick. The diffuse flux at the lower boundary and the cloud albedo were tabulated. The forward peak and other features in the single scattered phase function caused the radiance in many cases to be very different from that for Rayleigh scattering. The variation of the radiance with both the zenith or nadir angle and the azimuthal angle is more marked, and the relative limb darkening under very thick layers is greater, for haze than for Rayleigh scattering. The downward diffuse flux at the lower boundary for A = O is always greater and the cloud albedo is always less for haze than for Rayleigh layers.

Kattawar, G. W.↗

Earthshine as an Illumination Source at the Moon

Earthshine is the dominant source of natural illumination on the surface of the Moon during lunar night, and at some locations within permanently shadowed regions (PSRs) near the poles that never receive direct sunlight. As such, earthshine has the potential to enable the scientific investigation and exploration of conditions in areas of the Moon that are either temporarily or permanently hidden from the Sun. Earthshine has also been used to refer to Earthlight reflected from the lunar surface, but in this study we use it to refer specifically to Earthlight incident at the Moon. Under certain circumstances, the heat flux from earthshine could also influence the transport and cold-trapping of volatiles present in the very coldest areas within PSRs. In this study, Earth's spectral irradiance, as it would appear at the Moon in the solar reflectance band (0.3–3.0 µm) and at thermal emission wavelengths (3–50 µm), is examined with a suite of model image cubes and whole-disk spectra created using the Virtual Planetary Laboratory (VPL) three-dimensional (latitude, longitude and altitude) modeling capability. At the Moon, the broadband, hemispherical irradiance from Earth at full-phase is approximately 0.15 W m−2 with comparable contributions from solar reflectance and thermal emission; for context, this about 0.01% that of solar irradiance and has an equivalent temperature of around 40 K. Over the simulated timeframe, spanning two lunations, Earth's thermal irradiance shows very little net change (less than a few mW m−2 resulting from cloud variability and the south-to-north motion of the sub-observer latitude on Earth). In the solar band, Earth's diurnally averaged light curve at phase angles g ≤ 60° is well-fit using a Henyey–Greenstein integral phase function. At wavelengths longward of about 0.7 µm, near the well-known vegetation “red edge”, Earth's reflected solar radiance shows significant diurnal modulation as a result of the broad maximum in projected landmass at terrestrial longitudes between 60°W and 0°, as well as from the distribution of clouds. A simple formulation with adjustable coefficients is presented, condensed from the VPL model grid, for estimating Earth's hemispherical irradiance at the Moon as a function of wavelength, phase angle and sub-observer coordinates (terrestrial latitude and longitude). Uncertainties in any one prediction are estimated to be 10–12% at 0.3 µm, rising to >25% near 2.5 µm as a result of the increasing relative brightness and unpredictable influence of clouds. Although coefficient values are derived from a suite of spring equinox models, the approximation appears to be valid for all seasons, to within the stated uncertainties. It is demonstrated that earthshine is sufficiently bright to serve as a natural illumination source for optical measurements on a robotic lander/rover, allowing the identification of water ice mixed with regolith at the percent-level of fractional area.

Surface↗

Effect of gravity on liquid phase sintering

The unique conditions offered by the earth orbital environment in material processing involving both solid and liquid phases, such as liquid phase sintering, were studied. An experimental development program involving both test and theoretical work was initiated. Experimental work using material combinations selected such that maximum information about the effect of gravity can be derived has been conducted. Wetting of the solid phase by the liquid during sintering is an important phenomenon in liquid phase sintering, and gravity has influence on both capillary phenomenon and density segregation; hence, material combinations were selected such that these two effects can be suitably studies. The experimental work is meant to form the basis for similar comparative work done under low-g conditions. The part of the model dealing with the capillary phenomenon, as related to liquid phase sintering and the effect of gravity on it, suggest that gravity will have negligible effect on the Bond number and that the cohesive force is dependent on both the amount of liquid phase and the angle of contact.

Mookherji, T.↗

Phase-transition-induced thermal hysteresis in Type-II weyl semimetals M o T e 2 and M o 1 − x W x T e 2

The resistivity versus temperature measurement is commonly used for identifying temperature-induced phase-change and the resulting hysteresis loop. However, while the resistance is influenced by both the density of states and the carrier lifetimes, the Seebeck coefficient is influenced predominantly by the density of states, and hence is a better probe of the phase of the material. Here, 1T'-T d temperature-induced phase transition in MoTe 2 is studied using temperature-dependent X-ray diffraction, resistivity, and Seebeck coefficient measurements. A more distinct hysteresis is observed when measuring the Seebeck coefficient which is consistent with direct measurements of the crystallographic angle using the temperature-dependent X-ray diffraction. The Seebeck and electrical resistivity measurements indicate a competing contribution of the electrons and holes. The contribution of electron pockets becomes more dominant when molybdenum atoms are replaced by tungsten. In MoTe 2 , a topologically induced enhancement of the Nernst coefficient is observed at low temperatures and a relatively large phase-transition induced Thomson coefficient of 111 μV∙K -1 is measured at 254 K.

36 MATERIALS SCIENCE↗

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗

Strange Metal in Magic-Angle Graphene with near Planckian Dissipation

Recent experiments on magic-angle twisted bilayer graphene have discovered correlated insulating behavior and superconductivity at a fractional filling of an isolated narrow band. Here we show that magic-angle bilayer graphene exhibits another hallmark of strongly correlated systems—a broad regime of $\textit{T{$-linear resistivity above a small density-dependent crossover temperature—for a range of fillings near the correlated insulator. This behavior is reminiscent of similar behavior in other strongly correlated systems, often denoted “strange metals,” such as cuprates, iron pnictides, ruthenates, and cobaltates, where the observations are at odds with expectations in a weakly interacting Fermi liquid. We also extract a transport “scattering rate,” which satisfies a near Planckian form that is universally related to the ratio of $(k_{\text{B}}T/ℏ)$. In this work, our results establish magic-angle bilayer graphene as a highly tunable platform to investigate strange metal behavior, which could shed light on this mysterious ubiquitous phase of correlated matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Matrix operator theory of radiative transfer. II - Scattering from maritime haze.

Matrix operator theory is used to calculate the reflected and transmitted radiance of photons that have interacted with plane-parallel maritime haze layers. The results are presented for three solar zenith angles, three values of the surface albedo, and a range of optical thicknesses from very thin to very thick. The diffuse flux at the lower boundary and the cloud albedo are tabulated. The forward peak and other features in the single-scattered phase function cause the radiance in many cases to be very different from that for Rayleigh scattering. In particular, the variation of the radiance with both the zenith or nadir angle and the azimuthal angle is more marked and the relative limb darkening under very thick layers is greater for haze M than for Rayleigh scattering. The downward diffuse flux at the lower boundary for A = 0 is always greater and the cloud albedo is always less for haze M than for Rayleigh layers.

Kattawar, G. W.↗

High-pressure infrared sepctra of alpha-quartz, coesite, stishovite and silica glass

High-pressure infrared absorption spectra of alpha-quatz, coesite, stishovite, and SiO2 glass are consistent with the primary compression mechanism of the initially tetrahedrally bonded phases being the bending of the Si-O-Si angle at pressures less than 10-20 GPa. At higher pressures, up to 40 GPa, we observe a decline in the intensity of the infrared SiO4 asymmetric-stretching vibrations of all three phases, with an increase in the relative amplitude between 700 and 900/cm. This change in intensities is attributed to an increase in the average coordination number of silicon through extreme distortion of tetrahedra. At pressures above approximately 20 GPa, the low-pressure crystalline polymorphs gradually become amorphous, and the infrared spectra provide evidence for an increase in silicon coordination in these high-density amorphous phases. The pressure-amorphized samples prepared from quartz and coesite differ structurally both from each other and from silica glass that has been compressed, and the high pressure spectra indicate that these materials are considerably more disordered than stishovite under comparable pressure conditions. Average mode Grueneisen parameters calculated for quartz, stishovite and fused silica from both infrared and Raman spectra are compatible with the corresponding thermodynamic value of the Grueneisen parameter, however, that of coesite is significantly discrepant.

Williams, Q.↗