Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “failure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Modelling a nickel-cadmium battery as a heterogeneous device

A common failure mode for nickel-cadmium batteries involves a nonuniform degradation of the performance of individual cells. This behavior is manifested as a divergence in voltage during cyclic operation. An empirically based performance model is presented which simulates a battery composed of heterogeneous cells. Battery performance is studied as a function of the distribution of cell performance characteristics. The technique has successfully predicted cell voltage divergence behavior observed during cycle life testing. It is noted that this work enhances understanding of life-test results and may improve life-test design and data analysis.

Glueck, Peter R.↗
NASA NTRS2002

A Nonlinear Viscoelastic Model for Ceramics at High Temperatures

High-temperature creep behavior of ceramics is characterized by nonlinear time-dependent responses, asymmetric behavior in tension and compression, and nucleation and coalescence of voids leading to creep rupture. Moreover, creep rupture experiments show considerable scatter or randomness in fatigue lives of nominally equal specimens. To capture the nonlinear, asymmetric time-dependent behavior, the standard linear viscoelastic solid model is modified. Nonlinearity and asymmetry are introduced in the volumetric components by using a nonlinear function similar to a hyperbolic sine function but modified to model asymmetry. The nonlinear viscoelastic model is implemented in an ABAQUS user material subroutine. To model the random formation and coalescence of voids, each element is assigned a failure strain sampled from a lognormal distribution. An element is deleted when its volumetric strain exceeds its failure strain. Element deletion has been implemented within ABAQUS. Temporal increases in strains produce a sequential loss of elements (a model for void nucleation and growth), which in turn leads to failure. Nonlinear viscoelastic model parameters are determined from uniaxial tensile and compressive creep experiments on silicon nitride. The model is then used to predict the deformation of four-point bending and ball-on-ring specimens. Simulation is used to predict statistical moments of creep rupture lives. Numerical simulation results compare well with results of experiments of four-point bending specimens. The analytical model is intended to be used to predict the creep rupture lives of ceramic parts in arbitrary stress conditions.

Powers, Lynn M.↗
NASA NTRS2010

Simulation Modeling for Off-Nominal Conditions - Where Are We Today?

The modeling of aircraft flight characteris4cs in off-nominal or otherwise adverse conditions has become increasingly important for simulation in the loss-of-control arena. Adverse conditions include environmentally-induced upsets such as wind shear or wake vortex encounters; off-nominal flight conditions, such as stall or departure; on-board systems failures; and structural failures or aircraft damage. Spirited discussions in the research community are taking place as to the fidelity and data requirements for adequate representation of vehicle dynamics under such conditions for a host of research areas, including recovery training, flight controls development, trajectory guidance/planning, and envelope limiting. The increasing need for multiple sources of data (empirical, computational, experimental) for modeling across a larger flight envelope leads to challenges in developing methods of appropriately applying or combining such data, particularly in a dynamic flight environment with a physically and/or aerodynamically asymmetric vehicle. Traditional simplifications and symmetry assumptions in current modeling methodology may no longer be valid. Furthermore, once modeled, challenges abound in the validation of flight dynamics characteristics in adverse flight regimes

Shah, Gautam H.↗
NASA NTRS2015

Real-Time Prognostics of a Rotary Valve Actuator

Valves are used in many domains and often have system-critical functions. As such, it is important to monitor the health of valves and their actuators and predict remaining useful life. In this work, we develop a model-based prognostics approach for a rotary valve actuator. Due to limited observability of the component with multiple failure modes, a lumped damage approach is proposed for estimation and prediction of damage progression. In order to support the goal of real-time prognostics, an approach to prediction is developed that does not require online simulation to compute remaining life, rather, a function mapping the damage state to remaining useful life is found offline so that predictions can be made quickly online with a single function evaluation. Simulation results demonstrate the overall methodology, validating the lumped damage approach and demonstrating real-time prognostics.

actuator↗
NASA NTRS2018

Quantification of Error Associated with Using Misaligned Meshes in Continuum Damage Mechanics Material Models for Matrix Crack Growth Predictions in Composites

The ability of a material model to capture in-plane matrix mode I and mode II crack growth is an essential component for modeling ply level damage evolution in composite structures. Previous studies using a continuum damage mechanics (CDM) approach have shown success in satisfying benchmark solutions for mode I and II crack growth. However, success was shown using a fiber-aligned meshing strategy, which encourages matrix cracks to propagate in a single band of elements, along the fiber direction. Generating a fiber-aligned mesh becomes a highly involved process for laminates including off-axis (non 0° or 90°) plies. The objective of this study is to quantify the effect of non-fiber aligned mesh discretization on predictions of in-plane matrix crack propagation. The approach taken incrementally varies the mesh orientation angle relative to the fiber orientation; more specifically, misaligned meshes are used to quantify the effect of element angle orientation relative to the initial crack orientation on the energy released during matrix crack propagation simulations using a CDM method. CDM solutions obtained with the misaligned meshes are evaluated against known benchmarks for mode I and II matrix crack growth. The CDM solutions reveal a near-polynomial trend of increased predicted failure stress with increased mesh misalignment angle; hence implying a potential relationship between element orientation angle and apparent fracture toughness.

Justusson, Brian↗
DOE OSTI2023

RELAP5-3D Solutions to Exercise 1 of the OECD-NEA HTTF Benchmark

The high-temperature gas-cooled reactor (HTGR) is an advanced reactor concept that has received considerable attention over the last 60 years. HTGRs boast high coolant outlet temperatures, passive safety, and large margins to fuel failure. Recent years have seen increasing interest in HTGRs, whether large-scale HTGRs for electricity or process heat applications or as microreactors serving remote communities. HTGRs have considerable potential to fulfill these applications, but modeling and simulation tools typically used for reactor safety analysis have not been validated for HTGR modeling. The High-Temperature Test Facility (HTTF) was constructed at Oregon State University (OSU) to provide validation data for HTGR modeling. Recently, Idaho National Laboratory (INL), OSU, Argonne National Laboratory, Canadian Nuclear Laboratories, and the University of Tennessee Knoxville have collaborated to develop an HTGR thermal hydraulics benchmark based on HTTF experiments. This summary presents results from a portion of that benchmark using RELAP5-3D

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2011

A Method of Strengthening Composite/Metal Joints

The term tape setback method denotes a method of designing and fabricating bonded joints between (1) box beams or other structural members made of laminated composite (matrix/ fiber) materials and (2) metal end fittings used to fasten these structural members to other structural members. The basic idea of the tape setback method is to mask the bonded interface between the metallic end fitting and composite member such that the bond does not extend out to the free edges of the composite member. The purpose served by the tape setback method is to strengthen the joints by decoupling stress concentrations from edge defects, which can cause premature failures. A related prior method that serves a similar purpose, involving the use of tapered adherends at the joints, can be too difficult and costly to be acceptable in some applications. The tape setback method offers an easier, less costly alternative. The structural members to which the method was originally applied were box beams in the form of composite tubes having flat faces with rounded corners. The end fittings were plugs made of a low-thermal- expansion nickel/iron alloy (see figure). In computational-simulation studies of tensile and compressive loading of members without tape setback, stresses were found to be concentrated at the free end edges of the composite tubes, and inspection of members that had been subjected to real tension and compression tests showed that cracks started at the free end edges. As applied to these members, the tape setback method makes them less vulnerable to initiation of failure at edge defects produced during fabrication. In real tension tests of comparable members without and with tape setback, the average mean tensile strength of the members with tape setback was found to be 1.9 times that of the members without tape setback.

Polis, Daniel L.↗
NASA NTRS2010

Random vs. Combinatorial Methods for Discrete Event Simulation of a Grid Computer Network

This study compared random and t-way combinatorial inputs of a network simulator, to determine if these two approaches produce significantly different deadlock detection for varying network configurations. Modeling deadlock detection is important for analyzing configuration changes that could inadvertently degrade network operations, or to determine modifications that could be made by attackers to deliberately induce deadlock. Discrete event simulation of a network may be conducted using random generation, of inputs. In this study, we compare random with combinatorial generation of inputs. Combinatorial (or t-way) testing requires every combination of any t parameter values to be covered by at least one test. Combinatorial methods can be highly effective because empirical data suggest that nearly all failures involve the interaction of a small number of parameters (1 to 6). Thus, for example, if all deadlocks involve at most 5-way interactions between n parameters, then exhaustive testing of all n-way interactions adds no additional information that would not be obtained by testing all 5-way interactions. While the maximum degree of interaction between parameters involved in the deadlocks clearly cannot be known in advance, covering all t-way interactions may be more efficient than using random generation of inputs. In this study we tested this hypothesis for t = 2, 3, and 4 for deadlock detection in a network simulation. Achieving the same degree of coverage provided by 4-way tests would have required approximately 3.2 times as many random tests; thus combinatorial methods were more efficient for detecting deadlocks involving a higher degree of interactions. The paper reviews explanations for these results and implications for modeling and simulation.

Kuhn, D. Richard↗
NASA NTRS2011

First Stage Solid Propellant Multiply Debris Thermal Analysis

Destruction of a solid rocket stage of a launch vehicle can create a thermal radiation hazard for an aborting crew module. This hazard was assessed for the Constellation Program (Cx) crew and launch vehicle concept. For this concept, if an abort was initiated in first stage flight, the Crew Module (CM) will separate and be pulled away from the malfunctioning launch vehicle via a Launch Abort System (LAS). Having aborted the mission, the launch vehicle will likely be destroyed via a Flight Termination System (FTS) in order to prevent it from errantly traversing back over land and posing a risk to the public. The resulting launch vehicle debris field, composed primarily of first stage solid propellant, poses a threat to the CM. The harsh radiative thermal environment, caused by surrounding burning propellant debris, may lead to CM parachute failure. A methodology, detailed herein, has been developed to address this concern and to quantify the risk of first stage propellant debris leading to the thermal demise of the CM parachutes. Utilizing basic thermal radiation principles, a software program was developed to calculate parachute temperature as a function of time for a given abort trajectory and debris piece trajectory set. Two test cases, considered worst case aborts with regard to launch vehicle debris environments, were analyzed using the simulation: an abort declared at Mach 1 and an abort declared at maximum dynamic pressure (Max Q). For both cases, the resulting temperature profiles indicated that thermal limits for the parachutes were not exceeded. However, short duration close encounters by single debris pieces did have a significant effect on parachute temperature. Therefore while these two test cases did not indicate exceedance of thermal limits, in order to quantify the risk of parachute failure due to radiative effects from the abort environment, a more thorough probability-based analysis using the methodology demonstrated herein must be performed.

Toleman, Benjamin M.↗
DOE OSTI2021

New Physics in the Hot Electron Furnace (Final Report)

The primary purpose of this award was experiments and simulations related to the idea of a ‘Hot Electron Furnace’ where nonlocal heat transport would allow unique new experimental opportunities at x-ray free electron lasers. The period of this award included beamtime at the MEC endstation of the LCLS, and the execution of those experiments was a primary goal. Unfortunately, the beamtime had limited success due to (1) inability of the LCLS to deliver the expected two-color operations mode, and (2) failure of the MEC endstation’s own software to save the cavity timing data that would have allowed inference of time delays for the few hours when two-color mode briefly operated. The beamtime was then switched to operate in single-color (one pulse) mode, and studies were performed of x-ray heating and probe, simultaneously. Two interesting manuscripts are finally near submission, having been delayed by student graduations, need for collaborator simulations of x-ray heating processes, and then COVID. I am hopeful that they will soon be submitted (especially the paper on x-ray heating of MgO) because one of the principal students involved has returned from industry to take a postdoc position at SLAC. Below, we first list all publications and other work products (including the two that are nearing submission) and then summarize several outcomes from this effort. During the term of this award three UW physics graduate students were at least partially supported, and all three received their Ph.D.’s with significant components of their supported work in the dissertation.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

Micromechanical simulation of damage progression in carbon phenolic composites

Carbon/phenolic composites are used extensively as ablative insulating materials in the nozzle region of solid rocket motors. The current solid rocket motor (RSRM) on the space shuttle is fabricated from woven rayon cloth which is carbonized and then impregnated with the phenolic resin. These plies are layed up in the desired configuration and cured to form the finished part. During firing, the surface of the carbon/phenolic insulation is exposed to 5000 F gases from the rocket exhaust. The resin pyrolizes and the material chars to a depth which progresses with time. The rate of charring and erosion are generally predictable, and the insulation depth is designed to allow adequate safety margins over the firing time of the motor. However, anomalies in the properties and response of the carbon/phenolic materials can lead to severe material damage which may decrease safety margins to unacceptable levels. Three macro damage modes which were observed in fired nozzles are: ply lift, 'wedge out', and pocketing erosion. Ply lift occurs in materials with plies oriented nearly parallel to the surface. The damage occurs in a region below the charred material where material temperatures are relatively low - about 500 F. Wedge out occurs at the intersection of nozzle components whose plies are oriented at about 45 deg. The corner of the block of material breaks off along a ply interface. Pocketing erosion occurs in material with plies oriented normal to the surface. Thermal expansion is restrained in two directions resulting in large tensile strains and material failure normal to the surface. When a large section of material is removed as a result of damage, the insulation thickness is reduced which may lead to failure of the nozzle due to excessive heating of critical components. If these damage events cannot be prevented with certainty, the designer must increase the thickness of the insulator thus adding to both weight and cost. One of the difficulties in developing a full understanding of these macro damage mechanisms is that the loading environment and the material response to that environment are extremely complex. These types of damage are usually only observed in actual motor firings. Therefore, it is difficult and expensive to evaluate the reliability of new materials. Standard material tests which measure mechanical and thermal properties of test specimens can only provide a partial picture of how the material will respond in the service environment. The development of the ANALOG test procedure which can combine high heating rates and mechanical loads on a specimen will improve the understanding of the interactive effects of the various loads on the system. But a mechanistic model of material response which can account for the heterogeneity of the material, the progression of various micromechanical damage mechanisms, and the interaction of mechanical and thermal stresses on the material is required to accurately correlate material tests with response to service environments. A model based on fundamental damage mechanisms which is calibrated and verified under a variety of loading conditions will provide a general tool for predicting the response of rocket nozzles. The development of a micromechanical simulation technique was initiated and demonstrated to be effective for studying across-ply tensile failure of carbon/phenolic composites.

Slattery, Kerry T.↗
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Residual-strength tests of L-1011 vertical fin components after 10 and 20 years of simulated flight service

Part of the NASA/ACEE Program was to determine the effect of long-term durability testing on the residual strength of graphite-epoxy cover panel and spar components of the Lockheed L-1011 aircraft vertical stabilizer. The results of these residual strength tests are presented herein. The structural behavior and failure mode of both cover panel and spar components were addressed, and the test results obtained were compared with the static test results generated by Lockheed. The effect of damage on one of the spar specimens was described.

Lopez, O. F.↗
DOE OSTI2023

Temperature Rise on Liner Surfaces of Fuel Cell Electric Vehicle Tanks during Fueling Process

Herein, the temperature of the inner tank walls, or plastic liners, of composite pressure tanks in fuel cell electric vehicles during fueling using three‐dimensional computational fluid dynamics (CFD) models is evaluated. The liner materials are limited to 85.0 °C to prevent thermal stress causing material failure that would result in a hydrogen leak. Therefore, the temperature of hydrogen gas must be limited to below 85.0 °C during the fueling process as dictated by the current fueling protocol. However, there are limited experimental or simulation data confirming that the temperature changes do not exceed the threshold. Herein, the liner temperatures with CFD tank models for two sizes of type IV tanks representative of the upper and lower system bounds that are close to the SAE J2601 fueling protocol (36.0 and 244.0 L) are evaluated. First, each model's reliability is validated with experimental data and then analyzed, and the data are used to evaluate the maximum hydrogen and liner temperatures under real‐world fueling conditions. The evaluation shows that the maximum liner surface temperature of each tank model is at least 7 °C lower than that of the hydrogen. Additionally, there is at least 12 °C difference found between the upper limit and actual liner temperatures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Beyond the Stokes approximation: shallow visco-elastic ice-sheet models

Abstract The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

58 GEOSCIENCES↗
DOE OSTI2022

The influence of increasing atmospheric CO 2 , temperature, and vapor pressure deficit on seawater-induced tree mortality

We report increasing seawater exposure is killing coastal trees globally, with expectations of accelerating mortality with rising sea levels. However, the impact of concomitant changes in atmospheric CO 2 concentration, temperature, and vapor pressure deficit (VPD) on seawater-induced tree mortality is uncertain. We examined the mechanisms of seawater-induced mortality under varying climate scenarios using a photosynthetic gain and hydraulic cost optimization model validated against observations in a mature stand of Sitka-spruce (Picea sitchensis) trees in the Pacific Northwest, USA, that were dying from recent seawater exposure. The simulations matched well with observations of photosynthesis, transpiration, nonstructural carbohydrates concentrations, leaf water potential, the percentage loss of xylem conductivity, and stand-level mortality rates. The simulations suggest that seawater-induced mortality could decrease by ~16.7% with increasing atmospheric CO 2 levels due to reduced risk of carbon starvation. Conversely, rising VPD could increase mortality by ~5.6% because of increasing risk of hydraulic failure. Across all scenarios, seawater-induced mortality was driven by hydraulic failure in the first two years after seawater exposure began, with carbon starvation becoming more important in subsequent years. Changing CO 2 and climate appear unlikely to have a significant impact on coastal tree mortality under rising sea levels.

54 ENVIRONMENTAL SCIENCES↗