Search NASA⌕ Search

SEARCH · Search NASA

Results for “high order methods”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Mechanically Tuned Superconducting Main Injector Cavity

Radio Frequency (RF) superconductivity has been a mainstay of accelerator science for decades. However, its benefits have yet to be applied to proton synchrotrons with demanding tuning requirements. For example, the Main Injector (MI), Fermilab's high-energy proton synchrotron, currently utilizes 20+ ferrite-loaded cavities for a targeted 1.2 s acceleration cycle. Harnessing the extremely high gradients associated with superconductivity, the required number of cavities could be reduced by an order of magnitude, dramatically lowering operational power requirements even with cryogenic considerations. Additionally, the current plans for the Fermilab Accelerator Complex Evolution (ACE) initiative involve almost doubling the number of cavities in MI if the same designs are to be used, further highlighting the potential benefits of superconductivity. These advantages are attractive, but to date, no tunable superconducting cavity suitable for MI has been proposed due to the incompatibility of conventional broadband tuning methods with superconductivity. Here, we present a tunable superconducting cavity concept capable of record-breaking performance. Tuning will be accomplished by using high-speed linear actuators to vary the insertion depth of metallic plungers into the cavity volume. This tuning concept is theoretically viable with currently available technology and will be fully compatible with a superconducting cavity.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Prediction of s ± -Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La 4 Ni 3 O 10 under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La 4 ⁢Ni 3 ⁢O 10 under pressure, here we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d 3z 2 –r 2 orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p z orbitals inducing a “trimer” lattice in La 4⁢ Ni 3 ⁢O 10 , analogous to the dimers of La 3 ⁢Ni 2 ⁢O 7 . The Fermi surface consists of three electron sheets with mixed e g orbitals, and a hole and an electron pocket made up of the d 3⁢z 2 –r 2 orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s ± -wave channel due to partial nesting between the M = (π,π) centered pockets and portions of the Fermi surface centered at the Γ = (0,0) point. With changing electronic density n, the s ± instability remains leading and its pairing strength shows a domelike behavior with a maximum around n = 4.2 ( ~6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La 3 ⁢Ni 2 ⁢O 7 , correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T c superconductivity of La 3⁢ Ni 2 ⁢O 7 does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La 4 ⁢Ni 3 ⁢O 10 with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗
DOE OSTI2026

Towards the Ultimate Strength of Medium‐Entropy Alloys Through Pulsed Lasers

The tensile strength of metals at extreme strain rates is a key predictor of their performance in ballistic and structural impact applications. An important experimental method to reach these extreme strain rates is the use of high-amplitude, short-duration pulsed lasers. The Jupiter Laser Facility at the Lawrence Livermore National Laboratory enabled probing for the first time the mechanical response of several promising High Entropy Alloys at times on the order of nanoseconds (strain rates of ∼10 7 and ∼10 9 s −1 ). The measured strength is in the range of 6 to 10 GPa, ten times the quasistatic value. The mechanisms of plastic deformation and failure were identified and quantified through analysis and molecular dynamics simulation. The reflected wave amplitudes, obtained by VISAR, were used to determine the tensile (spall) stress. The high tensile strength obtained is due to two factors: the strain-rate dependence of plastic flow and the kinetics of void nucleation, growth, and coalescence. The experimental results are compared with an analytical prediction considering both grain-interior and grain-boundary void initiation. Molecular dynamics simulations, conducted at strain rates of 10 8 and 10 9 s −1 , rationalize the experimental results. In conclusion, they provide valuable information about the process of failure evolution, and reveal that grain boundary separation plays a pivotal role in spalling.

Materials science↗
DOE OSTI2024

Growth of ultra-flat ultra-thin alkali antimonide photocathode films

Ultra-flat, ultra-thin alkali antimonide photocathodes with high crystallinity can exhibit high quantum efficiency and low mean transverse energy of outgoing electrons, which are essential requirements for a variety of applications for photocathode materials. Here, we investigate the growth of Cs 3 Sb on graphene-coated 4H–SiC (Gr/4H–SiC), 3C–SiC, and Si 3 N 4 substrates. Sb is deposited using pulsed laser deposition, while Cs is deposited thermally and simultaneously. We demonstrate, employing x-ray analysis and quantum efficiency measurements, that this growth method yields atomically smooth Cs 3 Sb photocathodes with a high quantum efficiency (>10%), even in the ultra-thin limit (<30 nm). For the Si 3 N 4 substrate, film growth is shown to be polycrystalline, while films grown on Gr/4H–SiC show a high degree of ordering with signs of epitaxy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Detection of local-scale changes in greenhouse gas emissions in urban environments using micrometeorological methods and comparison to a high-resolution inventory

We used the Monin–Obukhov similarity theory (MOST) flux-variance relationship to estimate greenhouse gas (GHG) fluxes from high-precision mole fraction measurements at 3 instrumented urban communication towers over several years, demonstrating the ability of this method to detect and quantify changes in emissions. Depending on data availability, we used carbon dioxide (CO 2 ) and carbon monoxide (CO) measurements and/or tracer ratios to estimate fluxes at 1 urban site (Site 3) and 1 suburban site (Site 7) in Indianapolis, IN, USA, and 1 urban site (COM) in Los Angeles, CA, USA. We also compared the estimated fluxes of CO 2 from fossil fuel sources (CO 2 ff) at Sites 3 and 7 and the total CO 2 fluxes at Site 3 to 20 m, hourly resolution subdomains of the high-resolution CO 2 emissions inventory, Hestia, for the year 2020, introducing a new way to evaluate emissions inventories at small spatial and temporal scales. Using the flux-variance relationship, we detected and quantified abrupt decreases in CO and CO 2 fluxes at Site 3 and COM in April 2020, coinciding with the stay-at-home order due to COVID-19 pandemic, as well as abrupt decreases in CO and CO 2 fluxes at Site 3 in July 2018 coinciding with a highway closure next to the site. The Hestia emissions inventory detected a decrease in emissions in April 2020 at Sites 3 and 7, but this decrease differed in magnitude from those detected in the atmospheric estimates. Seasonal trends in emissions are similar between Hestia and the atmospheric estimates at Site 7. We use differences in seasonal and spatial trends between the flux estimation methods to identify potential sources of uncertainty in both the atmospheric and inventory methods. The results from this study show that the flux-variance estimation method is a useful tool to monitor local-scale emissions and evaluate high-resolution emissions inventories.

Carbon dioxide↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2024

A deep learning approach to fast analysis of collective Thomson scattering spectra

Fast analysis of collective Thomson scattering ion acoustic wave features using a deep convolutional neural network model is presented. The network was trained from spectra to predict the plasma parameters, including ion velocities, population fractions, and ion and electron temperatures. A fully kinetic particle-in-cell simulation was used to model a laboratory astrophysics experiment and simulate a diagnostic image of the ion acoustic wave feature. Network predictions were compared with Bayesian inference of the plasma model parameters for both the simulated and experimentally measured images. Both approaches were fairly accurate predicting the simulated image and the network predictions matched a good portion of the Bayesian results for the experimentally measured image. The Bayesian approach is more robust to noise and motivates future work to train deep learning models with realistic noise. The advantage of the deep learning model is making thousands of predictions in a few hundred milliseconds, compared to a few seconds to minutes per prediction for the optimization and Bayesian approaches presented here. The results demonstrate promising capabilities of deep learning models to analyze Thomson data orders of magnitude faster than conventional methods when using the neural network for standalone analysis. If more rigorous analysis is needed, neural network predictions can be used to quickly initialize other optimization methods and increase chances of success. This is especially useful when the dataset becomes very large or highly dimensional and manually refining initial conditions for the entire dataset are no longer tractable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermodynamics of the dipole-octupole pyrochlore magnet Ce 2 Hf 2 O 7 in applied magnetic fields

The recently discovered dipole-octupole pyrochlore magnet Ce 2 Hf 2 O 7 is a promising three-dimensional quantum spin liquid candidate which shows no signs of ordering at low temperature. Here we investigate the thermodynamic response to magnetic fields applied along the global [110] direction using specific heat measurements and fits using numerical methods, and solve the corresponding magnetic structure using neutron diffraction. Specific heat data in moderate fields are reproduced well, however, at high fields the agreement is not satisfactory. We especially observe a two-step release of entropy, a finding that demands a review of both theory and experiment. We address it within the framework of three possible scenarios, including an analysis of the crystal field Hamiltonian not restricted to the two-dimensional single-ion doublet subspace. We conclusively rule out two of these scenarios and find qualitative agreement with a simple model of field misalignment with respect to the crystalline direction. As a result, we discuss the implications of our findings for [111] applied fields and for future experiments on Ce 2 Hf 2 O 7 and its sister compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Elucidating the factors controlling the methanol synthesis from CO 2 on ZnPd-ZnO/TiO 2 catalysts prepared by direct/inverse impregnation

The direct hydrogenation of CO 2 to methanol over supported ZnPd intermetallic catalysts remains a challenge due to the complex nature of the active sites, which are highly dependent on both the dispersion of ZnPd intermetallics and their proximity to ZnO. This study investigates the effect of the Pd and Zn incorporation order on methanol synthesis over TiO 2 -supported ZnPd catalysts. Catalysts were prepared via direct (Pd/Zn), inverse (Zn/Pd), and co-impregnation (Pd+Zn) methods using anatase TiO 2 as the support. Comprehensive characterization by XRD, HRTEM, CO-DRIFTS and XPS revealed that the formation of small and crystalline intermetallic β-ZnPd particles and their interaction with ZnO/ZnO x depends on the sequence of Pd and Zn incorporation. It was found that, under the reaction conditions, the ZnO x species plays a critical role in facilitating the formation of the intermetallic ZnPd through the interactions between Pd and ZnO x . In addition, ZnO x species have also been shown to modulate key catalytic properties, such as CO 2 activation and methanol selectivity. Specifically, the ZnPd/Pd surface ratio, determined by Pd dispersion and proximity to ZnO/ZnO x , was shown to determine the methanol production from CO 2 directly or via a CO intermediate. The Pd/Zn catalyst showed superior selectivity due to the in-situ formation of intermetallic species facilitated by Pd-ZnO x interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Monolayer Films of Binary Metal Oxide Nanoparticles for Carbon Nanotube Synthesis

Binary metal oxide nanoparticles (biMO-NPs) combining transition and main-group metals with well-organized atomic architectures have unique potential as robust and economical heterogeneous catalysts. Herein, metal oxide nanoparticles (CoAlxOy) using cobalt and aluminum were synthesized, and their structure, morphology, and chemical composition were investigated using various characterization techniques. The biMO-NPs, which had a Gaussian-type size distribution (∼6 nm), were assembled on plain silicon substrates modified chemically with linker molecules bearing suitable end groups. Surface analysis revealed an ordered, uniform, close-packed arrangement of biMO-NPs forming a monolayer architecture with nanoscale thickness (∼12 nm), high surface uniformity, and negligible defects. The nanoparticle monolayer film (NP-MF) was employed as a catalyst to grow carbon nanotube (CNT) forests via the thermal chemical vapor deposition (CVD) method, similar to those grown from catalyst thin films deposited by physical vapor deposition. Structural characterization confirmed the growth of a dense, uniform, high-quality vertically aligned carbon nanotube (VA-CNT) forest with a length of ∼125 µm, which is comparable to the VA-CNTs grown from expensive thin films. Thus, CNT growth was successfully catalyzed by the biMO-NPs, highlighting a simple and inexpensive alternative for catalyst deposition. In this innovative wet-chemistry approach, the catalyst and catalyst support are combined within a single nanoparticle before NP-MF assembly. Interestingly, this approach provides a scalable and cost-effective pathway for CNT growth, eliminating the need for expensive thin deposition techniques.

Regmi, Bishow [University of Cincinnati]↗
DOE OSTI2026

Evaluating the Limits of QAOA Parameter Transfer at High-Rounds on Sparse Ising Models With Geometrically Local Cubic Terms

The emergent practical applicability of the Quantum Approximate Optimization Algorithm (QAOA) for approximate combinatorial optimization is a subject of considerable interest. One of the primary limitations of QAOA is the task of finding a set of good parameters, which is usually done using a variational optimization loop. Parameter transfer, or parameter concentration, is a phenomenon where QAOA angles trained on problem instances that are self-similar tend to perform well for other problem instances from that similar class. This suggests a potentially highly efficient and scalable non-variational learning method for QAOA angle finding. In this work, we systematically study QAOA parameter transferability from small problem sizes (16 and 27 decision variables) onto large problem instances (up to 156 qubits) for heavy-hex graph Ising models with geometrically local higher order terms using the Julia based QAOA simulation tool \texttt{JuliQAOA} to perform classical angle finding for up to $49$ QAOA layers ($p$). Parameter transfer of the fixed angles is validated using a combination of full statevector, Projected Entangled Pair States (PEPS), Matrix Product State (MPS), and LOWESA numerical simulations. We find that the QAOA parameter transfer from single instances applied to other (unseen) problem instances does not in general provide monotonically improving performance as a function of $p$ - there are many cases where the performance temporarily decreases as a function of $p$ - but despite this the transferred angles have a general trend of improved expectation value as the QAOA depth increases, in many cases converging close to the true ground-state energy of the $100+$ qubit instances. We also sample the hardware-compatible Ising models using the ensemble of transfer-learned QAOA parameters on several superconducting qubit IBM Quantum processors with 127, 133, and 156 qubits. We find continuous solution quality improvement of the hardware-compatible QAOA circuits run on the IBM NISQ processors up to $p=5$ on \texttt{ibm\_fez}, up to $p=9$ on \texttt{ibm\_torino}, and up to $p=10$ on \texttt{ibm\_pittsburgh}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Flash healing of laser-induced graphene

The advancement of laser-induced graphene (LIG) technology has streamlined the fabrications of flexible graphene devices. However, the ultrafast kinetics triggered by laser irradiation generates intrinsic amorphous characteristics, leading to high resistivity and compromised performance in electronic devices. Healing graphene defects in specific patterns is technologically challenging by conventional methods. Herein, we report the rapid rectification of LIG’s topological defects by flash Joule heating in milliseconds (referred to as F-LIG), whilst preserving its overall structure and porosity. The F-LIG exhibits a decreased I D /I G ratio from 0.84 – 0.33 and increased crystalline domain from Raman analysis, coupled with a 5-fold surge in conductivity. Pair distribution function and atomic-resolution imaging delineate a broader-range order of F-LIG with a shorter C-C bond of 1.425 Å. The improved crystallinity and conductivity of F-LIG with excellent flexibility enables its utilization in high-performance soft electronics and low-voltage disinfections. Notably, our F-LIG/polydimethylsiloxane strain sensor exhibits a gauge factor of 129.3 within 10% strain, which outperforms pristine LIG by 800%, showcasing significant potential for human-machine interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

d -Electron Heavy Fermion Behavior in a Near-Room-Temperature Polar Metallic Ferrimagnet: A Case of Mn 5 SiC

Polycrystalline Mn 5 SiC was synthesized by using a high-temperature solid-state method. Mn 5 SiC adopts a polar space group (Cmc2 1 ) with six crystallographic Mn sites confirmed by X-ray and neutron diffraction, transmission electron microscopy, and second harmonic generation experiments. The complex crystal structure features edge-sharing trigonal prisms and icosahedra, as well as face/edge-sharing pentagonal prisms. Magnetic measurements indicate ferrimagnetic ordering with a transition temperature of 284 K. The ferrimagnetic structure (magnetic space group Cm’c’2 1 ) was further identified by powder neutron diffraction, where collinear Mn spins align along the crystallographic c-axis. The refined magnetic moment for each crystallographic Mn site at 4 K is 1.8(2), −2.42(9), −1.72(8), 0.51(6), 0.50(4), and 1.7(2) μ B . Density functional theory calculations confirm both the metallic behavior and the ferrimagnetic structure observed experimentally and further provide insight into the observed Mn moment dependence across crystallographic sites. The resistivity and specific heat measurements and density functional theory calculations reveal a substantially large Kadowaki–Woods ratio of 5 × 10 –5 μΩ·cm/(mJ/mol) 2 and a many-body renormalization factor of 5.5, indicating the unusual heavy Fermion behavior in such an itinerant magnetic metal.

Messegee, Zachary T. [George Mason University, Fai↗
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↗