Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modified separator”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Microstructure and Hard Magnetic Properties of Sm 1-x Zr x (Fe,Co) 11.3-y Ti 0.7 By Ingots and Thick Melt-Spun Ribbons

Permanent magnets made from Sm(Fe,Co) 12 -based compounds are being actively pursued through nanostructuring and powder metallurgy. This study was aimed at the development of hard magnetic properties in bulk as-cast alloys and in melt-spun alloys for very low wheel speeds. Slower solidification rates and alloying with Zr promote the tetragonal ThMn 12 -type crystal structure, whereas higher solidification rates and alloying with B replace the ThMn 12 structure type with the TbCu 7 structure type. When introduced simultaneously, Zr and B dramatically reduce the alloy solidification rates required for both the refinement of the 1:12 crystallites and their replacement with the 1:7 phase. In bulk arc-melted alloys, this allowed for a microstructure of separated 1:12 crystallites 1–3 μm in size, although, because of the ferromagnetic nature of a minority phase, the coercivity of these fine-grained alloys reached only 0.73 kOe. A moderately accelerated solidification further refined the 1:12 crystallites and increased the coercivity; a Sm 0.7 Zr 0.4 (Fe,Co) 10.8 Ti 0.7 B 0.5 alloy exhibited a coercivity of 1.5 kOe and a maximum energy product of 3.4 MGOe when it was melt-spun into a 0.26-mm-thick ribbon. A more rapid solidification suppressed the 1:12 phase and after annealing at 800–850 °C, the alloys modified with Zr and B developed reasonably high coercivity and maximum energy product even when melt-spun at a wheel speed of 6 m/s. For the above-mentioned alloy, these values were 4.1 kOe and 7.8 MGOe, respectively. Further, a similarly processed very-Sm-lean Sm 0.5 Zr 0.6 (Fe,Co) 10.6 Ti 0.7 B 0.7 alloy exhibited a remanence of 8.8 kG and an energy product of 7.4 MGOe.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Turbulence modeling of shock separated boundary-layer flows

Computations of transonic and hypersonic shock-separated boundary-layer flows using zero-equation (algebraic), one-equation (kinetic energy), and two-equation (kinetic energy plus length scale) turbulence eddy viscosity models are described and compared with measurements. The computations make use of a new Navier-Stokes computer algorithm that has reduced computing times by one to two orders of magnitude. The algorithm, and how the turbulence models are incorporated into it, are described. Results for the transonic flow show that the unmodified one-equation model is superior to the zero-equation model in skin-friction predictions. For the hypersonic flow, a highly modified one-equation model that accurately predicts surface pressure and heat transfer is described. Preliminary two-equation model results are also presented.

Coakley, T. J.↗
NASA NTRS1994

Design and evaluation of a Stochastic Optimal Feed-forward and Feedback Technology (SOFFT) flight control architecture

This paper describes the design and evaluation of a stochastic optimal feed-forward and feedback technology (SOFFT) control architecture with emphasis on the feed-forward controller design. The SOFFT approach allows the designer to independently design the feed-forward and feedback controllers to meet separate objectives and then integrate the two controllers. The feed-forward controller has been integrated with an existing high-angle-of-attack (high-alpha) feedback controller. The feed-forward controller includes a variable command model with parameters selected to satisfy level 1 flying qualities with a high-alpha adjustment to achieve desired agility guidelines, a nonlinear interpolation approach that scales entire matrices for approximation of the plant model, and equations for calculating feed-forward gains developed for perfect plant-model tracking. The SOFFT design was applied to a nonlinear batch simulation model of an F/A-18 aircraft modified for thrust vectoring. Simulation results show that agility guidelines are met and that the SOFFT controller filters undesired pilot-induced frequencies more effectively during a tracking task than a flight controller that has the same feedback control law but does not have the SOFFT feed-forward control.

Ostroff, Aaron J.↗
NASA NTRS2008

A Numerical Study of Anti-Vortex Film Cooling Designs at High Blowing Ratio

A concept for mitigating the adverse effects of jet vorticity and liftoff at high blowing ratios for turbine film cooling flows has been developed and studied at NASA Glenn Research Center. This "anti-vortex" film cooling concept proposes the addition of two branched holes from each primary hole in order to produce a vorticity counter to the detrimental kidney vortices from the main jet. These vortices typically entrain hot freestream gas and are associated with jet separation from the turbine blade surface. The anti-vortex design is unique in that it requires only easily machinable round holes, unlike shaped film cooling holes and other advanced concepts. The anti-vortex film cooling hole concept has been modeled computationally for a single row of 30deg angled holes on a flat surface using the 3D Navier-Stokes solver Glenn-HT. A modification of the anti-vortex concept whereby the branched holes exit adjacent to the main hole has been studied computationally for blowing ratios of 1.0 and 2.0 and at density ratios of 1.0 and 2.0. This modified concept was selected because it has shown the most promise in recent experimental studies. The computational results show that the modified design improves the film cooling effectiveness relative to the round hole baseline and previous anti-vortex cases, in confirmation of the experimental studies.

Heidmann, James D.↗
DOE OSTI2021

Solar cell contacts: quantifying the impact of interfacial layers on selectivity, recombination, charge transfer, and V oc

Interfacial layers (IFLs) are ubiquitous in solar cells, but their precise impact on carrier transfer rates and the relation of these rates to performance metrics and the concepts of selectivity and recombination is lacking. We report the use of a well-defined interdigitated back-contact (IBC) silicon solar cell to determine the precise role of such IFLs. We characterize the action of the common IFL spiro-OMeTAD by making it a third contact to the IBC cell. This architecture creates three solar cells within a single structure that, with numerical simulation, provide the exchange current densities (i.e., charge transfer rates) for electrons ( J 0n ) and holes ( J 0p ) and the quasi Fermi-level splitting in the absorber, which measures the balance of generation and recombination. Further, we describe the relation of V oc to contact recombination, the asymmetry in electron/hole collection rates at a single contact (contact selectivity), and the asymmetry in collection rates of the same carrier at separate contacts (carrier selectivity). Relative to bare gold, neat spiro-OMeTAD reduces J 0n and J 0p (their geometric mean ( J 0n J 0p ) 0.5 decreases by 10 4 ), decreasing contact recombination. Addition of the common dopant Li-TFSI and air increase J 0p / J 0n by 10 6 with little effect on ( J 0n J 0p ) 0.5 , increasing contact hole selectivity. The significant increase in V oc observed by introducing spiro-OMeTAD/Li-TFSI IFLs into the cells studied, however, is due to an increase in the carrier selectivity rather than the contact selectivity or recombination of the spiro-OMeTAD-modified contact. Operando measurements further show voltage-dependent changes in the J 0 s, demonstrating that spiro-OMeTAD contributes to hysteresis. Broader context : Solar cells operate by photogenerating excess charge carriers in an absorber material and, in competition with recombination, asymmetrically extracting them at so-called carrier-selective contacts, one that ideally collects only electrons and the other, only holes. Particularly in emerging photovoltaics such as perovskites, thin layers of organic semiconductors or related materials are introduced between the absorber and contact to improve power conversion efficiency. In terms of interfacial charge transfer, a prevailing view is that such interfacial layers improve performance by helping block the collection of the undesired carrier, considered a form of recombination. We use a novel platform to study the simultaneous impact of spiro-OMeTAD, a common interfacial layer, on the collection of electrons, the collection of holes, and on the recombination of electrons and holes. We quantitatively demonstrate that spiro-OMeTAD layers indeed passivate the contact toward recombination, but that larger improvements in the open-circuit voltage, a key cell metric, can come not only from this but also from increasing the asymmetry of the collection of electrons in the system. Further, operando measurements show transient changes in the properties of spiro-OMeTAD which suggest that it contributes to hysteresis phenomena commonly observed in perovskite and other solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2000

Calculating Free Energies Using Scaled-Force Molecular Dynamics Algorithm

One common objective of molecular simulations in chemistry and biology is to calculate the free energy difference between different states of the system of interest. Examples of problems that have such an objective are calculations of receptor-ligand or protein-drug interactions, associations of molecules in response to hydrophobic, and electrostatic interactions or partition of molecules between immiscible liquids. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), 'native' state. Perhaps the best example of such a problem is folding of proteins or short RNA molecules. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to 'quasi non-ergodicity', whereby a part of phase space is inaccessible on timescales of the simulation. A host of strategies has been developed to improve efficiency of sampling the phase space. For example, some Monte Carlo techniques involve large steps which move the system between low-energy regions in phase space without the need for sampling the configurations corresponding to energy barriers (J-walking). Most strategies, however, rely on modifying probabilities of sampling low and high-energy regions in phase space such that transitions between states of interest are encouraged. Perhaps the simplest implementation of this strategy is to increase the temperature of the system. This approach was successfully used to identify denaturation pathways in several proteins, but it is clearly not applicable to protein folding. It is also not a successful method for determining free energy differences. Finally, the approach is likely to fail for systems with co-existing phases, such as water-membrane systems, because it may lead to spontaneous mixing. A similar difficulty may be encountered in any method relying on global modifications of phase space.

Darve, Eric↗
NASA NTRS1991

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗
DOE OSTI2020

Effects of transition metals on the evolution of polymer-derived SiOC ceramics

This work focuses on catalytic graphitization from different metals, iron (Fe), cobalt (Co), and nickel (Ni), for inducing carbon growth in silicon oxycarbides (SiOCs) during the pyrolysis. Fe, Co, and Ni-modified SiOCs were synthesized by pyrolysis to 900, 1100, and 1300 °C respectively in Ar. The transition metals induced the formation of the corresponding metal silicides, β-SiC, and graphitic carbon with the catalytic activity in the order of Fe > Co > Ni, in agreement with the activation energy calculation based on the carbon types and amounts. Lateral growth of turbostratic carbon followed a 2D grain growth process and its point-like defect density decreased based on the catalytic order of the transition metals, with SiOC/Fe having the lowest defect density. Quantitative analysis of the XPS results with Rietveld refinement illustrated that the phase separation of SiOC is more dominant than local carbothermal reduction between SiO 2 and C in the SiOC/M (M = Ni, Co, or Fe) systems. Understanding the catalytic graphitization effect of Ni, Co, and Fe on polymer derived ceramics offers new strategies in increasing high-temperature phase amounts and thus creating novel materials for heterogeneous catalysis, magnetic, and other applications.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Recent Developments in Gun Operating Techniques at the NASA Ames Ballistic Ranges

This paper describes recent developments in gun operating techniques at the Ames ballistic range complex. This range complex has been in operation since the early 1960s. Behavior of sabots during separation and projectile-target impact phenomena have long been observed by means of short-duration flash X-rays: new versions allow operation in the lower-energy ("soft") X-ray range and have been found to be more effective than the earlier designs. The dynamics of sabot separation is investigated in some depth from X-ray photographs of sabots launched in the Ames 1.0 in and 1.5 in guns; the sabot separation dynamics appears to be in reasonably good agreement with standard aerodynamic theory. Certain sabot packages appear to suffer no erosion or plastic deformation on traversing the gun barrel, contrary to what would be expected. Gun erosion data from the Ames 0.5 in, 1.0 in, and 1.5 in guns is examined in detail and can be correlated with a particular non- dimensionalized powder mass parameter. The gun erosion increases very rapidly as this parameter is increased. Representative shapes of eroded gun barrels are given. Guided by a computational fluid dynamics (CFD) code, the operating conditions of the Ames 0.5 in and 1.5 in guns were modified. These changes involved: (1) reduction in the piston mass, powder mass and hydrogen fill pressure and (2) reduction in pump tube volume, while maintaining hydrogen mass. These changes resulted in muzzle velocity increases of 0.5-0.8 km/sec, achieved simultaneously with 30-50 percent reductions in gun erosion.

Bogdanoff, D. W.↗
DOE OSTI2024

Organic Matter Concentration and Composition in November 2021 and April 2022 from 12 Streams Impacted by the 2020 Holiday Farm Fire (v2)

This dataset represents results from a field study aiming to understand storm induced transport of pyrogenic materials to streams impacted by varying degrees of burn severity. Time series samples were collected at 5 sites within the McKenzie River Watershed (Oregon, USA) whose catchment were each completely engulfed by the 2020 Holiday Farm Fire. An additional 7 sites were sampled once during the storm. The samples were collected during storm events in November 2020, January 2021, November 2021, and April 2022. Samples were characterized for benezenepolycarboxylic acids (BPCA), ultra-high resolution mass spectrometry, dissolved organic carbon and optics (absorbance and fluorescence). Fourier-transform ion cyclotron resonance mass spectrometry (FTICR) and dissolved organic carbon data from the November 2020 (referred to as “EWEB_2020”) sampling can be found in a separate data package (doi: 10.15485/1869708). NOTE: The 2020 samples were run on FTICR-MS in two unique instances. The first run can be found in the previous data package (EWEB_2020). The second run is included in this data package. These samples were run for a second time so that the data were more directly interoperable with the other samples in this data package. We have not done any investigation into the differences/similarities between these datasets and the previously ran/published data in the other data package. This data package was originally published in November 2024. It was updated in April 2025 (v2; new and modified files). See the change history section below for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset contains (1) file-level metadata; (2) data dictionary; (3) data package readme; (4) metadata; (5) methods information; (6) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data; (7) excitation emission matrix (EEM) methods; and (8) a sub-folder with processed EEM data (9) benzene polycarboxylic acid (BPCA) concentration data; (10) Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) methods; and (11) folder of high-resolution characterization of organic matter via 12 Tesla FTICR-MS generated through the Environmental Molecular Sciences Laboratory (EMSL; https://www.pnnl.gov/environmental-molecular-sciences-laboratory). The EEMs sub-folder contains two additional folders; the Absorbance and Fluorescence folders which contain the processed EEMs absorbance and fluorescence data respectively. This package contains the following file types: csv, xml, pdf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elucidating Photoinduced Processes of Photosystem I Via Multidimensional Electronic and Vibrational Spectroscopies

This project was motivated by an overarching goal to elucidate the mechanism of energy and electron transfer that governs the efficient charge separation in photosystem I (PSI) complexes. PSI is a natural light harvesting complex that drives oxygenic photosynthesis in plants, algae, and cyanobacteria. It uses ~300 tightly packed chlorophylls (Chls) to absorb photons, transfer the excitation energy to the reaction center (RC), and generate a charge separated state with near unity quantum efficiency (QE). A better understanding of the mechanism of energy transfer and charge separation in PSI is required for understanding the high QE of natural light harvesting complexes, and it could lead to the further development of artificial photosynthetic systems for solar energy conversion and modification of light harvesting complexes to improve crop yields. We applied two-dimensional optical spectroscopies to different cyanobacterial photosystem I complexes, including PSI complexes that contain Chl f molecules, to map energy transfer pathways and gain insight into the efficient light harvesting of PSI. We used two-dimensional electronic spectroscopies (2DES) to map energy transfer in Chl a and Chl f containing PSI complexes. To investigate the Chl f PSI complexes, we modified our spectrometer to probe the lower energy states associated with Chl f molecules. We interpreted the 2DES spectra through global analysis procedures to generate maps of energy transfer. We also constructed a two-dimensional electronic vibrational (2DEV) spectrometer that will be used to investigate charge transfer transitions and dynamics within PSI complexes. Measurements were performed on model systems to establish general data analysis procedures for interpreting 2D spectra and gain insight into protein cofactor interactions.

14 SOLAR ENERGY↗
NASA NTRS2002

Method and system for training dynamic nonlinear adaptive filters which have embedded memory

Described herein is a method and system for training nonlinear adaptive filters (or neural networks) which have embedded memory. Such memory can arise in a multi-layer finite impulse response (FIR) architecture, or an infinite impulse response (IIR) architecture. We focus on filter architectures with separate linear dynamic components and static nonlinear components. Such filters can be structured so as to restrict their degrees of computational freedom based on a priori knowledge about the dynamic operation to be emulated. The method is detailed for an FIR architecture which consists of linear FIR filters together with nonlinear generalized single layer subnets. For the IIR case, we extend the methodology to a general nonlinear architecture which uses feedback. For these dynamic architectures, we describe how one can apply optimization techniques which make updates closer to the Newton direction than those of a steepest descent method, such as backpropagation. We detail a novel adaptive modified Gauss-Newton optimization technique, which uses an adaptive learning rate to determine both the magnitude and direction of update steps. For a wide range of adaptive filtering applications, the new training algorithm converges faster and to a smaller value of cost than both steepest-descent methods such as backpropagation-through-time, and standard quasi-Newton methods. We apply the algorithm to modeling the inverse of a nonlinear dynamic tracking system 5, as well as a nonlinear amplifier 6.

Rabinowitz, Matthew↗
NASA NTRS2012

HD 172555: Detection of 63 micrometers [OI] Emission in a Debris Disc

Context. HD 172555 is a young A7 star belonging to the Beta Pictoris Moving Group that harbours a debris disc. The Spitzer IRS spectrum of the source showed mid-IR features such as silicates and glassy silica species, indicating the presence of a warm dust component with small grains, which places HD 172555 among the small group of debris discs with such properties. The IRS spectrum also shows a possible emission of SiO gas. Aims. We aim to study the dust distribution in the circumstellar disc of HD 172555 and to asses the presence of gas in the debris disc. Methods. As part of the GASPS Open Time Key Programme, we obtained Herschel-PACS photometric and spectroscopic observations of the source. We analysed PACS observations of HD 172555 and modelled the Spectral Energy Distribution (SED) with a modified blackbody and the gas emission with a two-level population model with no collisional de-excitation. Results. We report for the first time the detection of [OI] atomic gas emission at 63.18 micrometers in the HD 172555 circumstellar disc.We detect excesses due to circumstellar dust toward HD 172555 in the three photometric bands of PACS (70, 100, and 160 m). We derive a large dust particle mass of (4.8 plus-minus 0.6)x10(exp -4) Mass compared to Earth and an atomic oxygen mass of 2.5x10(exp -2)R(exp 2) Mass compared to Earth, where R in AU is the separation between the star and the inner disc. Thus, most of the detected mass of the disc is in the gaseous phase.

Riviere-Marichalar, P.↗
NASA NTRS2012

The Nucleus of Comet 9P-Tempel 1: Shape and Geology from Two Flybys

The nucleus of comet Tempel 1 has been investigated at close range during two spacecraft missions separated by one comet orbit of the Sun, 5 1/2 years. The combined imaging covers 70% of the surface of this object which has a mean radius of 2.83 +/- 0.1 km. The surface can be divided into two terrain types: rough, pitted terrain and smoother regions of varying local topography. The rough surface has round depressions from resolution limits (10 m/pixel) up to 1 km across, spanning forms from crisp steep-walled pits, to subtle albedo rings, to topographic rings, with all ranges of morphologic gradation. Three gravitationally low regions of the comet have smoother terrain, parts of which appear to be deposits from minimally modified flows, with other parts likely to be heavily eroded portions of multiple layer piles. Changes observed between the two missions are primarily due to backwasting of scarps bounding one of these probable flow deposits. This style of erosion is also suggested by remnant mesa forms in other areas of smoother terrain. The two distinct terrains suggest either an evolutionary change in processes, topographically- controlled processes, or a continuing interaction of erosion and deposition.

geological↗
NASA NTRS2015

Estimates of Power Plant NOx Emissions and Lifetimes from OMI NO2 Satellite Retrievals

Isolated power plants with well characterized emissions serve as an ideal test case of methods to estimate emissions using satellite data. In this study we evaluate the Exponentially-Modified Gaussian (EMG) method and the box model method based on mass balance for estimating known NOx emissions from satellite retrievals made by the Ozone Monitoring Instrument (OMI). We consider 29 power plants in the USA which have large NOx plumes that do not overlap with other sources and which have emissions data from the Continuous Emission Monitoring System (CEMS). This enables us to identify constraints required by the methods, such as which wind data to use and how to calculate background values. We found that the lifetimes estimated by the methods are too short to be representative of the chemical lifetime. Instead, we introduce a separate lifetime parameter to account for the discrepancy between estimates using real data and those that theory would predict. In terms of emissions, the EMG method required averages from multiple years to give accurate results, whereas the box model method gave accurate results for individual ozone seasons.

Emissions inventory↗
DOE OSTI2026

Hydration size-dependent transport of ions across nanoporous graphene membranes

Nanoporous atomically-thin membranes are promising candidates for metal ion separations due to their chemical stability and high permeance and selectivity, but experimental evidence of the mechanisms responsible for ion-ion selectivity is sparse. Here, in this work, we measured the simultaneous diffusion of a dilute mixture of ten different Group I, Group II, and rare earth cations in a salt background across nanoporous graphene (NPG) membranes with sub-nanometer pores. The membranes exhibited ion-ion selectivity, including between similarly-charged ions. Cation transport was governed primarily by the hydrated ion size that was consistent with continuum models of ion diffusion, with additional influence of the ionic charge. Selectivity enhancement was achieved by modifying ion sizes using an ion-selective complexing agent. Our study provides evidence of the importance of size-sieving and electrostatic mechanisms governing ion transport across NPG, and allows for quantitative prediction of transport rates to guide future development of ion-selective atomically-thin membranes.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Modelling Greenland Outlet Glaciers

The objective of this project was to develop simple yet realistic models of Greenland outlet glaciers to better understand ongoing changes and to identify possible causes for these changes. Several approaches can be taken to evaluate the interaction between climate forcing and ice dynamics, and the consequent ice-sheet response, which may involve changes in flow style. To evaluate the icesheet response to mass-balance forcing, Van der Veen (Journal of Geophysical Research, in press) makes the assumption that this response can be considered a perturbation on the reference state and may be evaluated separately from how this reference state evolves over time. Mass-balance forcing has an immediate effect on the ice sheet. Initially, the rate of thickness change as compared to the reference state equals the perturbation in snowfall or ablation. If the forcing persists, the ice sheet responds dynamically, adjusting the rate at which ice is evacuated from the interior to the margins, to achieve a new equilibrium. For large ice sheets, this dynamic adjustment may last for thousands of years, with the magnitude of change decreasing steadily over time as a new equilibrium is approached. This response can be described using kinematic wave theory. This theory, modified to pertain to Greenland drainage basins, was used to evaluate possible ice-sheet responses to perturbations in surface mass balance. The reference state is defined based on measurements along the central flowline of Petermann Glacier in north-west Greenland, and perturbations on this state considered. The advantage of this approach is that the particulars of the dynamical flow regime need not be explicitly known but are incorporated through the parameterization of the reference ice flux or longitudinal velocity profile. The results of the kinematic wave model indicate that significant rates of thickness change can occur immediately after the prescribed change in surface mass balance but adjustments in flow rapidly diminish these rates to a few cm/yr at most. The time scale for adjustment is of the order of a thousand years or so.

vanderVeen, Cornelis↗