Search NASA⌕ Search

SEARCH · Search NASA

Results for “EVAPORATION MODEL”

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 325 records · Page 18

Energy deficit of pulsed-laser field-ionized and field-emitted ions from non-metallic nano-tips

The energy deficit of pulsed-laser field-evaporated ions and field-ionized atoms of an inert gas from the surface of a non-metallic nano-metric tip is reported as a function of the laser intensity, ion current, and temperature. A new model is proposed to explain these results, taking into account the resistive properties of non-metallic nano-tips. A good agreement between the theoretical predictions and the experimental results is obtained for all parameters investigated experimentally. This model is also used to discuss the evaporation behavior of oxides analyzed in laser-assisted atom probe tomography. New insight into the contribution of the electrostatic field and the laser illumination on the evaporation process of non-metallic materials is given.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The atmospheric hydrologic cycle in the ACME v0.3 model

We examine the global water cycle characteristics in the Accelerated Climate Modeling for Energy v0.3 model (a close relative to version 5.3 of the Community Atmosphere Model) in atmosphere-only simulations spanning the years 1980–2005. We evaluate the simulations using a broad range of observational and reanalysis datasets, examine how the simulations change when the horizontal resolution is increased from 1° to 0.25, and compare the simulations against models participating in the the Atmosphere Model Intercomparison Project of the 5th Coupled Model Intercomparison Project (CMIP5). Particular effort has been made to evaluate the model using the best available observational estimates and verifying model biases with additional datasets when differences are known to exist among the observations. Regardless of resolution, the model exhibits several biases: global-mean precipitation, evaporation, and precipitable water are too high, light precipitation occurs too frequently, and the atmospheric residence time of water is too short. Many of these biases are shared by the multi-model mean climate of models participating in CMIP5. The reasons behind regional biases in precipitation are discussed by examining how different fields, such as local evaporation and transport of water vapor, contribute to the bias. Although increasing the horizontal resolution does not drastically change the water cycle, it does lead to a few differences: an increase in global mean precipitation rate, an increase in the fraction of total precipitation that falls over land, more frequent heavy precipitation (>30 mm day -1 ), and a decrease in precipitable water. One of the most notable changes is the shift of precipitation produced by the convective parameterization to that produced by the large-scale microphysics parameterization. We analyze how changes in moisture and circulation with resolution contribute to this shift in the precipitation partitioning. Because changing horizontal resolution requires some re-tuning, the effect of that tuning was evaluated by performing an additional simulation at 1 but using the tunings from the 0.25 simulation. In conclusion, the evaluation shows that the more frequent heavy precipitation, the decrease in precipitable water, and the shift from convective to large-scale precipitation are predominantly due to resolution changes, while tuning changes have a major influence on the global mean precipitation and the land/ ocean partitioning of precipitation.

54 ENVIRONMENTAL SCIENCES↗

ON THE BRIGHT LOOP TOP EMISSION IN POST-ERUPTION ARCADES

Observations of post-eruption arcades (PEAs) in X-rays and EUV reveal strong localized brightenings at the loop top regions. The origins of these brightenings and their dynamics are not well understood to date. Here, we study the dynamics of PEAs using one-dimensional hydrodynamic modeling, focusing on understanding the formation of localized brightening. Our findings suggest that these brightenings are the result of collisions between the counter-streaming chromospheric evaporation from both the footpoints. We perform forward modeling of the emission observed in simulated results in various spectral lines observed by the Extreme-Ultraviolet Imaging Telescope on board Hinode . The forward-modeled intensities in various spectral lines are in close agreement with a flare observed on 2006 December 17 by EIS.

79 ASTRONOMY AND ASTROPHYSICS↗

The Effect of an Equatorial Continent on the Tropical Rain Belt. Part 1: Annual Mean Changes in the ITCZ

The TRACMIP (Tropical Rain Belts with an Annual Cycle and Continent Model Intercomparison Project) ensemble includes slab-ocean aquaplanet controls and experiments with a highly idealized tropical continent, characterized by modified aquaplanet grid cells with increased evaporative resistance, increased albedo, reduced heat capacity, and no ocean heat transport (zero Q-flux). In the annual mean, an equatorial cold tongue develops west of the continent and induces dry anomalies and a split in the oceanic intertropical convergence zone (ITCZ). Ocean cooling is initiated by advection of cold, dry air from the winter portion of the continent; warm, humid anomalies in the summer portion are restricted to the continent by anomalous surface convergence. The surface energy budget suggests that ocean cooling persists and intensifies because of a positive feedback between a colder surface, drier and colder air, reduced downwelling longwave (LW) flux, and enhanced net surface LW cooling (LW feedback). A feedback between wind, evaporation, and SST (so-called WES feedback) also contributes to the establishment and maintenance of the cold tongue. Simulations with a gray-radiation model and simulations that diverge from protocol (with negligible winter cooling) confirm the importance of moist-radiative feedbacks and of rectification effects on the seasonal cycle. This mechanism coupling the continental and oceanic climate might be relevant to the double ITCZ bias. The key role of the LW feedback suggests that the study of interactions between monsoons and oceanic ITCZs requires full-physics models and a hierarchy of land models that considers evaporative processes alongside heat capacity as a defining characteristic of land.

54 ENVIRONMENTAL SCIENCES↗

Moisture sources for precipitation variability over the Arabian Peninsula

Abstract We apply the Lagrangian-based moisture back trajectory method to two reanalysis datasets to determine the moisture sources for wet season precipitation over the Arabian Peninsula, defined as land on the Asian continent to the south of the Turkish border and west of Iran. To accomplish this, we make use of the evaporative source region between 65°W–120°E and 30°S–60°N, which is divided into twelve sub-regions. Our comparison of reanalyses and multiple observations allows us to validate datasets and highlight broad-scale similarities in characteristics, notwithstanding some inconsistencies in the southwest AP. The results indicate north-to-south spatiotemporal heterogeneity in the characteristics of dominant moisture sources. In the north, moisture for precipitation is mainly sourced from midlatitude land and water bodies, such as the Mediterranean and Caspian Seas. Areas further south are dependent on moisture transport from the Western Indian Ocean and parts of the African continent. The El Niño-Southern Oscillation (ENSO) exhibits an overall positive but sub-seasonally varying influence on the precipitation variability over the region, with noticeable moisture anomalies from all major source regions. A significant drying trend exists over parts of the Peninsula, which both reanalyses partially attribute to anomalies in the moisture advection from the Congo Basin and South Atlantic Ocean. However, considerable uncertainty in evaporation trends over the terrestrial evaporative sources in observations warrants additional modeling studies to further our understanding of key processes contributing to the negative trends.

54 ENVIRONMENTAL SCIENCES↗

Climate Conundrum: A Wet or Dry European and Northern African Climate During the Middle Miocene

Abstract End of 21st‐century hydroclimate projections suggest an expansion of subtropical dry zones, with Mediterranean and Sahel regions becoming much drier. However, paleobotanical assemblage evidence from the middle Miocene (17‐12 Ma), suggests both regions were instead humid environments. Here we show that by modifying regional sea surface temperatures (SST) in an Earth System Model (CESM1.2) simulation of the middle Miocene, the increased ocean evaporation and integrated water vapor flux overrides any drying effects associated with warming‐induced land‐surface evaporation driven by atmospheric CO 2 concentrations. These modifications markedly reduce the bias in the model‐data comparison for this period. A vegetation model (BIOME4) forced with simulated climatologies predicts both regions were dominated by mixed forest, which is largely consistent with the paleobotanical record. This study unveils the potential for wetter subtropical Mediterranean climates associated with warming, presenting an alternative scenario from future drying projections with localized SST warming governing regional climate change.

Acosta, R. P.↗

PHOTOPHORETIC SEPARATION OF METALS AND SILICATES: THE FORMATION OF MERCURY-LIKE PLANETS AND METAL DEPLETION IN CHONDRITES

Mercury's high uncompressed mass density suggests that the planet is largely composed of iron, either bound within metal (mainly Fe-Ni) or iron sulfide. Recent results from the MESSENGER mission to Mercury imply a low temperature history of the planet which questions the standard formation models of impact mantle stripping or evaporation to explain the high metal content. Like Mercury, the two smallest extrasolar rocky planets with mass and size determination, CoRoT-7b and Kepler-10b, were found to be of high density. As they orbit close to their host stars, this indicates that iron-rich inner planets might not be a nuisance of the solar system but be part of a general scheme of planet formation. From undifferentiated chondrites, it is also known that the metal to silicate ratio is highly variable, which must be ascribed to preplanetary fractionation processes. Due to this fractionation, most chondritic parent bodies-most of them originated in the asteroid belt-are depleted in iron relative to average solar system abundances. The astrophysical processes leading to metal silicate fractionation in the solar nebula are essentially unknown. Here, we consider photophoretic forces. As these forces particularly act on irradiated solids, they might play a significant role in the composition of planetesimals forming at the inner edge of protoplanetary disks. Photophoresis can separate high thermal conductivity materials (iron) from lower thermal conductivity solids (silicate). We suggest that the silicates are preferentially pushed into the optically thick disk. Subsequent planetesimal formation at the edge moving outward leads to metal-rich planetesimals close to the star and metal depleted planetesimals farther out in the nebula.

79 ASTRONOMY AND ASTROPHYSICS↗

Atomic layer deposited aluminum nitride for interface protection of lithium germanium thiophosphate solid electrolytes

Lithium germanium thiophosphate (LGPS) is an attractive solid-state electrolyte material due to its exceptionally high ionic conductivity (~1.2 × 10 –2 S cm –1 ), comparable to many organic liquid electrolytes commonly used in batteries. Despite the high conductivity of LGPS, the susceptibility of LGPS to deleterious degradation reactions has impeded its commercial adoption into solid-state batteries. In particular, the poor voltage stability of LGPS with high-voltage cathode or lithium metal potentials often results in dramatically increasing cell impedance during galvanic cycling. Here, we use aluminum nitride (AlN) as a protection layer for Li metal anode, applied directly to the LGPS at 250 °C using plasma-enhanced atomic layer deposition, to enhance cell performance by preventing LGPS-Li metal reactions. We compare the surface chemistry and electrochemical cycling performance of atomic layer-deposited AlN grown using both plasma N 2 and NH 3 precursors. Galvanostatic cycling and electrochemical impedance spectroscopy show that AlN-coated LGPS cells perform better than bare LGPS cells in contact with Li metal anodes, with the AlN able to improve cycling longevity by over a factor of 3 in certain cases. Finally, we utilize x-ray photoelectron spectroscopy (XPS) line scans to highlight the slow room-temperature reactivity between AlN and evaporated lithium metal, and a computational model is built to aid further XPS analysis.

25 ENERGY STORAGE↗

Measurement of isospin mixing in {sup 80}Zr{sup *} at finite temperature

Isospin mixing induced by Coulomb interaction has been measured in the compound nucleus Z = N = 40 {sup 80}Zr at T{approx}2 MeV produced in a fusion-evaporation reaction. The observable sensitive to the isospin purity of the compound nucleus is the Giant Dipole Resonance {gamma} decay. The degree of mixing of the compound nucleus has been obtained via Statistical Model analysis of the measured {gamma} spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-F-T-217: A Comprehensive Monte-Carlo Study of Out-Of-Field Secondary Neutron Spectra in a Scanned-Beam Proton Therapy Treatment Room

Purpose: To simulate secondary neutron radiation-fields produced at different positions during phantom irradiation inside a scanning proton therapy gantry treatment room. Further, to identify origin, energy distribution and angular emission as function of proton beam energy. Methods: GEANT4 and FLUKA Monte-Carlo codes were used to model the relevant parts of the treatment room in a gantry-equipped pencil beam scanning proton therapy facility including walls, floor, metallic gantry-components, patient table and the homogeneous PMMA target. The proton beams were modeled based on experimental beam ranges in water and spot shapes in air. Neutron energy spectra were simulated at 0°, 45°, 90° and 135° relative to the beam axis at 2m distance from isocenter, as well as 11×11 cm2 fields for 75MeV, 140MeV, 200MeV and for 118MeV with 5cm PMMA range-shifter. The total neutron energy distribution was recorded for these four positions and proton energies. Additionally, the room-components generating secondary neutrons in the room and their contributions to the total spectrum were identified and quantified. Results: FLUKA and GEANT4 simulated neutron spectra showed good general agreement in the whole energy range of 10{sup −}9 to 10{sup 2} MeV. Comparison of measured spectra with the simulated contributions of the various room components helped to limit the complexity of the room model, by identifying the dominant contributions to the secondary neutron spectrum. The iron of the bending magnet and counterweight were identified as sources of secondary evaporation-neutrons, which were lacking in simplified room models. Conclusion: Thorough Monte-Carlo simulations have been performed to complement Bonner-sphere spectrometry measurements of secondary neutrons in a clinical proton therapy treatment room. Such calculations helped disentangling the origin of secondary neutrons and their dominant contributions to measured spectra, besides providing a useful validation of widely used Monte-Carlo packages in comparison to experimental data. Cluster of Excellence of the German Research Foundation (DFG) “Munich-Centre for Advanced Photonics (MAP)”.

60 APPLIED LIFE SCIENCES↗

Clumped isotope constraints on changes in latest Pleistocene hydroclimate in the northwestern Great Basin: Lake Surprise, California

Abstract During the Last Glacial Maximum (LGM) and subsequent deglaciation, the Great Basin in the southwestern United States was covered by numerous extensive closed-basin lakes, in stark contrast with the predominately arid climate observed today. This transition from lakes in the Late Pleistocene to modern aridity implies large changes in the regional water balance. Whether these changes were driven by increased precipitation rates due to changes in atmospheric dynamics, decreased evaporation rates resulting from temperature depression and summer insolation changes, or some combination of the two remains uncertain. The factors contributing to these large-scale changes in hydroclimate are critical to resolve, given that this region is poised to undergo future anthropogenic-forced climate changes with large uncertainties in model simulations for the 21st century. Furthermore, there are ambiguous constraints on the magnitude and even the sign of changes in key hydroclimate variables between the Last Glacial Maximum and the present day in both proxy reconstructions and climate model analyses of the region. Here we report thermodynamically derived estimates of changes in temperature, precipitation, and evaporation rates, as well as the isotopic composition of lake water, using clumped isotope data from an ancient lake in the northwestern Great Basin, Lake Surprise (California). Compared to modern climate, mean annual air temperature at Lake Surprise was 4.7 °C lower during the Last Glacial Maximum, with decreased evaporation rates and similar precipitation rates to modern. During the mid-deglacial period, the growth of Lake Surprise implied that the lake hydrologic budget briefly departed from steady state. Our reconstructions indicate that this growth took place rapidly, while the subsequent lake regression took place over several thousand years. Using models for precipitation and evaporation constrained from clumped isotope results, we determine that the disappearance of Lake Surprise coincided with a moderate increase in lake temperature, along with increasing evaporation rates outpacing increasing precipitation rates. Concomitant analysis of proxy data and climate model simulations for the Last Glacial Maximum are used to provide a robust means to understand past climate change, and by extension, predict how current hydroclimates may respond to expected future climate forcings. We suggest that an expansion of this analysis to more basins across a larger spatial scale could provide valuable insight into proposed climate forcings, and aid in climate model process depiction. Ultimately, our analysis highlights the importance of temperature-driven evaporation as a mechanism for lake growth and retreat in this region.

Geology↗

Combustion of petroleum-based transportation fuels and their blends with biofuels: a new approach for developing surrogates and understanding the effects of blending

Accurate and versatile reaction mechanisms are necessary to simulate combustion of transportation fuels blended with biofuels for improving fuel efficiency, reducing fuel consumption and mitigating the formation and emissions of combustion products harmful to the atmosphere. This project involved a combined effort of simulation and experimentation of liquid fuel burning using a combustion configuration amenable to numerical modeling with a level of detail not previously achieved for blends of biofuels with gasoline certification and surrogate fuels. The configuration of an isolated droplet burning under conditions where gas transport arises solely from fuel evaporation was selected as the platform for experiments and numerical modeling. The spherical symmetry and one-dimensional gas transport that results enabled simulating the droplet burning process with an precedented level of detail. Processes associated with unsteady gas and liquid transport, formation of particulate and gaseous products, radiation, multicomponent phase equilibrium at the droplet surface and moving boundary effects from droplet evaporation, were incorporated in a single numerical framework. The numerical model was based on the open source code OPENSmoke++ (OS) adapted to incorporate the effects noted above. In addition to simulation, an experimental design was developed for the isolated droplet configuration to acquire the data needed for validating the numerical model. Several broad accomplishments of the project were the following: demonstrating generally excellent agreement between measured and simulated combustion parameters for the fuel blends examined that included mixtures of heptane and isobutanol as a model system, and surrogates comprised of up to seven miscible components mixed with ethanol or isobutanol as representative biofuel additives; using a new approach to validate reaction mechanisms of biofuel blends which incorporated fuel evaporation into the process along with developing an experimental design for acquiring data to compare with simulations; and showing that a wealth of information could be obtained on the combustion physics of biofuels from experiments requiring volumes on the order of only nanoliters at a time thus opening the way to evaluating biofuels synthesized by new processes early in development.

02 PETROLEUM↗

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

Drizzle and Turbulence below closed cellular marine stratocumulus clouds

Drizzle is ubiquitous in marine boundary layer stratocumulus clouds with much of it evaporating before reaching the surface. Ten days of observations made at the Atmospheric Radiation Measurement (ARM)’s Eastern North Atlantic (ENA) site during closed cellular stratocumulus cloud conditions are used to characterize drizzle below the cloud base and its impact on the boundary layer turbulence. Cloud and drizzle micro- and macro-physical properties were retrieved by combining the data from vertically pointing Doppler cloud radar, ceilometer, and microwave radiometer. On average the drizzle shafts were 28.14 km wide, with cloud base rain rate and modal diameter of 1.26 mm day-1 and 145.26 µm respectively. The rain rate at the surface was negligible yielding an average diabatic cooling of -30.60 Wm-2 in the sub-cloud layer. The LWP and turbulence within the boundary layer increased with an increase in the cloud top radiative cooling, however none of these variables exhibited any relationship with cloud base rain rate. For a similar amount of radiative cooling at the cloud top, the average variance of vertical velocity in the sub-cloud layer was about 10% lower during strongly precipitating conditions as compared to lightly precipitating conditions. The reduction in the variance of vertical velocity due to drizzle evaporation was primarily confined to the upper half of the sub-cloud layer and was due reduction in the strengths of the downdrafts. Collectively our results show substantial impact of drizzle evaporation on boundary layer turbulence, necessitating its accurate representation in the Earth System Models.

drizzle↗

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of the photomolecular effect

It is well-known that water in liquid and vapor phases exhibits weak visible-light absorption. Recent experiments, however, show that at the liquid-air interface, absorption drastically increases, accelerating evaporation beyond thermal limits by 2–5 times. Strikingly, evaporation peaks at green wavelengths, despite no corresponding absorptance peak. The underlying mechanism of this observation, termed the photomolecular effect, remains puzzling, particularly as water molecules do not exhibit resonance peaks in the visible spectrum. Here, we present a theoretical model explaining the effect. We show that surface-bound water clusters undergo non-resonant photon-driven evaporation, with green light being particularly effective. Crucially, we do not expect green light to couple more strongly than other wavelengths to the molecules at the surface, rather the energy of green light is more effectively used to vaporize water. Our model accounts for why the evaporation peak does not align with absorptance and provides a quantum mechanical explanation of how a single photon can vaporize an entire molecular cluster. This model challenges conventional views on water-light interactions, revealing a fundamentally non-thermal evaporation mechanism. Furthermore, our findings have implications for water purification, energy-efficient drying, and climate modeling, opening new pathways for optimizing evaporation-based technologies.

Interfacial evaporation↗

A RECONNECTION-DRIVEN MODEL OF THE HARD X-RAY LOOP-TOP SOURCE FROM FLARE 2004 FEBRUARY 26

A compact X-class flare on 2004 February 26 showed a concentrated source of hard X-rays at the tops of the flare’s loops. This was analyzed in previous work and interpreted as plasma heated and compressed by slow magnetosonic shocks (SMSs) generated during post-reconnection retraction of the flux. That work used analytic expressions from a thin flux tube (TFT) model, which neglected many potentially important factors such as thermal conduction and chromospheric evaporation. Here we use a numerical solution of the TFT equations to produce a more comprehensive and accurate model of the same flare, including those effects previously omitted. These simulations corroborate the prior hypothesis that slow-mode shocks persist well after the retraction has ended, thus producing a compact, loop-top source instead of an elongated jet, as steady reconnection models predict. Thermal conduction leads to densities higher than analytic estimates had predicted, and evaporation enhances the density still higher, but at lower temperatures. X-ray light curves and spectra are synthesized by convolving the results from a single TFT simulation with the rate at which flux is reconnected, as measured through motion of flare ribbons, for example. These agree well with light curves observed by RHESSI and GOES and spectra from RHESSI . An image created from a superposition of TFT model runs resembles one produced from RHESSI observations. This suggests that the HXR loop-top source, at least the one observed in this flare, could be the result of SMSs produced in fast reconnection models like Petschek’s.

79 ASTRONOMY AND ASTROPHYSICS↗

Human tear-production rate from closed-eye Schirmer-strip capillary dynamics

A Schirmer tear test (STT) is commonly used to gauge human tear production, especially when dry-eye symptoms present. In an STT, the rounded tip of a standardized paper strip is inserted into the lower fornix of the eye, and the wetted length extending out from the lower lid is recorded after 5 min of eye closure. Longer wetted lengths suggest higher tear production rates. To date, however, there is no methodology to transform STT transient wetting lengths into basal tear- production rates. We develop a physical model to elucidate wetting kinetics in a Schirmer strip. Tear evaporation from the exposed portion of the strip and gravity are accounted for. Careful consideration of the initial depletion of tear in the closed-eye tear prism reveals an initial fast increase in wetted length followed by slower growth. Excellent agreement of the proposed model is achieved with experimental observation. When evaporation is negligible, the slow-growth regime exhibits a linear increase of wetted length in time. The linear-length-growth time regime permits simple calculation of quantitative tear-production rates. In this paper, we suggest measuring several dynamic wetting lengths along a sheathed Schirmer strip and near the 5-min insertion duration followed by fitting to a straight line. The slope of the length-versus-time data gives the basal lacrimal-supply rate.

59 BASIC BIOLOGICAL SCIENCES↗