SEARCH · Engineering Papers

Results for “standard normal variate”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Amplitude calibration of spaceborne synthetic aperture radars

Problems encountered during attempts to calibrate SAR imagery, recent successful experiments conducted with SEASAT SAR data, and a proposed program for the calibration and validation of the radar imagery from the forthcoming SIR-B SAR are discussed. The SEASAT SAR data for 10 passes over Death Valley, California, were processed with a modified digital correlator. The procedure included a preliminary screening of the data to check for raw data saturation, compensation of waveforms and estimation of the amplitude of the pilot tone. All data was normalized to this pilot tone signal to reduce the effects of variable gains in the data links and ground receivers. The digital correlation algorithm generated image data. Evaluation of 6 passes results in a maximum pass to pass gain variation of only 1.1 dB and a standard deviation amongst the passes of 0.35 dB. previously announced in STAR as N83-26215

Held, D. N.↗

DOE OSTI2021

Clutter Distributions for Tomographic Image Standardization in Ground-Penetrating Radar

Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce 3-D distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude and vary across different arrays, which makes a direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or log-normal distributions. Based on the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter- to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Lastly, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development that needs to be robust across multiple different arrays.

42 ENGINEERING↗

NASA NTRS1991

Signal processing schemes for Doppler global velocimetry

Two schemes for processing signals obtained from the Doppler global velocimeter are described. The analog approach is a simple, real-time method for obtaining an RS-170 video signal containing the normalized intensity image. Pseudocolors are added using a monochromatic frame grabber producing a standard NTSC video signal that can be monitored and/or recorded. The digital approach is more complicated, but maintains the full resolution of the acquisition cameras with the ability to correct the signal image for pixel sensitivity variations and to remove background light. Prototype circuits for each scheme are described, and example results from the investigation of the vortical flow field above a 75-deg delta wing are presented.

Meyers, James F.↗

NASA NTRS1991

Signal Processing Schemes for Doppler Global Velocimetry

Two schemes for processing signals obtained from the Doppler global velocimeter are described. The analog approach is a simple, real time method for obtaining an RS-170 video signal containing the normalized intensity image. Pseudo colors are added using a monochromatic frame grabber producing a standard NTSC video signal that can be monitored and/or recorded. The digital approach is more complicated, but maintains the full resolution of the acquisition cameras with the capabilities to correct the signal image for pixel sensitivity variations and to remove of background light. Prototype circuits for each scheme are described and example results from the investigation of the vortical flow field above a 75-degree delta wing presented.

Meyers, James F.↗

DOE OSTI2020

Performance of Photovoltaic Systems Recorded by Open Solar Performance and Reliability Clearinghouse (oSPARC)

This article analyzes data from 2,200 photovoltaic systems collected through the Open Solar Performance and Reliability Clearinghouse (oSPARC, 2018) to draw conclusions about the measured performance of the systems, how it compares to expectations, and how wide the variation is. The systems range in size from 2 kW to 1,200 kW, with an average size of 229 kW, and were activated in oSPARC between March 13, 2007 and Dec. 16, 2016. Data regarding co-incident insolation and ambient temperature are provided through the oSPARC platform for the population of systems and summarized here. Results indicate that - after correcting for measured insolation, derating for temperature, and correcting for balance-of-system efficiency and degradation as a result of system age - the performance ratio of these 2,200 systems averages 91.7%. Adjusting for known parameters such as age and efficiency, the ratio would ideally be 100%. Thus, the 'underperformance' of this system sample is measured at 9.3% and suggests that as much as 9.3% could be recovered if the systems were caused to perform as expected through optimal operation and maintenance. Efforts to correlate the performance ratio with environmental conditions resulted in very low coefficients of regression, but the method already normalizes for actual site insolation and temperature, which are responsible for much of the site-specific variation. An exemption is air-quality nonattainment areas, where the measured performance of plants averages 86.9% and is consistently lower than the average of 95.9% in areas that are in compliance with air quality standards, likely because of particulate matter such as diesel soot exacerbating the soiling rate.

14 SOLAR ENERGY↗

NASA NTRS1977

Veterinary medical considerations for the use of nonhuman primates in space research

The validity of biomedical research using animal subjects is highly dependent on the use of 'normal' and healthy animals. The current costs of research programs dictate that a minimum number of animals and test replicates be used to obtain the desired data. The use of healthy and standardized animals increases the probability of obtaining valid data while also permitting greater economy by reducing the between-individual variation, thus allowing the use of fewer animals. Areas of concern when planning animal payloads include constraints of the flight on candidate species selection, screening for physiological and psychological normalcy, procedures for routine care and quarantine of new animals and those returning from space, ground-based studies to determine experimental protocol, selection of instrumentation, stress during transportation for flight operations, housing and care facilities at launch and recovery sites, and the overall veterinary program.

Simmonds, R. C.↗

NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets such as sex or age of the model organism used. In the present study, NASA GeneLab-hosted RNAseq datasets from rodent liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC, to determine statistical differences between datasets before and after correction, Principal Component Analysis, to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the standard approach. Thus, the most robust standard correction will be implemented in the GeneLab Visualization 2.0 platform when datasets are combined.

GeneLab, RNA-seq, Batch Correction↗

NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the respective standard approach. Of the methods tested, standard ComBat and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

GeneLab↗

NASA NTRS2000

Evaluation of Non-Ozone-Depleting-Chemical Cleaning Methods for Space Mechanisms Using a Vacuum Spiral Orbit Rolling Contact Tribometer

Because CFC 113, an ozone depleting chemical (ODC), can no longer be produced, alternative bearing cleaning methods must be studied. The objective of this work was to study the effect of the new cleaning methods on lubricant lifetime using a vacuum bearing simulator (spiral orbit rolling contact tribometer). Four alternative cleaning methods were studied: ultra-violet (UV) ozone, aqueous levigated alumina slurry (ALAS), super critical fluid (SCF) CO2 and aqueous Brulin 815GD. Baseline tests were done using CFC 113. Test conditions were the following: a vacuum of at least 1.3 x 10(exp -6) Pa, 440C steel components, a rotational speed of 10 RPM, a lubricant charge of between 60-75 micrograms, a perfluoropolyalkylether lubricant (Z-25), and a load of 200N (44.6 lbs., a mean Hertzian stress of 1.5 GPa). Normalized lubricant lifetime was determined by dividing the total number of ball orbits by the amount of lubricant. The failure condition was a friction coefficient of 0.38. Post-test XPS analysis was also performed, showing slight variations in post-cleaning surface chemistry. Statistical analysis of the resultant data was conducted and it was determined that the data sets were most directly comparable when subjected to a natural log transformation. The natural log life (NL-Life) data for each cleaning method were reasonably normally (statistically) distributed and yielded standard deviations that were not significantly different among the five cleaning methods investigated. This made comparison of their NL-Life means very straightforward using a Bonferroni multiple comparison of means procedure. This procedure showed that the ALAS, UV-ozone and CFC 113 methods were not statistically significantly different from one another with respect to mean NL-Life. It also found that the SCF CO2 method yielded a significantly higher mean NL-Life than the mean NL-Lives of the ALAS, UV-ozone and CFC 113 methods. It also determined that the aqueous Brulin 815GD method yielded a mean NL-Life that was statistically significantly higher than the mean NL-Lives of each of the other four methods. Baseline tests using CFC 113 cleaned parts yielded a mean NL-Life 3.62 orbits/micro-g. ALAS and UV-ozone yielded similar mean NL-Life (3.31 orbits/mg and 3.33 orbits/micro-g, respectively). SCF CO2, gave a mean NL-Life of 4.08 orbits/mg and aqueous Brulin 8l5GD data yielded the longest mean NL-Life (4.66 orbits/micro-g).

Jansen, Mark J.↗

NASA NTRS2012

Discontinuous Mode Power Supply

A document discusses the changes made to a standard push-pull inverter circuit to avoid saturation effects in the main inverter power supply. Typically, in a standard push-pull arrangement, the unsymmetrical primary excitation causes variations in the volt second integral of each half of the excitation cycle that could lead to the establishment of DC flux density in the magnetic core, which could eventually cause saturation of the main inverter transformer. The relocation of the filter reactor normally placed across the output of the power supply solves this problem. The filter reactor was placed in series with the primary circuit of the main inverter transformer, and is presented as impedance against the sudden changes on the input current. The reactor averaged the input current in the primary circuit, avoiding saturation of the main inverter transformer. Since the implementation of the described change, the above problem has not reoccurred, and failures in the main power transistors have been avoided.

Lagadinos, John↗

NASA NTRSDate not supplied

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗

DOE OSTI2025

Estimation of intensity, footprint, and capacity of surface urban heat islands using a direction-enhanced adaptive synchronous extraction (DEASE) method

Here, the surface urban heat island (SUHI) effect, assessed through remotely sensed land surface temperature (LST), remains a focal point in urban climate research. Conventional indicators like SUHI intensity (SUHII) and footprint (SUHIF) capture peak values and spatial extent but fail to account for the cumulative thermal load—a critical dimension reflecting the total heat exposure imposed by spatially continuous warming, which directly limits a holistic assessment of ecological and societal impacts of the SUHI effect. Therefore, this study introduces an indicator termed SUHI capacity (SUHIC), designed to quantify the aggregated SUHI effect by integrating the magnitude of the warming signal across all affected areas, thereby enabling a more comprehensive evaluation of urban thermal environments. Furthermore, a direction-enhanced adaptive synchronous extraction (DEASE) method is proposed for the quantification of SUHIC. This method can dynamically identify the optimal background reference area based on the urban-rural LST gradients in various directions within the city, without relying on predefined mathematical models as previously. The results from 102 European cities first confirm that the directional variations in urban-rural LST gradients, and the DEASE method can effectively capture these distinctions for the simultaneous estimation of SUHII, SUHIF, and SUHIC. Secondly, the spatial patterns of absolute SUHIC values show strong associations with those of SUHIF (R2>0.86), while its relative values (normalized by the area of urban) align more closely with SUHII (R2 > 0.64). More importantly, SUHIC can serve as a crucial reference for assessing the urban thermal signal when SUHII and SUHIF diverge. The proposed method and framework contribute to standardizing the quantification of the SUHI effect.

Indicator↗

NASA NTRS1950

Flight Calibration of four airspeed systems on a swept-wing airplane at Mach numbers up to 1.04 by the NACA radar-phototheodolite method

The calibrations of four airspeed systems installed in a North American F-86A airplane have been determined in flight at Mach numbers up to 1.04 by the NACA radar-phototheodolite method. The variation of the static-pressure error per unit indicated impact pressure is presented for three systems typical of those currently in use in flight research, a nose boom and two different wing-tip booms, and for the standard service system installed in the airplane. A limited amount of information on the effect of airplane normal-force coefficient on the static-pressure error is included. The results are compared with available theory and with results from wind-tunnel tests of the airspeed heads alone. Of the systems investigated, a nose-boom installation was found to be most suitable for research use at transonic and low supersonic speeds because it provided the greatest sensitivity of the indicated Mach number to a unit change in true Mach number at very high subsonic speeds, and because it was least sensitive to changes in airplane normal-force coefficient. The static-pressure error of the nose-boom system was small and constant above a Mach number of 1.03 after passage of the fuselage bow shock wave over the airspeed head.

Thompson, Jim Rogers↗

NASA NTRS2001

NASA/American Cancer Society High-Resolution Flow Cytometry Project-I

BACKGROUND: The NASA/American Cancer Society (ACS) flow cytometer can simultaneously analyze the electronic nuclear volume (ENV) and DNA content of cells. This study describes the schematics, resolution, reproducibility, and sensitivity of biological standards analyzed on this unit. METHODS: Calibrated beads and biological standards (lymphocytes, trout erythrocytes [TRBC], calf thymocytes, and tumor cells) were analyzed for ENV versus DNA content. Parallel data (forward scatter versus DNA) from a conventional flow cytometer were obtained. RESULTS: ENV linearity studies yielded an R value of 0.999. TRBC had a coefficient of variation (CV) of 1.18 +/- 0.13. DNA indexes as low as 1.02 were detectable. DNA content of lymphocytes from 42 females was 1.9% greater than that for 60 males, with a noninstrumental variability in total DNA content of 0.5%. The ENV/DNA ratio was constant in 15 normal human tissue samples, but differed in the four animal species tested. The ENV/DNA ratio for a hypodiploid breast carcinoma was 2.3 times greater than that for normal breast tissue. CONCLUSIONS: The high-resolution ENV versus DNA analyses are highly reliable, sensitive, and can be used for the detection of near-diploid tumor cells that are difficult to identify with conventional cytometers. ENV/DNA ratio may be a useful parameter for detection of aneuploid populations.

NASA Center JSC↗

DOE OSTI2022

Global centroid moment tensor solutions in a heterogeneous earth: the CMT3D catalogue

SUMMARY For over 40 yr, the global centroid-moment tensor (GCMT) project has determined location and source parameters for globally recorded earthquakes larger than magnitude 5.0. The GCMT database remains a trusted staple for the geophysical community. Its point-source moment-tensor solutions are the result of inversions that model long-period observed seismic waveforms via normal-mode summation for a 1-D reference earth model, augmented by path corrections to capture 3-D variations in surface wave phase speeds, and to account for crustal structure. While this methodology remains essentially unchanged for the ongoing GCMT catalogue, source inversions based on waveform modelling in low-resolution 3-D earth models have revealed small but persistent biases in the standard modelling approach. Keeping pace with the increased capacity and demands of global tomography requires a revised catalogue of centroid-moment tensors (CMT), automatically and reproducibly computed using Green's functions from a state-of-the-art 3-D earth model. In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters, centroid latitude, longitude, depth and centroid time of global earthquakes. We take the GCMT solutions as a point of departure but update them to account for the effects of a heterogeneous earth, using the global 3-D wave speed model GLAD-M25. We generate synthetic seismograms from Green's functions computed by the spectral-element method in the 3-D model, select observed seismic data and remove their instrument response, process synthetic and observed data, select segments of observed and synthetic data based on similarity, and invert for new model parameters of the earthquake’s centroid location, time and moment tensor. The events in our new, preliminary database containing 9382 global event solutions, called CMT3D for ‘3-D centroid-moment tensors’, are on average 4 km shallower, about 1 s earlier, about 5 per cent larger in scalar moment, and more double-couple in nature than in the GCMT catalogue. We discuss in detail the geographical and statistical distributions of the updated solutions, and place them in the context of earlier work. We plan to disseminate our CMT3D solutions via the online ShakeMovie platform.

58 GEOSCIENCES↗

NASA NTRS1984

A comparative review of radiation aerosol models

Aerosol characteristics to be introduced into a set of models within a framework of a Standard Radiation Atmosphere (SRA) are reviewed, and a reconciliation of the Toon and Pollack (TP) global model with the SRA models is attempted. Among the topics discussed are physical and radiative characteristics, radiation aerosol models, albedo for single scattering, and phase function and spectral variation of extinction. Size distribution models are also covered, inclusive the Junge power law, and the modified gamma-, log normal-, and zold-distribution. Recommendations include modifying the SRA profile VI to resemble the Elterman profile modified in the TP model, improving discontinuity at the tropopause level for all SRA models, and retaining the same log-normal size distribution for all SRA components.

Lenoble, J.↗

NASA NTRS2017

Detecting Inter-Annual Variations in the Phenology of Evergreen Conifers Using Long-Term MODIS Vegetation Index Time Series

Long-term observations of vegetation phenology can be used to monitor the response of terrestrial ecosystems to climate change. Satellite remote sensing provides the most efficient means to observe phenological events through time series analysis of vegetation indices such as the Normalized Difference Vegetation Index (NDVI). This study investigates the potential of a Photochemical Reflectance Index (PRI), which has been linked to vegetation light use efficiency, to improve the accuracy of MODIS-based estimates of phenology in an evergreen conifer forest. Timings of the start and end of the growing season (SGS and EGS) were derived from a 13-year-long time series of PRI and NDVI based on a MAIAC (multi-angle implementation of atmospheric correction) processed MODIS dataset and standard MODIS NDVI product data. The derived dates were validated with phenology estimates from ground-based flux tower measurements of ecosystem productivity. Significant correlations were found between the MAIAC time series and ground-estimated SGS (R (sup 2) equals 0.36-0.8), which is remarkable since previous studies have found it difficult to observe inter-annual phenological variations in evergreen vegetation from satellite data. The considerably noisier NDVI product could not accurately predict SGS, and EGS could not be derived successfully from any of the time series. While the strongest relationship overall was found between SGS derived from the ground data and PRI, MAIAC NDVI exhibited high correlations with SGS more consistently (R (sup 2) is greater than 0.6 in all cases). The results suggest that PRI can serve as an effective indicator of spring seasonal transitions, however, additional work is necessary to confirm the relationships observed and to further explore the usefulness of MODIS PRI for detecting phenology.

Ulsig, Laura↗

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↗