Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large-scale simulation”

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 1,009 records · Page 56

Large-Scale Transport Responses to Tropospheric Circulation Changes Using GEOS-5

The mean age since air was last at the Northern Hemisphere midlatitude surface is a fundamental property of tropospheric transport. Recent comparisons among chemistry climate models, however, reveal that there are large differences in the mean age among models and that these differences are most likely related to differences in tropical (parameterized) convection. Here we use aquaplanet simulations of the Goddard Earth Observing System Model Version 5 (GEOS-5) to explore the sensitivity of the mean age to changes in the tropical circulation. Tropical circulation changes are forced by prescribed localized off-equatorial warm sea surface temperature anomalies that (qualitatively) reproduce the convection and circulation differences among the comprehensive models. Idealized chemical species subject to prescribed OH loss are also integrated in parallel in order to illustrate the impact of tropical transport changes on interhemispheric constituent transport.

Orbe, Clara↗

Development of Experimental Icing Simulation Capability for Full-Scale Swept Wings: Hybrid Design Process, Years 1 and 2

This report presents the key results from the first two years of a program to develop experimental icing simulation capabilities for full-scale swept wings. This investigation was undertaken as a part of a larger collaborative research effort on ice accretion and aerodynamics for large-scale swept wings. Ice accretion and the resulting aerodynamic effect on large-scale swept wings presents a significant airplane design and certification challenge to air frame manufacturers, certification authorities, and research organizations alike. While the effect of ice accretion on straight wings has been studied in detail for many years, the available data on swept-wing icing are much more limited, especially for larger scales.

Aircraft icing↗

Two-thirds of Global Cropland Area Impacted by Climate Oscillations

The El Nino Southern Oscillation (ENSO) peaked strongly during the boreal winter 2015-2016, leading to food insecurity in many parts of Africa, Asia and Latin America. Besides ENSO, the Indian Ocean Dipole (IOD) and the North Atlantic Oscillation (NAO) are known to impact crop yields worldwide. Here we assess for the first time in a unified framework the relationship between ENSO, IOD and NAO and simulated crop productivity at the sub-country scale. Our findings reveal that during 1961-2010, crop productivity is significantly influenced by at least one large-scale climate oscillation in two-thirds of global cropland area. Besides observing new possible links, especially for NAO in Africa and the Middle East, our analyses confirm several known relationships between crop productivity and these oscillations. Our results improve the understanding of climatological crop productivity drivers, which is essential for enhancing food security in many of the most vulnerable places on the planet.

IOD↗

Measurements to Understand the Flow Mechanisms Contributing to Tandem-Rotor Outwash

Downwash and outwash characteristics of a model-scale tandem-rotor system in the presence of the ground were analyzed by identifying and understanding the physical mechanisms contributing to the observed flow field behavior. A building block approach was followed in simplifying the problem, separating the effects of the fuselage, effects of one rotor on the other, etc. Flow field velocities were acquired in a vertical plane at four aircraft azimuths of a small-scale tandem rotor system using the particle image velocimetry (PIV) technique for radial distances up to 4 times the rotor diameter. Results were compared against full-scale CH-47D measurements. Excellent correlation was found between the small- and full-scale mean flow fields (after appropriate normalization using rotor and wall jet parameters). Following the scalability analysis, the effect of rotor height on the outwash was also studied. Close to the aircraft, an increase in rotor height above ground decreased the outwash velocity at all aircraft azimuths. However, farther away, the longitudinal and lateral axes of the aircraft showed increasing and decreasing outwash velocities, respectively, with increasing rotor height. Measurements also indicated the presence of large-scale (of the size of the rotor height) shear-layer vortical structures along the ground that could be the source of low-frequency (approximately 1 Hz) flow variation observed in the full-scale measurements. Flow visualization studies and PIV measurements were also made on jets of different sizes to complement the observations made on rotors wherever possible. Baseline rotor measurements were made out-of-ground effect to understand the nature of inflow distribution for realistic rotor configurations and their modified characteristics in the presence of ground. Lastly, a feasibility study on applying high-fidelity CFD simulations for outwash study was conducted using Helios to model an isolated rotor configuration IGE at full scale Reynolds number. The results were encouraging and demonstrated the practical challenges associated with predicting rotor outwash.

Measurements↗

Measurements to Understand the Flow Mechanisms Contributing to Tandem-Rotor Outwash

Downwash and outwash characteristics of a model-scale tandem-rotor system in the presence of the ground were analyzed by identifying and understanding the physical mechanisms contributing to the observed flow field behavior. A building block approach was followed in simplifying the problem, separating the effects of the fuselage, effects of one rotor on the other, etc. Flow field velocities were acquired in a vertical plane at four aircraft azimuths of a small-scale tandem rotor system using the particle image velocimetry (PIV) technique for radial distances up to 4 times the rotor diameter. Results were compared against full-scale CH-47D measurements. Excellent correlation was found between the small- and full-scale mean flow fields (after appropriate normalization using rotor and wall jet parameters). Following the scalability analysis, the effect of rotor height on the outwash was also studied. Close to the aircraft, an increase in rotor height above ground decreased the outwash velocity at all aircraft azimuths. However, farther away, the longitudinal and lateral axes of the aircraft showed increasing and decreasing outwash velocities, respectively, with increasing rotor height. Measurements also indicated the presence of large-scale (of the size of the rotor height) shear-layer vortical structures along the ground that could be the source of low-frequency (approximately 1 Hz) flow variation observed in the full-scale measurements. Flow visualization studies and PIV measurements were also made on jets of different sizes to complement the observations made on rotors wherever possible. Baseline rotor measurements were made out-of-ground effect to understand the nature of inflow distribution for realistic rotor configurations and their modified characteristics in the presence of ground. Lastly, a feasibility study on applying high-fidelity CFD simulations for outwash study was conducted using Helios to model an isolated rotor configuration IGE at full-scale Reynolds number. The results were encouraging and demonstrated the practical challenges associated with predicting rotor outwash.

Ramasamy, Manikandan↗

A Decade of MWC 758 Disk Images: Where are the Spiral-Arm-Driving Planets?

Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and corotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10 year timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 L'-band observations in 2017. MWC 758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of 0.°6(sup +3.°3, sub -0.°6)/yr at 3σ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best-fit orbit is at 89 au (0."59), with a 3σ lower limit of 30 au (0."20). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.

MWC 758↗

An Evaluation of Structural Analysis Methodologies for Space Deployable Structures

Benchmarks are introduced for evaluating the performance of numerical simulations of space deployable structures. These benchmarks embody the key challenges of interest to future large space deployable structures, including large angle motion, contact between flexible bodies, and the presence of both soft and stiff mechanical components. The benchmarks were used in companion studies to evaluate the ADAMS multibody dynamics code, the LS-Dyna nonlinear finite element code, and the Sierra large-scale parallel nonlinear finite element code. In the past, only multibody codes would have been considered for this application. This study found that all three codes could be used for these benchmarks, a finding that may lead to larger scale, higher fidelity simulations in the future.

Mobrem, Mehran↗

Northern Late Winter Planetary Waves: MRO/MARCI Observations and Mars Climate Model Simulations

As does Earth, Mars presents pronounced global atmospheric circulation patterns. Solar differential heating drives mean meridional overturning (Hadley) circulations which are deep and intense, are hemispherically asymmetric, and where a cross-equatorial single cell dominates. Within middle and high latitudes, thermally indirect eddy-driven (Ferrel) circulation cells have been indicated. Differently, however, large-amplitude orography on planetary and continental scales on Mars can force very non-Earth-like hemispheric circulation patterns. Recent observations from the Mars Reconnaissance Orbiter, "Mars Color Imager" (MARCI) instrument are utilized that emphasize water ice clouds in ultra-violet (UV) wavelengths, and these measurements have been binned into "daily global maps" (DGMs) of water-ice cloud optical depth. The presence of large-scale, extratropical quasi-stationary atmospheric wave disturbances in middle and late winter of the northern hemisphere have been found to be present in such DGMs. In combination with such observations, a full-physics Mars global climate model (NASA ARC marsgcm 2.1) is applied to place the observations into context. During late northern winter, it is found that strong, forced Rossby modes (i.e., planetary waves) exist, and with direct correlation to columnintegrated cloud opacity undulating spatial patterns. At this season, zonal wavenumber s = 2 dominates (in contrast to wavenumber s = 1), consistent with MGS/TES analyses at this particular season (Banfield et al., 2003). Large-scale, planetary waves dictate the "coherence" of the northern polar vortex. Fundamentally, such forced planetary waves influence the polar vortex's impermeability (wave-induced) to tracer transport (e.g., dust and water-ice aerosol) and temporal mean water vapor spatial variations. The large-scale dynamical features of such planetary waves will be highlighted and discussed.

Northern polar vortex↗

Development of Large-Scale Magnetic Calorimeter Arrays

We describe performance of large-scale arrays of metallic magnetic calorimeters (MMCs) we are developing to meet requirements of the Lynx X-ray Microcalorimeter (LXM) instrument in the astrophysics mission concept Lynx. We have fabricated prototypes with 55,800 x-ray pixels thermally connected to 5,688 MMC sensors. Subarrays demonstrate three types of pixels, which have different energy and spatial resolution goals in LXM. Pixel pitch is 50 or 25 µm. For two subarrays, use of position sensitive detectors, in a "hydra" configuration with 25 x-ray absorbers per sensor, helps achieve a large focal plane by increasing pixel count relative to readout channels. Since each absorber has a different thermal link to its sensor, it generates a different pulse shape and enables discrimination of pixel position. Superconducting wiring from all sensors was brought out to the perimeter of the overall array using multiple buried metal layers planarized by Chemical Mechanical Planarization to achieve high critical current, low inductance, and high fabrication yield. An automated, algorithmic approach was used to layout the complex wiring pattern. For readout with existing small arrays of Superconducting Quantum Interference (SQUID) amplifiers, 112 selected pixels were connected to wire bond pads. In hydras, we successfully identified 25 pulse shapes using rise-time and pulse height. Pulse shapes were similar to simulations. We measured noise and responsivity to substrate temperature and applied x-ray pulses. While we designed sensors to have high inductance to relative to stray inductance of wiring, we had only mismatched (lower input inductance) SQUIDs available in our current apparatus, giving non-optimized resolution. We compared measured and calculated energy resolutions. We have the required sensitivity to reach < 3 eV FWHM resolution for 6 keV x-rays, required for the LXM Main Array, and the ability to discriminate the 25 different absorber pixels down to energies < 300 eV.

microcalorimeters↗

Fully-Coupled Fluid-Structure Interaction Simulations of a Supersonic Parachute

A validated computational fluid-structure interaction method for simulating the complex interaction between the large deformation of very thin, highly deformable structures and compressible flows is extended to consider large-scale problems in supersonic flows using parallel computing. The coupled fluid-structure interaction system is solved in a partitioned, or weakly-coupled, manner. The foundations of the applied fluid-structure interaction method are a higher-order, block-structured Cartesian, sharp immersed boundary method for the compressible Navier-Stokes equations and a computational structural dynamics solver employing a geometrically nonlinear 3-node shell element based on the mixed interpolation of tensorial components formulation. The method is applied to large deformation fluid-structure interaction validation cases before being applied to the inflation of a supersonic parachute in the upper Martian atmosphere where the goal is to demonstrate the capabilities of the solver when considering large-scale problems in supersonic flows.

Boustani, Jonathan↗

Influence of Sea-Ice Anomalies on Antarctic Precipitation Using Source Attribution in the Community Earth System Model

We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50∘ S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr−1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr−1). The contrast in contribution from the Southern Ocean, 102 Gt yr−1, is even more significant compared to the interannual variability of 35 Gt yr−1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.

Wang, Hailong↗

Fit2Fly: A Proof of Concept for Testing the Commercial Feasibility of Unmanned Aerial System Operations

As autonomous Unmanned Aerial Systems (UAS) become more prevalent, improvements in airworthiness assessments, data security, establishing vehicle trust, and more will become necessary to make fleet operations of these systems routine. Currently, these are verified for every UAS vehicle through a time-intensive manual inspection performed by human individuals. This process becomes infeasible when introduced into large-scale operations, so the Fit2Fly project is attempting to address these concerns. Fit2Fly is a multi-year project with a large scope, so we chose to narrow our focus to building a proof-of-concept feasibility simulation in tandem with physical UAS demonstration. This work resulted in a software simulation capable of running asynchronous commercial operations with persistent data modeling in coordination with several autonomous drones using integrated radio and GPS configurations.

Ryan Bonk↗

Model Design and Pre-Test CFD Analysis for a Supersonic Retropropulsion Wind Tunnel Test

Future Mars human lander missions will require using powered descent beginning at supersonic conditions, something which has never been done before at Mars. Significant aerosciences challenges exist due to interactions between the retrorocket exhaust plumes and the freestream flow that alter the aerodynamic behavior of the powered descent vehicle. Historically, wind tunnel tests have been used to study supersonic retropropulsion with inert gas exhaust simulants. Also, SpaceX has successfully decelerated the Falcon 9 first stage numerous times by using powered flight at supersonic conditions. On the computational side, Reynolds-Averaged Navier-Stokes flowfield simulations are regularly completed at full-scale conditions. However, the available ground and flight data do not provide a basis for calibrating the computational uncertainties for aerodynamic interference forces and moments on proposed Mars descent vehicles, either because of insufficient data or dissimilar vehicles geometries and/or conditions. Thus, additional ground testing is needed to continue addressing the aerosciences risks for large-scale human Mars landers. To that end, a retropropulsion wind tunnel test will be conducted in the NASA Langley Unitary Plan Wind Tunnel in 2020. The test is designed to improve upon similar past tests, both in terms of model design and measured data. The test campaign will use subscale model geometries derived from the two current NASA powered descent reference vehicles: a blunt low lift-to-drag vehicle and a more slender geometry. This paper covers the reference vehicles, test objectives, model design, scaling parameters, test matrix, and computational analysis to date.

Karl T Edquist↗

Lightweight Thrust Chamber Assemblies using Multi-Alloy Additive Manufacturing and Composite Overwrap

Additive Manufacturing (AM) has brought significant design and fabrication opportunities for complex components with internal features such as liquid rocket engine thrust chambers not previously possible. This technology allows for significant cost savings and schedule reductions in addition to new performance optimization through weight reduction and increased margins. Specific to regeneratively-cooled combustion chambers and nozzles for liquid rocket engines, additive manufacturing offers the ability to form the complex internal coolant channels and the closeout of the channels to contain the high pressure liquid propellants with a single operation. Much of additive manufacturing development has focused on monolithic alloys using Laser Powder Bed Fusion (L-PBF), which do not allow for complete optimization of the structure. The National Aeronautics and Space Administration (NASA) completed feasibility of an AM bimetallic L-PBF GRCop-84 copper-alloy combustion chamber with an AM electron beam freeform Inconel 625 structural jacket under the Low Cost Upper Stage Propulsion (LCUSP) Project. A follow-on project called Rapid Analysis and Manufacturing Propulsion Technology (RAMPT) is under development to further expand large-scale multi-alloy thrust chambers while maturing composite overwrap technology for significant weight savings opportunities. The RAMPT project has three primary objectives: 1) Advancing blown powder Directed Energy Deposition (DED) to fabricate integral-channel large scale nozzles, 2) Develop composite overwrap technology to reduce weight and provide structural capability for thrust chamber assemblies, and 3) Develop bimetallic and multi-metallic additively manufactured radial and axial joints to optimize material performance. In addition to these primary manufacturing developments, analytical modeling efforts compliment the process development to simulate the AM processes to reduce build failures and distortions. The RAMPT project is also maturing the supply chain for various manufacturing processes described above in addition to L-PBF of GRCop-42. This paper will present an overview of the RAMPT project, the process development and hardware progress to date, material and hot-fire testing results, along with future developments.

Additive Manufacturing↗

Impact of Dust-Cloud-Radiation-Precipitation Dynamical Feedback on Subseasonal-to-Seasonal Variability of the Asian Summer Monsoon in Global Variable-Resolution Simulations With MPAS-CAM5

In this study, we investigate the effects of increased dust emission from the Middle East deserts on subseasonal-to-seasonal (S2S) variability of the Asian summer monsoon (ASM). Numerical experiments are performed using the Model for Prediction Across Scales (MPAS) coupled with the Community Atmosphere Model (CAM5) physics, with regional refinement at 30 km grid resolution over South Asia and the surrounding regions. Result shows that increased dust emission and transport from the Middle East/West Asia region induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a colder surface over the desert regions and the western Tibetan Plateau, but warmer and moister troposphere with enhanced cloudiness and precipitation over the Pakistan/Northwest India (PNWI) region. The latter changes are amplified by the dust aerosol induced Elevated Heat Pump (EHP) mechanism along the West Himalayas/Iranian Plateau foothill regions, most pronounced during May- June. During July-August, cloud radiation feedback further enhances the warming of the upper troposphere, and cooling of the land surface over the PNWI and adjacent regions overWest Asia. The upper tropospheric heating and increased precipitation over PNWI spur a large-scale anomalous Rossby wavetrain and a northward displacement of the subtropical jetstream, manifesting in a contraction of the South Asian High and westward shift of the Western Pacific Subtropical High. As a result, the entire ASM precipitation-cloud system is displaced westward. Precipitation and cloudiness are intensified over northwest and western India and west Asia, but suppressed over southern and central East Asia. Analyses of the S2S variability of the upper level vorticity balance suggests that heating by Middle East dust plays an important role in exciting, and anchoring a teleconnection pattern through interactions among dust cloud radiation, precipitation heating over the PNWI, the development of an upper level Rossby wavetrain, and the northward shift of the boreal summer jetstream over Eurasia.

Asian summer monsoon↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Comparison of SAR and CYNSS surface water extent metrics over the Yucatan Lake wetland site

Wetlands have a major role in the carbon cycle, outgassing large quantities of carbon dioxide and methane through processes that are directly and strongly influenced by the duration and timing of inundation. Therefore, understanding the seasonal pattern of inundation can be a component for regional to global scale carbon models. Measurement of inundation extent also establishes a benchmark for the current status of wetland areas, useful in assessing the future impacts of climate change. The incorporation of frequent measurements of inundation extent into large-scale hydrological models would permit the evaluation of more detailed seasonal and longer-term floodplain dynamics and their associated management implications.The Cyclone GNSS (CYGNSS) constellation of satellites launched in 2016, and carries receivers capable of receiving data from L-band GNSS reflections. Delay Doppler maps (DDM) are generated on board and telemetered to the ground, along with a small number of raw data takes that can be used for special processing on the ground for evaluation purposes. It has been previously shown that these data can be sensitive to inundation.The NASA ISRO Synthetic Aperture Radar, currently planned for launch in January 2023, has both an L-band and S-band SAR for earth imaging. The L-band SAR, which will image the Earth's land mass twice every 12 days, has a requirement for measuring wetland inundation extent at the 1 ha scale. One of the sites that will be used to validate this requirement is Yucatan Lake, Louisiana. This oxbow lake and surrounding area located adjacent to the Mississippi river experiences periodic and extensive flooding in the surrounding forest areas.In 2019, NASA's UAVSAR fully polarimetric airborne L-band SAR conducted a flight campaign to image a dozen sites in the SE USA at approximately 12-day intervals and both in the morning and evening, to simulate the type of data NISAR will obtain. One site imaged during this campaign was the Yucatan Lake area, spanning water stages from low to high flood conditions.It has been demonstrated previously that GNSS reflectometry such as that measured by CYGNSS may be used to characterize surface inundation. It has also been known for decades that L-band SAR may be used to characterize not only the presence of open water, but also the presence of subcanopy inundation in forested areas. In this paper we will present results comparing data from these two types of instruments.

Lavalle, Marco↗

On the Settling Depth of Meltwater Escaping From Beneath Antarctic Ice Shelves

Antarctic glacial meltwater is thought to play an important role in determining large-scale Southern Ocean climate trends, yet recent modeling efforts have proceeded without a good understanding of how its vertical distribution in the water column is set. To rectify this, here we conduct new large-eddy simulations of the ascent of a buoyant meltwater plume after its escape from beneath an Antarctic ice shelf. We find that the meltwater’s settling depth is primarily a function of the buoyancy forcing per unit width of the source and the ambient stratification, consistent with the classical theory of turbulent buoyant plumes and in contrast to previous work that suggested an important role for centrifugal instability. Our results further highlight the significant role played by localized variability in stratification; this helps explain observed interannual variability in the vertical meltwater distribution near Pine Island Glacier. Because of the vast heterogeneity in mass loss rates and ambient conditions at different Antarctic ice shelves, a dynamic parameterization of meltwater settling depth may be crucial for accurately simulating high-latitude climate in a warming world; we discuss how this may be developed following this work, and where the remaining challenges lie.

Antarctic glacial meltwater↗