Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WIND 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 19 records

Vertical profiles of the 3-D wind velocity retrieved from multiple wind lidars performing triple range-height-indicator scans

Abstract. Vertical profiles of 3-D wind velocity are retrieved from triple range-height-indicator (RHI) scans performed with multiple simultaneous scanning Doppler wind lidars. This test is part of the eXperimental Planetary boundary layer Instrumentation Assessment (XPIA) campaign carried out at the Boulder Atmospheric Observatory. The three wind velocity components are retrieved and then compared with the data acquired through various profiling wind lidars and high-frequency wind data obtained from sonic anemometers installed on a 300 m meteorological tower. The results show that the magnitude of the horizontal wind velocity and the wind direction obtained from the triple RHI scans are generally retrieved with good accuracy. However, poor accuracy is obtained for the evaluation of the vertical velocity, which is mainly due to its typically smaller magnitude and to the error propagation connected with the data retrieval procedure and accuracy in the experimental setup.

47 OTHER INSTRUMENTATION↗

A STATISTICAL ANALYSIS OF 50 YEARS OF DAILY SOLAR WIND VELOCITY DATA

Daily mean value of solar wind velocity measured by various spacecraft near the Earth’s orbit from 1963 November 27 to 2013 December 31 given by OMNIWeb is utilized to investigate its characteristics through statistical analyses. The percent probability distribution of solar wind velocity can be described well by the Γ distribution function with the most probable velocity to be 373 km s{sup -1}. It is found that solar wind could be statistically classified into three groups: (1) the low-velocity wind, v<450 km s{sup -1}, which positively responds to, and in the cycle phase lags, the solar activity cycle; (2) the high-velocity wind, 450≤v<725 km s{sup -1}, which negatively responds to, and in the cycle phase leads, the solar activity cycle, but in a short timescale (one-rotation-period) lags the solar activity cycle; and (3) the extreme-high-velocity wind, v≥725 km s{sup -1}, which positively responds to, and in cycle phase leads, the solar activity cycle. A period of about 27 days is determined for solar wind in the first two groups, so that solar wind in the groups is modulated by the solar rotation and related with solar long-life magnetic structures. Solar wind of extreme high velocity appears mainly at the descending phases of the solar cycles.

79 ASTRONOMY AND ASTROPHYSICS↗

A physics-informed data-driven low order model for the wind velocity deficit at the wake of isolated buildings

Detailed simulations of the Atmospheric boundary layer (ABL) provides information on wind velocity, direction and turbulence for optimizing the design of wind turbines and wind farms and their operations. Distributed wind turbines have been considered for use in urban or suburban areas for covering energy needs using a sustainable energy solution. Simulations of high Reynolds number turbulent ABL flows past obstacles can be used for the prediction of the velocity deficit and turbulence characteristics in the leeward side of building-like obstacles, allowing for the estimation of power losses and fatigue loads in wind turbines. However, fully three-dimensional simulations of high Reynolds number ABL flows using Large Eddy Simulation (LES) or Reynolds Average Navier-Stokes Simulations (RANS) continue to require large computational/ execution times. As a result, these 3D simulations may not be ideal tools for the rapid assessment of wind energy potential. In the present study, 3D numerical results from wall resolved RANS simulations were used for the development of a low order model that allows the rapid assessment of wind velocity and turbulence in the leeward side of an isolated building.

17 WIND ENERGY↗

Deep Learning Parameterization of Vertical Wind Velocity Variability via Constrained Adversarial Training

Atmospheric models with typical resolution in the tenths of kilometers cannot resolve the dynamics of air parcel ascent, which varies on scales ranging from tens to hundreds of meters. Small-scale wind fluctuations are thus characterized by a subgrid distribution of vertical wind velocity W with standard deviation σ W . The parameterization of σ W is fundamental to the representation of aerosol–cloud interactions, yet it is poorly constrained. Using a novel deep learning technique, this work develops a new parameterization for σ W merging data from global storm-resolving model simulations, high-frequency retrievals of W , and climate reanalysis products. The parameterization reproduces the observed statistics of σ W and leverages learned physical relations from the model simulations to guide extrapolation beyond the observed domain. Incorporating observational data during the training phase was found to be critical for its performance. The parameterization can be applied online within large-scale atmospheric models, or offline using output from weather forecasting and reanalysis products.

54 ENVIRONMENTAL SCIENCES↗

Site E - Virtual Tower High-Frequency Wind Velocity / Derived Data

This dataset contains Dual-Doppler wind reconstruction from stare lidar scans two rotor diameters south of the King Plains G02 turbine. The line-of-sight velocity data can be found in sa5.lidar.z03.b0 and sa1.lidar.z03.b0. The reconstruction is done neglecting the contribution from vertical velocity.

17 WIND ENERGY↗

Anorthite sputtering by H + and Ar q+ (q = 1-9) at solar wind velocities

Here, we report sputtering measurements of anorthite-like material, taken to be representative of soils found in the lunar highlands, impacted by singly and multicharged ions representative of the solar wind. The ions investigated include protons, as well as singly and multicharged Ar ions (as proxies for the nonreactive heavy solar wind constituents), in the charge state range +1 to +9, at fixed solar wind-relevant impact velocities of 165 and 310 km/s (0.25 keV/amu and 0.5 keV/amu). A quartz microbalance approach (QCM) for determination of total sputtering yields was used. The goal of the measurements was to determine the sputtering contribution of the heavy, multicharged minority solar wind constituents in comparison to that due to the dominant H + fraction. The QCM results show a yield increase of a factor of about 80 for Ar + versus H + sputtering and an enhancement by a factor of 1.67 between Ar 9+ and Ar + , which is a clear indication of a potential sputtering effect.

36 MATERIALS SCIENCE↗

Electronic excitation and charge transfer processes in collisions of H{sup +}, H{sub 2}{sup +}, and H{sub 3}{sup +} ions with carbon monoxide at typical solar-wind velocities

Luminescence in the 200-580 nm spectral region was observed in the collisions of H{sup +}, H{sub 2}{sup +}, and H{sub 3}{sup +} with CO in the 50-1000 eV projectile energy range. Using computer simulations, we have identified emission of the following products in the observed spectra: the CO{sup +}(A-X) comet-tail system, CO{sup +}(B-X) first negative system, CO{sup +}(B-A) Baldet-Johnson system, and CO(b-a) third positive system. Also, an emission from atomic hydrogen (H{sub β} line at 486nm) has been observed. From the analysis of the experimental spectra, we have determined the absolute emission cross-sections for the formation of the observed products. Computer simulations gave the excited-product population distributions over vibrational and rotational energy levels. The vibrational level distribution from the CO{sup +}(A-X) comet-tail system is compared with the data for CO excited by 100 eV electrons and extreme ultraviolet radiation (XUV) photons. We have used these data to analyze the excitation conditions in the comet Humason (1961e). From the vibrational population distributions observed in the comet, we found that this distribution can be reproduced if electrons produce 25%, protons 70%, and XUV photons produce 5% of the emitting molecules. We find that the ratio of the CO{sup +}(B-X) emission to the sum of two main emissions (CO{sup +}(A-X)+CO{sup +}(B-X)) is velocity dependent and does not depend on the projectile ion type. For small velocities (below 100 km s{sup –1}) the ratio is about 5%, while for higher velocities it increases to 30%. For these data, we have found an empirical formula that satisfactorily describes the experimental data: R = R {sub max}(1 – v {sub th}/v), (where R {sub max} = 33%, v {sub th} = 87 km s{sup –1}). This could be used to infer the velocity of ions producing the observed emission of CO{sup +} products.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of Cloud Base Height, Wind Velocity, and Short-Range Cloud Structure Using Multiple Sky Imagers Field Campaign Report

Clouds are a central focus of the U.S. Department of Energy (DOE)’s Atmospheric System Research (ASR) program and Atmospheric Radiation Measurement (ARM) Climate Research Facility, and more broadly are the subject of much investigation because of their important effects on atmospheric radiation and, through feedbacks, on climate sensitivity. Significant progress has been made by moving from a vertically pointing (“soda-straw”) to a three-dimensional (3D) view of clouds by investing in scanning cloud radars through the American Recovery and Reinvestment Act of 2009. Yet, because of the physical nature of radars, there are key gaps in ARM's cloud observational capabilities. For example, cloud radars often fail to detect small shallow cumulus and thin cirrus clouds that are nonetheless radiatively important. Furthermore, it takes five to twenty minutes for a cloud radar to complete a 3D volume scan and clouds can evolve substantially during this period. Ground-based stereo-imaging is a promising technique to complement existing ARM cloud observation capabilities. It enables the estimation of cloud coverage, height, horizontal motion, morphology, and spatial arrangement over an extended area of up to 30 by 30 km at refresh rates greater than 1 Hz (Peng et al. 2015). With fine spatial and temporal resolution of modern sky cameras, the stereo-imaging technique allows for the tracking of a small cumulus cloud or a thin cirrus cloud that cannot be detected by a cloud radar. With support from the DOE SunShot Initiative, the Principal Investigator (PI)’s team at Brookhaven National Laboratory (BNL) has developed some initial capability for cloud tracking using multiple distinctly located hemispheric cameras (Peng et al. 2015). To validate the ground-based cloud stereo-imaging technique, the cloud stereo-imaging field campaign was conducted at the ARM Facility’s Southern Great Plains (SGP) site in Oklahoma from July 15 to December 24. As shown in Figure 1, the cloud stereo-imaging system consisted of two inexpensive high-definition (HD) hemispheric cameras (each cost less than $1,500) and ARM’s Total Sky Imager (TSI). Together with other co-located ARM instrumentation, the campaign provides a promising opportunity to validate stereo-imaging-based cloud base height and, more importantly, to examine the feasibility of cloud thickness retrieval for low-view-angle clouds.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K2O, LaBr3, KBr, Cu, Zr, F, Ga, Ga2O3, NaNO2, Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated assuming 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradioactive chemical

54 ENVIRONMENTAL SCIENCES↗

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradradioactive chemicals modeled. And finally, the calculations assume all 12 tests will be performed at the same place at both locations, and all 11 radioactive material types will be used for each test. This is conservative because it is anticipated that no more than two material types will be used per test.

99 GENERAL AND MISCELLANEOUS↗

Introduction to and comparison of deep learning and optimization approaches to analytical wake modeling of a tilted wind turbine

This paper introduces innovative optimization and deep learning techniques to enhance the prediction of complex wake dynamics in the downstream wind velocity of tilted wind turbines. Traditional methods for calibrating the Bastankhah wake model often lead to increased errors in wind velocity distribution due to overfitting of the local wake characteristics. To address this issue, we propose an additional global optimization step to reduce errors in wind velocity predictions with respect to various wake parameters. Despite this improvement, the Bastankhah model's axisymmetric Gaussian wake shape limits its accuracy for complex wake structures. Therefore, we also propose a deep learning approach, which demonstrates promising results by accurately modeling complex wake shapes across a broader range of tilt angles with minimal computational cost. The deep learning approach achieves near-identical predictions to high-fidelity large-eddy simulations, representing a promising advancement in wake modeling.

17 WIND ENERGY↗

Machine Learning for Daily Forecasts of Arctic Sea Ice Motion: An Attribution Assessment of Model Predictive Skill

Physics-based simulations of Arctic sea ice are highly complex, involving transport between different phases, length scales, and time scales. Resultantly, numerical simulations of sea ice dynamics have a high computational cost and model uncertainty. We employ data-driven machine learning (ML) to make predictions of sea ice motion. The ML models are built to predict present-day sea ice velocity given present-day wind velocity and previous-day sea ice concentration and velocity. Models are trained using reanalysis winds and satellite-derived sea ice properties. We compare the predictions of three different models: persistence (PS), linear regression (LR), and a convolutional neural network (CNN). We quantify the spatiotemporal variability of the correlation between observations and the statistical model predictions. Additionally, we analyze model performance in comparison to variability in properties related to ice motion (wind velocity, ice velocity, ice concentration, distance from coast, bathymetric depth) to understand the processes related to decreases in model performance. Results indicate that a CNN makes skillful predictions of daily sea ice velocity with a correlation up to 0.81 between predicted and observed sea ice velocity, while the LR and PS implementations exhibit correlations of 0.78 and 0.69, respectively. The correlation varies spatially and seasonally: lower values occur in shallow coastal regions and during times of minimum sea ice extent. LR parameter analysis indicates that wind velocity plays the largest role in predicting sea ice velocity on 1-day time scales, particularly in the central Arctic. Regions where wind velocity has the largest LR parameter are regions where the CNN has higher predictive skill than the LR.

54 ENVIRONMENTAL SCIENCES↗

Online Learning of Effective Turbine Wind Speed in Wind Farms

To develop better wind farm controllers that can meet more complex objectives, methods of modeling the wind turbine wakes at low computational expense are needed. Gaussian process (GP) regression offers a computationally inexpensive framework for learning complex functions from noisy measurements with very few datapoints. In this work, an online learning approach is presented to learn the rotor-averaged wind velocity at downstream wind turbines with GPs, using the available datastream of wind field measurements and wind turbine control set-points. This framework can readily be integrated into model-based controls methods because the model a) is updated online at low computational expense, b) assumes a mathematically favorable Gaussian form, and c) explicitly quantifies the stochastic nature of the wake field so that the trade-off between exploration and exploitation, and the uncertainty in the prediction, can be utilized. We show that a GP-learned model can match true values with errors within 0.5% on average, with as few as 5 training data points.

Gaussian process↗

Starburst-driven galactic superbubbles radiating to 10 K

Our three-dimensional hydrodynamical simulations of starbursts examine the formation of superbubbles over a range of driving luminosities and mass loadings that determine superbubble growth and wind velocity. From this we determine the relationship between the velocity of a galactic wind (GW) and the power of the starburst. We find a threshold for the formation of a wind, above which the speed of the wind is not affected by grid resolution or the temperature floor of our radiative cooling. We investigate the effect that two different temperature floors in our radiative cooling prescription have on wind kinematics and content. We find that cooling to 10 K instead of to 10{sup 4} K increases the mass fraction of cold neutral and hot X-ray gas in the GW, while halving that in warm Hα. Our simulations show that the mass of cold gas transported into the lower halo does not depend on the starburst strength. Optically bright filaments form at the edge of merging superbubbles, or where a cold dense cloud has been disrupted by the wind. Filaments formed by merging superbubbles will persist and grow to >400 pc in length if anchored to a star forming complex. Filaments embedded in the hot GW contain warm and cold gas that moves 300−1200 km s{sup −1} slower than the surrounding wind, with the coldest gas hardly moving with respect to the Galaxy. Warm and cold matter in the GW show asymmetric absorption profiles consistent with observations, with a thin tail up to the wind velocity.

79 ASTRONOMY AND ASTROPHYSICS↗

Consideration of the Impact of the Environmental Conditions when Designing Heat-Receiving Systems of the Solar Cylindrical Parabolic Modules

The aerohydrodynamic and heat exchange characteristics in the cylindrical heat receiver of a solar cylindrical parabolic concentrating module are investigated. The mathematical model is based on the Navier- Stokes equation system for a laminar flow in a channel. The model is implemented using the control volume method. The created computer algorithm allows carrying out parametric studies and revealing the most rational geometric, dynamic, and power parameters of the process. The proposed model describes the general physical outlines of the heat exchange in a cylindrical heat receiver which allows estimating the influence of the environment—wind velocity, ambient temperature, etc.—on it.

14 SOLAR ENERGY↗

Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

79 ASTRONOMY AND ASTROPHYSICS↗

Extended main sequence turnoffs in intermediate-age star clusters: a correlation between turnoff width and early escape velocity

We present a color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1-2 Gyr old) star clusters in the Magellanic Clouds, including eight clusters for which new data were obtained. We find that all star clusters in our sample feature extended main-sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simple stellar population (including unresolved binary stars). FWHM widths of the MSTOs indicate age spreads of 200-550 Myr. We evaluate the dynamical evolution of clusters with and without initial mass segregation. Our main results are (1) the fraction of red clump (RC) stars in secondary RCs in eMSTO clusters scales with the fraction of MSTO stars having pseudo-ages of ≲1.35 Gyr; (2) the width of the pseudo-age distributions of eMSTO clusters is correlated with their central escape velocity v {sub esc}, both currently and at an age of 10 Myr. We find that these two results are unlikely to be reproduced by the effects of interactive binary stars or a range of stellar rotation velocities. We therefore argue that the eMSTO phenomenon is mainly caused by extended star formation within the clusters; and (3) we find that v {sub esc} ≥ 15 km s{sup –1} out to ages of at least 100 Myr for all clusters featuring eMSTOs, and v {sub esc} ≤ 12 km s{sup –1} at all ages for two lower-mass clusters in the same age range that do not show eMSTOs. We argue that eMSTOs only occur for clusters whose early escape velocities are higher than the wind velocities of stars that provide material from which second-generation stars can form. The threshold of 12-15 km s{sup –1} is consistent with wind velocities of intermediate-mass asymptotic giant branch stars and massive binary stars in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗