Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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,657 records · Page 92

Relationship between the blue response and open-circuit voltage of high performance silicon solar cells

The relationship between the response at blue wavelengths and the open-circuit voltage of high performance silicon solar cells with a virtually 'transparent' emitter is analyzed. It is shown that a one-dimensional cell model cannot simultaneously model the optical collection properties and operating characteristics of such cells. The analysis highlights the importance of surface conditions and shows that previous conclusions regarding bounds on Auger coefficients imposed by correlating these parameters must be treated with caution.

Green, M. A.↗

Transient water stress in a vegetation canopy - Simulations and measurements

Consideration is given to observational and modeling evidence of transient water stress, the effects of the transpiration plateau on the canopy radiometric temperature, and the factors responsible for the onset of the transpiration plateau, such as soil moisture. Attention is also given to the point at which the transient stress can be detected by remote measurement of surface temperature.

Carlson, Toby N.↗

Aviation Fuel Tracer Simulation: Model Intercomparison and Implications

An upper limit for aircraft-produced perturbations to aerosols and gaseous exhaust products in the upper troposphere and lower stratosphere (UT/LS) is derived using the 1992 aviation fuel tracer simulation performed by eleven global atmospheric models. Key Endings are that subsonic aircraft emissions: (1) have not be responsible for the observed water vapor trends at 40 deg N; (2) could be a significant source of soot mass near 12 km, but not at 20 km; (3) might cause a noticeable increase in the background sulfate aerosol surface area and number densities (but not mass density) near the northern mid-latitude tropopause; and (4) could provide a global, annual mean top of the atmosphere radiative forcing up to +0.006 W/sq m and -0.013 W/sq m due to emitted soot and sulfur, respectively.

Danilin, M. Y.↗

An Approach to Solving Enclosure Radiation Problems in A Multi-Physics Context

Thermal protection system analysis of complex features or damage sites can sometimes require modeling of high temperature enclosures. Implementing efficient and accurate view-factor algorithms required to model such problems is complex. The current work leverages the Non-equilibrium Radiation (NERO) software, which solves the radiation transport equation in a finite-volume scheme, to alleviating challenges often faced with view-factor calculations. By assuming heat transfer occurs only between grey bodies and that the medium is non-participating, computational cost of the method is significantly reduced. The enclosure physics are modeled through emitting and reflecting boundary conditions in NERO. The emitted radiative flux is dependent on the wall temperature which is a solution to the material response, obtained from Icarus, in this context. The Ares framework manages the time-advancement and exchange of the necessary data between the solvers. The surface energy balance is modified to account for the enclosure terms within the material response boundary condition. The methodology was verified against analytical solutions including radiating parallel plates, a hollow cylinder (shown in Fig. 1), and a hemisphere. Application of the methodology to inform the design of components of the Dragonfly system will be shown.

Ablation↗

Fault Characterization and Diagnostics Supporting Condition-Based Operation and Maintenance of Gas Turbine Engines

Condition-Based Operation and Maintenance (CBOM) is the state-of-the-art in maintenance approaches for gas turbine engines. CBOM applies engine sensor information to the optimization of future operation and maintenance procedures; this technique reduces costs for engine operators by minimizing unplanned outages and catastrophic engine degradation. However, due to the complexities of gas turbine operation and the lack of available engine sensors, further development is required to fully realize the benefits of CBOM. In the turbine section, rotating components interact with high temperature flows, which creates intense thermal and mechanical stresses. As a result, there are numerous mechanisms of component degradation in the turbine section. Furthermore, turbine components – like stator vanes and rotor blades – are among the most expensive in the engine because they are complex to design and manufacture. For these reasons, this dissertation addresses two main questions: (i) which parameters or faults within the turbine section are most important to monitor, and (ii) how can these parameters or faults be monitored in an engine-relevant environment? Although many turbine parameters and faults have been investigated in the open literature, there are some faults that are still not well understood. Rotor-casing eccentricity, which causes a non-constant blade tip clearance around the annulus, has not been investigated in terms of its effects on turbine efficiency. Therefore, the first study in this dissertation quantifies overall and local turbine efficiency for varying levels of rotor-casing eccentricity. Results showed negligible variations to overall turbine efficiency, meaning rotor-casing eccentricity only becomes relevant to CBOM when its severity causes rotordynamic issues. Purge flow is critical to turbine hardware longevity because it prevents ingestion of hot main gas path (MGP) flow into the under-platform regions. Despite its importance, there are currently no methods for monitoring purge flow performance in an engine environment. Therefore, the second study in this dissertation develops a predictive model for sealing effectiveness using inputs from two fast-response pressure sensors. Results exhibited low prediction errors across a full range of purge flow rates, which supports the viability of the modelling approach for CBOM. The final two studies in this dissertation address blade coolant flow monitoring. This cooling flow is responsible for protecting the turbine blades from the MGP flow, which exits the combustor at temperatures greater than the blade melting point. These studies showed that temperature measurements on the blade surface can be used to accurately predict blade coolant flow rate, and that defining the candidate features relative to the coolant trajectory is important for maintaining accuracy as coolant flow rate degradation occurs. This work enables blade coolant flow monitoring, which is currently not possible through existing condition monitoring techniques.

condition-based, gas turbines, diagnostics,↗

Evaluation of the Large-Scale and Regional Climatic Response Across North Africa to Natural Variability in Oceanic Modes and Terrestrial Vegetation Among the CMIP5 Models

Hydrologic variability is a serious threat to the poverty-stricken regions of North Africa. Meanwhile, the scientific community struggles to attribute these extreme climatic episodes to specific oceanic and land drivers. Prior modeling studies have not assessed simulated feedbacks over North Africa against an established observational benchmark, so their results are considered as untested, model-specific findings. The Coupled Model Intercomparison Project Phase Five (CMIP5) archive represents the state-of-the-art in climate projections, with most CMIP5 models now containing interactive vegetation phenology; however, there have been no studies to date of their simulated vegetation feedbacks. The capability of CMIP5 models at accurately simulating the forcing of both sea-surface temperature (SST) and regional leaf area index (LAI) anomalies on North African climate needs to be a key consideration in assessing the models’ overall credibility and determining appropriate model weighting for developing climate projections. The Generalized Equilibrium Feedback Assessment (GEFA) is a promising statistical method, which can be applied to either model output or observations, for isolating the local and remote impacts of individual oceanic or terrestrial forcings on regional climate. We performed a combined observational and modeling assessment of land-ocean-atmosphere interactions across the distinct ecological and moisture gradients of North Africa. We evaluated and demonstrated the reliability of the GEFA statistical method over North Africa using the Community Earth System Model (CESM), applied GEFA to observational data to quantify the observed forcing of ocean basin SST anomalies and regional LAI anomalies on North African climate, evaluated the CMIP5 models’ performance in terms of representing these key observed feedbacks, and formulated CMIP5 feedback performance metrics for weighting North African climate projections. Our study represented the first attempt to separate the observed roles of oceanic and vegetation feedbacks across North Africa, the first systematic assessment and intercomparison of land-ocean-atmosphere feedbacks in CMIP5, and the first exploration of vegetation feedbacks among CMIP5 models. The following work was accomplished. (1) The team generated 14 scientific publications and 30 presentations based on the DOE-funded research. (2) A stepwise version of the GEFA statistical method was developed in which unimportant forcings were dropped to increase the reliability of the results. (3) SGEFA was successfully validated through experiments with CESM in terms of its ability to isolate the atmospheric responses to individual oceanic or land forcings. (4) Observational evidence was revealed for the Sahel’s positive vegetation-rainfall feedback on the seasonal to interannual time scale, and it was attributed to a moisture recycling mechanism rather than an albedo mechanism. (5) An approach with SGEFA was developed in which the individual contributions of soil moisture versus vegetation forcings could be separated, revealing that the former forcing outweighed the latter for sub-Saharan Africa. (6) Tropical ocean temperatures were found to be key regulators of pan-tropical vegetation variability, especially for arid and semi-arid regions, including sub-Saharan Africa. (7) The CMIP5 models largely underestimated the importance of land feedbacks across the Sahel. (8) The general consensus among CMIP5 models indicates, for the late 21st century, a diminished seasonal predictability of sub-Saharan African regional climate and an elevated role of the land surface compared to oceanic drivers in regulating regional climate variability. (9) Land-ocean-atmosphere interactions were demonstrated to be key contributors to the seasonal predictability of African wildfire activity. (10) The El Djouf was determined to be of greater importance than the Bodélé depression in terms of providing trans-Atlantic dust transport to the Americas.

54 ENVIRONMENTAL SCIENCES↗

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Buffet characteristics of the F-8 supercritical wing airplane

The buffet characteristics of the F-8 supercritical wing airplane were investigated. Wing structural response was used to determine the buffet characteristics of the wing and these characteristics are compared with wind tunnel model data and the wing flow characteristics at transonic speeds. The wingtip accelerometer was used to determine the buffet onset boundary and to measure the buffet intensity characteristics of the airplane. The effects of moderate trailing edge flap deflections on the buffet onset boundary are presented. The supercritical wing flow characteristics were determined from wind tunnel and flight static pressure measurements and from a dynamic pressure sensor mounted on the flight test airplane in the vicinity of the shock wave that formed on the upper surface of the wing at transonic speeds. The comparison of the airplane's structural response data to the supercritical flow characteristics includes the effects of a leading edge vortex generator.

Deangelis, V. M.↗

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Effect of Configuration Pitching Motion on Twin Tail Buffet Response

The effect of dynamic pitch-up motion of delta wing on twin-tail buffet response is investigated. The computational model consists of a delta wing-twin tail configuration. The computations are carried out on a dynamic multi-block grid structure. This multidisciplinary problem is solved using three sets of equations which consists of the unsteady Navier-Stokes equations, the aeroelastic equations, and the grid displacement equations. The configuration is pitched-up from zero up to 60 deg. angle of attack, and the freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. With the twin tail fixed as rigid surfaces and with no-forced pitch-up motion, the problem is solved for the initial flow conditions. Next, the problem is solved for the twin-tail response for uncoupled bending and torsional vibrations due to the unsteady loads on the twin tail and due to the forced pitch-up motion. The dynamic pitch-up problem is also solved for the flow response with the twin tail kept rigid. The configuration is investigated for inboard position of the twin tail which corresponds to a separation distance between the twin tail of 33% wing chord. The computed results are compared with the available experimental data.

Sheta, Essam F.↗

Understanding Measured Water Rotational Temperatures and Column Densities in the Very Innermost Coma of Comet 73P/Schwassmann-Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann-Wachmann 3 obtained with a very high spatial resolution of approx. 30 km using IRCS/Subaru in May 2006. For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10-20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

comet↗

Climate Response to the Increase in Tropospheric Ozone since Preindustrial Times: A Comparison between Ozone and Equivalent CO2 Forcings

We examine the characteristics of the climate response to anthropogenic changes in tropospheric ozone. Using a general circulation model, we have carried out a pair of equilibrium climate simulations with realistic present-day and preindustrial ozone distributions. We find that the instantaneous radiative forcing of 0.49 W m(sup -2) due to the increase in tropospheric ozone since preindustrial times results in an increase in global mean surface temperature of 0.28 C. The increase is nearly 0.4 C in the Northern Hemisphere and about 0.2 C in the Southern Hemisphere. The largest increases (greater than 0.8 C) are downwind of Europe and Asia and over the North American interior in summer. In the lower stratosphere, global mean temperatures decrease by about 0.2 C due to the diminished upward flux of radiation at 9.6 micrometers. The largest stratospheric cooling, up to 1.0 C, occurs over high northern latitudes in winter, with possibly important implications for the formation of polar stratospheric clouds. To identify the characteristics of climate forcing unique to tropospheric ozone, we have conducted two additional climate equilibrium simulations: one in which preindustrial tropospheric ozone concentrations were increased everywhere by 18 ppb, producing the same global radiative forcing as present-day ozone but without the heterogeneity; and one in which CO2 was decreased by 25 ppm relative to present day, with ozone at present-day values, to again produce the same global radiative forcing but with the spectral signature of CO2 rather than ozone. In the first simulation (uniform increase of ozone), the global mean surface temperature increases by 0.25 C, with an interhemispheric difference of only 0.03 C, as compared with nearly 0.2 C for the heterogeneous ozone increase. In the second simulation (equivalent CO2), the global mean surface temperature increases by 0.36 C, 30% higher than the increase from tropospheric ozone. The stronger surface warming from CO2 is in part because CO2 forcing (obscured by water vapor) is shifted relatively poleward where the positive ice-albedo feedback amplifies the climate response and in part because the magnitude of the CO2 forcing in the mid-troposphere is double that of ozone. However, we find that CO2 is far less effective than tropospheric ozone in driving lower stratospheric cooling at high northern latitudes in winter.

Mickley L. J.↗

4-D Model of CO2 Deposition at North and South of Mars from HEND/Odyssey and MOLA/MGS

The first 1.5 year of neutron mapping measurements onboard Mars Odyssey spacecraft are presented based on High Energy Neutron Detector (HEND) observations. HEND instrument is a part of GRS suite responsible for registration of epithermal and fast neutrons originating in Mars subsurface layer. The scattering of fast neutrons in Mars surface caused by primary cosmic rays is strongly sensitive to presence of hydrogen atoms. Even several percents of subsurface water significantly depress epithermal and fast neutron flux. It turns orbit neutron spectroscopy into one of most efficient methods for finding distribution of subsurface water. The Mars Odyssey observations revealed huge water- ice regions above 60N and 60S latitudes. It was founded that distribution of subsurface water has layered structure at these regions. It is thought that more than 50% wt water ice covered by relatively dry layer with different thickness.

M L Litvak↗

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

A bulk flow model of a brush seal system

Fibers can be readily fabricated into a variety of seal configurations that are compliant and responsive to high speed or lightly loaded systems. A linear, circular, or contoured brush seal system is a contact seal consisting of the bristle pattern and hardened interface. When compared to a labyrinth seal, the brush seal system is superior and features low leakage, dynamic stability, and permits compliant structures. But in turn, the system usually requires a hardened smooth interface and permits only limited pressure drops. Wear life and wear debris for operations with static or dynamic excitation are largely undetermined. A seal system involves control of fluid within specific boundaries. The brush and rub ring (or rub surface) form a seal system. Design similitudes, a bulk flow model, and rub ring (interface) coatings are discussed. The bulk flow model calculations are based on flows in porous media and filters. The coatings work is based on experience and expanded to include current practice.

Hendricks, R. C.↗

Responses of Precipitation and Hydrologic Processes to Tropical SST

The goal of the research is to identify the mechanisms in the response of tropical precipitation and atmospheric hydrologic cycle to sea surface temperature (SST) variability at seasonal-to-interannual time scales, and to utilize the knowledge for better understanding of climate feedback processes relevant to global change. As a first step to achieve the goal, we characterize the inter-relationship among convective/stratiform rain, ice/water clouds water vapor, and SST using TRMM satellite data and a cloud-resolving model. We examined the daily hydrologic variables [column water vapor (PW), cloud liquid water (CW), rainfall rates (RR)] as a function of SST using high-resolution data (0.25 x 0.25, daily) derived from TRMM satellite measurements. Comparing the winter of 97/98 (El Nino condition) against the winter of 99/00 (La Nina condition), area-mean values of all four hydrologic variables in cloudy areas within the tropical Pacific are higher in the El Nino winter than in the La Nina winter. This is consistent with previous observational analyses and SST warming experiments (idealized or ENSO-like) that showed the interaction between hydrologic cyclic and radiation at the seasonal to interannual time scales leads to intensified tropical circulation and hydrologic cycle. However, there is evidence that the enhanced hydrologic cycle over the warm pool is accompanied by an expansion of radiatively -driven subsidence in response to a stronger SST gradient between warm pool and surrounding cold pool. The expanding subsidence effectively reduces cloud amounts over the warm pool. Our analysis of daily variability further indicates a more vigorous water cycle characterized by higher PW, CW, and RR in response to overall warming. This is expected from the Clausius Clapeyron relation as a thermodynamic response to warming. However cloudy areas decrease in response to overall warming. This may be due to factors that are fundamentally different. One possibility is that in a warmer climate, there may be alteration in the microphysical processes in clouds, e.g., increased coalescence of cloud size droplets and increased entrainment of dry air from above due to the invigorated convective motions in response to overall warming. As a result, cloud amounts are less. Another interpretation is that, in colder air masses, more ice particles will be produced and, because of their lower densities, are less likely to precipitate, hence the higher cloud amount. This suggests that, when averaged over a sufficiently large area covering both the rising and sinking branches of the tropical circulation, there maybe an increase in the precipitation efficiency in stratiform clouds, in conjunction with increased subsidence in an enhanced water cycle. This is supported by the residence time at different rainrates.

Sui, C.-H.↗

The effects of free surfaces on deformation twinning in HCP metals

Deformation twinning is a predominant mode of plastic deformation for hexagonal close packed metals, like Mg and Ti. The heterogenous microstructure and the local stresses associated with twinning play a key role in their mechanical response and fracture. Surface analyses, like electron microscopy, are frequently employed to spatially map microstructure and micromechanical fields in order to study twinning behavior. However, these measurements are inherently influenced by the vicinity of the free surface. In this work, an elasto-visco-plastic fast-Fourier-transform (EVP-FFT) polycrystal modeling approach is employed to investigate the effects of free surfaces on twin development before and after loading. We compare calculated micromechanical fields on free surfaces with those calculated inside the bulk and, in some cases, experimental surface measurements. The results indicate that the creation of free surfaces can promote twin propagation and growth and can influence twin morphology by causing a twin lamella to become larger, more blunted and irregular. The structure along the twin boundaries are also affected, due to the higher driving stresses that extend prismatic-basal and basal-prismatic facets. Furthermore, free surfaces invoke different slip activities in the twin and the surrounding parent crystal by enhancing basal, prismatic and pyramidal slip in some localized regions, while reducing slip in others. We demonstrate that the simulated free-surface effects lead to better qualitative and quantitative agreement with experimental measurements from scanning electron microscopy and digital image correlation.

36 MATERIALS SCIENCE↗

Biological mechanisms may contribute to soil carbon saturation patterns

Increasing soil organic carbon (SOC) storage is a key strategy to mitigate rising atmospheric CO 2 , yet SOC pools often appear to saturate, or increase at a declining rate, as carbon (C) inputs increase. Soil C saturation is commonly hypothesized to result from the finite amount of reactive mineral surface area available for retaining SOC, and is accordingly represented in SOC models as a physicochemically determined SOC upper limit. However, mineral-associated SOC is largely microbially generated. In this perspective, we present the hypothesis that apparent SOC saturation patterns could emerge as a result of ecological constraints on microbial biomass—for example, via competition or predation—leading to reduced C flow through microbes and a reduced rate of mineral-associated SOC formation as soil C inputs increase. Microbially explicit SOC models offer an opportunity to explore this hypothesis, yet most of these models predict linear microbial biomass increases with C inputs and insensitivity of SOC to input rates. Synthesis of 54 C addition studies revealed constraints on microbial biomass as C inputs increase. Different hypotheses limiting microbial density were embedded in a three-pool SOC model without explicit limits on mineral surface area. As inputs increased, the model demonstrated either no change, linear, or apparently saturating increases in mineral-associated and particulate SOC pools. Taken together, our results suggest that microbial constraints are common and could lead to reduced mineral-associated SOC formation as input rates increase. We conclude that SOC responses to altered C inputs—or any environmental change—are influenced by the ecological factors that limit microbial populations, allowing for a wider range of potential SOC responses to stimuli. Understanding how biotic versus abiotic factors contribute to these patterns will better enable us to predict and manage soil C dynamics.

54 ENVIRONMENTAL SCIENCES↗