Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatially explicit modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Grand challenge scientific questions in coupled modeling

Most convective field experiments in the past (e.g., SESAME, CCOPE, CINDE) have attempted to resolve only the immediate scales of moist convection using network arrays that spanned two or three atmospheric scales at most. Furthermore, these scales have been defined more on practical considerations (cost, manpower, etc.) than on a clear understanding of their theoretical significance. Unfortunately, this has precluded a description of the entire life cycle of MCS's and their interaction with larger scale systems, the land surface, and trace species. Fortunately, the following factors now make it possible to attempt to simulate scale contraction processes from the synoptic scale down to the cloud scale, as well as interactions between complex meteorological, land surface, precipitation, chemical, and hydrologic processes with coupled, multiscale models: the availability of new technology to sample meteorological fields at high temporal and spatial resolution over a broad region made possible by the weather observing modernization program; increased computer power and improved numerical approaches to run limited area models with nonhydrostatic precipitation physics so as to explicitly resolve MCS (Mesoscale Convective System) processes; and four dimensional assimilation of non-conventional data to provide dynamically consistent datasets for diagnostic analysis of nonlinear scale-interactive dynamics. Several examples of scale-interactive processes which present grand challenges for coupled, multiscale modeling were presented.

Koch, Steven↗
NASA NTRS2009

Sensitivity of a Cloud-Resolving Model to Bulk and Explicit Bin Microphysical Schemes: Cloud Microphysics and Storm Dynamics Interactions - Part 2

Part I of this paper compares two simulations, one using a bulk and the other a detailed bin microphysical scheme, of a long-lasting, continental mesoscale convective system with leading convection and trailing stratiform region. Diagnostic studies and sensitivity tests are carried out in Part II to explain the simulated contrasts in the spatial and temporal variations by the two microphysical schemes and to understand the interactions between cloud microphysics and storm dynamics. It is found that the fixed raindrop size distribution in the bulk scheme artificially enhances rain evaporation rate and produces a stronger near surface cool pool compared with the bin simulation. In the bulk simulation, cool pool circulation dominates the near-surface environmental wind shear in contrast to the near-balance between cool pool and wind shear in the bin simulation. This is the main reason for the contrasting quasi-steady states simulated in Part I. Sensitivity tests also show that large amounts of fast-falling hail produced in the original bulk scheme not only result in a narrow trailing stratiform region but also act to further exacerbate the strong cool pool simulated in the bulk parameterization. An empirical formula for a correction factor, r(q(sub r)) = 0.11q(sub r)(exp -1.27) + 0.98, is developed to correct the overestimation of rain evaporation in the bulk model, where r is the ratio of the rain evaporation rate between the bulk and bin simulations and q(sub r)(g per kilogram) is the rain mixing ratio. This formula offers a practical fix for the simple bulk scheme in rain evaporation parameterization.

Li, Xiaowen↗
NASA NTRS1984

Aberrations of varied line-space grazing incidence gratings in converging light beams

Analyses of the imaging properties of several designs for varied-line space gratings in converging beams of light in grazing-incidence spectrometers are presented. An explicit model is defined for the case of a plane-reflection grating intercepting light that converges and is reflected to a stigmatic point associated with the zero-order image of the grating. Smooth spatial variation of the grating constant then permits aberration correction. The aberrations are expressed as polynomials in the grating lens coordinates using power series expansions. Application of the model is illustrated in terms of aberrations experienced with the short wavelength spectrometer on the EUVE satellite. Attention is given to straight and parallel in-plane grooves, curved groove in-plane designs and off-plane grooves. Aberrations due to dispersions and misalignment are also considered.

Hettrick, M. C.↗
NASA NTRS1987

Efficiency of a model human image code

Hypothetical schemes for neural representation of visual information can be expressed as explicit image codes. Here, a code modeled on the simple cells of the primate striate cortex is explored. The Cortex transform maps a digital image into a set of subimages (layers) that are bandpass in spatial frequency and orientation. The layers are sampled so as to minimize the number of samples and still avoid aliasing. Samples are quantized in a manner that exploits the bandpass contrast-masking properties of human vision. The entropy of the samples is computed to provide a lower bound on the code size. Finally, the image is reconstructed from the code. Psychophysical methods are derived for comparing the original and reconstructed images to evaluate the sufficiency of the code. When each resolution is coded at the threshold for detection artifacts, the image-code size is about 1 bit/pixel.

Watson, Andrew B.↗
NASA NTRS2000

A Comparison of Three Programming Models for Adaptive Applications

We study the performance and programming effort for two major classes of adaptive applications under three leading parallel programming models. We find that all three models can achieve scalable performance on the state-of-the-art multiprocessor machines. The basic parallel algorithms needed for different programming models to deliver their best performance are similar, but the implementations differ greatly, far beyond the fact of using explicit messages versus implicit loads/stores. Compared with MPI and SHMEM, CC-SAS (cache-coherent shared address space) provides substantial ease of programming at the conceptual and program orchestration level, which often leads to the performance gain. However it may also suffer from the poor spatial locality of physically distributed shared data on large number of processors. Our CC-SAS implementation of the PARMETIS partitioner itself runs faster than in the other two programming models, and generates more balanced result for our application.

Shan, Hong-Zhang↗
NASA NTRS1986

On the nature of models in remote sensing

An explicit framework can provide a better understanding of remote sensing models and their interrelationships. This framework distinguishes between the scene, which is real and exists on the ground, and the image, which is a collection of spatially arranged masurements drawn from the scene. The scene model generalizes and parameterizes the essential qualities of the scene. Scene models may be discrete, in which the scene model consists of discrete elements with boundaries, or continuous, in which matter and energy flows are taken to be continuous and there are no clear or sharp boundaries in the scene. In the discrete case, there are two possibilities for models: H- and L-resolution. In the H-resolution case, the resolution cells of the image are smaller than the elements, and thus the elements may be individually resolved. In the L-resolution case, the resolution cells are larger than the elements and cannot be resolved. Most canopy models are L-resolution, deterministic, and noninvertible in nature; image processing models, however, tend to be H-resolution, empirical, and invertible. This taxonomy helps add insight to the development of remote sensing theory and point the way to new, productive areas of research.

Strahler, A. H.↗
NASA NTRS2012

Understanding Climate Policy Data Needs. NASA Carbon Monitoring System Briefing: Characterizing Flux Uncertainty, Washington D.C., 11 January 2012

Climate policy in the United States is currently guided by public-private partnerships and actions at the local and state levels. This mitigation strategy is made up of programs that focus on energy efficiency, renewable energy, agricultural practices and implementation of technologies to reduce greenhouse gases. How will policy makers know if these strategies are working, particularly at the scales at which they are being implemented? The NASA Carbon Monitoring System (CMS) will provide information on carbon dioxide fluxes derived from observations of earth's land, ocean and atmosphere used in state of the art models describing their interactions. This new modeling system could be used to assess the impact of specific policy interventions on CO2 reductions, enabling an iterative, results-oriented policy process. In January of 2012, the CMS team held a meeting with carbon policy and decision makers in Washington DC to describe the developing modeling system to policy makers. The NASA CMS will develop pilot studies to provide information across a range of spatial scales, consider carbon storage in biomass, and improve measures of the atmospheric distribution of carbon dioxide. The pilot involves multiple institutions (four NASA centers as well as several universities) and over 20 scientists in its work. This pilot study will generate CO2 flux maps for two years using observational constraints in NASA's state-of -the-art models. Bottom-up surface flux estimates will be computed using data-constrained land and ocean models; comparison of the different techniques will provide some knowledge of uncertainty in these estimates. Ensembles of atmospheric carbon distributions will be computed using an atmospheric general circulation model (GEOS-5), with perturbations to the surface fluxes and to transport. Top-down flux estimates will be computed from observed atmospheric CO2 distributions (ACOS/GOSAT retrievals) alongside the forward-model fields, in conjunction with an inverse approach based on the CO2 model of GEOS ]Chem. The forward model ensembles will be used to build understanding of relationships among surface flux perturbations, transport uncertainty and atmospheric carbon concentration. This will help construct uncertainty estimates and information on the true spatial resolution of the top-down flux calculations. The relationship between the top-down and bottom-up flux distributions will be documented. Because the goal of NASA CMS is to be policy relevant, the scientists involved in the flux modeling pilot need to understand and be focused on the needs of the climate policy and decision making community. If policy makers are to use CMS products, they must be aware of the modeling effort and begin to design policies that can be evaluated with information. Improving estimates of carbon sequestered in forests, for example, will require information on the spatial variability of forest biomass that is far more explicit than is presently possible using only ground observations. Carbon mitigation policies being implemented by cities around the United States could be designed with the CMS data in mind, enabling sequential evaluation and subsequent improvements in incentives, structures and programs. The success of climate mitigation programs being implemented in the United States today will hang on the depth of the relationship between scientists and their policy and decision making counterparts. Ensuring that there is two-way communication between data providers and users is important for the success both of the policies and the scientific products meant to support them..

Brown, Molly E.↗
NASA NTRS2024

Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations

Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.

Matthew D Duez↗
NASA NTRS1993

A time domain, weighted residual formulation of Maxwell's equations

A finite element model is developed and used to simulate two-dimensional electromagnetic wave propagation and scattering. The spatial discretization of the time-domain electrodynamic equations is accomplished by a Galerkin approach. The semi-discrete equations are solved explicitly using a second-order Runge-Kutta scheme. Both the electric and magnetic fields are discretized using a single grid, with the divergence-free conditions satisfied through a correction approach. Examples depicting the scattering of plane waves in 2D geometries are given to demonstrate the validity of the methodology.

Young, Jeffrey L.↗
NASA NTRS2018

Estimation of Terrestrial Global Gross Primary Production (GPP) with Satellite Data-Driven Models and Eddy Covariance Flux Data

We estimate global terrestrial gross primary production (GPP) based on models that use satellite data within a simplified light-use efficiency framework that does not rely upon other meteorological inputs. Satellite-based geometry-adjusted reflectances are from the MODerate-resolution Imaging Spectroradiometer (MODIS) and provide information about vegetation structure and chlorophyll content at both high temporal (daily to monthly) and spatial (1 km) resolution. We use satellite-derived solar-induced fluorescence (SIF) to identify regions of high productivity crops and also evaluate the use of downscaled SIF to estimate GPP. We calibrate a set of our satellite-based models with GPP estimates from a subset of distributed eddy covariance flux towers (FLUXNET 2015). The results of the trained models are evaluated using an independent subset of FLUXNET 2015 GPP data. We show that variations in light-use efficiency (LUE) with incident PAR are important and can be easily incorporated into the models. Unlike many LUE-based models, our satellite-based GPP estimates do not use an explicit parameterization of LUE that reduces its value from the potential maximum under limiting conditions such as temperature and water stress. Even without the parameterized downward regulation, our simplified models are shown to perform as well as or better than state-of-the-art satellite data-driven products that incorporate such parameterizations. A significant fraction of both spatial and temporal variability in GPP across plant functional types can be accounted for using our satellite-based models. Our results provide an annual GPP value of 140 Pg C year 1 for 2007 that is within the range of a compilation of observation-based, model, and hybrid results, but is higher than some previous satellite observation-based estimates

CO2↗
NASA NTRS1976

Wavelength dependence in radio-wave scattering and specular-point theory

Radio-wave scattering from natural surfaces contains a strong quasispecular component that at fixed wavelengths is consistent with specular-point theory, but often has a strong wavelength dependence that is not predicted by physical optics calculations under the usual limitations of specular-point models. Wavelength dependence can be introduced by a physical approximation that preserves the specular-point assumptions with respect to the radii of curvature of a fictitious, effective scattering surface obtained by smoothing the actual surface. A uniform low-pass filter model of the scattering process yields explicit results for the effective surface roughness versus wavelength. Interpretation of experimental results from planetary surfaces indicates that the asymptotic surface height spectral densities fall at least as fast as an inverse cube of spatial frequency. Asymptotic spectral densities for Mars and portions of the lunar surface evidently decrease more rapidly.

Tyler, G. L.↗
NASA NTRS2006

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗
NASA NTRS2024

Explicitly Resolving Lightning and Electrification Processes from the 10-12 April 2019 Thundersnow Outbreak

The 10-12 April 2019 thundersnow (i.e., lightning within snowfall) outbreak was examined via ground- and space-based lightning observations and was simulated using a numerical weather prediction model with an explicit electrification parameterization. When compared to observations, the simulation propagated the synoptic snowband two to six hours faster while also exaggerating the 3-D reflectivity structure. Throughout the event, the simulation produced 1,733 thundersnow flashes which was less than what was observed by ground- and space-based lightning sensors. In general, simulated thundersnow flashes were spatially offset from the largest reflectivities within the synoptic snowband and tended to occur within elevated convection that traversed isentropically along the top of mid-level frontogenesis. These simulated thundersnow flashes were associated with a tripole charge structure with ice/snow hydrometeors contributing most to the main negative charge region. Both simulated and observed thundersnow flashes initiated in conditionally unstable environments. Lastly, a conceptual model was developed to explain the spatial separation between the largest reflectivities in the snowband and the occurrence of thundersnow. It is hypothesized that the spatial offset of thundersnow initiation from the reflectivity cores within the synoptic snowband arose from a thermal circulation – induced by mid-level frontogenesis – that advects positively charged ice/snow hydrometeors towards the surface and creates a nearly homogeneous vertical charge structure.

lightning↗
NASA NTRSDate not supplied

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗
NASA NTRS1976

Some factors influencing radiation of sound from flow interaction with edges of finite surfaces

Edges of surfaces which are exposed to unsteady flow cause both strictly acoustic effects and hydrodynamic effects, in the form of generation of new hydrodynamic sources in the immediate vicinity of the edge. An analytical model is presented which develops the explicit sound-generation role of the velocity and Mach number of the eddy convection past the edge, and the importance of relative scale lengths of the turbulence, as well as the relative intensity of pressure fluctuations. The Mach number (velocity) effects show that the important paramater is the convection Mach number of the eddies. The effects of turbulence scale lengths, isotropy, and spatial density (separation) are shown to be important in determining the level and spectrum of edge sound radiated for the edge dipole mechanism. Experimental data is presented which provides support for the dipole edge noise model in terms of Mach number (velocity) scaling, parametric dependence on flow field parameter, directivity, and edge diffraction effects.

Hayden, R. E.↗
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗
NASA NTRS1994

Evolution of solitary density waves in stellar winds of early-type stars: A simple explanation of discrete absorption component behavior

We model the evolution of a density shell propagating through the stellar wind of an early-type star, in order to investigate the effects of such shells on UV P Cygni line profiles. Unlike previous treatments, we solve the mass, momentum, and energy conservation equations, using an explicit time-differencing scheme, and present a parametric study of the density, velocity, and temperature response. Under the assumed conditions, relatively large spatial scale, large-amplitude density shells propagate as stable waves through the supersonic portion of the wind. Their dynamical behavior appears to mimic propagating 'solitary waves,' and they are found to accelerate at the same rate as the underlying steady state stellar wind (i.e., the shell rides the wind). These hydrodynamically stable structures quantitatively reproduce the anomalous 'discrete absorption component' (DAC) behavior observed in the winds of luminous early-type stars, as illustrated by comparisons of model predictions to an extensive International Ultraviolet Explorer (IUE) time series of spectra of zeta Puppis (O4f). From these comparisons, we find no conclusive evidence indicative of DACs accelerating at a significantly slower rate than the underlying stellar wind, contrary to earlier reports. In addition, these density shells are found to be consistent within the constraints set by the IR observations. We conclude that the concept of propagating density shells should be seriously reconsidered as a possible explanation of the DAC phenomenon in early-type stars.

Waldron, Wayne L.↗