Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vessel 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 487 records · Page 27

Stability of tropical forest tree carbon-water relations in a rainfall exclusion treatment through shifts in effective water uptake depth

Increasing severity and frequency of drought is predicted for large portions of the terrestrial biosphere, with major impacts already documented in wet tropical forests. Using a 4-year rainfall exclusion experiment in the Daintree Rainforest in northeast Australia, we examined canopy tree responses to reduced precipitation and soil water availability by quantifying seasonal changes in plant hydraulic and carbon traits for 11 tree species between control and drought treatments. Even with reduced soil volumetric water content in the upper 1 m of soil in the drought treatment, we found no significant difference between treatments for predawn and midday leaf water potential, photosynthesis, stomatal conductance, foliar stable carbon isotope composition, leaf mass per area, turgor loss point, xylem vessel anatomy, or leaf and stem nonstructural carbohydrates. While empirical measurements of aboveground traits revealed homeostatic maintenance of plant water status and traits in response to reduced soil moisture, modeled belowground dynamics revealed that trees in the drought treatment shifted the depth from which water was acquired to deeper soil layers. Finally, these findings reveal that belowground acclimation of tree water uptake depth may buffer tropical rainforests from more severe droughts that may arise in future with climate change.

54 ENVIRONMENTAL SCIENCES↗

Attached Jet Velocity Profiles in Mixing Tanks

Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.

jet, attached jet, offset jet↗

Optimized Tandem Catalyst Patterning for CO 2 Reduction Flow Reactors

Tandem catalysis involves two or more catalysts arranged in proximity within a single reaction vessel, with the aim of synergistically aligning the catalysts’ reaction pathways to maximize overall system performance. This study presents a proof of concept showing the integration of continuum transport modeling with design optimization in a simplified two-dimensional flow reactor setup for electrochemical CO 2 reduction. Ag catalysts provide the CO 2 ⟶ CO reaction capability, and Cu catalysts provide the CO ⟶ high-value products reaction capability. Given a set of input parameters, the optimization algorithm uses adjoint methods to modify the Ag/Cu surface patterning in order to maximize the current density toward high-value products, such as ethylene. The optimized designs yield significant performance enhancement especially at more negative applied voltages (i.e., stronger surface reactions) and for larger numbers of patterning sections. For an applied voltage of −1.7 V vs. SHE, the 12-section optimized design increases the current density toward ethylene by up to 65% compared to the unoptimized 2-section design. For the optimized cases, observed differences in the production and consumption of CO (the key intermediate species) and minimized zones of low CO reactant surface concentration on Cu sections explain the improved reactor performance.

CO2 reduction↗

MELCOR Accident Progression and Source Term Demonstration Calculations for a HTGR

MELCOR is an integrated thermal hydraulics, accident progression, and source term code for reactor safety analysis that has been developed at Sandia National Laboratories for the United States Nuclear Regulatory Commission (NRC) since the early 1980s. Though MELCOR originated as a light water reactor (LWR) code, development and modernization efforts over the past decades have expanded its application scope to include non-LWR reactor concepts. Current MELCOR development efforts include providing the NRC with the analytical capabilities to support regulatory readiness for licensing non-LWR technologies under Strategy 2 of the NRC's near-term Implementation Action Plans. Beginning with the Next Generation Nuclear Project (NGNP), MELCOR ha s undergone a range of enhancements to provide analytical capabilities for modeling the spectrum of advanced non-LWR concepts. This report describes the generic plant model developed to demonstrate MELCOR capabilities to perform high-temperature gas reactor (HTGR) safety evaluations. The generic plant model is based on publicly available PMBR-400 design information. For plant aspects (e.g., reactor building size and leak rate) that are not described in the PBMR-400 references, the analysts made assumptions needed to construct a MELCOR full-plant model. The HTGR model uses a TRi-structural ISOtropic (TRISO)-particle fuel pebble-bed reactor with a primary system rejecting heat to a recuperative heat exchange r. Surrounding the reactor vessel is a reactor cavity contained within a confinement room cooled by the Reactor Cavity Cooling System (RCCS). Example calculations are performed to show the plant response and MELCOR capabilities to characterize a range of accident conditions. The accidents selected for evaluation consider a range of degraded and failed modes of operation for key safety functions providing reactivity control, primary system heat removal and reactor vessel decay heat removal, and confinement cooling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Alterations in Skeletal Muscle Microcirculation of Head-Down Tilted Rats

In this study we assessed the function of microscopic blood vessels in skeletal muscle (cremaster muscle) for alterations which may contribute to the observed elevation of blood pressure associated with head-down tilted whole body suspension (HDT/WBS), a model of weightlessness. Arteriolar baseline diameters, vasoconstrictor responses to norepinephrine (NE) and vasodilation to nitroprusside (NP) were assessed in control rats, rats suspended for 7 or 14 day HDT/WBS rats, and rats allowed to recover for 1 day after 7 days HDT/WBS. Neither baseline diameters nor ability to dilate were influenced by HDT/WBS. Maximum vasoconstriction to norepinephrine was significantly greater in arterioles of hypertensive 14 day HDT/WBS rats. This first study of the intact microvasculature in skeletal muscle indicates that an elevated contractility of arterioles to norepinephrine in suspended rats, and suggests an elevated peripheral resistance in striated muscle may contribute to the increase in blood pressures among animals subjected to HDT/WBS.

Musacchia, X. J.↗

The fluid dynamic and shear environment in the NASA/JSC rotating-wall perfused-vessel bioreactor

The rotating-wall perfused-vessel (RWPV) bioreactor, used for both microgravity and Earth-based cell science experiments, is characterized in terms of the fluid dynamic and fluid shear stress environment. A numerical model of the flow field is developed and verified with laser Doppler velocimeter measurements. The effects of changes in operating conditions, including rotation rates and fluid perfusion rates, are investigated with the numerical model. The operating conditions typically used for ground-based experiments (equal rotation of the inner and outer cylinders) leads to flow patterns with relatively poor mass distribution characteristics. Approximately 50% of the inlet-perfused fluid bypasses the bulk of the fluid volume and flows to the perfusion exit. For operating conditions typical in microgravity, small differential rotation rates between the inner and outer cylinders lead to greatly improved flow distribution patterns and very low fluid shear stress levels over a large percentage of the fluid volume. Differences in flow patterns for the different operating conditions are explored. Large differences in the hydrodynamic environments for operating conditions typical of true microgravity and ground-based "microgravity simulations" are demonstrated.

Bioreactors↗

Femoral Blood Flow and Cardiac Output During Blood Flow Restricted Leg Press Exercise

Low load blood flow restricted resistance exercise (LBFR) causes muscle hypertrophy that may be stimulated by the local ischemic environment created by the cuff pressure. However, local blood flow (BF) during such exercise is not well understood. PURPOSE: To characterize femoral artery BF and cardiac output (CO) during leg press exercise (LP) performed at a high load (HL) and low load (LL) with different levels of cuff pressure. METHODS: Eleven subjects (men/women 4/7, age 31.4+/-12.8 y, weight 68.9+/-13.2 kg, mean+/-SD) performed 3 sets of supine left LP to fatigue with 90 s of rest in 4 conditions: HL (%1-RM/cuff pressure: 80%/0); LL (20%/0); LBFR(sub DBP) (20%/1.3 x diastolic blood pressure, BP); LBFR(sub SBP) (20%/1.3 x supine systolic BP). The cuff remained inflated throughout the LBFR exercise sessions. Artery diameter, velocity time integral (VTI), and stroke volume (SV) were measured using Doppler ultrasound at rest and immediately after each set of exercise. Heart rate (HR) was monitored using a 3-lead ECG. BF was calculated as VTI x vessel cross-sectional area. CO was calculated as HR x SV. The data obtained after each set of exercise were averaged and used for analyses. Multi-level modeling was used to determine the effect of exercise condition on dependent variables. Statistical significance was set a priori at p< 0.05. RESULTS: Artery diameter did not change from baseline. BF increased (p<0.05) after exercise in each condition except LBFR(sub SBP) in the order of HL (12.73+/-1.42 cm3,mean+/-SE) > LL (9.92+/-0.82 cm3) > LBFR(sub dBP)(6.47+/-0.79 cm3) > LBFR(sub SBP) (3.51+/-0.59 cm3). Blunted exercise induced increases occurred in HR, SV, and CO after LBFR compared to HL and LL. HR increased 45% after HL and LL and 28% after LBFR (p<0.05), but SV increased (p<0.05) only after HL. Consequently, the increase (p<0.05) in CO was greater in HL and LL (approximately 3 L/min) than in LBFR (approximately 1 L/min). CONCLUSION: BF during LBFR(sub SBP) was 1/3 of that observed in LL, which supports the hypothesis that local ischemia stimulates the LBFR hypertrophic response. As the cuff did not compress the artery, the ischemia may have occurred because of the blunted rise in CO or because arterial BP cannot overcome the cuff pressure. As LBFR(sub DBP) effectively reduced BF and CO with cuff pressures less than systolic BP, future studies should investigate the hypertrophic potential of LBFR at even lower cuff pressures.

Everett, M. E.↗

Evaluation of Human Spaceflight-Related Tissue Weight Relief Using Whole Body Finite Element Model Simulations

Tissue Weight Relief (TWR) is a physiological condition observed in human spaceflight. It not only impacts the injury biomechanics of soft tissue but also the physiological responses of the cardiovascular system both due to fluid redistribution and the effect of tissue-related transmural pressure on the large venous blood vessels. Understanding the effects of tissue weight relief is especially important because of the role it may play in understanding the cause of Space Associated Neuro-Ocular Syndrome (SANS). SANS can be characterized by a number of ocular changes which reduce visual acuity and SANS related symptoms occur in up to 51% of astronauts. A prevailing theory for the causation of SANS is that of headward (cephalad) fluid shift and a prolonged increase of Intracranial Pressure (ICP) similar to intracranial hypertension, which is not fully supported by the experimental data or astronaut symptom reporting. However, it is still believed that SANS is caused by a pressure change in the eye and the surrounding tissues. It has been proposed that TWR plays a substantial role in affecting internal pressures and fluid shifts in microgravity. In this effort, two whole-body Finite Element (FE) models – Elemance and THUMS – are used to ascertain the microgravity-associated TWR of the musculature surrounding the lower body veins. Elemance and THUMS are physics-based computational models that have been validated and verified for several automotive and domestic applications, and as such, can simulate the relief of soft tissue weight due to changes in the gravitation vector. Specifically, the current effort modeled the transition of the gravitational vector from 1G to 0G, applied across the whole-body model in a supine position. For each 1G to 0G transition simulation, the lower body vein’s transmural pressure-time profile was extracted and averaged around the anterior portion of the thigh muscle. The ascertained transmural pressure changes from 1G to 0G transition are given in Figure 1 for the Elemance and the THUMS FE models. The transmural pressure changes of 10 mmHg and 21 mmHg are in the same order of magnitude as Lu’s value of 44 mmHg. It is to be noted that Lu implemented a 0D to 1D lumped parameter model and the Elemance and THUMS are 3D higher order computational models. This proof-of-concept approach demonstrates that TWR pressure can be adequately estimated with in silico techniques however, further in silico investigations need to be conducted to address the unique contributions to the transmural pressure from each of the computational models.

Finite element modeling↗

Orion Meteoroid and Debris Analysis with a Single Kevlar® Layer Enhancement to the Propulsion Tank Multi-Layer Insulation: Comparison of the Observed Data and the Adjusted Ballistic Models

In a continual mass optimization effort for the Orion vehicle by NASA and its prime contractor Lockheed Martin Space, along with the European Space Agency (ESA) and its contractors Airbus and Thales Alenia Space-Italy (TAS-I), a series of shield characterization shots have been performed to quantify the ballistic enhancement of a single-layer of Kevlar® in multi-layer insulation (MLI). This enhanced MLI is part of Orion’s propulsion subsystem, and it is planned to replace a double-layer Kevlar® configuration starting at the Artemis IV mission. This effort has included fifty-two shield characterization shots for five different shield configurations over the four Orion propellant storage vessels. These storage vessels are located in the European Service Module (ESM). This shield characterization has been managed by the Hypervelocity Impact Technology (HVIT) group in NASA Astromaterials Research and Exploration Sciences at Johnson Space Center (JSC) for the NASA Engineering and Safety Center (NESC) and the Multi-Purpose Crew Vehicle (MPCV) program office in coordination with Lockheed Martin Space and Thales Alenia Space-Italy (TAS-I). The shield characterization has used the two-stage, light-gas-gun at the Remote Hypervelocity Test Laboratory (RHTL) of NASA JSC White Sands Test Facility (WSTF) for acceleration of projectiles to representative orbital impact speeds. From this shield characterization effort, a revised ballistic limiting equation (BLE) has been developed for each of the five representative shielding configurations of the ESM pressure vessels. This document provides a brief description of the shield configurations considered for the Orion program and the shield characterization findings. These results are compared to the assumed BLE for each configuration going into the Artemis missions that utilize this version of ballistic enhanced MLI. The models are integrated into a fourteen day High Lunar Orbit (HLO) mission as a representative reliability prediction impact for Orion, and it has been found that the predicted risk of a loss-of-crew (LOC) event is actually a little less than half of the assumed value that launched this design modification opportunity.

Joshua E Miller↗

Report on Data Quality Required for Modeling an MSR

Calculations were performed to assess how uncertainties in molten salt property values impact the results of behavior models of molten salt reactor systems under steady-state conditions. The thermophysical and thermochemical properties of salt mixtures being considered for use in molten salt reactors are not yet well characterized and methods used to measure property values to molten salts are not yet well established. The limited property data for salts of interest that are available in the literature are inconsistent and unreliable. Many methods used to measure property values remain developmental in the sense that not all variables affecting measured values have been identified or controlled during the measurement and not all aspects of the measurements affecting data quality are calibrated. A prior report described the repeatability of measurements made at Argonne and the reproducibility of measurements made with the same salts at several institutions is being evaluated. There is significant uncertainty in measured property values due to instrumental limitations and the effects of contaminants in the salts used to make the measurements that impacts predictions of salt behavior in a reactor. Most salts of interest are multicomponent eutectic or near-eutectic mixtures and both uncertainty in the salt composition and uncertainty in measured values. Furthermore, the salt composition evolves over time due to the uptake of impurities from the atmosphere and corrosion of containment vessels. More significantly, fission products will accumulate during operation that change the salt compositions and property values. MSR developers need to understand how this uncertainty in the basic property data of the fuel will impact modeling predictions for normal reactor operation as well was in transient scenarios that will impact the safety basis. The calculations and analyses in this report address the impact of uncertainty in measured property values under steady-state conditions through a generalized parametric study. The pooled uncertainty from contributions by a number of factors were considered and quantified, including varying composition due to impurities and fission products, imprecision in the measured values, and uncontrolled aspects of the measurement. In this way, we can begin to gain an understanding of the magnitude of the impact uncertainty will have in modeling MSR systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Simulation of the NuScale SMR and Investigation of the Effect of Load-Following on Component Lifetimes

The NuScale SMR has been modelled using the Virtual Environment for Reactor Applications (VERA) multiphysics environment and the results compared with the publicly reported data in the Design Certification Application (DCA). The results show an excellent agreement for the compared axial and radial power distributions, temperature coefficients of reactivity, boron and control rod worths, and fast neutron flux. This NuScale model is then used to investigate the effect of different operational modes on reactor components to determine how flexible load-following operation may affect control rod and reactor pressure vessel (RPV) lifetimes. The control rod degradation is confirmed to primarily affect the silver-indium-cadmium (AIC) rod tip. The degradation rate is observed to follow a non-linear function of core power level where the increase in degradation decreases with insertion depth. For the variation in core power levels expected with current load-following schemes, the total control rod degradation is found to be mild, at 5-10% of usable life per cycle for a reactor operating at <80% power. Nonetheless, this enables load following strategies to be confirmed and/or modified to ensure that control rods do not need to be replaced during the 60+ year life of the reactor. The RPV degradation was found to be almost directly proportional to the core power level and was not overly sensitive to flux shape perturbations. Future work is planned using these damage functions to optimize operation over multiple NuScale SMR units and develop strategies for prognostics and health management.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Elastic-Plastic Thermal Stress Analysis of a High-Pressure Cryogenic Storage Tank

The thermal stresses on a cryogenic storage tank contribute strongly to the state of stress of the tank material and its ability to withstand operational stresses. These thermal stresses also affect the growth of any surface damage that might occur in the tank walls. These stresses are particularly of concern during the initial cooldown period for a new tank placed into service, and during any subsequent thermal cycles. A previous preliminary elastic analysis showed that the thermal stress on the inner wall would reach approximately 1,000MPa (145,000 psi). This stress far exceeds the ASTM specified room temperature values for both yield (170MPa) and ultimate (485 MPa) strength for 304L stainless steel. The present analysis determines the thermal stresses using an elastic-plastic model. The commercial software application ANSYS was used to determine the transient spatial temperature profile and the associated spatial thermal stress profiles in a segment of a thick-walled vessel during a typical cooldown process. A strictly elastic analysis using standard material properties for 304L stainless steel showed that the maximum thermal stress on the inner and outer walls was approximately 960 MPa (tensile) and - 270 MPa (compressive) respectively. These values occurred early in the cooldown process, but at different times, An elastic-plastic analysis showed significantly reducing stress, as expected due to the plastic deformation of the material. The maximum stress for the inner wall was approximately 225 MPa (tensile), while the maximum stress for the outer wall was approximately - 130 MPa (compressive).

Barker, J. Mark↗

Reliability Considerations for Composite Overwrapped Pressure Vessels on Spacecraft

Composite Overwrapped Pressure Vessels (COPVs) are used to store gases under high pressure onboard spacecraft. These are used for a variety of purposes such as propelling liquid fuel etc, Kevlar, glass, Carbon and other more recent fibers have all been in use to overwrap the vessels. COPVs usually have a thin metallic liner with the primary purpose of containing the gases and prevent any leakage. The liner is overwrapped with filament wound composite such as Kevlar, Carbon or Glass fiber. Although the liner is required to perform in the leak before break mode making the failure a relatively benign mode, the overwrap can fail catastrophically under sustained load due to stress rupture. It is this failure mode that is of major concern as the stored energy of such vessels is often great enough ta cause loss of crew and vehicle. The present paper addresses some of the reliability concerns associated specifically with Kevlar Composite Overwrapped Pressure Vessels. The primary focus of the paper is on how reliability of COPV's are established for the purpose of deciding in general their flight worthiness and continued use. Analytical models based on existing design data will be presented showing how to achieve the required reliability metric to the end of a specific period of performance. Uncertainties in the design parameters and how they affect reliability and confidence intervals will be addressed as well. Some trade studies showing how reliability changes with time during a program period will be presented.

Murthy, Pappu L. N.↗

Nondestructive Evaluation and Monitoring Results from COPV Accelerated Stress Rupture Testing, NASA White Sands Test Facility (WSTF)

Develop and demonstrate NDE techniques for real-time characterization of CPVs and, where possible, identification of NDE capable of assessing stress rupture related strength degradation and/or making vessel life predictions (structural health monitoring or periodic inspection modes). Secondary: Provide the COPV user and materials community with quality carbon/epoxy (C/Ep) COPV stress rupture progression rate data. Aid in modeling, manufacturing, and application of COPVs for NASA spacecraft.

Saulsberry Regor↗

Computational Fluid Dynamics Simulation of Methane Slosh and Drain Experiments

NASA possesses a wealth of historical cryogenic experiments that provide valuable insights into the design, troubleshooting, and understanding involved in the complex fluid and thermodynamics of managing cryogenic propellants. One approach to leveraging these historical datasets is by simulating these experiments to validate the accuracy of simulation environments and constituent models. This study focused on simulating a selection of the K site test series for both static and sloshing pressurized liquid methane draining experiments conducted at NASA in the 1970’s, utilizing computational fluid dynamics. The simulations were performed in the ANSYS FLUENT environment using the Volume of Fluid (VOF) numerical approach. A k-omega turbulence model was used with interfacial turbulence damping, and accurately predicted the amount of pressurant needed to maintain the required tank pressure throughout the static drain. The static simulation predicted the temperature stratification in the ullage observed at the end of the drain. During the methane expulsion with sloshing test, many features were successfully captured using the k-omega turbulence model with interfacial turbulence damping included. The rate of phase change and liquid temperature was overpredicted compared to the experimental measurements. The overprediction may be attributed to uncertainties in the vessel geometry, methane pressurant temperature and composition, and methodological differences in how the sloshing frequency was adjusted during the expulsion.

Cryogenic Propellants↗

Effect of ATF Cr-coated-Zircaloy on BWR In-vessel Accident Progression during a Station Blackout

The deposition of protective coatings on nuclear fuel cladding has been considered as a near-term Accident Tolerant Fuel (ATF) concept that will reduce the high-temperature oxidation rate and enhance accident tolerance of the cladding while providing additional benefits during normal and transient. In this study, the performance of the proposed ATF concept of Cr-coated-Zircaloy is assessed using a generic Boiling Water Reactor MELCOR plant model considering a Short-term Station Blackout (STSBO) scenario. Simulation results indicate that the use of Cr-coated-Zircaloy as cladding and canister material mitigates the core degradation process as compared to the traditional Zircaloy cladding and canister design. The onset of fuel degradation and collapse is delayed by over thirty minutes, and the extent of fuel degradation is reduced. Specifically, the gross in-vessel hydrogen generation decreased by almost a factor of three. Although the eutectic reaction between Cr-coating and Zircaloy could cause an early failure of the coating, the improvement in the delay of fuel degradation is still notable. In addition, a thicker coating is found helpful to obtain additional coping time and to decrease hydrogen generation. In addition to the eutectic formation that may compromise Cr-coated Zr, a different failure mode is identified for the Cr-coated-Zr when compared to Zircaloy; i.e., a complete melt of base material leads to component collapse before the coating is oxidized and consumed. These findings can help the industry focus on productive areas of research and development for accident-tolerant fuel concepts and enhancement of core safety margins.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Graphite Degradation Modeling and Analysis

A graphite component in a nuclear reactor core is subjected to variety of stresses and can experience degradation during normal and off-normal operation. Understanding how a graphite component will behave in service is essential to ensuring core structural stability and safe reactor operation. This report summarizes a graphite modeling tool currently under development at Idaho National Laboratory. The model incorporates several of the anticipated stresses during service and includes the effects of oxidation and irradiation prior to turnaround. This tool is intended to be used to help assess the design of graphite components by utilizing design code rules found in Section III, Division 5 of the American Society of Mechanical Engineering Boiler and Pressure Vessel Code. Specifically, the tool uses the methodologies found within the Full and Simplified assessments from Article HHA-3000 to verify that a graphite component has an acceptably low probability of failure.

36 MATERIALS SCIENCE↗

Localizing Clinical Patterns of Blast Traumatic Brain Injury Through Computational Modeling and Simulation

Blast traumatic brain injury is ubiquitous in modern military conflict with significant morbidity and mortality. Yet the mechanism by which blast overpressure waves cause specific intracranial injury in humans remains unclear. Reviewing of both the clinical experience of neurointensivists and neurosurgeons who treated service members exposed to blast have revealed a pattern of injury to cerebral blood vessels, manifested as subarachnoid hemorrhage, pseudoaneurysm, and early diffuse cerebral edema. Additionally, a seminal neuropathologic case series of victims of blast traumatic brain injury (TBI) showed unique astroglial scarring patterns at the following tissue interfaces: subpial glial plate, perivascular, periventricular, and cerebral gray-white interface. The uniting feature of both the clinical and neuropathologic findings in blast TBI is the co-location of injury to material interfaces, be it solid-fluid or solid-solid interface. This motivates the hypothesis that blast TBI is an injury at the intracranial mechanical interfaces. In order to investigate the intracranial interface dynamics, we performed a novel set of computational simulations using a model human head simplified but containing models of gyri, sulci, cerebrospinal fluid (CSF), ventricles, and vasculature with high spatial resolution of the mechanical interfaces. Simulations were performed within a hybrid Eulerian—Lagrangian simulation suite (CTH coupled via Zapotec to Sierra Mechanics). Because of the large computational meshes, simulations required high performance computing resources. Twenty simulations were performed across multiple exposure scenarios—overpressures of 150, 250, and 500 kPa with 1 ms overpressure durations—for multiple blast exposures (front blast, side blast, and wall blast) across large variations in material model parameters (brain shear properties, skull elastic moduli). All simulations predict fluid cavitation within CSF (where intracerebral vasculature reside) with cavitation occurring deep and diffusely into cerebral sulci. These cavitation events are adjacent to high interface strain rates at the subpial glial plate. Larger overpressure simulations (250 and 500kPa) demonstrated intraventricular cavitation—also associated with adjacent high periventricular strain rates. Additionally, models of embedded intraparenchymal vascular structures—with diameters as small as 0.6 mm—predicted intravascular cavitation with adjacent high perivascular strain rates. The co-location of local maxima of strain rates near several of the regions that appear to be preferentially damaged in blast TBI (vascular structures, subpial glial plate, perivascular regions, and periventricular regions) suggest that intracranial interface dynamics may be important in understanding how blast overpressures leads to intracranial injury.

59 BASIC BIOLOGICAL SCIENCES↗