Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “injection rate”

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,423 records · Page 79

Experimental studies of shock-wave/wall-jet interaction in hypersonic flow

Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer -- the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.

Holden, Michael S.↗

Experimental studies of shock-wave/wall-jet interaction in hypersonic flow, part A

Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident-shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer, the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.

Holden, Michael S.↗

In-flight imaging of transverse gas jets injected into transonic and supersonic crossflows: Design and development

The design and development of an airborne flight-test experiment to study nonreacting gas jets injected transversely into transonic and supersonic crossflows is presented. Free-stream/crossflow Mach numbers range from 0.8 to 2.0. Planar laser-induced fluorescence (PLIF) of an iodine-seeded nitrogen jet is used to visualize the jet flow. Time-dependent images are obtained with a high-speed intensified video camera synchronized to the laser pulse rate. The entire experimental assembly is configured compactly inside a unique flight-test-fixture (FTF) mounted under the fuselage of the F-104G research aircraft, which serves as a 'flying wind tunnel' at NASA Dryden Flight Research Center. The aircraft is flown at predetermined speeds and altitudes to permit a perfectly expanded (or slightly underexpanded) gas jet to form just outside the FTF at each free-stream Mach number. Recorded gas jet images are then digitized to allow analysis of jet trajectory, spreading, and mixing characteristics. Comparisons will be made with analytical and numerical predictions. This study shows the viability of applying highly sophisticated groundbased flow diagnostic techniques to flight-test vehicle platforms that can achieve a wide range of thermo/fluid dynamic conditions. Realistic flow environments, high enthalpies, unconstrained flowfields, and moderate operating costs are also realized, in contrast to traditional wind-tunnel testing.

Wang, Kon-Sheng Charles↗

An Assessment of the Ozone Loss During the 1999-2000 SOLVE Arctic Campaign

Ozone observations from ozonesondes, the lidars aboard the DC-8, in situ ozone measurements from the ER-2, and satellite ozone measurements from Polar Ozone and Aerosol Measurement III (POAM) were used to assess ozone loss during the Sage III Ozone Loss and Validation Experiment (SOLVE) 1999-2000 Arctic campaign. Two methods of analysis were used. In the first method a simple regression analysis is performed on the ozonesonde and POAM measurements within the vortex. In the second method, the ozone measurements from all available ozone data were injected into a free running diabatic trajectory model and carried forward in time from December 1 to March 15. Vortex ozone loss was then estimated by comparing the ozone values of those parcels initiated early in the campaign with those parcels injected later in the campaign. Despite the variety of observational techniques used during SOLVE, the measurements provide a fairly consistent picture. Over the whole vortex, the largest ozone loss occurs between 550 and 400 K potential temperatures (approximately 23-16 km) with over 1.5 ppmv lost by March 15, the end of the SOLVE mission period. An ozone loss rate of 0.04-0.05 ppmv/day was computed for March 15. Ozonesondes launched after March 15 suggest that an additional 0.5 ppmv or more ozone was lost between March 15 and April 1. The small disagreement between ozonesonde and POAM analysis of January ozone loss is found to be due to biases in vortex sampling. POAM makes most of its solar occultation measurements at the vortex edge during January 2000 which bias samples toward air parcels that have been exposed to sunlight and likely do experience ozone loss. Ozonesonde measurements and the trajectory technique use observations that are more distributed within the interior of the vortex. Thus the regression analysis of the POAM measurements tends to overestimate mid-winter vortex ozone loss. Finally, our loss calculations are broadly consistent with other loss computations using ER-2 tracer data and MLS satellite data, but we find no evidence for the 1992 high mid-January loss reported using the Match technique.

Schoeberl, Mark R.↗

Wide-Temperature-Range Integrated Operational Amplifier

A document discusses a silicon-on-insulator (SOI) complementary metal oxide/semiconductor (CMOS) integrated- circuit operational amplifier to be replicated and incorporated into sensor and actuator systems of Mars-explorer robots. This amplifier is designed to function at a supply potential less than or equal to 5.5 V, at any temperature from -180 to +120 C. The design is implemented on a commercial radiation-hard SOI CMOS process rated for a supply potential of less than or equal to 3.6 V and temperatures from -55 to +110 C. The design incorporates several innovations to achieve this, the main ones being the following: NMOS transistor channel lengths below 1 m are generally not used because research showed that this change could reduce the adverse effect of hot carrier injection on the lifetimes of transistors at low temperatures. To enable the amplifier to withstand the 5.5-V supply potential, a circuit topology including cascade devices, clamping devices, and dynamic voltage biasing was adopted so that no individual transistor would be exposed to more than 3.6 V. To minimize undesired variations in performance over the temperature range, the transistors in the amplifier are biased by circuitry that maintains a constant inversion coefficient over the temperature range.

Mojarradi, Mohammad↗

Performance of the Fermilab Linac Injector

The Fermilab Linac injection line consists of a 35 keV magnetron-type H- ion source, two-solenoid Low Energy Beam Transport (LEBT), 201 MHz 4-rod 750 keV Radio Frequency Quadrupole (RFQ), and a Medium Energy Transport (MEBT) containing 4 quadrupoles and a bunching cavity. The injector delivers 25 mA, 48 s pulses to drift-tube Linac at a repetition rate of 15Hz. The transmission efficiency has been lower than expected since commissioning. Recent beam current measurements suggest that the beam is primarily lost upstream of the RFQ exit. Numerical simulations indicate that ions passing through the non-linear field region of the solenoids could produce a beam with an increased emittance resulting in up to 50 % of the LEBT beam current failing to meet the RFQ acceptance. An aperture restriction was installed upstream of the first solenoid to remove these ions. This report describes the results of measurements and simulations as well as the LEBT tuning.

Jones, Daniel C.↗

Results of a real-time irradiation of lithium P/N and conventional N/P silicon solar cells.

Eight types of lithium-diffused P/N and three types of conventional 10 ohm-cm N/P silicon solar cells were irradiated at four different temperatures with a strontium-90 radioisotope at a rate typical of that expected in earth orbit. The six-month irradiation confirmed earlier accelerator results, showed that certain cell types outperform others at the various temperatures, and, in general, verified the recent improvements and potential usefulness of lithium solar cells. The experimental approach and statistical methods and analyses employed yielded increased confidence in the validity of the results. Injection level effects were observed to be significant.

Reynard, D. L.↗

Orion Heat Shield Manufacturing Producibility Improvements for the EM-1 Flight Test Program

This paper describes how the Orion program is incorporating improvements in the heat shield design and manufacturing processes reducing programmatic risk and ensuring crew safety in support of NASA’s Exploration Missions. The Orion program successfully completed the EFT-1 flight test in 2014 and is currently developing the EM-1 spacecraft to meet the test objectives of an orbital mission to the moon and return to earth in 2019. Lessons learned from the EFT-1 manufacturing and flight test experience are being incorporated into a wide variety of vehicle systems and manufacturing processes to reduce risk to the Orion missions and flight crew. A critical contributor to crew safety is the heat shield that protects the crew capsule during re-entry through the earth’s atmosphere for return from deep space. The first flight test vehicle, EFT-1, was manufactured and tested in the Neil Armstrong Operations and Checkout (O&C) facility at KSC to demonstrate early risk reduction including the functionality of the Thermal Protection Systems (TPS) for capsule safe return to earth. The approach for the EFT-1 heat shield utilized a low risk approach using Apollo heritage design and manufacturing processes using an Avcoat TPS ablator with a honeycomb substrate to provide a one piece heat shield to meet the mission re-entry heating environments. The manufacturing processes used honeycomb cell injection guns which were redeveloped from the Apollo Lunar Program processes to build the EFT-1 heat shield. The completed heat shield was transported across the country by aircraft to the O&C at KSC for installation onto the capsule. The EFT-1 heat shield successfully performed its mission and experienced ~80% of the re-entry velocity (50% heating rate) for a lunar return for an Exploration Mission. The second flight test vehicle is the EM-1 mission which will have additional flight systems installed to fly to the moon and return. Heat shield design and producibility improvements have been incorporated in the EM-1 vehicle to meet deep space mission and programmatic requirements. The design continues to use the Avcoat material, but in a “block” configuration to enable improvements in the application processes as well as additional improvements in the carrier structure design and manufacturing operations. Incorporating flight test results and producibility improvements from EFT-1 for the heat shield system design and processes have improved the thermal protection capability, improved the producibility, and cost for EM-1 flight test.

Producibility↗

Performance of the Fermilab Linac Injector

The Fermilab Linac injection line consists of a 35 keV magnetron-type H- ion source, two-solenoid Low Energy Beam Transport (LEBT), 201 MHz 4-rod 750 keV Radio Frequency Quadrupole (RFQ), and a Medium Energy Transport (MEBT) containing 4 quadrupoles and a bunching cavity. The injector delivers 25 mA, 48 s pulses to drift-tube Linac at a repetition rate of 15Hz. The transmission efficiency has been lower than expected since commissioning. Recent beam current measurements suggest that the beam is primarily lost upstream of the RFQ exit. Numerical simulations indicate that ions passing through the non-linear field region of the solenoids could produce a beam with an increased emittance resulting in up to 50 % of the LEBT beam current failing to meet the RFQ acceptance. An aperture restriction was installed upstream of the first solenoid to remove these ions. This report describes the results of measurements and simulations as well as the LEBT tuning.

Jones, Daniel C.↗

Effect of hydrogen isotope on plasma impedance and thermal pinching induced by single aerosol droplets injected into an inductively coupled plasma

The mass of the hydrogen isotope in micrometer-sized aqueous droplets introduced into an ICP has been found to affect the degree of resulting plasma perturbation. Studied plasma perturbations include local plasma cooling, plasma shrinkage due to thermal pinching, and changes in plasma impedance. The local cooling effect caused by a vaporizing aerosol droplet was found to be similar between H 2 O and D 2 O. This finding is likely a result of the similar heat capacities and vaporization rates of H 2 O and D 2 O droplets and of the comparable thermal conductivities and dissociation kinetics of H 2 O and D 2 O vapor. In contrast, a clear isotope effect was observed for plasma shrinkage and impedance change, with H 2 O causing a considerably stronger effect. Plasma shrinkage, gauged by the reduction in Ar emission at the edge of the plasma, was detected at all observation heights within the load-coil region, with D 2 O droplets causing as much as 30% less of an effect than H 2 O. Likewise, D 2 O droplets induce significantly less change in plasma impedance than H 2 O droplets. Depending on experimental conditions and the probing method, changes in plasma impedance caused by monodisperse D 2 O droplets were only 52%–80% that of H 2 O. The lower influence on plasma shrinkage and impedance change for D 2 O droplets is believed to be related to the much lower thermal conductivity of atomic deuterium, around 70% that of atomic hydrogen (protium). The time required for the plasma impedance to return to its original steady-state value was comparatively long (~10 ms) and similar for H 2 O and D 2 O droplets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of a latent fault detector in a digital system

Methods of modeling the detection time or latency period of a hardware fault in a digital system are proposed that explain how a computer detects faults in a computational mode. The objectives were to study how software reacts to a fault, to account for as many variables as possible affecting detection and to forecast a given program's detecting ability prior to computation. A series of experiments were conducted on a small emulated microprocessor with fault injection capability. Results indicate that the detecting capability of a program largely depends on the instruction subset used during computation and the frequency of its use and has little direct dependence on such variables as fault mode, number set, degree of branching and program length. A model is discussed which employs an analog with balls in an urn to explain the rate of which subsequent repetitions of an instruction or instruction set detect a given fault.

Nagel, P. M.↗

Parametric and time resolved studies of autoignition and flameholding in a clean-air supersonic combustor

A model scramjet combustor was tested using transverse sonic injection of gaseous hydrogen behind a rearward facing step. It was found that ignition in the constant area duct initiates well downstream and propagates upstream with velocities ranging from 300 to nearly 1,000 m/s in the laboratory frame. The downstream ignition develops into a detonation wave prior to reaching the region of observation. A number of parameters including fuel flow rate, air stagnation temperature and pressure, and the operating conditions of the downstream air ejector affected ignition conditions.

Whitehurst, R. B.↗

The Impact of Geoengineering Aerosols on Stratospheric Temperature and Ozone

Anthropogenic greenhouse gas emissions are warming the global climate at an unprecedented rate. Significant emission reductions will be required soon to avoid a rapid temperature rise. As a potential interim measure to avoid extreme temperature increase, it has been suggested that Earth's albedo be increased by artificially enhancing stratospheric sulfate aerosols. We use a 3D chemistry climate model, fed by aerosol size distributions from a zonal mean aerosol model, to simulate continuous injection of 1-10 Mt/a into the lower tropical stratosphere. In contrast to the case for all previous work, the particles are predicted to grow to larger sizes than are observed after volcanic eruptions. The reason is the continuous supply of sulfuric acid and hence freshly formed small aerosol particles, which enhance the formation of large aerosol particles by coagulation and, to a lesser extent, by condensation. Owing to their large size, these particles have a reduced albedo. Furthermore, their sedimentation results in a non-linear relationship between stratospheric aerosol burden and annual injection, leading to a reduction of the targeted cooling. More importantly, the sedimenting particles heat the tropical cold point tropopause and, hence, the stratospheric entry mixing ratio of H2O increases. Therefore, geoengineering by means of sulfate aerosols is predicted to accelerate the hydroxyl catalyzed ozone destruction cycles and cause a significant depletion of the ozone layer even though future halogen concentrations will be significantly reduced.

Heckendorn, P.↗

The Impact of Geoengineering Aerosols on Stratospheric Temperature and Ozone

Anthropogenic greenhouse gas emissions are warming the global climate at an unprecedented rate. Significant emission reductions will be required soon to avoid a rapid temperature rise. As a potential interim measure to avoid extreme temperature increase, it has been suggested that Earth's albedo be increased by artificially enhancing stratospheric sulfate aerosols. We use a 3D chemistry climate model, fed by aerosol size distributions from a zonal mean aerosol model. to simulate continuous injection of 1-10 Mt/a into the lower tropical stratosphere. In contrast to the case for all previous work, the particles are predicted to grow to larger sizes than are observed after volcanic eruptions. The reason is the continuous supply of sulfuric acid and hence freshly formed small aerosol particles, which enhance the formation of large aerosol particles by coagulation and, to a lesser extent, by condensation. Owing to their large size, these particles have a reduced albedo. Furthermore, their sedimentation results in a non-linear relationship between stratospheric aerosol burden and annual injection, leading to a reduction of the targeted cooling. More importantly, the sedimenting particles heat the tropical cold point tropopause and, hence, the stratospheric entry mixing ratio of H2O increases. Therefore, geoengineering by means of sulfate aerosols is predicted to accelerate the hydroxyl catalyzed ozone destruction cycles and cause a significant depletion of the ozone layer even though future halogen concentrations will he significantly reduced.

Heckendorn, P.↗

GUI for Computational Simulation of a Propellant Mixer

Control Panel is a computer program that generates a graphical user interface (GUI) for computational simulation of a rocket-test-stand propellant mixer in which gaseous hydrogen (GH2) is injected into flowing liquid hydrogen (LH2) to obtain a combined flow having desired thermodynamic properties. The GUI is used in conjunction with software that models the mixer as a system having three inputs (the positions of the GH2 and LH2 inlet valves and an outlet valve) and three outputs (the pressure inside the mixer and the outlet flow temperature and flow rate). The user can specify valve characteristics and thermodynamic properties of the input fluids via userfriendly dialog boxes. The user can enter temporally varying input values or temporally varying desired output values. The GUI provides (1) a set-point calculator function for determining fixed valve positions that yield desired output values and (2) simulation functions that predict the response of the mixer to variations in the properties of the LH2 and GH2 and manual- or feedback-control variations in valve positions. The GUI enables scheduling of a sequence of operations that includes switching from manual to feedback control when a certain event occurs.

Figueroa, Fernando↗

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: III. Novel sub-cyclotron modes

Dedicated DIII-D experiments studied instabilities above and below the ion cyclotron frequency by injecting thirteen different neutral beam populations into L-mode plasmas with different thermal compositions of H, D, and 3 He at five different values of toroidal field. Two companion papers describe the results of these experiments in detail. This paper focuses on data and analysis of two instabilities that occur at ${\sim}0.6$ of the deuterium and hydrogen cyclotron frequencies $f_\mathrm {cD}$ and $f_\mathrm {cH}$ in mixed hydrogen and deuterium plasmas. The modes below $f_\mathrm {cD}$ are often driven unstable by deuterium beams but are rarely destabilized by hydrogen beams. The modes below $f_\mathrm {CD}$ are reported here for the first time in a fusion plasma. These $\mathrm{sub-}f_\mathrm {cH}$ modes are only driven by deuterium beams but, unlike the $\mathrm{sub-}f_\mathrm {cD}$ modes, they are only unstable in plasmas with deuterium concentrations that are ${\lesssim}20\%$; they also have a smaller amplitude than the sub-deuterium modes. Addition of fully-ionized 3 He to the plasma has little effect on the sub-cyclotron deuterium modes but is stabilizing for the sub-cyclotron hydrogen modes. The modes are excited in the core and have toroidal mode numbers around 20. These trends are compared with calculations of the growth rate in a simplified model. The theory successfully predicts global Alfvén eigenmode instability with frequencies and mode numbers close to experimental values, the stronger amplitude of the sub-deuterium modes, and the species-mix dependence of the sub-hydrogen modes. Increased damping associated with the ion-ion hybrid, cutoff, and crossover frequencies plausibly accounts for the stabilization of the $\mathrm{sub-}f_\mathrm {cH}$ modes when 3 He is puffed. However, the relative stability of $\mathrm{sub-}f_\mathrm {cH}$ modes during hydrogen injection remains unexplained.

electromagnetic ion cyclotron instability↗

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗