SEARCH · Engineering Papers

Results for “surface complexation modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Broken Ergodicity in MHD Turbulence in a Spherical Domain

Broken ergodicity (BE) occurs in Fourier method numerical simulations of ideal, homogeneous, incompressible magnetohydrodynamic (MHD) turbulence. Although naive statistical theory predicts that Fourier coefficients of fluid velocity and magnetic field are zero-mean random variables, numerical simulations clearly show that low-wave-number coefficients have non-zero mean values that can be very large compared to the associated standard deviation. In other words, large-scale coherent structure (i.e., broken ergodicity) in homogeneous MHD turbulence can spontaneously grow out of random initial conditions. Eigenanalysis of the modal covariance matrices in the probability density functions of ideal statistical theory leads to a theoretical explanation of observed BE in homogeneous MHD turbulence. Since dissipation is minimal at the largest scales, BE is also relevant for resistive magnetofluids, as evidenced in numerical simulations. Here, we move beyond model magnetofluids confined by periodic boxes to examine BE in rotating magnetofluids in spherical domains using spherical harmonic expansions along with suitable boundary conditions. We present theoretical results for 3-D and 2-D spherical models and also present computational results from dynamical simulations of 2-D MHD turbulence on a rotating spherical surface. MHD turbulence on a 2-D sphere is affected by Coriolus forces, while MHD turbulence on a 2-D plane is not, so that 2-D spherical models are a useful (and simpler) intermediate stage on the path to understanding the much more complex 3-D spherical case.

Shebalin, John V.↗

NASA NTRS2012

Complex Protostellar Chemistry

Two decades ago, our understanding chemistry in protostars was simple -- matter either fell into the central star or was trapped in planetary-scale objects. Some minor chemical changes might occur as the dust and gas fell inward, but such effects were overwhelmed by the much larger-scale processes that occurred even in bodies as small as asteroids. The chemistry that did occur in the nebula was relatively easy to model because the fall from the cold molecular cloud into the growing star was a one-way trip down a well-known temperature pressure gradient; the only free variable was time. However, just over 10 years ago it was suggested that some material could be processed in the inner nebula, flow outward, and become incorporated into comets. This outward flow was confirmed when the Stardust mission returned crystalline mineral fragments from Comet Wild 2 that must have been processed close to the Sun before they were incorporated into the comet. In this week's Science Express, Ciesla and Sandford demonstrate that even the outermost regions of the solar nebula can be a chemically active environment. Their finding could have consequences for the rest of the nebula. Our understanding of the chemistry in protostellar systems has made enormous progress over the last few decades, fueled by an increased awareness of the complex dynamics of these evolving energetic nebulae. We can no longer consider just the simple local environment to explain the composition of a planet, asteroid, or comet as was done in the past, but must now consider chemical processes that might take place within the nebula as a whole as well as the probability of transport and mixing the products of such reactions throughout the system. just as we now find it impossible to explain the complex chemistry of the terrestrial atmosphere without reference to detailed transport models that interconnect highly dissimilar chemical environments, so chemical models of protostars and of the solar nebula must eventually treat these environments as tightly coupled, interactive systems. The demonstration that the chemistry on the surfaces of outward-flowing, dynamically mixing icy grain surfaces both mimics the chemistry in cold cloud cores and strikes at the central assumption of the photochemical self-shielding model for oxygen isotopes in solar system solids only adds emphasis to this conclusion.

Nuth, Joseph A., III↗

NASA NTRS2018

Nozzle Plume/Shock Interaction Experimental and Computational Sonic Boom Analyses from the NASA Ames 9- by 7-Foot Supersonic Wind Tunnel

A wind tunnel test and a computational study were conducted to investigate the complex interactions between a supersonic nozzle plume and shock waves of differing strengths generated from various aft surfaces typical of supersonic aircraft. These analytically-defined aft surfaces were representative of horizontal tails of various sizes, and an aft deck. CFD simulations of many proposed model configurations allowed for assessments of the detailed flow interactions of components in close proximity to the nozzle, as well as assessments of the nozzle jet flow itself. The evaluation of the computational results for many candidate configurations guided the design of model components. The interactions of the waveforms from these surfaces with the jet exhaust plume can have significant adverse effects on the loudness of the sonic boom if the surfaces are not carefully integrated into an aircraft design. The greatest discrepancy in estimating sonic boom loudness for low-boom flight vehicles is currently in predicting the signatures from the aft part of an aircraft, including the interactions with the plume flow. The objectives of this test were to gain a better understanding of these interactions, and to provide a detailed experimental database from multiple sources for use as validation cases for CFD tool development. The subject test was run in the NASA Ames 9- by 7-Ft Supersonic Wind Tunnel in February 2016 at Mach numbers of 1.6 and 2.0, and was funded by the NASA Commercial Supersonics Technology (CST) Project. The nozzle flow was provided by high-pressure air (HPA) pumped through the model, and pressure signature data were acquired with the NASA 14-inch sonic boom pressure rail. The rail measured the locations of the shocks and expansions at various distances and off-track angles from the model. This enabled the impact of the nozzle plume/shock interactions on the near- and mid-field sonic boom pressure waveforms to be quantified. Schlieren images of the flow field around and behind the model were obtained with an RBOS (Retroreflective Background-Oriented Schlieren) technique to determine the origins of the shock and expansion waves, to identify the shape and boundaries of the plume, and to determine the changes in incoming and exiting waveforms within the plume. A total pressure rake was positioned closely behind the model nozzle in order to measure the total pressure profiles of the flow above, within, and below the nozzle exhaust. Model angles and positions in the tunnel were measured by photogrammetry using two cameras since the lack of a model force balance prevented the measurement of model deflections under load.Navier-Stokes computations using two different CFD codes were compared to the experimental sonic boom pressure signature data, and the rake total pressure data in the plume. A computational schlieren technique was used to compare the computed flow field with the RBOS images. The computational results were also used to complement the test data with flow field quantities that could not be measured, such as Mach number and pressure distributions to distinguish shock waves and expansion waves.

sonic boom↗

NASA NTRS2019

Satellite‐Detected Ocean Ecosystem Response to Volcanic Eruptions in the Subarctic Northeast Pacific Ocean

Volcanic eruptions in the Aleutian Archipelago during the summer of 2008 deposited large quantities of iron‐laden ash to the subarctic North Pacific Ocean. The surface ocean ecosystem response has been previously characterized using limited in situ and autonomous field measurements and numerical modeling, but basin‐scale satellite remote sensing has been limited to a simple description of chlorophyll variability. Here, we show that the ecosystem response is a complex combination of phytoplankton biomass and physiology that can be described with satellite ocean color diagnostics such as the chlorophyll‐to‐carbon biomass ratio (Chl:Cphyto) and chlorophyll fluorescence yield. Together, these quantities outline a more complete picture of ecological responses spanning unique signals of iron stress (and relief from), photo acclimation, changes in phytoplankton growth rate, increases in biomass, and timescales of decay for these processes.

T.K. Westberry↗

NASA NTRSDate not supplied

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗

NASA NTRS2003

Subgrid Scale Modeling in Solar Convection Simulations using the ASH Code

The turbulent solar convection zone has remained one of the most challenging and important subjects in physics. Understanding the complex dynamics in the solar con- vection zone is crucial for gaining insight into the solar dynamo problem. Many solar observatories have generated revealing data with great details of large scale motions in the solar convection zone. For example, a strong di erential rotation is observed: the angular rotation is observed to be faster at the equator than near the poles not only near the solar surface, but also deep in the convection zone. On the other hand, due to the wide range of dynamical scales of turbulence in the solar convection zone, both theory and simulation have limited success. Thus, cutting edge solar models and numerical simulations of the solar convection zone have focused more narrowly on a few key features of the solar convection zone, such as the time-averaged di erential rotation. For example, Brun & Toomre (2002) report computational finding of differential rotation in an anelastic model for solar convection. A critical shortcoming in this model is that the viscous dissipation is based on application of mixing length theory to stellar dynamics with some ad hoc parameter tuning. The goal of our work is to implement the subgrid scale model developed at CTR into the solar simulation code and examine how the differential rotation will be a affected as a result. Specifically, we implement a Smagorinsky-Lilly subgrid scale model into the ASH (anelastic spherical harmonic) code developed over the years by various authors. This paper is organized as follows. In x2 we briefly formulate the anelastic system that describes the solar convection. In x3 we formulate the Smagorinsky-Lilly subgrid scale model for unstably stratifed convection. We then present some preliminary results in x4, where we also provide some conclusions and future directions.

Young, Y.-N.↗

NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗

NASA NTRS1979

A study of the blown flap/jet flap analogy

A study of the blown flap/jet flap analogy has been undertaken. Analytical predictions were made using both improved lifting line and optimized vortex lattice models for the jet flap. Results were compared with experimental data for three propulsive lift systems; the jet augmented flap, the externally blown flap, and the upper surface blown flap. Force increments due to changes in geometry and jet parameters were well approximated in most cases, although the absolute values of the aerodynamic forces were usually underestimated. The relatively simple jet-flap models gave performance predictions of accuracy comparable to more complex analyses.

Hough, G. R.↗

NASA NTRS1979

NASCAP modelling of high-voltage power system interactions with space charged-particle environments

The NASA Charging Analyzer Program (NASCAP), an engineering tool capable of analyzing the impact of the charged particle environment on spacecraft surfaces and systems, is described. NASCAP is a quasi-static computational program which analyzes the charging of a 3-dimensional complex body as a function of time and system-generated voltages for given space environmental conditions. The material properties of the surfaces are taken into account; the surface potentials, low energy sheath, potential distribution in space and particle trajectories are calculated. An application of NASCAP to a simple space solar power station consisting of two 6 m by 18 m solar array wings surrounding a central body is presented. Each solar array wing is considered to be divided into three regions operating at 2000 volts. Results of NASCAP analysis of the system for a normal environment and a moderate geomagnetic substorm environment are discussed.

Stevens, N. J.↗

NASA NTRS1990

An ocean scatter propagation model for aeronautical satellite communication applications

In this paper an ocean scattering propagation model, developed for aircraft-to-satellite (aeronautical) applications, is described. The purpose of the propagation model is to characterize the behavior of sea reflected multipath as a function of physical propagation path parameters. An accurate validation against the theoretical far field solution for a perfectly conducting sinusoidal surface is provided. Simulation results for typical L band aeronautical applications with low complexity antennas are presented.

Moreland, K. W.↗

NASA NTRS1993

A revised catalog of CfA galaxy groups in the Virgo/Great Attractor flow field

A new identification of groups and clusters in the CfAl Catalog of Huchra, et al. (1983) is presented, using a percolation algorithm to identify density enhancements. The procedure differs from that of the original Geller and Huchra (1983; GH) catalog in several important respects; galaxy distances are calculated from the Virgo-Great Attractor flow model of Faber and Burnstein (1988), the adopted distance linkage criteria is only approx. 1/4 as large as in the Geller and Huchra catalog, the sky link relation is taken from Nolthenius and White (1987), correction for interstellar extinction is included, and 'by-hand' adjustments to group memberships are made in the complex regions of Virgo/Coma I/Ursa Major and Coma/A1367 (to allow for varying group velocity dispersions and to trim unphysical 'spider arms'). Since flow model distances are poorly determined in these same regions, available distances from the IR Tully-Fisher planetary nebula luminosity function and surface brightness resolution methods are adopted if possible.

Nolthenius, Richard↗

NASA NTRS1997

The solid state environment orchestrates embryonic development and tissue remodeling

Cell interactions with extracellular matrix and with other cells play critical roles in morphogenesis during development and in tissue homeostasis and remodeling throughout life. Extracellular matrix is information-rich, not only because it is comprised of multifunctional structural ligands for cell surface adhesion receptors, but also because it contains peptide signaling factors, and proteinases and their inhibitors. The functions of these groups of molecules are extensively interrelated. In this review, three primary cell culture models are described that focus on adhesion receptors and their roles in complex aspects of morphogenesis and remodeling: the regulation of proteinase expression by fibronectin and integrins in synovial fibroblasts; the regulation of osteoblast differentiation and survival by fibronectin, and the regulation of trophoblast differentiation and invasion by integrins, cadherins and immunoglobulin family adhesion receptors.

NASA Discipline Musculoskeletal↗

NASA NTRS2009

Computer Code for Nanostructure Simulation

Due to their small size, nanostructures can have stress and thermal gradients that are larger than any macroscopic analogue. These gradients can lead to specific regions that are susceptible to failure via processes such as plastic deformation by dislocation emission, chemical debonding, and interfacial alloying. A program has been developed that rigorously simulates and predicts optoelectronic properties of nanostructures of virtually any geometrical complexity and material composition. It can be used in simulations of energy level structure, wave functions, density of states of spatially configured phonon-coupled electrons, excitons in quantum dots, quantum rings, quantum ring complexes, and more. The code can be used to calculate stress distributions and thermal transport properties for a variety of nanostructures and interfaces, transport and scattering at nanoscale interfaces and surfaces under various stress states, and alloy compositional gradients. The code allows users to perform modeling of charge transport processes through quantum-dot (QD) arrays as functions of inter-dot distance, array order versus disorder, QD orientation, shape, size, and chemical composition for applications in photovoltaics and physical properties of QD-based biochemical sensors. The code can be used to study the hot exciton formation/relation dynamics in arrays of QDs of different shapes and sizes at different temperatures. It also can be used to understand the relation among the deposition parameters and inherent stresses, strain deformation, heat flow, and failure of nanostructures.

Filikhin, Igor↗

NASA NTRS2010

Magnetic Topology of Coronal Hole Linkages

In recent work, Antiochos and coworkers argued that the boundary between the open and closed field regions on the Sun can be extremely complex with narrow corridors of open ux connecting seemingly disconnected coronal holes from the main polar holes, and that these corridors may be the sources of the slow solar wind. We examine, in detail, the topology of such magnetic configurations using an analytical source surface model that allows for analysis of the eld with arbitrary resolution. Our analysis reveals three important new results: First, a coronal hole boundary can join stably to the separatrix boundary of a parasitic polarity region. Second, a single parasitic polarity region can produce multiple null points in the corona and, more important, separator lines connecting these points. Such topologies are extremely favorable for magnetic reconnection, because it can now occur over the entire length of the separators rather than being con ned to a small region around the nulls. Finally, the coronal holes are not connected by an open- eld corridor of finite width, but instead are linked by a singular line that coincides with the separatrix footprint of the parasitic polarity. We investigate how the topological features described above evolve in response to motion of the parasitic polarity region. The implications of our results for the sources of the slow solar wind and for coronal and heliospheric observations are discussed.

Titov, V. S.↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗

NASA NTRS2019

Satellite-Detected Ocean Ecosystem Response to Volcanic Eruption in the Subarctic Northeast Pacific Ocean

Volcanic eruptions in the Aleutian archipelago during the summer of 2008 deposited large quantities of iron-laden ash to the Subarctic North Pacific Ocean. The surface ocean ecosystem response has been previously characterized using limited in situ and autonomous field measurements and numerical modeling, but basin-scale satellite remote sensing has been limited to a simple description of chlorophyll variability. Here, we show that the ecosystem response is a complex combination of phytoplankton biomass and physiology that can be described with satellite ocean color diagnostics such as the chlorophyll to carbon biomass ratio (Chl:Cphyto) and chlorophyll fluorescence yield. Together, these quantities outline a more complete picture of ecological response spanning unique signals of iron stress (and relief from), photo acclimation, changes in phytoplankton growth rate, increases in biomass, and timescales of decay for these processes.

Ocean Ecosystem↗

NASA NTRS2022

Constraining Arctic Climate Projections of Wintertime Warming With Surface Turbulent Flux Observations and Representation of Surface-Atmosphere Coupling

The drivers of rapid Arctic climate change—record sea ice loss, warming SSTs, and a lengthening of the sea ice melt season—compel us to understand how this complex system operates and use this knowledge to enhance Arctic predictability. Changing energy flows sparked by sea ice decline, spotlight atmosphere-surface coupling processes as central to Arctic system function and its climate change response. Despite this, the representation of surface turbulent flux parameterizations in models has not kept pace with our understanding. The large uncertainty in Arctic climate change projections, the central role of atmosphere-surface coupling, and the large discrepancy in model representation of surface turbulent fluxes indicates that these processes may serve as useful observational constraints on projected Arctic climate change. This possibility requires an evaluation of surface turbulent fluxes and their sensitivity to controlling factors (surface-air temperature and moisture differences, sea ice, and winds) within contemporary climate models (here Coupled Model Intercomparison Project 6). The influence of individual controlling factors and their interactions is diagnosed using a multi-linear regression approach. This evaluation is done for four sea ice loss regimes, determined from observational sea ice loss trends, to control for the confounding effects of natural variability between models and observations. The comparisons between satellite- and model-derived surface turbulent fluxes illustrate that while models capture the general sensitivity of surface turbulent fluxes to declining sea ice and to surface-air gradients of temperature and moisture, substantial mean state biases exist. Specifically, the central Arctic is too weak of a heat sink to the winter atmosphere compared to observations, with implications to the simulated atmospheric circulation variability and thermodynamic profiles. Models were found to be about 50% more efficient at turning an air-sea temperature gradient anomaly into a sensible heat flux anomaly relative to observations. Further, the influence of sea ice concentration on the sensible heat flux is underestimated in models compared to observations. The opposite is found for the latent heat flux variability in models; where the latent heat flux is too sensitive to a sea ice concentration anomaly. Lastly, the results suggest that present-day trends in sea ice retreat regions may serve as suitable observational constraints of projected Arctic warming.

turbulent fluxes↗