Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interactive materials 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 253 records · Page 14

Source of the optical red-slope in iron-rich meteorites

The relationship between ordinary chondrites and S-type asteroids is an unresolved issue in meteorite science. S-type asteroids exhibit a positively red-sloped spectrum that is interpreted to indicate the presence of elemental iron on the surfaces. The characteristic red-sloped spectrum of iron-rich meteorites is produced by only the specular component of the reflectance. Complex metallic surfaces can be modeled as linear mixtures of specular and nonspecular components. It is the geometry of the metal on a surface and its interaction with surrounding material, rather than the absolute amount of metal, that determine the redness of resulting spectra. In order to distinguish between ordinary chondrite and differentiated parent bodies it is important to understand how regolith processes affect the nature and form of metal on asteroid surfaces.

Britt, Daniel T.↗

Explosive mafic volcanism on Earth and Mars

Deposits within Amazonia Planitia, Mars, have been interpreted as ignimbrite plains on the basis of their erosional characteristics. The western flank of Hecates Tholus appears to be mantled by an airfall deposit, which was produced through magma-water interactions or exsolution of magmatic volatiles. Morphologic studies, along with numerical and analytical modeling of Martian plinian columns and pyroclastic flows, suggest that shield materials of Tyrrhena and Hadriaca paterae are composed of welded pyroclastic flows. Terrestrial pyroclastic flows, ignimbrites, and airfall deposits are typically associated with silicic volcanism. Because it is unlikely that large volumes of silicic lavas have been produced on Mars, we seek terrestrial analogs of explosives, mafic volcanism. Plinian basaltic airfall deposits have been well-documented at Masaya, Nicaragua, and basaltic ignimbrite and surge deposits also have been recognized there. Ambrym and Yasour, both in Vanuatu, are mafic stratovolcanioes with large central calderas, and are composed of interbedded basaltic pyrocalstic deposits and lava flows. Zavaritzki, a mafic stratovolcano in the Kurile Islands, may have also produced pyroclastic deposits, although the exact nature of these deposits in unknown. Masaya, Ambrym and Yasour are known to be located above tensional zones. Hadriaca and Tyrrhena Paterae may also be located above zones of tension, resulting from the formation and evolution of Hellas basin, and, thus, may be directly analogous to these terrestrial mafic, explosive volcanoes.

Gregg, Tracy K. P.↗

Containerless, Low-Gravity Undercooling of Ti-Ce Alloys in the MSFC Drop Tube

Previous tests of the classical nucleation theory as applied to liquid-liquid gap miscibility systems found a discrepancy between experiment and theory in the ability to undercool one of the liquids before the L1-L2 separation occurs. To model the initial separation process in a two-phase liquid mixture, different theoretical approaches, such as free-energy gradient and density gradient theories, have been put forth. If there is a large enough interaction between the critical liquid and the crucible, both models predict a wetting temperature (T(sub w)) above which the minority liquid perfectly wets and layers the crucible interface, but only on one side of the immiscibility dome. Materials with compositions on the other side of the dome will have simple surface adsorption by the minority liquid before bulk separation occurs when the coexistence (i.e., binoidal) line in reached. If the interaction between the critical liquid and the crucible were to decrease, T(sub w) would increase, eventually approaching the critical consolute temperature (T(sub cc)). If this situation occurs, then there could be large regions of the miscibility gap where non-perfect wetting conditions prevail resulting in droplets of L1 liquid at the surface having a non-zero contact angle. The resulting bulk structure will then depend on what happens on the surface and the subsequent processing conditions. In the past several decades, many experiments in space have been performed on liquid metal binary immiscible systems for the purpose of determining the effects that different crucibles may have on the wetting and separation process of the liquids. Potard performed experiments that showed different crucible materials could cause the majority phase to preferentially wet the container and thus produce a dispersed microstructure of the minority phase. Several other studies have been performed on immiscibles in a semi-container environment using an emulsion technique. Only one previous study was performed using completely containerless processing of immiscible metals and the results of that investigation are similar to some of the emulsion studies. In all the studies, surface wetting was attributed as the cause for the similar microstructures or the asymmetry in the ability to undercool the liquid below the binoidal on one side of the immiscibility dome. By removing the container completely from the separation process, it was proposed that the loss of the crucible/liquid interaction would produce a large shift in T(sub w) and thus change the wetting characteristics at the surface. By investigating various compositions across the miscibility gap, a change in the type and amount of liquid wetting at the surface of a containerless droplet should change the surface nucleating behavior of the droplet - whether it be the liquid-liquid wetting or the liquid-to-solid transition. Undercooling of the liquid into the metastable region should produce significant differences in the separation process and the microstructure upon solidification. In this study, we attempt to measure these transitions by monitoring the temperature of the sample by optical pyrometry. Microstructural analysis will be made to correlate with the degree of undercooling and the separation mechanisms involved.

Robinson, M. B.↗

Event Generators for Simulating Heavy Ion Interactions of Interest in Evaluating Risks in Human Spaceflight

Simulating the Space Radiation environment with Monte Carlo Codes, such as FLUKA, requires the ability to model the interactions of heavy ions as they penetrate spacecraft and crew member's bodies. Monte-Carlo-type transport codes use total interaction cross sections to determine probabilistically when a particular type of interaction has occurred. Then, at that point, a distinct event generator is employed to determine separately the results of that interaction. The space radiation environment contains a full spectrum of radiation types, including relativistic nuclei, which are the most important component for the evaluation of crew doses. Interactions between incident protons with target nuclei in the spacecraft materials and crew member's bodies are well understood. However, the situation is substantially less comfortable for incident heavier nuclei (heavy ions). We have been engaged in developing several related heavy ion interaction models based on a Quantum Molecular Dynamics-type approach for energies up through about 5 GeV per nucleon (GeV/A) as part of a NASA Consortium that includes a parallel program of cross section measurements to guide and verify this code development.

Wilson, Thomas L.↗

Probabilistic constitutive relationships for material strength degradation models

In the present probabilistic methodology for the strength of aerospace propulsion system structural components subjected to such environmentally-induced primitive variables as loading stresses, high temperature, chemical corrosion, and radiation, time is encompassed as an interacting element, allowing the projection of creep and fatigue effects. A probabilistic constitutive equation is postulated to account for the degradation of strength due to these primitive variables which may be calibrated by an appropriately curve-fitted least-squares multiple regression of experimental data. The resulting probabilistic constitutive equation is embodied in the PROMISS code for aerospace propulsion component random strength determination.

Boyce, L.↗

Recent spectroscopic findings concerning clay/water interactions at low humidity: Possible applications to models of Martian surface reactivity

A feasibility study assessing the utility of the adaptation of near infrared correlation spectroscopy to quantifying iron and adsorbed water in some clay-based Mars soil analog materials (MarSAM's). The work was intended to constitute Phase 1 of an approach to identifying optical analytical wavelength regions, not only for important mineral classes, but for chemically active centers within them. Many of these centers are common to unrelated mineral classes and of disproportionate influence relative to the mineral structure as a whole in determining the surface reactivity of mineral surfaces. We previously reported linearity between reflectance and total iron and total moisture over a large range of both key variables. We also discovered interesting relationships between the intensity of iron bands and the relative humidity of the systems. These relationships were confirmed. We also show that, in the low humidity range, reflectance is linearly dependent on a different kind of water from that best representing the full humidity range (the kind of water associated, in clays, with surface acidity). These relationships and the sensitivity and capability of quantitation of near infrared data indicate high promise with the production of reactive surface intermediates of products of surface reactions.

Coyne, L.↗

Nonlinear response of unidirectional boron/aluminum

Experimental results obtained for unidirectional boron/aluminum subjected to combined loading using off-axis tension, compression and Iosipescu shear specimens are correlated with a nonlinear micromechanics model. It is illustrated that the nonlinear response in the principal material directions is markedly influenced by the different loading modes and different ratios of the applied stress components. The observed nonlinear response under pure and combined loading is discussed in terms of initial yielding, subsequent hardening, stress-interaction effects and unloading-reloading characteristics. The micromechanics model is based on the concept of a repeating unit cell representative of the composite-at-large and employs the unified theory of Bodner and Partom to model the inelastic response of the matrix. It is shown that the employed micromechanics model is sufficiently general to predict the observed nonlinear response of unidirectional boron/aluminum with good accuracy.

Pindera, M.-J.↗

Ablative Rocket Deflector Testing and Computational Modeling

A deflector risk mitigation program was recently conducted at the NASA Stennis Space Center. The primary objective was to develop a database that characterizes the behavior of industry-grade refractory materials subjected to rocket plume impingement conditions commonly experienced on static test stands. The program consisted of short and long duration engine tests where the supersonic exhaust flow from the engine impinged on an ablative panel. Quasi time-dependent erosion depths and patterns generated by the plume impingement were recorded for a variety of different ablative materials. The erosion behavior was found to be highly dependent on the material s composition and corresponding thermal properties. For example, in the case of the HP CAST 93Z ablative material, the erosion rate actually decreased under continued thermal heating conditions due to the formation of a low thermal conductivity "crystallization" layer. The "crystallization" layer produced near the surface of the material provided an effective insulation from the hot rocket exhaust plume. To gain further insight into the complex interaction of the plume with the ablative deflector, computational fluid dynamic modeling was performed in parallel to the ablative panel testing. The results from the current study demonstrated that locally high heating occurred due to shock reflections. These localized regions of shock-induced heat flux resulted in non-uniform erosion of the ablative panels. In turn, it was observed that the non-uniform erosion exacerbated the localized shock heating causing eventual plume separation and reversed flow for long duration tests under certain conditions. Overall, the flow simulations compared very well with the available experimental data obtained during this project.

Allgood, Daniel C.↗

Neptune's inner magnetosphere and aurora: Energetic particle constraints

A dramatic and peculiar dropout of greater than 500-keV ions (but not electrons) was observed within Neptune's inner magnetosphere near 2 R(sub N) as the Voyager 2 spacecraft approached the planet. Unlike a number of other energetic particle features this feature could not be accounted for by known material bodies in the context of the most utilized magnetic field models (neither the offset tilted dipole models nor the spehrical harmonic model 'O8'). However, the configuration of Neptune's inner magnetosphere is highly uncertain. By applying a novel technique, utilizing energetic particle measurements, to constrain the magnetic field configuration of the inner regions, we show that appeals to unobserved materials within Neptune's system are unnecessary, and that the ion dropout feature was, in all likelihood, the result of ion interactions with maximum L excursions of the ring 1989N1R. The constraints also favor the se of the M2 magnetic field model (Selesnick, 1992) over the previous models. An electron feature was probably absent because the electron interactions with the ring occurred substantially before the ion interactions (about 2 hours for the electrons versus a few minutes for the ions). Pitch-angle scattering apparently eliminated the electron signature. Minimum scattering rates determined based on this premise yield enough electron precipitation power to explain the brightest component of Neptune's aurora. We propose that this bright component is analogous to the Earth's diffuse aurora.

Mauk, B. H.↗

Automated technologies needed to prevent radioactive materials from reentering the atmosphere

Project SIREN (Search, Intercept, Retrieve, Expulsion Nuclear) was created to identify and evaluate the technologies and operational strategies needed to rendezvous with and capture aerospace radioactive materials (e.g., a distressed or spent space reactor core) before such materials can reenter the terrestrial atmosphere and then to safely move these captured materials to an acceptable space destination for proper disposal. A major component of the current project SIREN effort is the development of an interactive technology model (including a computerized data base) that explores, in building-block fashion, the interaction of the technologies and procedures needed to successfully accomplish a SIREN mission. The SIREN model will include appropriate national and international technology elements--both contemporary and projected into the next century. To obtain maximum flexibility and use, the SIREN technology data base is being programmed for use on 286-class PC's. The major technical elements for a successful SIREN mission include: ground and space-based tracking, launch vehicles of needed payload capacity, telerobotic systems, sensors, capture technologies, and space transport and disposal. However, Project SIREN also will impose specialized requirements including the use of dextrous aerospace systems capable of properly functioning in intense radiation and thermal environments.

Buden, David↗

Automated technologies needed to prevent radioactive materials from reentering the atmosphere

Project SIREN (Search, Intercept, Retrieve, Expulsion Nuclear) has been created to identify and evaluate the technologies and operational strategies needed to rendezvous with and capture aerospace radioactive materials (e.g., a distressed or spent space reactor core) before such materials can reenter the terrestrial atmosphere and then to safely move these captured materials to an acceptable space destination for proper disposal. A major component of the current Project SIREN effort is the development of an interactive technology model (including a computerized data base) that explores in building block fashion the interaction of the technologies and procedures needed to successfully accomplish a SIREN mission. This SIREN model will include appropriate national and international technology elements-both contemporary and projected into the next century. To permit maximum flexibility and use, the SIREN technology data base is being programmed for use on 386-class PC's.

Buden, David↗

Coupled electromechanical response of composite beams with embedded piezoelectric sensors and actuators

Unified mechanics are developed with the capability to model both sensory and active composite laminates with embedded piezoelectric layers. A discrete-layer formulation enables analysis of both global and local electromechanical response. The mechanics include the contributions from elastic, piezoelectric, and dielectric components. The incorporation of electric potential into the state variables permits representation of general electromechanical boundary conditions. Approximate finite element solutions for the static and free-vibration analysis of beams are presented. Applications on composite beams demonstrate the capability to represent either sensory or active structures and to model the complicated stress-strain fields, the interactions between passive/active layers, interfacial phenomena between sensors and composite plies, and critical damage modes in the material. The capability to predict the dynamic characteristics under various electrical boundary conditions is also demonstrated.

Saravanos, D. A.↗

Thermoviscoplastic analysis of engine cowl leading edge subjected to oscillating shock-shock interaction

A finite element thermoviscoplastic analysis method, which employs a unified constitutive model proposed by Bodner and Partom, is used to predict rate-dependent nonlinear structural behavior. The method is evaluated by predicting stress-strain behavior of a uniaxially loaded bar of nickel-based superalloy (B1900 + Hf) material. The method is used to predict the time-dependent thermoviscoplastic response of a B1900 + Hf leading edge subjected to oscillating shock-shock interaction loading. Viscoplastic analysis shows that the leading edge experiences significant plastic straining. The plastic region increases with cyclic loading in the high heat flux area.

Pandey, Ajay K.↗

Modeling of Substitutional Site Preference in Ordered Intermetallic Alloys

We investigate the site substitution scheme of specific alloying elements in ordered compounds and the dependence of site occupancy on compound stoichiometry, alloy concentration. This basic knowledge, and the interactions with other alloying additions are necessary in order to predict and understand the effect of various alloying schemes on the physical properties of a material, its response to various temperature treatments, and the resulting mechanical properties. Many theoretical methods can provide useful but limited insight in this area, since most techniques suffer from constraints in the type of elements and the crystallographic structures that can be modeled. With this in mind, the Bozzolo-Ferrante-Smith (BFS) method for alloys was designed to overcome these limitations, with the intent of providing an useful tool for the theoretical prediction of fundamental properties and structure of complex systems. After a brief description of the BFS method, its use for the determination of site substitution schemes for individual as well as collective alloying additions to intermetallic systems is described, including results for the concentration dependence of the lattice parameter. Focusing on B2 NiAl, FeAl and CoAl alloys, the energetics of Si, Ti, V, Cr, Fe, Co, Ni, Cu, Zr, Nb, Mo, Ru, Hf, Ta and W alloying additions are surveyed. The effect of single additions as well as the result of two simultaneous additions, discussing the interaction between additions and their influence on site preference schemes is considered. Finally, the BFS analysis is extended to ternary L1(sub 2) (Heusler phase) alloys. A comparison between experimental and theoretical results for the limited number of cases for which experimental data is available is also included.

Bozzolo, Guillermo↗

Characterizing Wheel-Soil Interaction Loads Using Meshfree Finite Element Methods: A Sensitivity Analysis for Design Trade Studies

A wheel experiencing sinkage and slippage events poses a high risk to planetary rover missions as evidenced by the mobility challenges endured by the Mars Exploration Rover (MER) project. Current wheel design practice utilizes loads derived from a series of events in the life cycle of the rover which do not include (1) failure metrics related to wheel sinkage and slippage and (2) performance trade-offs based on grouser placement/orientation. Wheel designs are rigorously tested experimentally through a variety of drive scenarios and simulated soil environments; however, a robust simulation capability is still in development due to myriad of complex interaction phenomena that contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree nite element approaches enable simulations that capture su cient detail of wheel-soil interaction while remaining computationally feasible. This study implements the JPL wheel-soil benchmark problem in the commercial code environment utilizing the large deformation modeling capability of Smooth Particle Hydrodynamics (SPH) meshfree methods. The nominal, benchmark wheel-soil interaction model that produces numerically stable and physically realistic results is presented and simulations are shown for both wheel traverse and wheel sinkage cases. A sensitivity analysis developing the capability and framework for future ight applications is conducted to illustrate the importance of perturbations to critical material properties and parameters. Implementation of the proposed soil-wheel interaction simulation capability and associated sensitivity framework has the potential to reduce experimentation cost and improve the early stage wheel design proce

terramechanics↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Tidal heating and mass loss in neutron star binaries - Implications for gamma-ray burst models

A neutron star in a close binary orbit around another neutron star (or stellar-mass black hole) spirals inward owing to gravitational radiation. We discuss the effects of tidal dissipation during this process. Tidal energy dissipated in the neutron star's core escapes mainly as neutrinos, but heating of the crust, and outward diffusion of photons, blows off the outer layers of the star. This photon-driven mass loss precedes the final coalescence. The presence of this eject material impedes the escape of gamma-rays created via neutrino interactions. If an e(+) - e(-) fireball, created in the late stages of coalescence, were loaded with (or surrounded by) material with the mean column density of the ejecta, it could not be an efficient source of gamma-rays. Models for cosmologically distant gamma-rays burst that involve neutron stars must therefore be anisotropic, so that the fireball expands preferentially in directions where the column density of previously blown-off material is far below the spherically averaged value which we have calculated. Some possible 'scenarios' along these lines are briefly discussed.

Meszaros, P.↗