Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow velocity”

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 919 records · Page 51

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗

Stokes-dependent droplet collection efficiency on a NACA 0012 airfoil from droplet-informed simulations with statistical overloading

Accurate modelling of ice accretion on aircraft wings requires analysing droplet impingement on the surface to optimize the design of ice-protection systems. We perform Euler–Lagrange simulations of a droplet-laden flow impinging on a NACA 0012 airfoil. Our study includes water droplets with eight discrete sizes ranging from 1 to 160 microns. We vary the free-stream velocity of the incoming airflow in the range 60 ≤ U ≤ 240 m s −1 and the chord length of the airfoil in the range 0.5 ≤ c ≤ 2 m. Due to the dilute nature of supercooled clouds, one-way coupling is used in the simulations. The effects of droplet breakup and collision are also neglected. To reduce the computational cost, we employ statistical overloading of droplets, allowing us to simulate millions of impinging droplets in a time span on the order of milliseconds. Our results show that the droplet collection efficiency, which measures the likelihood of droplet impingement on the airfoil surface, increases with droplet size and free-stream velocity but decreases with airfoil size. We demonstrate that collection efficiency, impingement velocity and impingement angle are primarily dictated by a single non-dimensional parameter, the droplet Stokes number. We also identify a critical stagnation-streamline Stokes number below which impingements do not occur and use it to estimate the minimum droplet size for impingement. In addition, we observe droplet behaviour to become Stokes number independent at large values of the Stokes number. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗

Integrated Turbine Component Cooling Designs Facilitated by Additive Manufacturing and Optimization

This research program involved the development of improved gas turbine cooling by optimizing the combined internal and film cooling configurations used for turbine blades and vanes. Improved cooling of turbine airfoils allows the gas turbines to be operated at higher temperatures which improves efficiency for gas turbine-based power generation. This also facilitates operation with hydrogen-based fuels which inherently generate higher operating temperatures. Designs of the new cooling configurations were based on the added geometric flexibility that is available when using metal additive manufacturing techniques. When making these designs, adjustments were made to account for imperfections occurring in the metal additive manufacturing process. New designs for film cooling holes were developed using computational adjoint gradient optimization techniques which generated a unique film cooling hole with 70% greater film cooling effectiveness compared to conventional shaped film cooling holes. A major factor in improving the film cooling effectiveness were small external protrusions on either side of the new hole which generated vortices that pushed coolant towards the surface and increased the lateral spreading of the coolant. Our studies also included optimization of film cooling holes that can be built with conventional manufacturing techniques and optimizing the internal coolant channels that supply coolant to the film cooling holes. Various internal coolant channel shapes were investigated and designs that maximized the cooling of the airfoil while minimizing pressure losses in the channel flows were identified. Testing and verification of the performance of the final configurations was done with a transonic wind tunnel facility which allowed testing at engine realistic high velocities. Final testing of the combined internal and film cooling configurations was done by incorporating the configurations in Penn State’s National Experimental Turbine (NExT) vane. These tests confirmed the enhanced overall cooling effectiveness provided by the new cooling configuration designs.

03 NATURAL GAS↗

The role of RF-induced E×B flows in the mitigation of scrape-off-layer convective transport during ion cyclotron resonance heating

Abstract While multiple experiments have reported a decrease in intermittent fluctuations in the far Scrape-Off-Layer (SOL) during ion cyclotron resonance heating (Antar et al 2010 Phys. Rev. Lett. 105 165001, Li et al 2022 Nucl. Eng. Technol. 54 207–19, Antar et al 2012 Nucl. Fusion 52 103005), the physical mechanism behind this observation has not been fully established yet. In this work, we demonstrate, for the first time, a direct correlation between the amplitude of RF-induced E×B flows and turbulence suppression in the far SOL. Using the Gas Puff Imaging (GPI) diagnostic on Alcator C-Mod, we show again that Ion Cyclotron Range of Frequencies can significantly alter the flow in the SOL and introduce a shear layer in regions magnetically connected close to the antenna (Cziegler et al 2012 Plasma Phys. Control. Fusion 54 105019). With the 4-strap field-aligned antenna operated in dipole phasing, the ratio of the power coupled by the central two straps to the power coupled by the outer two straps was varied. The resulting RF-induced radial electric field magnitude thus varied, and we show that the impact on the far SOL turbulence correlates with the modified E×B velocity. We then apply a newly-developed blob tracking algorithm (Han et al 2022 Sci. Rep. 12 18142) to higher-resolution GPI videos in order to directly observe the process of blob shearing by RF-induced E×B flows. We show that the radially sheared poloidal flows act as a transport barrier by stretching, stopping, and destroying filaments, which is consistent with the observed difference in turbulence statistics.

Diab, R. (ORCID:0000000342065868)↗

Numerical study to predict the effect of surface roughness on the thermal and hydraulic performance of additively manufactured heat exchangers

Heat exchangers are used in many industries to transfer thermal energy between two or more objects with an associated thermal and hydraulic performance indicated by thermal effectiveness and pressure drop respectively. Many investigations have aimed to enhance the thermal effectiveness of heat exchangers by introducing turbulent flow inducers that increase the heat transfer rate at the expense of increasing pressure drop. Flow inducers like surface roughness can be controlled with advanced manufacturing methods like additive manufacturing for heat exchanger production. Selecting and manufacturing an optimal surface roughness for heat exchangers requires a more comprehensive approach that accounts for the surface roughness effect on the thermal and hydraulic performance of a sample. Since fabricating samples with additive manufacturing to a desired surface roughness can be an extensive experimental process, this study will fundamentally investigate the effect of surface roughness on the thermal and velocity boundary layers for various surface roughness profiles with physics-based computational modeling. To accurately model the surface roughness effect on thermohydraulic performance, correlations like roughness ratio, friction factor, and Nusselt number are investigated. This study expects to indicate correlations for different flow regimes which can successfully predict the effect of surface roughness on thermohydraulic performance. Results generated from modeling surface roughness effect with unique correlations are then compared to experimental data obtained from literature for validation. The results of this study will provide insight to understand the thermal fluid physics associated with surface roughness and design guidelines for heat exchangers produced with additive manufacturing.

36 MATERIALS SCIENCE↗

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗

The vertical-velocity skewness in the atmospheric boundary layer without buoyancy and Coriolis effects

One of the main features of near-neutral atmospheric boundary layer (ABL) turbulence is the positive vertical velocity skewness $Sk_w$ above the roughness sublayer or the buffer region in smooth-walls. The $Sk_w$ variations are receiving renewed interest in many climate-related parameterizations of the ABL given their significance to cloud formation and to testing sub-grid schemes for Large Eddy Simulations (LES). The vertical variations of $Sk_w$ are explored here using wind tunnel and flume experiments collected above smooth, rough, and permeable-walls in the absence of buoyancy and Coriolis effects. These laboratory experiments form a necessary starting point to probe the canonical structure of $Sk_w$ as they deal with a key limiting case (i.e., near-neutral conditions). Diagnostic models based on cumulant expansions, realizability constraints, and constant mass flux approach routinely employed in the convective boundary layer as well as prognostic models based on third-order budgets are used to explain variations in $Sk_w$ for the idealized laboratory conditions. The failure of flux-gradient relations to model $Sk_w$ from the gradients of the vertical velocity variance σ$_w^2$ are explained and corrections based on models of energy transport offered. Novel links between the diagnostic and prognostic models are also featured, especially for the inertial term in the third-order budget of the vertical velocity fluctuation. The co-spectral properties of w′/σ w vs w′ 2 /σ$_w^2$ are also presented for the first time to assess the dominant scales governing $Sk_w$ in the inner and outer layers, where w′ is the fluctuating vertical velocity and σ w is the vertical velocity standard deviation.>

Boundary layer flow↗

Fluidization of Group A Glass Particles: Experiments and Preliminary Validation

This work is aimed at providing reliable high-quality data from fluidization experiments. Geldart Group A glass beads are used as bed material. The fluidizing medium is air at atmospheric pressure which enters the system at 40% relative humidity. The air inlet velocity is varied over a broad range from 2 U mf to 26 U mf , where the minimum fluidization velocity, U mf , is 0.608 cm/s. There is not a significant change in the mean values of differential pressure; however, there are noticeable trends in their standard deviation values. Preliminary validation studies are performed using the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). The results from Eulerian-Eulerian (MFiX-TFM, two-fluid model) and Eulerian-Lagrangian (MFiX-PIC, particle-in-cell) formulations are compared. The work presented in this report highlights the applicability of MFiX-PIC for large-scale applications, and reveals that PIC simulation offers a suitable alternative to TFM simulation when effects due to multiple components or polydispersity are significant. PIC modeling provides a readily extendible framework for such systems. However, PIC models do rely on an empirical closure for inter-particle stress. Some solution discrepancies were evident in MFiX-PIC predictions which require further analysis. As such, the study underlines the need for a systematic approach to quantify the sources of numerical uncertainty within the PIC model. However, even with known empiricism, MFiX-PIC has demonstrated considerable potential in analyzing gas-solid systems. The results obtained are encouraging and warrant further investigation to improve the predictive capability of MFiX-PIC.

36 MATERIALS SCIENCE↗

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗

Wind power production from very large offshore wind farms

In this work, we provide the first quantitative assessment of power production and wake generation from offshore wind energy lease areas along the U.S. east coast. Deploying 15-MW wind turbines, with spacing equal to the European average, yields electricity production of 116 TWh/year or 3% of current national supply. However, power production is reduced by one-third due to wakes caused by upwind wind turbines and wind farms. Under some flow conditions whole wind-farm wakes can extend up to 90 km downwind of the largest lease areas, and the frequency-weighted average area with a 5% velocity deficit is 2.6 times the footprint of the lease areas. Simulations including maritime corridors demonstrate reduction in the wake effects leading to power-efficiency gains and may offer contingent benefits. First-order scaling rules are developed that describe how “wake shadows” from large offshore wind farms scale with prevailing meteorology and wind turbine installed densities.

15 MW↗

Hydrodynamics of two-dimensional compressible fluid with broken parity: Variational principle and free surface dynamics in the absence of dissipation

In this paper, we consider an isotropic compressible nondissipative fluid with broken parity subject to free surface boundary conditions in two spatial dimensions. The hydrodynamic equations describing the bulk dynamics of the fluid and the free surface boundary conditions depend explicitly on the parity-breaking nondissipative odd viscosity term. We construct an effective action which gives both bulk hydrodynamic equations and free surface boundary conditions. The free surface boundary conditions require an additional boundary term in the action which resembles a 1+1D chiral boson field coupled to the background geometry. We solve the linearized hydrodynamic equations for the deep water case and derive the dispersion of chiral surface waves. We show that in the long-wavelength limit the flow profile exhibits an oscillating vortical boundary layer near the free surface. The layer thickness is controlled by the ratio between the odd viscosity (ν o ) and the sound velocity (c s ), δ ~ ν o /c s . In the incompressible limit, c s → ∞, the vortical boundary layer becomes singular with the vorticity within the layer diverging as ω ~ c s . The boundary layer is formed by odd viscosity coupling the divergence of velocity ∇ • v to vorticity ∇ × v. It results in nontrivial chiral free surface dynamics even in the absence of external forces. The structure of the odd-viscosity-induced boundary layer is very different from the conventional free surface boundary layer associated with dissipative shear viscosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Automated nuclear cloud feature extraction from film

Chemical, biological, radiological, nuclear, and explosives incidents require rapid detection and characterization for appropriate response. For a nuclear detonation, visible-light cameras may be used to locate the cloud and characterize fallout deposition when coupled with numerical models. Films from the United States’ nuclear testing era compose the only sizeable collection of imagery depicting high-yield detonations. These films offer unique insights into characteristics of flows involving scales that are difficult to replicate experimentally, and they are a valuable source of data for the validation of models for nuclear fallout transport, either as part of emergency response or forensic activities. In this work, we implement modern computer vision and machine learning techniques to identify and track the cloud automatically and subsequently determine the time dependence of some of its features. We trained a ResNet-18 image classifier on hundreds of images to categorize nuclear cloud morphology. Each category or cloud regime is determined by early cloud evolution and is associated to constitutive properties of the flow, such as distribution of vorticity. Next, we identified keypoint features using the KAZE algorithm and tracked these keypoints in the images, allowing us to determine the dimensions and velocities of the cloud across film frames. These measurements converted to real-world units provide valuable experimental data that can be used in the development and validation of nuclear cloud models. We compared the results of this method against manual cloud rise measurements from two different films. In one, our automated method accelerated the feature extraction process without sacrificing measurement accuracy.

Khristy, Joel [ORNL] (ORCID:0000000209963060)↗