Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sigma Models”

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 163 records · Page 9

X-ray sources in regions of star formation. 5: The low mass stars of the Upper Scorpius association

We report followup investigations of Einstein x-ray observations of the Upper Scorpius OB association. We identify 28 low mass pre-main sequence stars as counterparts of x-ray sources in the approximately = 7 square degrees of the OB association observed. Based on kinematics and lithium abundances, these stars are low mass members of the OB association. We use optical spectra and optical and near-IR photometry to determine the stellar luminosities, effective temperatures, masses, and ages. We show that the bolometric corrections and effective temperatures of the G and K stars are consistent with those of subgiants. The low mass stars have isochronal ages of 1-2 Myr, depending on the choice of evolutionary models, with very small dispersion (sigma approximately = 1 Myr). This age is significantly younger than the 5-6 Myr found for the more massive B stars. The small dispersion in stellar ages, less than 10% the sound-crossing time of the association, suggests that star formation was triggered. We present two scenarios for star formation in this association. In the two-episode scenario, formation of the low mass stars was triggered by a supernova explosion, and the low mass stars form quickly, with high efficiency. Alternatively, high and low mass star formation was all initiated at the same time, some 5-6 Myr ago, and the apparent systematic age difference is an artifact of how the isochrones are dated. The effect of the supernova is to terminate mass accretion and yield an apparently coeval population. We show that the incompleteness in the x-ray sampling is about 65%, and is strongly dependent on stellar mass. After correction for incompleteness, we estimate there are about 2000 low mass members (stellar mass less than 2 solar mass) of this association. The mass function in this association is indistinguishable from that of the field. The ratio of naked to classical T Tauri stars is much larger than in Tau-Aur, and may be attributable to the local environment. We also present observations of eight ROX sources associated with the rho Oph cloud, and observations of non-pre-main sequence (PMS) stars in our fields.

Walter, Frederick M.↗

Evidence of a Supermassive Black Hole in the Galaxy NGC 1023 From The Nuclear Stellar Dynamics

We analyze the nuclear stellar dynamics of the SBO galaxy NGC 1023, utilizing observational data both from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and from the ground. The stellar kinematics measured from these long-slit spectra show rapid rotation (V equals approx. 70 km/s at a distance of O.1 deg = 4.9 pc from the nucleus) and increasing velocity dispersion toward the nucleus (where sigma = 295 +/- 30 km/s). We model the observed stellar kinematics assuming an axisymmetric mass distribution with both two and three integrals of motion. Both modeling techniques point to the presence of a central dark compact mass (which presumably is a supermassive black hole) with confidence > 99%. The isotropic two-integral models yield a best-fitting black hole mass of (6.0 +/- 0.4) x 10(exp 7) solar masses and mass-to-light ratio (M/L(sub v)) of 5.38 +/- 0.08, and the goodness-of-fit (CHI(exp 2)) is insensitive to reasonable values for the galaxy's inclination. The three-integral models, which non-parametrically fit the observed line-of-sight velocity distribution as a function of position in the galaxy, suggest a black hole mass of (3.9 +/- 0.4) x 10(exp 7) solar masses and M/L(sub v) of 5.56 +/- 0.02 (internal errors), and the edge-on models are vastly superior fits over models at other inclinations. The internal dynamics in NGC 1023 as suggested by our best-fit three-integral model shows that the velocity distribution function at the nucleus is tangentially anisotropic, suggesting the presence of a nuclear stellar disk. The nuclear line of sight velocity distribution has enhanced wings at velocities >= 600 km/s from systemic, suggesting that perhaps we have detected a group of stars very close to the central dark mass.

Bower, G. A.↗

Distributed deformation and block rotation in 3D

The authors address how block rotation and complex distributed deformation in the Earth's shallow crust may be explained within a stationary regional stress field. Distributed deformation is characterized by domains of sub-parallel fault-bounded blocks. In response to the contemporaneous activity of neighboring domains some domains rotate, as suggested by both structural and paleomagnetic evidence. Rotations within domains are achieved through the contemporaneous slip and rotation of the faults and of the blocks they bound. Thus, in regions of distributed deformation, faults must remain active in spite of their poor orientation in the stress field. The authors developed a model that tracks the orientation of blocks and their bounding faults during rotation in a 3D stress field. In the model, the effective stress magnitudes of the principal stresses (sigma sub 1, sigma sub 2, and sigma sub 3) are controlled by the orientation of fault sets in each domain. Therefore, adjacent fault sets with differing orientations may be active and may display differing faulting styles, and a given set of faults may change its style of motion as it rotates within a stationary stress regime. The style of faulting predicted by the model depends on a dimensionless parameter phi = (sigma sub 2 - sigma sub 3)/(sigma sub 1 - sigma sub 3). Thus, the authors present a model for complex distributed deformation and complex offset history requiring neither geographical nor temporal changes in the stress regime. They apply the model to the Western Transverse Range domain of southern California. There, it is mechanically feasible for blocks and faults to have experienced up to 75 degrees of clockwise rotation in a phi = 0.1 strike-slip stress regime. The results of the model suggest that this domain may first have accommodated deformation along preexisting NNE-SSW faults, reactivated as normal faults. After rotation, these same faults became strike-slip in nature.

Scotti, Oona↗

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

Dependence of the Normalized Radar Cross Section of Water Waves on Bragg Wavelength-Wind Speed Sensitivity

Measurements of the normalized radar cross section (sigma(sup o)) made by the YSCAT ultrawideband scatterometer during an extended deployment on the Canada Centre for Inland Waters(CCIW) Research Tower located at Lake Ontario are analyzed and compared with anemometer wind measurements to study the sensitivity of (sigma(sup o)) to the wind speed as a function of the Bragg wavelength. This paper concentrates on upwind and downwind azimuth angles in the wind speed range of 4.5-12 m/s. While YSCAT collected measurements of sigma(sup o) at a variety of frequencies and incidence angles, this paper focuses on frequencies of 2.0, 3.05, 5.30, 10.02, and 14.0 GHz and incidence angles within the Bragg regime, 30-50 deg. Adopting a power law model to describe the relationship between sigma(sup o) and wind speed, both wind speed exponents and upwind/downwind (u/d) ratios of sigma(sup o) are found using least squares linear regression. The analysis of the wind speed exponents and u/d ratios show that shorter Bragg wavelengths (Lambda less than 4 cm) are the most sensitive to wind speed and direction. Additionally, vertical polarization (V-pol) sigma(sup o) is shown to be more sensitive to wind speed than horizontal polarization (H-pol) sigma(sup o), while the H-pol u/d ratio is larger than the V-pol u/d ratio.

Long, David G.↗

Improvements to a Response Surface Thermal Model for Orion Mated to the International Space Station

This study is an extension of previous work to evaluate the applicability of Design of Experiments (DOE)/Response Surface Methodology to on-orbit thermal analysis. The goal was to determine if the methodology could produce a Response Surface Equation (RSE) that predicted the thermal model temperature results within +/-10 F. An RSE is a polynomial expression that can then be used to predict temperatures for a defined range of factor combinations. Based on suggestions received from the previous work, this study used a model with simpler geometry, considered polynomials up to fifth order, and evaluated orbital temperature variations to establish a minimum and maximum temperature for each component. A simplified Outer Mold Line (OML) thermal model of the Orion spacecraft was used in this study. The factors chosen were the vehicle's Yaw, Pitch, and Roll (defining the on-orbit attitude), the Beta angle (restricted to positive beta angles from 0 to 75), and the environmental constants (varying from cold to hot). All factors were normalized from their native ranges to a non-dimensional range from -1.0 to 1.0. Twenty-three components from the OML were chosen and the minimum and maximum orbital temperatures were calculated for each to produce forty-six responses for the DOE model. A customized DOE case matrix of 145 analysis cases was developed which used analysis points at the factor corners, mid-points, and center. From this data set, RSE s were developed which consisted of cubic, quartic, and fifth order polynomials. The results presented are for the fifth order RSE. The RSE results were then evaluated for agreement with the analytical model predictions to produce a +/-3(sigma) error band. Forty of the 46 responses had a +/-3(sigma) value of 10 F or less. Encouraged by this initial success, two additional sets of verification cases were selected. One contained 20 cases, the other 50 cases. These cases were evaluated both with the fifth order RSE and with the analytical model. For the maximum temperature predictions, 12 of the 23 components had all predictions within +/-10 F and 17 were within +/-20 F. For the minimum temperature predictions, only 4 of the 23 components (the four radiator temperatures), were within the 10 F goal. The maximum temperature RSEs were then run through 59,049 screening cases. The RSE predictions were then filtered to find 55 cases that produced the hottest temperatures. These 55 cases were then analyzed using the thermal model and the results compared against the RSE predictions. As noted earlier, 12 of the 23 responses were within +/-10 F at 17 within +/-20 F. These results demonstrate that if properly formulated, an RSE can provide a reliable, fast temperature prediction. Despite this progress, additional work is needed to determine why the minimum temperatures responses and 6 of the hot temperature responses did not produce reliable RSEs. Recommend focus areas are the model itself (arithmetic vs. diffusion nodes) and seeking consultations with statistical application experts.

Miller, StephenW.↗

Ensemble Methodologies for Astronaut Cancer Risk Assessment in the face of Large Uncertainties

A new approach to NASA space radiation risk modeling has successfully extended the current NASA probabilistic cancer risk model to an ensemble framework able to consider sub-model parameter uncertainty (e.g. uncertainty in a radiation quality parameter) as well as model-form uncertainty associated with differing theoretical or empirical formalisms (e.g. combined dose-rate and radiation quality effects). Ensemble methodologies are already widely used in weather prediction, modeling of infectious disease outbreaks, and certain terrestrial radiation protection applications to better understand how uncertainty may influence risk decision-making. Applying ensemble methodologies to space radiation risk projections offers the potential to efficiently incorporate emerging research results, allow for the incorporation of future (including international) models, improve uncertainty quantification for underlying sub-models developed against sparse experimental data, and reduce the impact of subjective bias on risk projections. Moreover, risk forecasting across an ensemble of multiple predictive models can provide stakeholders additional information on risk acceptance if current health/medical standards cannot be met or the level of knowledge doesn’t permit a specific risk or exposure limit to be developed for future space exploration missions. In this work, ensemble risk projections implementing multiple sub-models of radiation quality, dose and dose-rate effectiveness factors, excess risk, and latency as ensemble members are presented. Initial consensus methods for ensemble model weights and correlations to account for individual model bias are discussed. In these analyses, the ensemble forecast compares well to results from NASA's current operational cancer risk projection model used to assess permissible exposure limits and permissible mission durations for astronauts. However, a large range of projected risk values are obtained at the upper 95th confidence level where models must extrapolate beyond available biological data sets; closer agreement is seen at the median + one sigma due to the inherent similarities in available models. Future work, including the addition of new models and methods for statistical correlation between predictive members are discussed to define alternate ways of thinking about risk and ‘acceptable’ uncertainty with respect to NASA’s current permissible exposure limits.

space radiation↗

Ensemble Cancer Risk Model for Astronaut Risk Assessment

A new approach to NASA space radiation risk modeling has successfully extended the current NASA probabilistic cancer risk model to an ensemble framework able to consider sub-model parameter uncertainty (e.g. uncertainty in a radiation quality parameter) as well as model-form uncertainty associated with differing theoretical or empirical formalisms (e.g. combined dose-rate and radiation quality effects). Ensemble methodologies are already widely used in weather prediction, modeling of infectious disease outbreaks, and certain terrestrial radiation protection applications to better understand how uncertainty may influence risk decision-making. Applying ensemble methodologies to space radiation risk projections offers the potential to efficiently incorporate emerging research results, allow for the incorporation of future (including international) models, improve uncertainty quantification for underlying sub-models developed against sparse experimental data, and reduce the impact of subjective bias on risk projections. Moreover, risk forecasting across an ensemble of multiple predictive models can provide stakeholders additional information on risk acceptance if current health/medical standards cannot be met or the level of knowledge doesn’t permit a specific risk or exposure limit to be developed for future space exploration missions. In this work, ensemble risk projections implementing multiple sub-models of radiation quality, dose and dose-rate effectiveness factors, excess risk, and latency as ensemble members are presented. Initial consensus methods for ensemble model weights and correlations to account for individual model bias are discussed. In these analyses, the ensemble forecast compares well to results from NASA's current operational cancer risk projection model used to assess permissible exposure limits and permissible mission durations for astronauts. However, a large range of projected risk values are obtained at the upper 95th confidence level where models must extrapolate beyond available biological data sets; closer agreement is seen at the median + one sigma due to the inherent similarities in available models. Future work, including the addition of new models and methods for statistical correlation between predictive members are discussed to define alternate ways of thinking about risk and ‘acceptable’ uncertainty with respect to NASA’s current permissible exposure limits.

Lisa C Simonsen↗

Experimental and Computational Study fo CH, CH*, and OH* in an Axisymmetric Laminar Diffusion Flame

In this study, we extend the results of previous combined numerical and experimental investigations of an axisymmetric laminar diffusion flame in which difference Raman spectroscopy, laser-induced fluorescence (LIF), and a multidimensional flame model were used to generate profiles of the temperature and major and minor species. A procedure is outlined by which the number densities of ground-state CH (X(sup 2)II) excited-state CH (A(sup 2)Delta, denoted CH*), and excited-state OH (A(sup 2)Sigma, denoted OH*) are measured and modeled. CH* and OH* number densities are deconvoluted from line-of-sight flame-emission measurements. Ground-state CH is measured using linear LIF. The computations are done with GRI Mech 2.11 as well as an alternate hydrocarbon mechanism. In both cases, additional reactions for the production and consumption of CH* and OH* are added from recent kinetic studies. Collisional quenching and spontaneous emission are responsible for the de-excitation of the excited-state radicals. As with our previous investigations, GRI Mech 2.11 continues to produce very good agreement with the overall flame length observed in the experiments, while significantly under predicting the flame lift-off height. The alternate kinetic scheme is much more accurate in predicting lift-off height but overpredicts the over-all flame length. Ground-state CH profiles predicted with GRI Mech 2.11 are in excellent agreement with the corresponding measurements, regarding both spatial distribution and absolute concentration (measured at 4 ppm) of the CH radical. Calculations of the excited-state species show reasonable agreement with the measurements as far as spatial distribution and overall characteristics are concerned. For OH*, the measured peak mole fraction, 1.3 x 10(exp -8), compared well with computed peaks, while the measured peak level for CH*, 2 x 10(exp -9), was severely underpredicted by both kinetic schemes, indicating that the formation and destruction kinetics associated with excited-state species in flames require further research.

Walsh, K. T.↗

PECHCV, PECHFV, PEFHCV and PEFHFV: A set of atmospheric, primitive equation forecast models for the Northern Hemisphere, volume 3

As part of the SEASAT program of NASA, a set of four hemispheric, atmospheric prediction models were developed. The models, which use a polar stereographic grid in the horizontal and a sigma coordinate in the vertical, are: (1) PECHCV - five sigma layers and a 63 x 63 horizontal grid, (2) PECHFV - ten sigma layers and a 63 x 63 horizontal grid, (3) PEFHCV - five sigma layers and a 187 x 187 horizontal grid, and (4) PEFHFV - ten sigma layers and a 187 x 187 horizontal grid. The models and associated computer programs are described.

Wellck, R. E.↗

Implicit model following and parameter identification of unstable aircraft

A transformation in the s-plane is described which has utility in implicit model-following optimal control design application and in estimation or parameter identification problems. The objective of the transformation is, for the control problem, to achieve an unstable closed-loop system, and, for the estimation problem, to alleviate algorithm convergence problems that may arise in identifying unstable systems. For the control problem, the transformation is a shift along the real (sigma) axis of the plant and model poles and zeros. This transformation is shown to be equivalent to a modified performance index but offers the advantage of compatibility with existing optimal control solution algorithms. For the estimation problem, the data are multiplied by an exponential function and the assumed measurement and process noise covariances are appropriately modified. Examples of both control and estimation applications are presented.

Lebacqz, J. V.↗

A mesoscale sixth-order numerical modelling system

A numerical simulation system is currently under development for NASA which is intended to improve the modeling of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development and atmospheric transport processes. The model utilizes a standard hydrostatic sigma-p coordinate primitive equation set, with x,y-space differencing accurate to eighth order. A three-step dynamic initialization procedure is employed between the analysis of real-time data and grid interpolation. Results of an 18-hour simulation during which synoptic scale cyclogenesis, subsynoptic scale jet streak adjustments, mesoscale convergence zones and tornadic storms were observed have shown the present model to have the potential for simulating the fine-scale structure of features associated with cyclogenesis and intense squall-line development. The mesoscale model was also found to produce less truncation than the NWS LFM model, although a frictionless version of the mesoscale model somewhat overdeepens and overaccelerates features.

Kaplan, M. L.↗

A preliminary C-band scatterometer model function for the ERS-1 AMI instrument

Backscatter cross section (sigma(sub 0)) measurements from the ERS-1 scatterometer are collocated with surface wind velocity data from operational Numerical Weather Prediction (NWP) analyses to yield a fully empirical C band model function. The empirical model does not obey a power law at speeds less than 5 m/s, and upwind/crosswind ratios are small for all incidence angles at low wind speeds and for small incidence angles at all wind speeds. Downwind sigma(sub 0) exceeds upwind values for incidence angles below approximately 28 degrees. The full ERS-1 scatterometer data set has been reprocessed using the NWP model function.

Freilich, M. H.↗

A nested-grid limited-area model for short term weather forecasting

The present investigation is concerned with a mesoscale atmospheric simulation system (MASS), incorporating the sigma-coordinate primitive equations. The present version of this model (MASS 3.0) has 14 vertical layers, with the upper boundary at 100 mb. There are 128 x 96 grid points in each layer. The earlier version of this model (MASS 2.0) has been described by Kaplan et al. (1982). The current investigation provides a summary of major revisions to that version and a description of the parameterization schemes which are presently included in the model. The planetary boundary layer (PBL) is considered, taking into account aspects of generalized similarity theory and free convection, the surface energy budget, the surface moisture budget, and prognostic equations for the depth h of the PBL. A cloud model is discussed, giving attention to stable precipitation, and cumulus convection.

Wong, V. C.↗

Retrieval of pine forest biomass using JPL AIRSAR data

The analysis of Jet Propulsion Laboratory (JPL) Airborne Synthetic Aperture Radar (AIRSAR) data over the Landes forest in South-West France revealed strong correlation between L- and especially P-band sigma degrees and the pine forest biomass. To explain the physical link of radar backscatter to biomass, a polarimetric backscattering model was developed and validated. Then the model was used in a simulation study to predict sigma degree sensitivity to undesired canopy and environmental parameters. Main results concerning the data analysis, modeling, and simulation at P-band are reported.

Beaudoin, A.↗

Sensitivity of an atmospheric photochemistry model to chlorine perturbations including consideration of uncertainty propagation

Models of stratospheric photochemistry are generally tested by comparing their predictions for the composition of the present atmosphere with measurements of species concentrations. These models are then used to make predictions of the atmospheric sensitivity to perturbations. Here the problem of the sensitivity of such a model to chlorine perturbations ranging from the present influx of chlorine-containing compounds to several times that influx is addressed. The effects of uncertainties in input parameters, including reaction rate coefficients, cross sections, solar fluxes, and boundary conditions, are evaluated using a Monte Carlo method in which the values of the input parameters are randomly selected. The results are probability distributions for present atmosheric concentrations and for calculated perturbations due to chlorine from fluorocarbons. For more than 300 Monte Carlo runs the calculated ozone perturbation for continued emission of fluorocarbons at today's rates had a mean value of -6.2 percent, with a 1-sigma width of 5.5 percent. Using the same runs but only allowing the cases in which the calculated present atmosphere values of NO, NO2, and ClO at 25 km altitude fell within the range of measurements yielded a mean ozone depletion of -3 percent, with a 1-sigma deviation of 2.2 percent. The model showed a nonlinear behavior as a function of added fluorocarbons. The mean of the Monte Carlo runs was less nonlinear than the model run using mean value of the input parameters.

Stolarski, R. S.↗

Estimation of biophysical properties of forest canopies using C-band microwave data

A scatterometer ranging experiment is described in which C-band data and boresight photography were collected using a helicopter so as to provide a capability to study scattering processes in forest canopies in the Superior National Forest in Minnesota. An inversion scheme is used to determine C-band volume extinction and scattering coefficients for high density aspen sites. Analysis of data through the season indicates that VV, HH, and VH volume extinction coefficients change during the year and are presumably affected by the emergence and senescence of leaves. A linear relationship was observed between sigma(0) (VV) and leaf area index for low and medium density aspen sites, but a large decrease occurred in sigma(0) for both high density sites. Calculations using the Fung disk model, which accounts only for the leaves, shows an underestimate of sigma(0) by a factor of two or three indicating that scattering by branches and soil background may be important at C-band.

Pitts, David E.↗