Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TOPOGRAPHY”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Information or resolution: Which is required from an SEM to study bulk inorganic materials?: Evaluate SEMs’ practical performance

Significant technological advances in scanning electron microscopy (SEM) have been achieved over the past years. Different SEMs can have significant differences in functionality and performance. This work presents the perspectives on selecting an SEM for research on bulk inorganic materials. Understanding materials demands quantitative composition and orientation information, and informative and interpretable images that reveal subtle differences in chemistry, orientation/structure, topography, and electronic structure. The capability to yield informative and interpretable images with high signal-to-noise ratios and spatial resolutions is an overall result of the SEM system as a whole, from the electron optical column to the detection system. The electron optical column determines probe performance. The roles of the detection system are to capture, filter or discriminate, and convert signal electrons to imaging information. The capability to control practical operating parameters including electron probe size and current, acceleration voltage or landing voltage, working distance, detector selection, and signal filtration is inherently determined by the SEM itself. As a platform for various accessories, e.g. an energydispersive spectrometer and an electron backscatter diffraction detector, the properties of the electron optical column, specimen chamber, and stage greatly affect the performance of accessories. Ease-of-use and ease-of-maintenance are of practical importance. It is practically important to select appropriate test specimens, design suitable imaging conditions, and analyze the specimen chamber geometry and dimensions to assess the overall functionality and performance of an SEM. Finally, for an SEM that is controlled/operated with a computer, the stable software and user-friendly interface significantly affect the usability of the SEM.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Role of Liquid Ink Transport in the Direct Placement of Quantum Dot Emitters onto Sub-Micrometer Antennas by Dip-Pen Nanolithography

Here, dip-pen nanolithography (DPN) is used to precisely position core/thick-shell (“giant”) quantum dots (gQDs; ≥10 nm in diameter) exclusively on top of silicon nanodisk antennas (≈500 nm diameter pillars with a height of ≈200 nm), resulting in periodic arrays of hybrid nanostructures and demonstrating a facile integration strategy toward next-generation quantum light sources. A three-step reading-inking-writing approach is employed, where atomic force microscopy (AFM) images of the pre-patterned substrate topography are used as maps to direct accurate placement of nanocrystals. The DPN “ink” comprises gQDs suspended in a non-aqueous carrier solvent, o-dichlorobenzene. Systematic analyses of factors influencing deposition rate for this non-conventional DPN ink are described for flat substrates and used to establish the conditions required to achieve small (sub-500 nm) feature sizes, namely: dwell time, ink-substrate contact angle and ink volume. Finally, it is shown that the rate of solvent transport controls the feature size in which gQDs are found on the substrate, but also that the number and consistency of nanocrystals deposited depends on the stability of the gQD suspension. In conclusion, the results lay the groundwork for expanded use of nanocrystal liquid inks and DPN for fabrication of multi-component nanostructures that are challenging to create using traditional lithographic techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Automated Experiments of Local Non–Linear Behavior in Ferroelectric Materials

An automated experiment in multimodal imaging to probe structural, chemical, and functional behaviors in complex materials and elucidate the dominant physical mechanisms that control device function is developed and implemented. Here, the emergence of non-linear electromechanical responses in piezoresponse force microscopy (PFM) is explored. Non-linear responses in PFM can originate from multiple mechanisms, including intrinsic material responses often controlled by domain structure, surface topography that affects the mechanical phenomena at the tip-surface junction, and the presence of surface contaminants. Using an automated experiment to probe the origins of non-linear behavior in ferroelectric lead titanate (PTO) and ferroelectric Al 0.93 B 0.07 N films, it is found that PTO shows asymmetric nonlinear behavior across a/c domain walls and a broadened high nonlinear response region around c/c domain walls. In contrast, for Al 0.93 B 0.07 N, well-poled regions show high linear piezoelectric responses, when paired with low non-linear responses regions that are multidomain show low linear responses and high nonlinear responses. It is shown that formulating dissimilar exploration strategies in deep kernel learning as alternative hypotheses allows for establishing the preponderant physical mechanisms behind the non-linear behaviors, suggesting that automated experiments can potentially discern between competing physical mechanisms. Furthermore, this technique can also be extended to electron, probe, and chemical imaging.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Weather radar in Nepal: opportunities and challenges in a mountainous region

Extreme rainfall is one of the major causes of natural hazards in the central Himalayan region, including Nepal. The performance of strategies to manage hazards and related risks relies on the accuracy of quantitative hydrometeorological prediction. Rain gauges have traditionally been used to measure the rainfall amount. However, point measurements with limited gauge coverage cannot accurately represent spatial precipitation variability in complex topography. Weather radar have shown potential for useful information on accurate areal rainfall estimates. The Department of Hydrology and Meteorology (DHM) in Nepal installed their first weather radar in 2019 in the western region of the country. Two more radars will be added to the planned radar network in the near future, in the country's central and eastern regions, respectively. We highlight both the opportunities and challenges with radar installation and observation in the mountainous regions. Radar-rainfall estimates across the Himalayas can be useful to inform decision-making in a broad range of infrastructure sectors, including water, energy, construction, transportation, and agriculture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

3D Cosmic Ray Muon Tomography from an Underground Tunnel

Here, we present an underground cosmic ray muon tomographic experiment imaging 3D density of overburden, part of a joint study with differential gravity. Muon data were acquired at four locations within a tunnel beneath Los Alamos, New Mexico, and used in a 3D tomographic inversion to recover the spatial variation in the overlying rock–air interface, and compared with a priori knowledge of the topography. Densities obtained exhibit good agreement with preliminary results of the gravity modeling, which will be presented elsewhere, and are compatible with values reported in the literature. The modeled rock–air interface matches that obtained from LIDAR within 4 m, our resolution, over much of the model volume. This experiment demonstrates the power of cosmic ray muons to image shallow geological targets using underground detectors, whose development as borehole devices will be an important new direction of passive geophysical imaging.

3-D inversion↗
DOE OSTI2022

Improving Simulations of Vegetation Dynamics over the Tibetan Plateau: Role of Atmospheric Forcing Data and Spatial Resolution

The efficacy of vegetation dynamics simulations in offline land surface models (LSMs) largely depends on the quality and spatial resolution of meteorological forcing data. In this study, the Princeton Global Meteorological Forcing Data (PMFD) and the high spatial resolution and upscaled China Meteorological Forcing Data (CMFD) were used to drive the Simplified Simple Biosphere model version 4/Top-down Representation of Interactive Foliage and Flora Including Dynamics (SSiB4/TRIFFID) and investigate how meteorological forcing datasets with different spatial resolutions affect simulations over the Tibetan Plateau (TP), a region with complex topography and sparse observations. By comparing the monthly Leaf Area Index (LAI) and Gross Primary Production (GPP) against observations, we found that SSiB4/TRIFFID driven by upscaled CMFD improved the performance in simulating the spatial distributions of LAI and GPP over the TP, reducing RMSEs by 24.3% and 20.5%, respectively. The multi-year averaged GPP decreased from 364.68 gC m –2 yr –1 to 241.21 gC m –2 yr –1 with the percentage bias dropping from 50.2% to –1.7%. When using the high spatial resolution CMFD, the RMSEs of the spatial distributions of LAI and GPP simulations were further reduced by 7.5% and 9.5%, respectively. This study highlights the importance of more realistic and high-resolution forcing data in simulating vegetation growth and carbon exchange between the atmosphere and biosphere over the TP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Investigating added value of regional climate modeling in North American winter storm track simulations

Extratropical Cyclone (EC) characteristics depend on a combination of large-scale factors and regional processes. However, the latter are considered to be poorly represented in global climate models (GCMs), partly because their resolution is too coarse. This work describes a framework using possibilities given by regional climate models (RCMs) to gain insight into storm activity during winter over North America (NA). Recent past climate period (1981–2005) is considered to assess EC activity over NA using the NCEP regional reanalysis (NARR) as a reference, along with the European reanalysis ERA-Interim (ERAI) and two CMIP5 GCMs used to drive the Canadian Regional Climate Model—version 5 (CRCM5) and the corresponding regional-scale simulations. While ERAI and GCM simulations show basic agreement with NARR in terms of climatological storm track patterns, detailed bias analyses show that, on the one hand, ERAI presents statistically significant positive biases in terms of EC genesis and therefore occurrence while capturing their intensity fairly well. On the other hand, GCMs present large negative intensity biases in the overall NA domain and particularly over NA eastern coast. In addition, storm occurrence over the northwestern topographic regions is highly overestimated. When the CRCM5 is driven by ERAI, no significant skill deterioration arises and, more importantly, all storm characteristics near areas with marked relief and over regions with large water masses are significantly improved with respect to ERAI. Conversely, in GCMdriven simulations, the added value contributed by CRCM5 is less prominent and systematic, except over western NA areas with high topography and over the Western Atlantic coastlines where the most frequent and intense ECs are located. Despite this significant added-value on seasonalmean characteristics, a caveat is raised on the RCM ability to handle storm temporal ‘seriality’, as a measure of their temporal variability at a given location. In fact, the driving models induce some significant footprints on the RCM skill to reproduce the intra-seasonal pattern of storm activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new spatially distributed added value index for regional climate models: the EURO-CORDEX and the CORDEX-CORE highest resolution ensembles

The added value of using regional climate models (RCMs) to downscale data from general circulation models (GCMs) has often been questioned and researched. Although several studies have used different methods to identify (and in some cases quantify) the added value, there is still a need to find a general metric that quantifies the added value of any variable. This paper builds on past studies to propose a new metric of added value in the simulation of present-day climate which measures the difference in the probability density functions (PDFs) at each grid-cell between a model and an observation source, and then compares the results of the RCM and GCM in order to spatially compute the added value index. The same method is also adapted to quantify the climate change downscaling signal in a way that is consistent with the present-day metric. These new metrics are tested on the daily precipitation output from the EURO-CORDEX and CORDEX-CORE projection ensembles and reveal an overall positive added value of RCMs, especially at the tail-end of the distribution. Higher added value is obtained in areas of complex topography and coast-lines, as well as in tropical regions. Areas with large added value in present-day climate are consistent with areas of significant climate change downscaling signal in the RCP 8.5 far future simulations, and when the analysis is repeated at a low-resolution. The use of different resolution observations shows that the added value tends to decrease when models are compared to low-resolution observation datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Drivers of coupled climate model biases in representing Labrador Sea convection

Abstract This study investigates the representation of ocean convection in the Labrador Sea in seven Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 5 and 6 datasets. The relative role of the oceanic and atmospheric biases in the subpolar North Atlantic gyre are explored using regional ocean simulations where the atmospheric forcing or the ocean initial and boundary conditions are replaced by reanalysis data in the absence of interactive air-sea coupling. Commonalities and differences among model behaviors are discussed with the objective of finding a pathway forward to improve the representation of the ocean mean state and variability in a region of fundamental importance for climate variability and change. Results highlight that an improved representation of ocean stratification in the North Atlantic subpolar gyre is urgently needed to constrain future climate change projections. While improving the ocean model resolution in the North Atlantic alone may contribute a better representation of both boundary currents and propagation of heat and freshwater anomalies into the Labrador Sea, it may not be sufficient. Addressing the atmospheric heat flux bias with better resolution in the atmosphere and land topography may allow for deep convection to occur in the Labrador Sea in some of the models that miss it entirely, but the greatest priority remains improving the representation of ocean stratification.

Liu, Guangpeng↗
DOE OSTI2018

Topographic controls on a phreatomagmatic maar-diatreme eruption: field and numerical results from the Holocene Dotsero volcano (Colorado, USA)

We analyzed the eruptive products of Holocene Dotsero volcano (Colorado, USA; ~4.15 ka), which recorded evidence of a progression from effusive magmatic activity through explosive phreatomagmatic maar-forming activity to a final explosive magmatic phase. The nature of the deposits suggests that the irregular and mountainous pre-eruptive topography strongly influenced the formation of directed pyroclastic density currents (PDCs) during the phreatomagmatic phase, where the topographically high northern crater rim acted as a barrier and promoted transport to the south. Furthermore, the crater shape strongly controlled the grain size of the final deposits, providing further evidence that deposit grain size can be a misleading proxy for fragmentation or explosion intensity. We test these hypotheses, as gleaned from fieldwork, by implementing a multiphase numerical model to simulate discrete explosions and their resulting PDCs at maar-diatreme volcanoes. We confirm that a crater shape where one rim is higher than the other promotes multiple pulses of particle transport in the direction of the lower rim. The initial PDCs contain both coarse and fine particles, whereas PDCs generated by later pulses contain, predominately, finer particles. This mechanism provides an alternative explanation to the formation of lapilli tuff/tuff couplets often found at maar-diatreme tephra rings (including Dotsero volcano). Additional modeling shows how the crater depth can control grain size of the resulting PDC deposits. In the case of a deep crater and a relatively weak explosion, it is possible that only the fine particles escape following collapse of the eruptive column, resulting in a fine-grained deposit. This work allows us to improve models developed to describe emplacement of PDCs for different crater geometries and informs on the dynamics of maar-diatreme eruptions, but also highlights the limitations of our knowledge of key parameters, especially those related to the processes of magma-water interaction and syn-eruptive hydrogeology, during such events.

58 GEOSCIENCES↗
DOE OSTI2017

Numerical Modeling Describing the Effects of Heterogeneous Distributions of Asperities on the Quasi-static Evolution of Frictional Slip

A better understanding of how slip accumulates along faults and its relation to the breakdown of shear stress is beneficial to many engineering disciplines, such as, hydraulic fracture and understanding induced seismicity (among others). Asperities forming along a preexisting fault resist the relative motion of the two sides of the interface and occur due to the interaction of the surface topographies. We employ a finite element model to simulate circular partial slip asperities along a nominally flat frictional interface. Shear behavior of our partial slip asperity model closely matched the theory described by Cattaneo. The asperity model was employed to simulate a small section of an experimental fault formed between two bodies of polymethyl methacrylate, which consisted of multiple asperities whose location and sizes were directly measured using a pressure sensitive film. The quasi-static shear behavior of the interface was modeled for cyclical loading conditions, and the frictional dissipation (hysteresis) was normal stress dependent. We further our understanding by synthetically modeling lognormal size distributions of asperities that were randomly distributed in space. Synthetic distributions conserved the real contact area and aspects of the size distributions from the experimental case, allowing us to compare the constitutive behaviors based solely on spacing effects. Traction-slip behavior of the experimental interface appears to be considerably affected by spatial clustering of asperities that was not present in the randomly spaced, synthetic asperity distributions. Finally, estimates of bulk interfacial shear stiffness were determined from the constitutive traction-slip behavior and were comparable to the theoretical estimates of multi-contact interfaces with non-interacting asperities.

58 GEOSCIENCES↗
DOE OSTI2014

Quantifying the sensitivity of ephemeral streams to land disturbance activities in arid ecosystems at the watershed scale

Large areas of public lands administered by the Bureau of Land Management and located in arid regions of the southwestern United States are being considered for the development of utility-scale solar energy facilities. Land-disturbing activities in these desert, alluvium-filled valleys have the potential to adversely affect the hydrologic and ecologic functions of ephemeral streams. Regulation and management of ephemeral streams typically falls under a spectrum of federal, state, and local programs, but scientifically based guidelines for protecting ephemeral streams with respect to land-development activities are largely nonexistent. This study developed an assessment approach for quantifying the sensitivity to land disturbance of ephemeral stream reaches located in proposed solar energy zones (SEZs). The ephemeral stream assessment approach used publicly-available geospatial data on hydrology, topography, surficial geology, and soil characteristics, as well as highresolution aerial imagery. These datasets were used to inform a professional judgment-based score index of potential land disturbance impacts on selected critical functions of ephemeral streams, including flow and sediment conveyance, ecological habitat value, and groundwater recharge. The total sensitivity scores (sum of scores for the critical stream functions of flow and sediment conveyance, ecological habitats, and groundwater recharge) were used to identify highly sensitive stream reaches to inform decisions on developable areas in SEZs. Total sensitivity scores typically reflected the scores of the individual stream functions; some exceptions pertain to groundwater recharge and ecological habitats. The primary limitations of this assessment approach were the lack of high-resolution identification of ephemeral stream channels in the existing National Hydrography Dataset, and the lack of mechanistic processes describing potential impacts on ephemeral stream functions at the watershed scale.The primary strength of this assessment approach is that it allows watershed-scale planning for low-impact development in arid ecosystems; the qualitative scoring of potential impacts can also be adjusted to accommodate new geospatial data, and to allow for expert and stakeholder input into decisions regarding the identification and potential avoidance of highly sensitive stream reaches.

Arid ecosystems↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗
DOE OSTI2020

Zonal-Based Emission Source Term Model for Predicting Particulate Emission Factors in Wildfire Simulations

A physics/chemistry-based numerical model for predicting the emission of fine particles from wildfires is proposed. This model implements the fundamental mechanisms of soot formation in a combustion environment: soot nucleation, surface growth, agglomeration, oxidation, and particle fragmentation. These mechanisms occur on a scale too fine for the discretization of most wildfire models, which need to simulate landscape-scale dynamics. As a result this model implements a zonal approach, where the computed soot particle distribution is partitioned into process zones within a single resolved grid cell. These process zones include: an inception zone (for nucleation), a heating zone (for coagulation, surface growth, and fragmentation), a reaction zone (for oxidation), and a quenched zone (for atmospheric processes). Governing mechanisms are applied to the appropriate zones to predict total particle growth and emission. Additionally, the proposed model is implemented into HIGRAD/FIRETEC, a physics-based wildfire simulation code which couples interactions between fire, fuels, atmosphere, and topography on a landscape scale. Fire simulations among grasslands and conifer forests are performed and compared against experimental data for emission factors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Influence of Ba Substitutions on Microstructure and Magnetic Properties of La–Ca–Co System Permanent Ferrite

La{sub 0.6}Ca{sub 0.4−x}Ba{sub x}O⋅n/2(Fe{sub 1−y}Co{sub y}){sub 2}O{sub 3}(x = 0, 0.05, 0.1, 0.15 and 0.2) permanent ferrite materials were prepared by ceramic technique. The topography and the structure were characterized by scanning electron microscopy (SEM) and x-ray diffraction (XRD), respectively. The magnetic properties were measured using vibrating sample magnetometer (VSM) and digital magnetometer. Single M-phase and flake-like grains were obtained for samples without Ba substitution(x = 0), and the magnetic properties of the pre-sintering La–Ca–Co ferrite were as follows: B{sub r} = 440 mT, H{sub cb} = 296 kA/m, H{sub cj} = 426 kA/m, (BH){sub max} = 36.6 kJ/m{sup 3} and H{sub k}/H{sub cj} = 0.83. When Ca was substituted by Ba, pre-sintered powder presents single M-phase structure and close three-dimensional size of grains, resulting in an improved rectangle ratio of sintered magnet. The rectangle ratio of sintered magnet reaches a maximum of 0.93 when x = 0.1. By optimizing process parameters, M-type La–Ca–Ba–Co sintered magnet was prepared with B{sub r} = 445 mT, H{sub cb} = 337 kA/m, H{sub cj} = 420 kA/m, (BH){sub max} = 38.6 kJ/m{sup 3} and H{sub k}/H{sub cj} = 0.95.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Advanced geophysical underground coal gasification monitoring

Underground Coal Gasification (UCG) produces less surface impact, atmospheric pollutants and greenhouse gas than traditional surface mining and combustion. Therefore, it may be useful in mitigating global change caused by anthropogenic activities. Careful monitoring of the UCG process is essential in minimizing environmental impact. In this paper, we first summarize monitoring methods that have been used in previous UCG field trials. We then discuss in more detail a number of promising advanced geophysical techniques. These methods – seismic, electromagnetic, and remote sensing techniques – may provide improved and cost-effective ways to image both the subsurface cavity growth and surface subsidence effects. Active and passive seismic data have the promise to monitor the burn front, cavity growth, and observe cavity collapse events. Electrical resistance tomography (ERT) produces near real time tomographic images autonomously, monitors the burn front and images the cavity using low-cost sensors, typically running within boreholes. Interferometric synthetic aperture radar (InSAR) is a remote sensing technique that has the capability to monitor surface subsidence over the wide area of a commercial-scale UCG operation at a low cost. It may be possible to infer cavity geometry from InSAR (or other surface topography) data using geomechanical modeling. The expected signals from these monitoring methods are described along with interpretive modeling for typical UCG cavities. They are illustrated using field results from UCG trials and other relevant subsurface operations.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2017

Simulations of Seismic Wave Propagation on Mars

In this paper, we present global and regional synthetic seismograms computed for 1D and 3D Mars models based on the spectral-element method. For global simulations, we implemented a radially-symmetric Mars model with a 110 km thick crust. For this 1D model, we successfully benchmarked the 3D seismic wave propagation solver SPECFEM3D_GLOBE against the 2D axisymmetric wave propagation solver AxiSEM at periods down to 10 s. We also present higher-resolution body-wave simulations with AxiSEM down to 1 s in a model with a more complex 1D crust, revealing wave propagation effects that would have been difficult to interpret based on ray theory. For 3D global simulations based on SPECFEM3D_GLOBE, we superimposed 3D crustal thickness variations capturing the distinct crustal dichotomy between Mars’ northern and southern hemispheres, as well as topography, ellipticity, gravity, and rotation. The global simulations clearly indicate that the 3D crust speeds up body waves compared to the reference 1D model, whereas it significantly changes surface waveforms and their dispersive character depending on its thickness. We also perform regional simulations with the solver SES3D based on 3D crustal models derived from surface composition, thereby addressing the effects of various distinct crustal features down to 2 s. The regional simulations confirm the strong effects of crustal variations on waveforms. Finally, we conclude that the numerical tools are ready for examining more scenarios, including various other seismic models and sources.

58 GEOSCIENCES↗
DOE OSTI2024

Enhancing the Range and Reliability of the Spacer Layer Imaging Method

The spacer layer imaging method (SLIM) is widely used to measure the thickness of additive and lubricant films, in lubricant development and evaluation, and for fundamental research into elastohydrodynamic lubrication and tribofilm formation mechanisms. The film thickness measurement, as implemented on several popular tribometers, provides powerful, non-destructive in-situ mapping of film topography with nanometre-scale height sensitivity. However, the results can be highly sensitive to experimental procedure, machine condition, and image analysis, in some cases reporting unphysical film thickness trends. The prevailing image analysis techniques make it challenging to interrogate these errors, often hiding their multivariate nonlinear behaviour from the user by spatial averaging. Herein, several common ‘silent errors’ in the SLIM measurement, including colour matching to incorrect fringe orders, and colour drift due to the optical properties of the system or film itself, are discussed, with examples. A robust suite of novel a priori and a posteriori methods to address these issues, and to improve the accuracy and reliability of the measurement, are also presented, including a novel, computationally inexpensive circle-finding algorithm for automated image processing. In combination, these methods allow reliable mapping of films up to at least 800 nm in thickness, representing a significant milestone for the utility of SLIM applied to elastohydrodynamic contact.

EHL film geometry↗