Search NASA⌕ Search

Engineering topics

Craig, Anthony P.

Publications and source records attributed to Craig, Anthony P..

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗
DOE OSTI2024

Sensitivities of the West Greenland Current to Greenland Ice Sheet Meltwater in a Mesoscale Ocean/Sea Ice Model

Abstract Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic.

58 GEOSCIENCES↗
DOE OSTI2024

Evaluation of dynamical downscaling in a fully coupled regional earth system model

A set of decadal simulations has been completed and evaluated for gains using the Regional Arctic System Model (RASM) to dynamically downscale data from a global Earth system model and two atmospheric reanalyses. RASM is a fully coupled atmosphere–land–ocean–sea ice regional Earth system model. Nudging to the forcing data is applied to approximately the top half of the atmospheric domain. RASM simulations were also completed with a modification to the atmospheric physics for evaluating changes to the modeling system. The results show that for the top half of the atmosphere, the RASM simulations follow closely to that of the forcing data, regardless of the forcing data. The results for the lower half of the atmosphere, as well as the surface, show a clustering of atmospheric state and surface fluxes based on the modeling system. At all levels of the atmosphere the imprint of the weather from the forcing data is present as indicated in the pattern of the annual means. Biases, in comparison to reanalyses, are evident in the Earth system model forced simulations for the top half of the atmosphere but are not present in the lower atmosphere. This suggests that bias correction is not needed for fully coupled dynamical downscaling simulations. While the RASM simulations tended to go to the same mean state for the lower atmosphere, there are a differences in the variability and changes of weather patterns across the ensemble of simulations. These differences in the weather result in variances in the sea ice and oceanic states.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data from: Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

This data collection represents all UH8to2 model data used for the publication "Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model." This includes ocean and sea ice characteristics from the 2017-2020 UH8to2 model run. Further details about the model setup, forcing, and our interpretation of these data can be found in the publication.

58 GEOSCIENCES↗
DOE OSTI2023

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗
DOE OSTI2023

Causes and evolution of winter polynyas north of Greenland

During the 42-year period (1979–2020) of satellite measurements, four major winter (December–March) polynyas have been observed north of Greenland: one in December 1986 and three in the last decade, i.e., February of 2011, 2017, and 2018. The 2018 polynya was unparalleled in its magnitude and duration compared to the three previous events. Given the apparent recent increase in the occurrence of these extreme events, this study aims to examine their evolution and causality, in terms of forced versus natural variability. The limited weather station and remotely sensed sea ice data are analyzed combining with output from the fully coupled Regional Arctic System Model (RASM), including one hindcast and two ensemble simulations. We found that neither the accompanying anomalous warm surface air intrusion nor the ocean below had an impact (i.e., no significant ice melting) on the evolution of the observed winter open-water episodes in the region. Instead, the extreme atmospheric wind forcing resulted in greater sea ice deformation and transport offshore, accounting for the majority of sea ice loss in all four polynyas. Our analysis suggests that strong southerly winds (i.e., northward wind with speeds greater than 10 m s –1 ) blowing persistently over the study region for at least 2 d or more were required over the study region to mechanically redistribute some of the thickest Arctic sea ice out of the region and thus to create open-water areas (i.e., a latent heat polynya). To assess the role of internal variability versus external forcing of such events, we carried out and examined results from the two RASM ensembles dynamically downscaled with output from the Community Earth System Model (CESM) Decadal Prediction Large Ensemble (DPLE) simulations. Out of 100 winters in each of the two ensembles (initialized 30 years apart: one in December 1985 and another in December 2015), 17 and 16 winter polynyas were produced north of Greenland, respectively. The frequency of polynya occurrence had no apparent sensitivity to the initial sea ice thickness in the study area pointing to internal variability of atmospheric forcing as a dominant cause of winter polynyas north of Greenland. We assert that dynamical downscaling using a high-resolution regional climate model offers a robust tool for process-level examination in space and time, synthesis with limited observations, and probabilistic forecasts of Arctic events, such as the ones being investigated here and elsewhere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An evaluation of the E3SMv1 Arctic ocean and sea-ice regionally refined model

The Energy Exascale Earth System Model (E3SM) is a state-of-the-science Earth system model (ESM) with the ability to focus horizontal resolution of its multiple components in specific areas. Regionally refined global ESMs are motivated by the need to explicitly resolve, rather than parameterize, relevant physics within the regions of refined resolution, while offering significant computational cost savings relative to the respective cost of configurations with high-resolution (HR) everywhere on the globe. In this paper, we document results from the first Arctic regionally refined E3SM configuration for the ocean and sea-ice components (E3SM-Arctic-OSI), while employing data-based atmosphere, land, and hydrology components. Our aim is an improved representation of the Arctic coupled ocean and sea-ice state, its variability and trends, and the exchanges of mass and property fluxes between the Arctic and the sub-Arctic. We find that E3SM-Arctic-OSI increases the realism of simulated Arctic ocean and sea-ice conditions compared to a similar low-resolution E3SM simulation without the Arctic regional refinement in ocean and sea-ice components (E3SM-LR-OSI). In particular, exchanges through the main Arctic gateways are greatly improved with respect to E3SM-LR-OSI. Other aspects, such as the Arctic freshwater content variability and sea-ice trends, are also satisfactorily simulated. Yet, other features, such as the upper-ocean stratification and the sea-ice thickness distribution, need further improvements, involving either more advanced parameterizations, model tuning, or additional grid refinements. Overall, E3SM-Arctic-OSI offers an improved representation of the Arctic system relative to E3SM-LR-OSI, at a fraction (15 %) of the computational cost of comparable global high-resolution configurations, while permitting exchanges with the lower-latitude oceans that cannot be directly accounted for in Arctic regional models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

An experimental low Reynolds number comparison of a Wortmann FX67-K170 airfoil, a NACA 0012 airfoil and a NACA 64-210 airfoil in simulated heavy rain

Wind tunnel experiments were conducted on Wortmann FX67-K170, NACA 0012, and NACA 64-210 airfoils at rain rates of 1000 mm/hr and Reynolds numbers of 310,000 to compare the aerodynamic performance degradation of the airfoils and to attempt to identify the various mechanisms which affect performance in heavy rain conditions. Lift and drag were measured in dry and wet conditions, a variety of flow visualization techniques were employed, and a computational code which predicted airfoil boundary layer behavior was used. At low angles of attack, the lift degradation in wet conditions varied significantly between the airfoils. The Wortmann section had the greatest overall lift degradation and the NACA 64-210 airfoil had the smallest. At high angles of attack, the NACA 64-210 and 0012 airfoils had improved aerodynamic performance in rain conditions due to an apparent reduction of the boundry layer separation. Performance degradation in heavy rain for all three airfoils at low angles of attack could be emulated by forced boundary layer transition near the leading edge. The secondary effect occurs at time scales consistent with top surface water runback times. The runback layer is thought to effectively alter the airfoil geometry. The severity of the performance degradation for the airfoils varied. The relative differences appeared to be related to the susceptibility of each airfoil to premature boundary layer transition.

Craig, Anthony P.↗
NASA NTRS1987

Comparative low Reynolds number tests of NACA 64-210, NACA 0012, and Wortmann FX67-K170 airfoils in heavy rain

The effect of simulated rain at 1000 mm/h on the aerodynamic performance of NACA 64-210, NACA 0012, and Wortmann FX67-K170 airfoils is investigated experimentally at Reynolds number 310,000, freestream velocity 70 mph, and angle of attack -10 to +20 deg in the 1 x 1-ft section of the MIT low-turbulence wind tunnel. The results are presented graphically and characterized in detail. Rain is found to produce lift degradation of 5-25 percent, an effect attributed to development of a premature transition at the leading edge followed (after 1-10 s) by alteration of the profile by surface water runback.

Hansman, R. John, Jr.↗