Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed surface roughness”

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 55 records · Page 3

Influence of Nut Design on the Load Distribution over the Turns of a Screw

Screw turns are known to be deformed as a result of the contact pliability of the nut and screw surfaces. The load distribution over the turns is investigated when fine thread is employed, when compression nuts are replaced by tension–compression and tension nuts, when the roughness of the working surfaces is increased, and when the load on the screw is increased.

42 ENGINEERING↗

Distinct roles of land cover in regulating spatial variabilities of temperature responses to radiative effects of aerosols and clouds

Surface temperature responses to radiative perturbations due to aerosols and clouds are complicated by the land surface properties. To disentangle these complexities, this study, from a terrestrial surface energy budget perspective, isolates the underlying biophysical processes from the instantaneous radiative effects of aerosols and clouds on surface temperature using the National Center for Atmospheric Research Community Earth System Model version 1.2.1. It is found that in comparison with the global heterogeneous distributions of instantaneous radiative perturbations at the surface induced by aerosols and clouds, the spatial variations of the corresponding surface temperature responses to aerosol direct radiative effects (DRE) during the daytime and cloud radiative effects (CRE) during the nighttime are amplified. It is because of the consistent global distribution of the local surface climate sensitivity (a function of land cover properties such as surface roughness and Bowen ratio) with those of daytime DRE and nighttime CRE. By applying identical anthropogenic aerosol and precursor emissions over eight major past, present and projected future anthropogenic aerosol emitting regions (i.e. Brazil, China, East Africa, India, Indonesia, South Africa, United States and Western Europe), surface temperature responses to aerosol radiative cooling in the daytime and cloud radiative warming in the nighttime over these regions positively regulated by local surface climate sensitivities are prominent.

54 ENVIRONMENTAL SCIENCES↗

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

36 MATERIALS SCIENCE↗

Designing a Convection‐Cloud Chamber for Collision‐Coalescence Using Large‐Eddy Simulation With Bin Microphysics

Abstract Collisional growth of cloud droplets is an essential yet uncertain process for drizzle and precipitation formation. To improve the quantitative understanding of this key component of cloud‐aerosol‐turbulence interactions, observational studies of collision‐coalescence in a controlled laboratory environment are needed. In an existing convection‐cloud chamber (the Pi Chamber), collisional growth is limited by low liquid water content and short droplet residence times. In this work, we use numerical simulations to explore various configurations of a convection‐cloud chamber that may intensify collision‐coalescence. We employ a large‐eddy simulation (LES) model with a size‐resolved (bin) cloud microphysics scheme to explore how cloud properties and the intensity of collision‐coalescence are affected by the chamber size and aspect ratio, surface roughness, side‐wall wetness, side‐wall temperature arrangement, and aerosol injection rate. Simulations without condensation and evaporation within the domain are first performed to explore the turbulence dynamics and wall fluxes. The LES wall fluxes are used to modify the Scalar Flux‐budget Model, which is then applied to demonstrate the need for non‐uniform side‐wall temperature (two side walls as warm as the bottom and the two others as cold as the top) to maintain high supersaturation in a tall chamber. The results of LES with full cloud microphysics reveal that collision‐coalescence is greatly enhanced by employing a taller chamber with saturated side walls, non‐uniform side‐wall temperature, and rough surfaces. For the conditions explored, although lowering the aerosol injection rate broadens the droplet size distribution, favoring collision‐coalescence, the reduced droplet number concentration decreases the frequency of collisions.

54 ENVIRONMENTAL SCIENCES↗

Experimental approach for development of a powder spreading metric in additive manufacturing

The powder spreading is a vital step of powder-based additive manufacturing (AM) processes. The quality of spread powder can considerably influence the properties of fabricated parts. Poorly packed powder beds with high surface roughness result in printed part layers with large porosity and low dimensional accuracy, leading to poor mechanical properties. Therefore, the powder spreadability and its dependence on process parameters and powder characteristics should be quantified to improve the efficiency of powder-based AM methods. This study proposes a novel dimensionless powder spreadability metric that can be commonly used in different powder-based AM processes. The quality of spread powder in terms of powder bed density and surface roughness was evaluated by adjusting the process parameters including recoating velocity and layer thickness, and powder characteristics including particle size distribution. In addition, the dynamic repose angle was proposed and examined as another powder spreadability metric. Furthermore, the results showed that these two proposed metrics were strongly correlated and lower recoating velocity and larger layer thickness led to higher spreadability and lower dynamic repose angle.

36 MATERIALS SCIENCE↗

Electric fields near undulating dielectric membranes

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. Here, we report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic “slab” geometry using a periodic summation of Green’s functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions.

Chemistry↗

Major improvements on the longwave radiative interactions between surface and clouds in the Polar Regions in atmospheric global circulation model (GCM). Final Project Report

1. Summary of Work and Accomplishments in Year 1 (01/15/2015-01/14/2016). We developed new ice cloud parameterizations based on state-of-the-art MODIS Collection 6 (MC6) ice particle model. Features of the model include the consideration of surface roughness of ice crystals and spectral consistency in cloud property retrievals. Using more than 14,000 particle size distributions (PSDs) from in-situ measurements and the latest ice particle single-scattering property library, our parameterizations provide band-averaged bulk-scattering properties. According to statistics of fittings, the parameterization schemes represent the band-averaged bulk-scattering properties very well. Compared to the ice particle model developed by Fu (1998), flux simulations with the MC6 particle model result in larger upward flux at the top of the atmosphere (TOA) and smaller at the surface for typical cirrus clouds. In the assessment of the effect of longwave light scattering, our results also support the conclusions of previous studies. The error in flux calculations reaches about 10 W/m 2 when a cirrus cloud is thin and its particle size is small. 2. Summary of Work and Accomplishments in Year 2 (01/15/2016-01/14/2017) Clouds are a major modulator of the global radiation budgets. However, representations of clouds in general circulation models (GCMs) neglect many important processes, particularly those related to ice clouds. One major uncertainty source is the radiative transfer scheme in GCMs. To reduce the computational burden in most GCMs, absorption is the only radiative transfer process considered in longwave spectral bands. Recently, a state-of-the-art MC6 ice model has been available, which provides spectral consistency in ice cloud retrievals. With an advanced light scattering library for ice particles, we performed sensitivity studies of the potential impacts on the climate including light scattering for MC6 model and other 11 ice cloud models. Specifically, we modified RRTMG_LW (the GCM version of Longwave Rapid Radiative Transfer Model), a radiative transfer module widely used in GCMs and numerical prediction models, to identify flux simulation differences with and without considering light scattering processes. The simulations with light scattering come from a rigorous radiative solver, DISORT (Discrete Ordinates Radiative Transfer Program for a Multi-Layered Plane-Parallel Medium), are used as benchmarks. The results show that the weighted annual mean biases of RRTMG_LW for the upward flux at the top of the atmosphere, the downward flux at the surface, and the net flux into the atmosphere are about 0.8 ± 0.3, -0.1 ± 0.05, and -0.7 ± 0.3 W m-2, respectively. According to the different the ice crystal shapes, the weighted annual mean bias for heating rate is about -0.006 ± 0.02 to 0.04 K/day. 3. Summary of Work and Accomplishments in Year 3 (01/15/2017-01/14/2018) Since absorption dominates optical properties of clouds in the longwave (LW) spectrum, most of general circulation models (GCMs) only take absorption properties of clouds into account in the radiative transfer modules in order to reduce computational costs. To quantify the biases of excluding LW scattering when clouds exist, we simulated fluxes and heating rates by using satellite data from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO), CloudSat, Clouds and the Earth’s Radiant Energy System (CERES) and Moderate Resolution Imaging Spectrometer (MODIS) merged products (CCCM) in 2010. In annual global average, neglecting LW scattering causes overestimation by 2.6 W/m 2 for the top-of-atmosphere (TOA) upward fluxes and underestimation by 1.2 W/m 2 for surface downward fluxes. In particular, regional extreme biases are approximately 12 W/m 2 for TOA upward fluxes and are approximately -3.6 W/m 2 for surface downward fluxes, or are approximately 5% of the global averaged outgoing longwave radiation (OLR) and 1% of the global averaged surface downward fluxes, respectively. In terms of heating rates, on average, LW scattering heats the whole atmosphere column by 0.0045 K/day, heats the cloud layers by 0.0420 K/day, heats the surface by 0.028 K/day, and cools the tropopause by 0.018 K/day.

58 GEOSCIENCES↗

Crystal truncation rods from miscut surfaces

Crystal truncation rods are used to study surface and interface structure. Since real surfaces are always somewhat miscut from a low index plane, it is important to study the effect of miscuts on crystal truncation rods. We develop a model that describes the truncation rod scattering from miscut surfaces that have steps and terraces. We show that nonuniform terrace widths and jagged step edges are both forms of roughness that decrease the intensity of the rods. Nonuniform terrace widths also result in a broad peak that overlaps the rods. We use our model to characterize the terrace width distribution and step edge jaggedness on three SrTiO 3 (001) samples, showing excellent agreement between the model and the data, confirmed by atomic force micrographs of the surface morphology. As a result, we expect our description of terrace roughness will apply to many surfaces, even those without obvious terracing.

36 MATERIALS SCIENCE↗

Atomic step disorder on polycrystalline surfaces leads to spatially inhomogeneous work functions

Structural disorder causes materials’ surface electronic properties, e.g., work function (ϕ), to vary spatially, yet it is challenging to prove exact causal relationships to underlying ensemble disorder, e.g., roughness or granularity. For polycrystalline Pt, nanoscale resolution photoemission threshold mapping reveals a spatially varying ϕ=5.70±0.03 eV over a distribution of (111) vicinal grain surfaces prepared by sputter deposition and annealing. With regard to field emission and related phenomena, e.g., vacuum arc initiation, a salient feature of the ϕ distribution is that it is skewed with a long tail to values down to 5.4 eV, i.e., far below the mean, which is exponentially impactful to field emission via the Fowler–Nordheim relation. We show that the ϕ spatial variation and distribution can be explained by ensemble variations of granular tilts and surface slopes via a Smoluchowski smoothing model wherein local ϕ variations result from spatially varying densities of electric dipole moments, intrinsic to atomic steps, that locally modify ϕ. Atomic step-terrace structure is confirmed with scanning tunneling microscopy (STM) at several locations on our surfaces, and prior works showed STM evidence for atomic step dipoles at various metal surfaces. From our model, we find an atomic step edge dipole μ=0.12 D/edge atom, which is comparable to values reported in studies that utilized other methods and materials. Our results elucidate a connection between macroscopic ϕ and the nanostructure that may contribute to the spread of reported ϕ for Pt and other surfaces and may be useful toward more complete descriptions of polycrystalline metals in the models of field emission and other related vacuum electronics phenomena, e.g., arc initiation.

36 MATERIALS SCIENCE↗

Characteristics of Shock-Induced Dust Flux from the Surface of Lead and Copper Samples under Pressure in the Shock Wave of 45–50 GPa with Unloading into Vacuum

The results of experimental studies of the shock-induced particle ejection (“dusting”) from a free rough (R{sub z} = 20) surface of lead and copper samples into an evacuated medium are presented. The experimental methods were based on different physical principles. To determine the size of particles by shadow laser-optical imaging more efficiently, a narrow strip was isolated in the center of the free surface of a sample from which an optically transparent stream of particles was ejected. The rest of the surface, substantially larger in size, ejected particles, parameters of which were reliably recorded using optical heterodyne interferometry detection, radiographic imaging, and piezoelectric techniques. This made it possible to obtain more reliable data on the particle size distribution, the velocity of the front of a particle flux, and the density (weight) distribution of the flow in the direction of its motion, necessary to refine the existing models and create more reliable models to describe the phenomenon. Using lead and copper samples, the material of which melts or does not melt under shock wave loading in selected close conditions (amplitude of the shock wave and roughness of the free surface), made it possible to clearly demonstrate the effect of melting on the qualitative pattern and quantitative characteristics of the shock-induced dusting process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Au assisted smooth ultrathin epitaxial ZnO film grown by pulsed laser deposition on sapphire(0001)

High quality ZnO has applications in photonics and electronics. Literature has shown that continuous, but rough ZnO film can be grown on sapphire substrate by pulsed laser deposition (PLD). In this work, ZnO and Au were co-deposited by PLD on sapphire(0001) substrates from a single ZnO target overlaid with an Au strip that serves as a catalyst. The reflection high-energy electron diffraction (RHEED) patterns show paired double stripes from sample of ZnO-Au where the difference of in-plane lattice parameters between ZnO and Au is 11.4 %. The two-dimensional (2D) reciprocal space maps constructed from azimuthal RHEED (ARHEED) patterns reveal parallel epitaxial relationship between ZnO and Au in the continuous ZnO-Au film on sapphire substrate. Transmission electron microscopy (TEM) cross section images and energy dispersive spectroscopy maps show a minute amount of Au nanocrystals distributed in the entire ultrathin ZnO film. Atomic force microscopy shows the root-mean-square roughness of the surface of ultrathin ZnO film with Au is in the sub-nm range. The ARHEED and TEM results show that a couple of atomic percentages of Au co-deposited with ZnO catalyzes the growth of smooth ultrathin epitaxial ZnO film.

36 MATERIALS SCIENCE↗

An electrochemical thermal transistor

Abstract The ability to actively regulate heat flow at the nanoscale could be a game changer for applications in thermal management and energy harvesting. Such a breakthrough could also enable the control of heat flow using thermal circuits, in a manner analogous to electronic circuits. Here we demonstrate switchable thermal transistors with an order of magnitude thermal on/off ratio, based on reversible electrochemical lithium intercalation in MoS 2 thin films. We use spatially-resolved time-domain thermoreflectance to map the lithium ion distribution during device operation, and atomic force microscopy to show that the lithiated state correlates with increased thickness and surface roughness. First principles calculations reveal that the thermal conductance modulation is due to phonon scattering by lithium rattler modes, c -axis strain, and stacking disorder. This study lays the foundation for electrochemically-driven nanoscale thermal regulators, and establishes thermal metrology as a useful probe of spatio-temporal intercalant dynamics in nanomaterials.

42 ENGINEERING↗

Sulfate Aerosol in the Arctic: Source Attribution and Radiative Forcing

Source attribution of Arctic sulfate and its radiative forcing due to aerosol-radiation interactions (RFari) for 2010–2014 are quantified in this study using the Community Earth System Model equipped with an explicit sulfur source-tagging technique. The model roughly reproduces the seasonal pattern of sulfate but has biases in simulating the magnitude of near-surface concentrations and vertical distribution. Regions that have high emissions and/or are near/within the Arctic present relatively large contributions to Arctic sulfate burden, with the largest contribution from sources in East Asia (27%). Seasonal variations of the contribution to Arctic sulfate burden from remote sources are strongly influenced by meteorology. The mean RFari of anthropogenic sulfate offsets one third of the positive top of the atmosphere (TOA) RFari from black carbon. A 20% global reduction in anthropogenic SO 2 emissions leads to a net Arctic TOA forcing increase of +0.019W m –2 . These results indicate that a joint reduction in BC and SO 2 emissions could prevent at least some of the Arctic warming from any future SO 2 emission reductions. Furthermore, sulfate RFari efficiency calculations suggest that source regions with short transport pathways and meteorology favoring longer lifetimes are more efficient in influencing the Arctic sulfate RFari. Based on Arctic climate sensitivity factors, about 0.19 K of the Arctic surface temperature cooling is attributed to anthropogenic sulfate, with –0.05 K of that from sources in East Asia, relative to preindustrial conditions.

54 ENVIRONMENTAL SCIENCES↗

Relationships between fault characteristics and seismic responses in a large lab-scale tri-axial injection test conducted on a faulted Castlegate Sandstone

Here to investigate mechanisms causing microseismicity (M w < 2) at a CO 2 injection site, a large-scale triaxial block experiment was carried out on a faulted (saw-cut) cubic-meter of Castlegate Sandstone. The experiment consisted of injection tests at varying differential stresses, while monitoring and recording pore pressure and acoustic emissions (AEs). During the experiment, ~33,000 AEs and ~14 mm of horizontal displacement/slip, like a strike-slip fault movement, occurred. To understand the AE responses and ascertain fault characteristics near the located AEs, we modeled the topography of the fault surface, fault aperture, and fault-gouge thickness using pre- and post-experiment laser scans of the fault surface on each half of the block. Additionally, we characterized surface roughness parallel and perpendicular to slip. Models show crushing and flattening of the fault surfaces can be linked to the spatiotemporal distribution of AEs within 50 mm of the fault surface. Approximately 65% of AEs were in areas with small aperture (≤300 μm); thicker fault gouge was observed in adjacent areas with wider aperture and shows a two-fold reduction in grain size relative to unaltered Castlegate Sandstone. This work provides a conceptual understanding on fault surface evolution, which can be applied towards modeling of seismic slip.

58 GEOSCIENCES↗

Quantitative Defect Analysis of Ablator Capsule Surfaces Using a Leica Confocal Microscope and a High-Density Atomic Force Microscope

High-density carbon (HDC) is being evaluated as an alternative to the current National Ignition Facility (NIF) point-design ablator material (glow discharge plasma, or GDP, plastic) due to its high density and optimal opacity, which leads to a higher implosion velocity. Chemical-vapor-deposition-coated HDC capsules have a near perfect surface figure but a microscopically rough surface. After polishing, the surface becomes smooth at nanometer scales but has numerous micron-sized surface pits, whose volumes, morphology, and distribution must be quantified to guide NIF target selection. Traditional metrology tools for GDP surface defects, such as the atomic force microscope (AFM) based Spheremapper and a phase-shifting differential interferometer, lack the resolution to characterize these localized features. In this paper, we describe how this metrology challenge is met by developing automated surface metrology solutions based on a high-density (HD) AFM and a Leica confocal microscope. Furthermore, these tools are complementary in nature. HD-AFM has a 0.1-μm spatial resolution and determines the overall shape distortion and pit statistics by tracing great circles on a capsule with high throughput.

36 MATERIALS SCIENCE↗

Simulations Supporting the Development of Northstar's Indirect Beam Parameters Monitoring System

NorthStar Medical Radioisotopes, LLC is planning to produce the important medical radioisotope molybdenum-99 (Mo-99), the parent of technetium-99m (Tc-99m), through photonuclear reactions in molybdenum-100 (Mo-100). In this approach, a target comprising multiple thin disks of enriched molybdenum metal is bombarded with a 40-MeV electron beam. Electrons impinged on the molybdenum target produce bremsstrahlung X-rays that cause the nuclear reaction. Because enriched Mo-100 is expensive, there is a desire to utilize as much beam power as possible to achieve maximum production yield and minimize the size of the target. This requirement leads to very high beam power density (and heat deposition in the target), which creates challenging requirements for the cooling of the target. The critical part of the target is the target window. It separates the high-pressure helium cooled target from the vacuum beamline and the subject of structural and thermal stress. The temperature of the target window is proportional to the energy density deposited by the beam, so it is critical to maintain the desired beam profile on the target window. The feasibility of indirectly monitoring the maximum energy density of the beam on the beam window through beam losses at the main collimator (Collimator) before the production target was verified. A model of the NorthStar beam transport line was constructed for this purpose using MAD-X and Tao/Bmad codes. Beam optics were computed for the standard operational scenario, followed by an investigation involving approximately 400 cases with parameter variations in the last tuning quadrupoles. This was done to assess the correlation between losses in the collimator and the peak energy density on the target. We developed a model to explore the potential application of Optical transition radiation (OTR) for controlling beam parameters in the NorthStar beam delivery system. This model was based on a generic formula derived from the fundamental solution of the inhomogeneous wave equation of the vector potential, and allowed us to consider various surfaces, even those with irregular or random features, using numerical integration. We applied the model to OTR generated by relativistic electrons impacting an Inconel® 718 beam window. We examined cases with different levels of the window’s surface roughness, ranging from 0.5 to 3.0 microns of root square mean (RMS) deviation. The results of the OTR simulations provided distributions of OTR photons that can be used to study the limitations of optical systems for controlling beam parameters.

43 PARTICLE ACCELERATORS↗

Computational Analysis of Water Film Thickness During Rain Events for Assessing Hydroplaning Risk Part 2: Rough Road Surfaces

Hydroplaning occurs when a water film forms on a road during a rain event and vehicles are traveling at a speed that does not provide sufficient time for the tires to push the water film out of the tire path. Under these conditions, the tire loses contact with the road and the driver may lose control over the vehicle. The present study is a part of a multiyear research effort aimed at gaining a better understanding of the factors that contribute to hydroplaning risk in order to minimize the occurrence of hydroplaning accidents. Three-dimensional computational fluid dynamics (CFD) analysis was used to investigate the water film thickness on multilane roads, one of the key parameters in evaluating hydroplaning risk. In the first part of the study, water films forming on wide roadways, i.e. multilane highways with 2, 4, and 6 lanes per side, were analyzed, with varying cross slope, longitudinal slope, and rainfall rate. Roadways with and without curbs and drainage were included in the analysis. The analysis in the first part was limited to nearly smooth asphalt or concrete surfaces because the maximum roughness height that can be specified in a CFD model using wall functions to determine shear stress at the road surface boundary is one half the thickness of computational cell layer adjacent to the boundary and the need to resolve the vertical velocity distribution in a thin water film requires cell layers that are a fraction of a millimeter thick. In this part of the study, the analysis was extended to consider the influence of pavement roughness over a range extending up 3.6 mm on the development of water film thickness. Various methods of modeling roughness were investigated and are presented: the roughness height model, a meshed-out geometry of the macrotexture of the pavement, and a porous region model. The computational results were analyzed and compared with experimental measurements performed by Gallaway et al.

42 ENGINEERING↗