Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortex flows”

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 163 records · Page 9

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗

Experimental and computational investigation of interactive flow induced oscillations of two tandem rough cylinders at 3×10 4 ≤Re≤1.2×10 5

The interactive flow induced oscillations (FIO) of two adjacent, elastically mounted, rigid, tandem, locally-rough cylinders in transverse flow direction are analyzed utilizing two-dimensional Unsteady Reynolds-Averaged Navier-Stokes (2-D URANS) equations and verified experimentally in the proximity-wake region. Three sets of experiments and simulations (K = 600N/m, ζ = 0.14; K = 600N/m, ζ = 0.26; K = 1,200N/m, ζ = 0.26; K is spring stiffness, ζ is damping ratio) of two tandem cylinders with turbulence stimulation are tested and simulated for Reynolds number in the range of 30,000 = Re ≤ 120,000. The reduced velocity range is U* = 2.48–14.22, the mass ratio is m* = 1.343, and the center-to-center in-flow spacing to diameter ratio is d/D = 2.57. The characteristics of amplitude response, frequency response, lift force, and interactive wake patterns are presented and discussed. The trends of the amplitude and frequency responses from numerical simulations are in good agreement with experimental results. The main conclusions of Reynolds number effect on interactive flow induced oscillations are: (1) Five significant flow patterns between two tandem locally-rough cylinders for different Reynolds numbers are observed through analyzing the complex but stable interactions between vortices and cylinders. (2) In the initial and upper VIV branches, the downstream cylinder's FIO is seriously interfered by the wake of the upstream cylinder. (3) The downstream cylinder is strongly impinged by the vortices shed from the upstream cylinder resulting in nearly 180O out-of-phase oscillations in transition from VIV to galloping and in-phase oscillations in galloping.

42 ENGINEERING↗

Visualization and analysis of coupling between plasmas self-organization and plasma-induced fluid circulation in 1 atm DC glows with liquid anode

The physical processes prevailing when plasma contacts liquid water are poorly understood; however, it is this very interaction that is the basis of new, exciting technologies for water purification, the treatment of cancer and disease and the production of new, high value chemical products. A chief impediment to widespread application of plasmas for the aforementioned applications is scale up and cost. Here we investigate plasma induced flows which facilitate transport of reactivity from the plasma to liquid water at rates much faster than diffusion. We find that this plasma induced motion varies depending on whether or not the plasma is self organized. The insight from this work paves the way to a better understanding of transport and ultimately informing plasma treatment technologies on how to best optimize the rate at which processes such as plasma based water purification occurs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lift Equivalence and Cancellation for Airfoil Surge–Pitch–Plunge Oscillations

A NACA 0018 airfoil in freestream velocity is oscillated in longitudinal, transverse, and angle-of-attack directions with respect to the freestream velocity, known as surge, plunge, and pitch. The lift-based equivalence method introduces phase shifts between these three motions to construct in-phase sinusoidal components for maximum lift, waveform construction. Lift cancellation is also determined with the exact negative pitch and plunge motion amplitudes found from the equivalence method to achieve out-of-phase wave destruction. Lift cancellation occurs when a combination of these motions is sought to obtain a constant lift magnitude throughout the oscillation cycle. To achieve both equivalence and cancellation of lift, a prescribed pure pitch amplitude through the Theodorsen theory equates the corresponding equivalent plunge amplitude and pitch–plunge phase shift. These Theodorsen, linear superposition findings of pitch–plunge are leveraged toward the Greenberg theory to determine a closed-form, surge–pitch–plunge solution through the addition of a surge–plunge phase shift and optimal surge amplitude for lift cancellation. The lift cancellation surge–pitch–plunge amplitudes define the equivalence amplitude investigated here and theoretically limit the experiment to combinations of the first lift harmonic of the Greenberg theory. The analytical results are then compared with experimental lift force measurements and dye visualization. Finally, the normalized lift differences due to unsteady wake and boundary-layer behavior are examined to explore the extents of the Greenberg theory for these cases of lift-based equivalence and cancellation.

42 ENGINEERING↗

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗

The Large-Scale Vorticity Balance of the Antarctic Continental Margin in a Fine-Resolution Global Simulation

The depth-integrated vorticity budget of a global, eddy-permitting ocean/sea ice simulation over the Antarctic continental margin (ACM) is diagnosed to understand the physical mechanisms implicated in meridional transport. The leading-order balance is between the torques due to lateral friction, nonlinear effects, and bottom vortex stretching, although details vary regionally. Maps of the time-averaged depth-integrated vorticity budget terms and time series of the spatially averaged, depth-integrated vorticity budget terms reveal that the flow in the Amundsen, Bellingshausen, and Weddell Seas and, to a lesser extent, in the western portion of East Antarctica, is closer to an approximate topographic Sverdrup balance (TSB) compared to other segments of the ACM. Correlation and coherence analyses further support these findings, and also show that inclusion of the vorticity tendency term in the response (the planetary vorticity advection and the bottom vortex stretching term) increases the correlation with the forcing (the vertical net stress curl), and also increases the coherence between forcing and response at high frequencies across the ACM, except for the West Antarctic Peninsula. Finally, these findings suggest that the surface stress curl, imparted by the wind and the sea ice, has the potential to contribute to the meridional, approximately cross-slope, transport to a greater extent in the Amundsen, Bellingshausen, Weddell, and part of the East Antarctic continental margin than elsewhere in the ACM.

58 GEOSCIENCES↗

Thermal Hydraulic Experimental Test Article (FY2020 Status Report)

The Thermal Hydraulic Experimental Test Article (THETA) is a facility that will be used to develop sodium components and instrumentation as well as acquire experimental data for validation of reactor thermal hydraulic and safety analysis codes. The facility will simulate nominal conditions as well as protected/unprotected loss of flow accidents in a sodium-cooled fast reactor (SFR). High fidelity distributed temperature profiles of the developed flow field will be acquired with Rayleigh backscatter based optical fiber temperature sensors. The facility is being designed in partnership with systems code experts to tailor the experiment to ensure the most relevant and highest quality data for code validation. THETA is comprised of a traditional primary and secondary system. The primary system is submerged in the pool of sodium and consists of a pump, electrically heated core, intermediate heat exchanger, and connected piping and thermal barriers (redan). The secondary system, located outside of the sodium pool, consists of a pump, sodium to air heat exchanger, and connected piping and valves. Figure 1 illustrates the main components of the primary system. THETA will be installed in the Mechanisms Engineering Test Loop (METL) with the primary system in the 28 inch Test Vessel #4, Figure 2. The design work for the facility has been completed. Since the FY19 THETA status report, the THETA primary system components have been assembled and tested using de-ionized water. Some results of this testing at select Richardson and Reynolds numbers are included herein. The custom THETA sodium flowmeter (permanent magnet submersible sodium flowmeter) has been built and calibrated using sodium in an auxiliary facility against a NIST traceable vortex shedder type flowmeter. The permanent magnet flowmeter showed excellent performance, following the theoretical equation for measured voltage quite well, with an error < 3.5% at flow rates ranging from 3-24 GPM at temperatures up to 400 °C. The secondary sodium system for heat rejection from intermediate heat exchanger has been analyzed according to ASME code. All secondary side piping, fittings and a U-stamped air-to-sodium heat exchanger have been procured.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experimental and numerical study of the shielding effect of two tandem rough cylinders in flow-induce oscillation

The shielding effect of the downstream cylinder in flow induced oscillation (FIO) of two cylinders arranged in tandem is studied experimentally and numerically at Reynolds number 30,000 to 120,000. Both cylinders are in one degree-of-freedom, transverse-oscillations, and have turbulence stimulation in the form of selective surface roughness to expand FIO beyond vortex-induced vibration (VIV) into galloping. Shielding of the downstream cylinder has a negative effect on harnessing hydrokinetic energy. Further, to study its effect and mechanics, selective cases are studied both numerically and experimentally and discussed to demonstrate the shielding effect on the downstream cylinder and understand its cause. The main conclusions are: (1) The shielding effect for the downstream cylinder shows a strong relation to the damping ratio. As the damping ratio increases, the shielding effect is mitigated. Additionally, the oscillation of the rear cylinder becomes stable and shows stable frequency. (2) In the VIV region, as the stiffness and natural frequency increase, the shielding effect decreases substantially. (3) In the VIV region, the vorticity of the vortices shedding from both the upper and the lower sides of the downstream cylinder does not accumulate enough due to the attraction by the vortices shed from the upstream cylinder, thus resulting in partial suppression of the oscillation on the downstream cylinder. (4) In the galloping region, the shielding effect for the downstream cylinder depends on whether the vorticity near the downstream cylinder is strengthened by the vortices generated by the shear layers of the upstream cylinder or weakened.

42 ENGINEERING↗

Three-Dimensional Aerodynamics and Vortex-Shedding Characteristics of Wind Turbine Airfoils over 360-Degree Angles of Attack

In this work, we present the first three-dimensional (3D) computational investigation of wind turbine airfoils over 360° angles of attack to predict unsteady aerodynamic loads and vortex-shedding characteristics. To this end, static–airfoil simulations are performed for the FFA-W3 airfoil family at a Reynolds number of 107 with the Improved Delayed Detached Eddy Simulation turbulence model. Aerodynamic forces reveal that the onset of boundary-layer instabilities and flow separation does not necessarily coincide with the onset of stall. In addition, a comparison with two-dimensional simulation data and flat plate theory extension of airfoil polars, suggest that, in the deep stall regime, 3D effects remain critical for predicting both the unsteady loads and the vortex-shedding dynamics. For all airfoils, the vortex-shedding frequencies are found to be inversely proportional to the wake width. In the case of slender airfoils, the frequencies are nearly independent of the airfoil thickness, and their corresponding Strouhal number St is approximately 0.15. Based on the calculated St, the potential for shedding frequencies to coincide with the natural frequencies of the International Energy Agency 15 MW reference wind turbine blades is investigated. The analysis shows that vortex-induced vibrations occur primarily at angles of attack of around ±90° for all airfoils.

17 WIND ENERGY↗

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Emergence and dynamics of unconfined self-organised vortices in active magnetic roller liquids

Actively driven colloids demonstrate complex out-of-equilibrium dynamics often rivaling self-organized patterns and collective behavior observed in living systems. Recent studies revealed the emergence of steady macroscopic states with multiple interacting vortices in an unconfined environment that emerge from the coupling between microscale particle rotation and translation. Yet, insights into the microscopic behavior during the vortex emergence, growth, and formation of a multi-vortical state remain lacking. Here, we investigate in experiments and simulations how the microscale magnetic roller behavior leads to the emergence of seed vortices, their aggregation or annihilation, and the formation of stable large-scale vortical structures. We reveal that the coupling of roller-induced hydrodynamic flows guides the local self-densifications and self-organization of the micro-rollers into seed vortices. Here, the resulting multi-vortical state is sensitive to the external magnetic field amplitude and allows tuning the rollers' number density in a vortex and its characteristic size.

36 MATERIALS SCIENCE↗

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

Experimental measurements of fluid flow in an 84-pin hexagonal rod bundle with spacer grid for a gas-cooled fast modular reactor

A 50 MW e helium-cooled fast modular reactor (FMR) is under development, and as part of the Department of Energy Integrated Research Project (IRP), Texas A&M University is conducting the thermal-hydraulic characterization of its baseline core configuration. Here, we experimentally investigated the axial and cross flow fields characteristics in the hexagonal fuel rod bundle composed of 84 rods, a central rod, and spacer grids at a Reynolds number of 12,000. A fully transparent experimental facility resembling of one unit of the fuel assembly was constructed. Time-resolved particle image velocimetry (TR-PIV) measurements were performed to characterize hydraulic behavior downstream the spacer grid. Velocity measurements were conducted in the axial and radial direction of the rod bundle. From the PIV velocity vector fields, the full-field flow statistics were computed for the mean velocity, vorticity, and Reynolds stresses. An analysis of the energy decay downstream of the spacer grid was conducted with calculations of secondary-flow intensities, turbulent kinetic energy, power spectrum, and spatial-temporal velocity cross correlations. The vorticity field was obtained and pairs of counter-rotating vortexes were identified using a 3D reconstruction of several measurement planes. The data from this experimental campaign will be used to validate numerical models based on Computational Fluid Dynamics and other codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

First Observations of an Ion Vortex in a Magnetic Hole in the Solar Wind by MMS

The coronal heating region is able to generate mirror mode structures by ion mirror instabilities. Linear magnetic holes are believed to be the remnants of mirror mode structures, thus they are believed to be messengers from the coronal heating region. They can be convected to ∼9 au with the solar wind flow, indicating that a stabilizing mechanism is necessary to make the magnetic holes survive for such a long time. Here, we investigate a magnetic hole with a size of ∼6.7 ρ {sub i} in the solar wind based on observations by the Magnetospheric Multiscale mission. The unprecedented high-resolution data enable us to reveal the existence of the ion vortex inside the structure for the first time. Such an ion vortex forms a ring-like current, which is consistent with the magnetic field depression. The self-consistent structure of the magnetic hole contributed by the ion vortex can help to further shed light on the mechanism of the long-term survival of magnetic holes in the astrophysical plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comprehensive review of nonlinear oscillators in hydrokinetic energy harnessing using flow-induced vibrations

A comprehensive review of hydrokinetic energy converters based on alternating lift technology (ALT) is provided. Emphasis is on nonlinear oscillators based on Flow Induced Vibration (FIV) or Oscillation (FIO). Due to strong coupling in Fluid-Structure Interaction (FSI), and in order to maximize the hydrokinetic harnessed energy, design of nonlinear oscillators and analysis by model tests or computational fluid dynamics dominates this area. Research confirmed that the nonlinear oscillator can harvest energy from a stochastic excitation modeled by a generic wide spectrum, and overcome the most severe oscillator limitations: specifically, the need for continuous frequency tuning due to the narrow bandwidth response, and low efficiency outside the narrow bandwidth oscillator response. This review covers the following aspects of nonlinear oscillators in ALT converters: (1) Geometric changes in oscillator cross-section; e.g., circular, square, rectangular, or trilateral shapes. (2) Passive turbulence control of FIV/FIO. (3) Position based nonlinear stiffness. (4) Multi-cylinder synergistic FIV/FIO. (5) Mechanically linked oscillators. (6) Velocity-based, nonlinear, adaptive harnessing damping.

42 ENGINEERING↗

Anomalous transverse resistance in the topological superconductor β-Bi2Pd

Abstract A supercurrent flowing in a superconductor meets no resistance. Yet an electric field may still be established within the superconductor in the presence of dissipative processes, such as vortex motion. Here we report the observation of a transverse voltage drop in superconducting β- Bi 2 Pd thin films. Unlike the Hall effect in general or in other superconductors, the sign of the observed transverse voltage does not depend on the external magnetic field. Instead, it is dictated by the broken inversion symmetry on the film interfaces. This anomalous transverse voltage, or transverse resistance, is indicative of a chirality that likely resonates with the topological surface states reported in β- Bi 2 Pd.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modelling of a Flow-Induced Oscillation, Two-Cylinder, Hydrokinetic Energy Converter Based on Experimental Data

The VIVACE Converter consists of cylindrical oscillators in tandem subjected to transverse flow-induced oscillations (FIOs) that can be improved by varying the system parameters for a given in-flow velocity: damping, stiffness, and in-flow center-to-center spacing. Compared to a single isolated cylinder, tandem cylinders can harness more hydrokinetic energy due to synergy in FIO. Experimental and numerical methods have been utilized to analyze the FIO and energy harnessing of VIVACE. A surrogate-based model of two tandem cylinders is developed to predict the power harvesting and corresponding efficiency by introducing a backpropagation neural network. It is then utilized to reduce excessive experimental or computational testing. The effects of spacing, damping, and stiffness on harvested power and efficiency of the established prediction-model are analyzed. At each selected flow velocity, optimization results of power harvesting using the prediction-model are calculated under different combinations of damping and stiffness. The main conclusions are: (1) The surrogate model, built on extensive experimental data for tandem cylinders, can predict the cylinder oscillatory response accurately. (2) Increasing the damping ratio range from 0–0.24 to 0–0.30 is beneficial for improving power efficiency, but has no significant effect on power harvesting. (3) In galloping, a spacing ratio of 1.57 has the highest optimal harnessed power and efficiency compared with other spacing values. (4) Two tandem cylinders can harness 2.01–4.67 times the optimal power of an isolated cylinder. In addition, the former can achieve 1.46–4.01 times the efficiency of the latter. (5) The surrogate model is an efficient predictive tool defining parameters of the Converter for improved energy acquisition.

backpropagation neural network↗

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗