Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theoretical description and modelling”

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

The Ghent Hybrid model in NuWro: a new neutrino single-pion production model in the GeV regime

Neutrino-induced single-pion production constitutes an essential interaction channel in modern neutrino oscillation experiments, with its products building up a significant fraction of the observable hadronic final states. Frameworks of oscillation analyses strongly rely on Monte Carlo neutrino event generators, which provide theoretical predictions of neutrino interactions on nuclear targets. Thus, it is crucial to integrate state-of-the-art single-pion production models with Monte Carlo simulations to prepare for the upcoming systematics-dominated landscape of neutrino measurements. In this work, we present the implementation of the Ghent Hybrid model for neutrino-induced single-pion production in the NuWro Monte Carlo event generator. The interaction dynamics includes coherently-added contributions from nucleon resonances and a non-resonant background, merged into the pythia branching predictions in the deep-inelastic regime, as instrumented by NuWro. This neutrino-nucleon interaction model is fully incorporated into the nuclear framework of the generator, allowing it to account for the influence of both initial- and final-state nuclear medium effects. We compare the predictions of this integrated implementation with recent pion production data from accelerator-based neutrino experiments. The results of the novel model show improved agreement of the generator predictions with the data and point to the significance of the refined treatment of the description of pion-production processes beyond the ∆ region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2016

Real-time Identification and Control of Satellite Signal Impairments Solution and Application of the Stratonovich Equation Part 1. Theoretical Development

As satellite communications systems become both more complex and reliant with respect to their operating environment, it has become imperative to be able to identify, during real-time operation, the onset of one or more impairments to the quality of overall communications system integrity. One of the most important aspects to monitor of a satellite link operating within the Earth's atmosphere is the signal fading due to the occurrence of rain and/or phase scintillations. This, of course, must be done in the presence of the associated measurement uncertainty or potentially faulty measurement equipment such as in the Advanced Communication Technology Satellite (ACTS) experiment. In the present work, an approach originally suggested in 1991, and apparently still considered iconoclastic, will be significantly developed and applied to the satellite communications link on which the deleterious composite signal fade is the result of one or many component fade mechanisms. Through the measurement (with the attendant uncertainty or 'error' in the measurement) of such a composite fading satellite signal, it is desired to extract the level of each of the individual fading mechanisms so they can be appropriately mitigated before they impact the overall performance of the communications network. Rather than employing simple-minded deterministic filtering to the real-time fading, the present approach is built around all the models and/or descriptions used to describe the individual fade components, including their dynamic evolution. The latter is usually given by a first-order Langevin equation. This circumstance allows the description of the associated temporal transition probability densities of each of the component processes. By using this description, along with the real-time measurements of the composite fade (along with the measurement errors), one can obtain statistical estimates of the levels of each of the component fading mechanisms as well as their predicted values into the future. This is all accomplished by the use of the well-known Stratonovich integro-differential equation that results from the model of the measured signal fade that is also tailored to adaptively adjust the values of the parameters used in the statistical models of the individual fade mechanisms. Three examples of increasing complexity are addressed and solved for the iterative determination of fade component levels from the measured composite signal fade in the presence of measurement error and, in the last case, with uncertainty in the model parameters.

Countermeasure↗
NASA NTRS1989

Interaction of a sodium ion with the water liquid-vapor interface

Molecular dynamics results are presented for the density profile of a sodium ion near the water liquid-vapor interface at 320 K. These results are compared with the predictions of a simple dielectric model for the interaction of a monovalent ion with this interface. The interfacial region described by the model profile is too narrow and the profile decreases too abruptly near the solution interface. Thus, the simple model does not provide a satisfactory description of the molecular dynamics results for ion positions within two molecular diameters from the solution interface where appreciable ion concentrations are observed. These results suggest that surfaces associated with dielectric models of ionic processes at aqueous solution interfaces should be located at least two molecular diameters inside the liquid phase. A free energy expense of about 2 kcal/mol is required to move the ion within two molecular layers of the free water liquid-vapor interface.

NASA Center ARC↗
DOE OSTI2022

Description of global EGAM in the maximum of local frequency during current ramp-up discharges in DIII-D

Energetic-particle-induced geodesic acoustic modes, EGAMs (Fu, Phys. Rev. Let., vol. 101, 2008, pp. 185002), driven by neutral beam injection (NBI), have been observed in many DIII-D tokamak experiments (Nazikian et al., Phys. Rev. Lett., vol. 101, 2008, pp. 185001). This mechanism has been theoretically investigated in (Qiu et al., Plasma Phys. Control. Fusion, vol. 52, 2010, pp. 095003), using a sharp energetic particle distribution function, and in (Qu et al., Plasma Phys. Control. Fusion, vol. 59, 2017, pp. 055018), where the dispersion relation and eigenmode behaviour were obtained for the situation of early beam scenario, that is, for times smaller than the beam slowing down time. In this work, we extend these studies determining the eigenmode for beyond the slowing down time, in a scenario with reverse safety factor q profile, where a small concentration of energetic ions can produce an off-axis maximum in the GAM dispersion relation. Here, the characteristics of EGAM are analytically studied with the drift kinetic equation together with the MHD code NOVA. The toroidal energetic ion transit frequency, coupled with the GAM frequency, produces the maximum in the dispersion relation where the eigenmode can be found. The quantitative correspondence of experimental results with the predictions of the proposed model is analysed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

Multi-Band Atmospheric Correction Algorithm for Ocean Color Retrievals

NASA's current Atmospheric Correction (AC) algorithm for ocean color utilizes two bands and their ratio in the Near Infrared (NIR) to estimate aerosol reflectance and aerosol type. The algorithm then extrapolates the spectral dependence of aerosol reflectance to the visible wavelengths based on modeled spectral dependence of the identified aerosol type. Future advanced ocean color sensors, such as the Ocean Color Instrument (OCI) that will be carried on the Plankton, Aerosol, Cloud, and ocean Ecosystem (PACE) satellite, will be capable of measuring the hyperspectral radiance from 340 to 890 nm at 5-nm spectral resolution and at 7 discrete Short-wave Infrared (SWIR) channels: 940, 1038,1250, 1378, 1615, 2130, and 2260 nm. To optimally employ this unprecedented instrument capability, we propose an improved AC algorithm that utilizes all atmospheric-window channels in the NIR to SWIR spectral range to reduce the uncertainty in the AC process. A theoretical uncertainty analysis of this, namely Multi-Band AC (MBAC), indicates that the algorithm can reduce the uncertainty in remote sensing reflectance (Rrs) retrievals of the ocean caused by sensor random noise. Furthermore, in optically complex waters, where the NIR signal is affected by contributions from highly-reflective turbid waters, the MBAC algorithm can be adaptively weighted to the strongly-absorbing SWIR channels to enable improved ocean color retrievals in coastal waters. We provide here a description of the algorithm and demonstrate the improved performance in ocean color retrievals, relative to the current NASA standard AC algorithm, through comparison with field measurements and assessment of propagated uncertainties in applying the MBAC algorithm to MODIS and simulated PACE OCI data.

PACE↗
DOE OSTI2021

Determining the Structure-Property Relationships of Quasi-Two-Dimensional Semiconductor Nanoplatelets.

We report a theoretical study of CdSe nanoplatelets aimed at identifying the main factors determining their photophysical properties. Using atomic configurations optimized with density functional theory calculations, we computed quasiparticle and exciton binding energies of nanoplatelets with two to seven monolayers. We employed many body perturbation theory at the GW level and solved the Bethe-Salpeter equation to obtain absorption spectra and excitonic properties. Our results, which agree well with recent experiments, were then used to design a model that allows us to disentangle the effects of quantum confinement, strain induced by passivating ligands, and dielectric environment on the electronic properties of nanoplatelets. We found that, for the model to accurately reproduce our first principle results, it is critical to account for surface stress and consider a finite potential barrier and energy-dependent effective masses when describing quantum confinement. Our findings call into question previous assumptions on the validity of an infinite barrier to describe carrier confinement in nanoplatelets, suggesting that it may be possible to optimize interfacial charge transfer and extraction by appropriately choosing passivating ligands. The model developed here is generalizable to core-shell platelets and enables the description of system sizes not yet directly treatable by first-principles calculations.

Greenwood, Arin R.↗
NASA NTRS2003

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The 2-7 May 1998 Storm

A complete description of a self-consistent model of magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves and back on waves are considered self-consistently by solving both equations on a global magnetospheric scale under nonsteady state conditions. The developed model is employed to simulate the entire 2-7 May 1998 storm period. First, the trapped number fluxes of the ring current protons are calculated and presented along with comparison with the data measured by the three- dimensional hot plasma instrument Polar/HYDRA. Incorporating in the model the wave-particle interaction leads to much better agreement between the experimental data and the model results. Second, examining of the wave (MLT, L shell) distributions produced by the model during the storm progress reveals an essential intensification of the wave emission about 2 days after the main phase of the storm. This result is well consistent with the earlier ground-based observations. Finally, the theoretical shapes and the occurrence rates of the wave power spectral densities are studied. It is found that about 2 days after the storm s main phase on 4 May, mainly non-Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗
DOE OSTI2026

Quantum Choreography of the Nucleus: Rotations, Vibrations, and Emergent Structure

Abstract Nuclei are complex many-body quantum systems where interactions of the neutrons and protons via the strong, the weak, and the electromagnetic forces lead to the emergence of simple patterns of energy states that have been described by various theoretical approaches. One of the goals of all the theoretical models is the development of a universal theory that can be applied across the entire chart of nuclides. Significant progress has been made by experiments as well as the increasing sophistication of models, but a universal theory has yet to be established. A recent review of nuclei in the Z = 50–82 region of the chart of nuclides has analyzed all the available compiled data from several decades of studies towards a clarification of the low-lying structure of nuclei. Other reviews have reported and explained the emergence of multiple different shapes in nuclei at somewhat higher excitation energies than the ground state. Some have challenged the interpretation of the first excited K π = 0 + band as a vibration of ground state. This work attempts to provide a guide to determining the nature of the first excited K π = 0 + band in nuclei by the combined use of nuclear lifetimes, energy level evolutions, dynamic moments of inertia, and intrinsic quadrupole moments extracted from transition probabilities. The result is that for a subset of the nuclei in this region, the K π = 0 + band is consistent with the traditional 𝛽-vibration description of an oscillation built on the ground state.

CHFB+5DCH↗
DOE OSTI2024

Wormholes without averaging

After averaging over fermion couplings, SYK has a collective field description that sometimes has “wormhole” solutions. We study the fate of these wormholes when the couplings are fixed. Working mainly in a simple model, we find that the wormhole saddles persist, but that new saddles also appear elsewhere in the integration space — “half-wormholes.” The wormhole contributions depend only weakly on the specific choice of couplings, while the half-wormhole contributions are strongly sensitive. The half-wormholes are crucial for factorization of decoupled systems with fixed couplings, but they vanish after averaging, leaving the non-factorizing wormhole behind.

1/N Expansion↗
DOE OSTI2023

Emergent geometry and duality in the carbon nucleus

Abstract The carbon atom provides the backbone for the complex organic chemistry composing the building blocks of life. The physics of the carbon nucleus in its predominant isotope, 12 C, is similarly full of multifaceted complexity. Here we provide a model-independent density map of the geometry of the nuclear states of 12 C using the ab initio framework of nuclear lattice effective field theory. We find that the well-known but enigmatic Hoyle state is composed of a “bent-arm” or obtuse triangular arrangement of alpha clusters. We identify all of the low-lying nuclear states of 12 C as having an intrinsic shape composed of three alpha clusters forming either an equilateral triangle or an obtuse triangle. The states with the equilateral triangle formation also have a dual description in terms of particle-hole excitations in the mean-field picture.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

A Finite Element Procedure for Calculating Fluid-Structure Interaction Using MSC/NASTRAN

This report is intended to serve two purposes. The first is to present a survey of the theoretical background of the dynamic interaction between a non-viscid, compressible fluid and an elastic structure is presented. Section one presents a short survey of the application of the finite element method (FEM) to the area of fluid-structure-interaction (FSI). Section two describes the mathematical foundation of the structure and fluid with special emphasis on the fluid. The main steps in establishing the finite element (FE) equations for the fluid structure coupling are discussed in section three. The second purpose is to demonstrate the application of MSC/NASTRAN to the solution of FSI problems. Some specific topics, such as fluid structure analogy, acoustic absorption, and acoustic contribution analysis are described in section four. Section five deals with the organization of the acoustic procedure flowchart. Section six includes the most important information that a user needs for applying the acoustic procedure to practical FSI problems. Beginning with some rules concerning the FE modeling of the coupled system, the NASTRAN USER DECKs for the different steps are described. The goal of section seven is to demonstrate the use of the acoustic procedure with some examples. This demonstration includes an analytic verification of selected FE results. The analytical description considers only some aspects of FSI and is not intended to be mathematically complete. Finally, section 8 presents an application of the acoustic procedure to vehicle interior acoustic analysis with selected results.

Chargin, Mladen↗
DOE OSTI2021

Mechanistic Understanding of Methane Combustion over Ni/CeO 2 : A Combined Experimental and Theoretical Approach

Catalytic oxidation of methane (CH 4 ) over nonprecious Ni/CeO 2 catalysts has received a lot of attention due to the large natural gas reserves found in North America and the prohibitive cost of palladium-based catalysts, commonly used for CH 4 oxidation. However, the catalytic mechanism of CH 4 oxidation over Ni/CeO 2 still remains unclear. Moreover, the parameters affecting the reaction rates, the interaction between nickel and CeO 2 , and the reaction intermediates are still not well understood. In this study, kinetic model fitting, CH 4 temperature-programmed reduction-mass spectroscopy (CH 4 TPR-MS), in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), and density functional theory (DFT) calculations were combined to elucidate the mechanism of complete oxidation of CH 4 over Ni/CeO 2 . CH 4 TPR-MS showed that the complete oxidation of CH 4 over Ni/CeO 2 requires 55–120 °C lower compared to bare CeO 2 or Ni/quartz sand; complete oxidation of CH 4 took place when the surface oxygen species were abundant, while partial oxidation products (CO, H 2 ) were formed when the oxygen species were depleted. In situ DRIFTS showed that CH 3 , CH 2 , CO, and CO 2 were formed after CH 4 activation over Ni/CeO 2 , while CH 3 O species were not observed. Combining those findings with kinetic model fitting, a redox Mars–van Krevelen (MvK) mechanism showed the best description of the experimental observations. The MvK mechanism involves the reaction of dissociated oxygen species with gas-phase CH 4 while water inhibits the reaction rate by adsorbing on the oxidized sites. Moreover, CH 4 activation leads to the reduction of the active sites and oxygen vacancy formation followed by reoxidation of the active sites by gas-phase O 2 . A CH 4 oxidation reaction pathway over Ni/CeO 2 is proposed by DFT calculations. In summary, the findings shown here suggest that CH 4 oxidation over Ni/CeO 2 follows a redox MvK mechanism and provides guidance for the rational design of non-precious-metal catalysts for CH 4 oxidation reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Model predictions and visualization of the particle flux on the surface of Mars

Model calculations of the particle flux on the surface of Mars due to the Galactic Cosmic Rays (GCR) can provide guidance on radiobiological research and shielding design studies in support of Mars exploration science objectives. Particle flux calculations for protons, helium ions, and heavy ions are reported for solar minimum and solar maximum conditions. These flux calculations include a description of the altitude variations on the Martian surface using the data obtained by the Mars Global Surveyor (MGS) mission with its Mars Orbiter Laser Altimeter (MOLA) instrument. These particle flux calculations are then used to estimate the average particle hits per cell at various organ depths of a human body in a conceptual shelter vehicle. The estimated particle hits by protons for an average location at skin depth on the Martian surface are about 10 to 100 particle-hits/cell/year and the particle hits by heavy ions are estimated to be 0.001 to 0.01 particle-hits/cell/year.

NASA Center JSC↗
NASA NTRS1984

The Shuttle Orbiter high resolution accelerometer package experiment - Preliminary flight results

A description of the High Resolution Accelerometer Package (HiRAP) experiment, designed to measure rarefied flow aerodynamic accelerations, is given. The ground test calibration factors, as well as post-flight data processing techniques to extract aerodynamic accelerations, are discussed and applied to the recorded reentry data of the recent STS-6 Shuttle Orbiter flight. The ratio of the measured normal-to-axial aerodynamic accelerations during reentry is used to obtain the first flight measurement of the lift-to-drag ratio (L/D) of a winged entry vehicle in the rarefied flow flight regime. The preliminary result for the free-molecule flow L/D is 0.10 + or - 0.03. The initial flight results on L/D are compared with current predictions which are based upon theoretical and empirical considerations. The measured free-molecule flow L/D value is higher than predicted by a factor of about three, indicating that surface reflection is not completely diffuse as currently assumed. In the rarefied flow transition regime, the Shuttle data book-bridging formula fits the flight data adequately. Upper altitude density profiles are also deduced from the measurements and presented. The density profiles show a wave phenomena with amplitude of about 60 percent, relative to a standard model.

Blanchard, R. C.↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Footprints of Kitaev spin liquid in the Fano lineshape of Raman-active optical phonons

We develop a theoretical description of the Raman spectroscopy in the spin-phonon-coupled Kitaev system and show that it can provide observable signatures of fractionalized excitations characteristic of the underlying spin-liquid phase. In particular, we obtain the explicit form of the phonon modes and construct the coupling Hamiltonian based on the D 3d symmetry. We then systematically compute the Raman intensity and show that the spin-phonon coupling renormalizes phonon propagators and generates the salient Fano lineshape. We find that the temperature evolution of the Fano lineshape displays two crossovers, and the low-temperature crossover shows pronounced magnetic-field dependence. We thus identify the observable effect of the Majorana fermions and the Z2 gauge fluxes encoded in the Fano lineshape. Our results are consistent with the phonon Raman scattering experiments in the candidate material α-RuCl 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Here, we use the public code EBHLIGHT to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Droplet-Surface Impingement Dynamics for Intelligent Spray Design

Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1) By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface; 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling; and 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions. Determining the splash/non-splash boundary is important for many practical applications. Coating and cooling processes would each benefit from near-term empirical relations and subsequent models. Such demonstrations can guide theoretical development by providing definitive testing of its predictive capabilities. Thus, empirical relations describing the boundary between splash and non-splash are given for drops impinging upon a dry solid surface and upon a thin fluid film covering a similar surface. Analytical simplification of the power laws describing the boundary between the splash and non-splash regions yields insight into the engineering parameters governing the splash and non-splash outcomes of the fluid droplets. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a dry solid surface. Splashing upon a dry surface is reasonably described by Ca greater than 0.85, reflecting the competing roles of surface tension and viscosity. The power law correlation is shown separating the splashing versus non-splashing regions as developed for droplets impinging upon a thin fluid film covering the solid surface. Splashing upon a thin fluid film, as described by v (pd/s) greater than 63, is governed by fluid density and surface tension, but is rather independent of viscosity. Finally, the data presented here suggests that a more direct dependence upon the surface tension and viscosity, given a better understanding of their interplay, would allow accurate description of the droplet-surface impacts for more complicated situations involving non-Newtonian fluids, specifically those exhibiting viscoelastic behavior.

Wal, Randy L. Vander↗