Search NASA⌕ Search

SEARCH · Search NASA

Results for “Statistical homogeneity”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Linearly Forced Isotropic Turbulence

Stationary isotropic turbulence is often studied numerically by adding a forcing term to the Navier-Stokes equation. This is usually done for the purpose of achieving higher Reynolds number and longer statistics than is possible for isotropic decaying turbulence. It is generally accepted that forcing the Navier-Stokes equation at low wave number does not influence the small scale statistics of the flow provided that there is wide separation between the largest and smallest scales. It will be shown, however, that the spectral width of the forcing has a noticeable effect on inertial range statistics. A case will be made here for using a broader form of forcing in order to compare computed isotropic stationary turbulence with (decaying) grid turbulence. It is shown that using a forcing function which is directly proportional to the velocity has physical meaning and gives results which are closer to both homogeneous and non-homogeneous turbulence. Section 1 presents a four part series of motivations for linear forcing. Section 2 puts linear forcing to a numerical test with a pseudospectral computation.

Lundgren, T. S.↗
NASA NTRS1985

Excitation of the Earth's Chandler wobble by southern oscillation/El Nino, 1900-1979

The southern oscillation/El Nino (ENSO) is the single most prominent interannual signal in global atmospheric/oceanic fluctuations. The following question is addressed: how important is the angular momentum carried by ENSO in exciting the Earth's Chandler wobble? The question is attacked through a statistical analysis of the coherence spectra (correlation as a function of frequency) between two data sets spanning 1900 to 1979-the southern oscillation index (SOI) time series and the excitation function psi (with x-component psi sub x and y-component psi sub y) of the Chandler wobble derived from the homogeneous ILS (International Latitude Service) polar motion data. The coherence power and phase in the Chandler frequency band (approx. 0.79 to 0.89 cpy) are studied. It is found that, during 1900 to 1979 the coherence between SOI and psi sub x is significant well over the 95% confidence threshold whereas that between SOI and psi sub y is practically nil. Quantitatively, the coherence study shows that ENSO provides some 20% of the observed Chandler wobble excitation power. Since earlier investigations have shown that the total atmospheric/oceanic variation can account for the Chandler wobble excitation at about 20% level, the implication is that ENSO maybe an important (interannual) part of the atmospheric/oceanic variation that is responsible for the Chandler wobble excitation during 1900 to 1979.

Chao, B. F.↗
NASA NTRS1988

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗
NASA NTRS2002

Broadband Noise Predictions Based on a New Aeroacoustic Formulation

A new analytic result in acoustics called 'Formulation 1B,' proposed by Farassat, is used to compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term. The formulation contains a far-field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is specified analytically from a result that is based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B, and to demonstrate its equivalence to Formulation 1A, of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. The predicted results also agree very well with those of Paterson and Amiet, who used a frequency-domain approach. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise.

Casper, J.↗
NASA NTRS2004

An Artificial Intelligence Classification Tool and Its Application to Gamma-Ray Bursts

Despite being the most energetic phenomenon in the known universe, the astrophysics of gamma-ray bursts (GRBs) has still proven difficult to understand. It has only been within the past five years that the GRB distance scale has been firmly established, on the basis of a few dozen bursts with x-ray, optical, and radio afterglows. The afterglows indicate source redshifts of z=1 to z=5, total energy outputs of roughly 10(exp 52) ergs, and energy confined to the far x-ray to near gamma-ray regime of the electromagnetic spectrum. The multi-wavelength afterglow observations have thus far provided more insight on the nature of the GRB mechanism than the GRB observations; far more papers have been written about the few observed gamma-ray burst afterglows in the past few years than about the thousands of detected gamma-ray bursts. One reason the GRB central engine is still so poorly understood is that GRBs have complex, overlapping characteristics that do not appear to be produced by one homogeneous process. At least two subclasses have been found on the basis of duration, spectral hardness, and fluence (time integrated flux); Class 1 bursts are softer, longer, and brighter than Class 2 bursts (with two second durations indicating a rough division). A third GRB subclass, overlapping the other two, has been identified using statistical clustering techniques; Class 3 bursts are intermediate between Class 1 and Class 2 bursts in brightness and duration, but are softer than Class 1 bursts. We are developing a tool to aid scientists in the study of GRB properties. In the process of developing this tool, we are building a large gamma-ray burst classification database. We are also scientifically analyzing some GRB data as we develop the tool. Tool development thus proceeds in tandem with the dataset for which it is being designed. The tool invokes a modified KDD (Knowledge Discovery in Databases) process, which is described as follows.

Hakkila, Jon↗
NASA NTRS2016

On the Importance of High Frequency Gravity Waves for Ice Nucleation in the Tropical Tropopause Layer

Recent investigations of the influence of atmospheric waves on ice nucleation in cirrus have identified a number of key processes and sensitivities: (1) ice concentrations produced by homogeneous freezing are strongly dependent on cooling rates, with gravity waves dominating upper tropospheric cooling rates; (2) rapid cooling driven by high-frequency waves are likely responsible for the rare occurrences of very high ice concentrations in cirrus; (3) sedimentation and entrainment tend to decrease ice concentrations as cirrus age; and (4) in some situations, changes in temperature tendency driven by high-frequency waves can quench ice nucleation events and limit ice concentrations. Here we use parcel-model simulations of ice nucleation driven by long-duration, constant-pressure balloon temperature time series, along with an extensive dataset of cold cirrus microphysical properties from the recent ATTREX high-altitude aircraft campaign, to statistically examine the importance of high-frequency waves as well as the consistency between our theoretical understanding of ice nucleation and observed ice concentrations. The parcel-model simulations indicate common occurrence of peak ice concentrations exceeding several hundred per liter. Sedimentation and entrainment would reduce ice concentrations as clouds age, but 1-D simulations using a wave parameterization (which underestimates rapid cooling events) still produce ice concentrations higher than indicated by observations. We find that quenching of nucleation events by high-frequency waves occurs infrequently and does not prevent occurrences of large ice concentrations in parcel simulations of homogeneous freezing. In fact, the high-frequency variability in the balloon temperature data is entirely responsible for production of these high ice concentrations in the simulations.

simulation↗
NASA NTRS1993

IUEAGN: A database of ultraviolet spectra of active galactic nuclei

In 13 years of operation, IUE has gathered approximately 5000 spectra of almost 600 Active Galactic Nuclei (AGN). In order to undertake AGN studies which require large amounts of data, we are consistently reducing this entire archive and creating a homogeneous, easy-to-use database. First, the spectra are extracted using the Optimal extraction algorithm. Continuum fluxes are then measured across predefined bands, and line fluxes are measured with a multi-component fit. These results, along with source information such as redshifts and positions, are placed in the IUEAGN relational database. Analysis algorithms, statistical tests, and plotting packages run within the structure, and this flexible database can accommodate future data when they are released. This archival approach has already been used to survey line and continuum variability in six bright Seyfert 1s and rapid continuum variability in 14 blazars. Among the results that could only be obtained using a large archival study is evidence that blazars show a positive correlation between degree of variability and apparent luminosity, while Seyfert 1s show an anti-correlation. This suggests that beaming dominates the ultraviolet properties for blazars, while thermal emission from an accretion disk dominates for Seyfert 1s. Our future plans include a survey of line ratios in Seyfert 1s, to be fitted with photoionization models to test the models and determine the range of temperatures, densities and ionization parameters. We will also include data from IRAS, Einstein, EXOSAT, and ground-based telescopes to measure multi-wavelength correlations and broadband spectral energy distributions.

Pike, G.↗
NASA NTRS1993

Application of the target decomposition theorem to a polarimetric random media model

With advances in polarimetric radar measurements of land surfaces, the need for understanding the underlying scattering mechanisms and dominant target features has become the focus of many studies. In particular, the maximum use of the polarimetric information to identify and/or separate parameters related to the surface features such as vegetation thickness, structure, water content, and soil surface characteristics will enhance the possibility of using polarimetric radars for monitoring the earth's surface from space. In this paper, Cloude's decomposition theorem is applied to a polarimetric random media model to simulate the radar measurements of vegetated canopies. The vegetated canopies are modeled as a three layer discrete random medium with leaves and branches in the first layer, tree trunks in the second layer, and a half space of homogeneous ground with rough interface as the bottom layer. The distorted born approximation (DBA) has been used to compute full Mueller matrix of the canopy, using canonical dielectric objects such as thin discs and cylinders as leaves, branches, and trunks, respectively. The Mueller matrix and the derived covariance matrix contain information on the second order statistics of radar signals at various polarizations from the canopy. To decompose the covariance matrix to its constituent targets, the eigenvalues and eigenvectors of the covariance matrix are computed in terms of the physical parameters of the canopy. In addition, each eigenvector explicitly shows the scattering mechanisms such as odd and even reflections in the canopy. Cloude's decomposition theorem is applied using the Pauli spin matrices as a basis and an expression for the degree of disorder or the entropy for the vegetated surface is found. Then, the physical parameters estimated from in situ measurements are used in the random media to obtain realistic covariance matrices. As a result, the sensitivity of the eigenvalue spectrums and the coefficients resulting from the target decomposition theorem to the physical parameters of the canopy are examined and the possible use of Cloude's theorem to estimate vegetation and soil parameters is discussed.

Saatchi, Sasan S.↗
NASA NTRS2021

Calibration of the DSCOVR EPIC Visible and NIR Channels using Multiple LEO Radiometers

The first five years of operation of the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) at the Lagrange one point have produced results that uniquely complement the data from currently operating low orbit Earth-observing instruments. In this paper we describe an updated unified approach to EPIC calibration. In this approach, calibration coefficients and their trends were obtained by comparing EPIC observations to the measurements from polar orbiting radiometers. In this study L1B reflectances from Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua and Terra satellites, Multi-angle Imaging Spectroradiometer (MISR) onboard Terra and Visible Infrared Imaging Radiometer (VIIRS) onboard the Suomi National Polar-orbiting Partnership (Suomi NPP) spacecraft were used to infer calibration coefficients for four EPIC visible and near-infrared channels: 443 nm, 551 nm, 680 nm, and 780 nm. EPIC Version three measurements made between June 2015 and August 2020 were used for comparison. The calibration procedure identifies the most homogeneous low Earth orbit radiometer scenes matching scattering angles that are temporarily and spatially collocated with EPIC observations. These scenes are used to determine reflectance to count (R/C) ratios in spectrally analogous channels. Seasonal average R/C ratios were analyzed to obtain EPIC calibration gains and trends. The trends for the full dataset period are not statistically significant except in the 443 nm channel. No significant changes in calibration were found after the instrument’s exit from safe hold in March 2020. The R/C ratios were also used to determine the differences in EPIC gains resulting from separate calibrations: against MODIS Aqua or Terra, as well as against forward or aftward MISR cameras. Statistical tests indicate that the differences between the two datasets are not significant except in the 780 nm channels where Aqua-derived coefficients may be around 2% lower compared to Terra. The dependence of EPIC calibration gains on the instrument scattering angle and on DSCOVR-Earth distance were investigated. Lastly, model Low Earth Orbit (LEO) reflectances calculated to match the EPIC viewing geometry were employed to study how EPIC calibration coefficients depend on EPIC-LEO viewing geometry differences. The effect of LEO and EPIC angular mismatch on calibration was shown to be small.

deep space climate observatory/earth polychromatic↗
DOE OSTI2021

Core–Multishell-Structured Digital-Gradient Cathode Materials with Enhanced Mechanical and Electrochemical Durability

Ni-rich cathode materials provide high energy density, but their structural and surface instability limits their cyclability and thermal stability. As one of the approaches to mitigate this problem, cathode materials comprising Ni-rich high-capacity core wrapped in Mn-rich multiple shells are produced successfully. In contrast to the conventional batch-type process for concentration-gradient materials, a digital-gradient cascade coprecipitation process described here achieves the improvements in productivity and quality consistency needed to move toward large-scale manufacturing. The core-multishell cathode materials produced in this manner not only have longer cycle life and improved rate performance compared to homogeneous Ni-rich cathode materials having the same overall composition, but also show remarkably enhanced thermal stability and low impedance growth characteristics. In a novel attempt to determine the correlation between the mechanical properties of the core-multishell cathode particles and their electrochemical cyclabilities, their breaking force and elasticity were successfully measured using a statistical approach, which indicates that a cathode particle with stable surface composition as well as high breaking force has improved capacity retention and durability. These results guide the realization of long life and high thermal stability in Ni-rich cathode materials through heterogeneous particle engineering.

25 ENERGY STORAGE↗
NASA NTRS1974

Thermal histories and crystal distributions in partly devitrified lunar glasses cooled by radiation

Calculations of the expected crystallization behavior of selected lunar compositions (60095 and 14259) are presented. Knowledge of this behavior combined with measurements of the state of crystallinity in partially devitrified lunar samples permits the determination of the thermal history of the samples. The state of crystallinity of a specimen is described by a statistical crystal distribution function psi. The heat flow problems of a plate of opaque glass quenched on an aluminum block and of a semi-transparent glass sphere cooling by radiation are analyzed. The results are combined with laboratory measurements of the crystal growth velocity and the viscosity, and with theoretical expressions for the nucleation frequency (homogeneous and heterogeneous) to obtain psi. It is found that the calculated psi is much too large in comparison with laboratory experiments on lunar samples. Possible reasons for the discrepancy are discussed, and the most likely cause is that the assumed nucleation barrier is too small. It also appears likely that these materials have relatively few nucleating heterogeneities.

Hopper, R. W.↗
NASA NTRS1990

Structure of turbulence at high shear rate

The structure of homogeneous turbulence subject to high shear rate has been investigated by using three-dimensional, time-dependent numerical simulations of the Navier-Stokes equations. This study indicates that high shear rate alone is sufficient for generation of the streaky structures, and that the presence of a solid boundary is not necessary. Evolution of the statistical correlations is examined to determine the effect of high shear rate on the development of anisotropy in turbulence. It is shown that the streamwise fluctuating motions are enhanced so profoundly that a highly anisotropic turbulence state with a 'one-component' velocity field and 'two-component' vorticity field develops asymptotically as total shear increases. Because of high-shear rate, rapid distortion theory predicts remarkably well the anisotropic behavior of the structural quantities.

Lee, Moon Joo↗
NASA NTRS2008

Gamma-Ray Burst Afterglows as Probes of Environment and Blast Wave Physics. II. The Distribution of rho and Structure of the Circumburst Medium

We constrain blast wave parameters and the circumburst media ofa subsample of 10 BeppoSAX gamma-ray bursts (GRBs). For this sample we derive the values of the injected electron energy distribution index, p, and the density structure index of the circumburst medium, k, from simultaneous spectral fits to their X-ray, optical, and NIR afterglow data. The spectral fits have been done in count space and include the effects ofmetallicity, and are compared with the previously reported optical and X-ray temporal behavior. Using the blast wave model and some assumptions which include on-axis viewing and standard jet structure, constant blast wave energy, and no evolution of the microphysical parameters, we find a mean value ofp for the sample as a whole of 9.... oa -0.003.0" 2 a_ statistical analysis of the distribution demonstrates that the p-values in this sample are inconsistent with a single universal value forp at the 3 _ level or greater, which has significant implications for particle acceleration models. This approach provides us with a measured distribution ofcircumburst density structures rather than considering only the cases of k ----0 (homogeneous) and k - 2 (windlike). We find five GRBs for which k can be well constrained, and in four of these cases the circumburst medium is clearly windlike. The fifth source has a value of 0 < k < 1, consistent with a homogeneous circumburst medium.

Starling, R. L. C.↗
NASA NTRS1981

Computing aerodynamic sound using advanced statistical turbulence theories

It is noted that the calculation of turbulence-generated aerodynamic sound requires knowledge of the spatial and temporal variation of Q sub ij (xi sub k, tau), the two-point, two-time turbulent velocity correlations. A technique is presented to obtain an approximate form of these correlations based on closure of the Reynolds stress equations by modeling of higher order terms. The governing equations for Q sub ij are first developed for a general flow. The case of homogeneous, stationary turbulence in a unidirectional constant shear mean flow is then assumed. The required closure form for Q sub ij is selected which is capable of qualitatively reproducing experimentally observed behavior. This form contains separation time dependent scale factors as parameters and depends explicitly on spatial separation. The approximate forms of Q sub ij are used in the differential equations and integral moments are taken over the spatial domain. The velocity correlations are used in the Lighthill theory of aerodynamic sound by assuming normal joint probability.

Hecht, A. M.↗
NASA NTRS2007

The Effects of Rainfall Inhomogeneity on Climate Variability of Rainfall Estimated from Passive Microwave Sensors

Passive microwave rainfall estimates that exploit the emission signal of raindrops in the atmosphere are sensitive to the inhomogeneity of rainfall within the satellite field of view (FOV). In particular, the concave nature of the brightness temperature (T(sub b)) versus rainfall relations at frequencies capable of detecting the blackbody emission of raindrops cause retrieval algorithms to systematically underestimate precipitation unless the rainfall is homogeneous within a radiometer FOV, or the inhomogeneity is accounted for explicitly. This problem has a long history in the passive microwave community and has been termed the beam-filling error. While not a true error, correcting for it requires a priori knowledge about the actual distribution of the rainfall within the satellite FOV, or at least a statistical representation of this inhomogeneity. This study first examines the magnitude of this beam-filling correction when slant-path radiative transfer calculations are used to account for the oblique incidence of current radiometers. Because of the horizontal averaging that occurs away from the nadir direction, the beam-filling error is found to be only a fraction of what has been reported previously in the literature based upon plane-parallel calculations. For a FOV representative of the 19-GHz radiometer channel (18 km X 28 km) aboard the Tropical Rainfall Measuring Mission (TRMM), the mean beam-filling correction computed in this study for tropical atmospheres is 1.26 instead of 1.52 computed from plane-parallel techniques. The slant-path solution is also less sensitive to finescale rainfall inhomogeneity and is, thus, able to make use of 4-km radar data from the TRMM Precipitation Radar (PR) in order to map regional and seasonal distributions of observed rainfall inhomogeneity in the Tropics. The data are examined to assess the expected errors introduced into climate rainfall records by unresolved changes in rainfall inhomogeneity. Results show that global mean monthly errors introduced by not explicitly accounting for rainfall inhomogeneity do not exceed 0.5% if the beam-filling error is allowed to be a function of rainfall rate and freezing level and does not exceed 2% if a universal beam-filling correction is applied that depends only upon the freezing level. Monthly regional errors can be significantly larger. Over the Indian Ocean, errors as large as 8% were found if the beam-filling correction is allowed to vary with rainfall rate and freezing level while errors of 15% were found if a universal correction is used.

Kummerow, Christian↗
NASA NTRS2017

Cancer Risk in Astronauts: A Constellation of Uncommon Consequences

Excess cancers resulting from external radiation exposures have been noted since the early 1950s, when a rise in leukemia rates was first reported in young atomic bomb survivors [1]. Further studies in atomic bomb survivors, cancer patients treated with radiotherapy, and nuclear power plant workers have confirmed that radiation exposure increases the risk of not only leukemia, but also a wide array of solid cancers [2,3]. NASA has long been aware of this risk and limits astronauts' risk of exposure-induced death (REID) from cancer by specifying permissible mission durations (PMD) for astronauts on an individual basis. While cancer is present among astronauts, current data does not suggest any excess of known radiation-induced cancers relative to a comparable population of U.S. adults; however, very uncommon cancers have been diagnosed in astronauts including nasopharyngeal cancer, lymphoma of the brain, and acral myxoinflammatory fibroblastic sarcoma. In order to study cancer risk in astronauts, a number of obstacles must be overcome. Firstly, several factors make the astronaut cohort considerably different from the cohorts that have previously been studied for effects resulting from radiation exposure. The high rate of accidents and the much healthier lifestyle of astronauts compared to the U.S. population make finding a suitable comparison population a problematic task. Space radiation differs substantially from terrestrial radiation exposures studied in the past; therefore, analyses of galactic cosmic radiation (GCR) in animal models must be conducted and correctly applied to the human experience. Secondly, a large enough population of exposed astronauts must exist in order to obtain the data necessary to see any potential statistically significant differences between the astronauts and the control population. Thirdly, confounders and effect modifiers, such as smoking, diet, and other space stressors, must be correctly identified and controlled for in those analyses. In order to begin work assessing the astronaut population, the earliest groups of astronauts (selection groups 1-7) provide a unique model. These astronauts were relatively homogenous, white males whose lifestyle characteristics were similar to an average U.S. citizen of the same birth cohort. This work reviews radiation exposure levels, age, and causes of mortality among these early NASA astronauts and discusses the benefits and limitations of assessing such a cohort for radiation-induced cancer risk.

Milder, Caitlin M.↗
DOE OSTI2022

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

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↗