Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical physics”

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 289 records · Page 16

A proposal for the intercomparison of the dynamical cores of atmospheric general circulation models

A benchmark calculation is proposed for evaluating the dynamical cores of atmospheric general circulation models (GCMs) independently of the physical parameterizations. The test focuses on the long-term statistical properties of a fully developed general circulation; thus, it is particularly appropriate for intercomparing the dynamics used in climate models. To illustrate the use of this benchmark, two very different atmospheric dynamical cores--one spectral, one finite difference--are compared. It is found that the long-term statistics produced by the two models are very similar. Selected results from these calculations are presented to initiate the intercomparison.

Held, Isaac M.↗

Outcome analysis of individualized vestibular rehabilitation protocols

OBJECTIVE: To determine the outcome of vestibular rehabilitation protocols in subjects with peripheral vestibular disorders compared with normal and abnormal control subjects. STUDY DESIGN: Prospective study using repeated measure, matched control design. Subjects were solicited consecutively according to these criteria: vestibular disorder subjects who had abnormal results of computerized dynamic posturography (CDP) sensory organization tests (SOTs) 5 and 6 and underwent rehabilitation; vestibular disorder subjects who had abnormal results of SOTs 5 and 6 and did not undergo rehabilitation; and normal subjects (normal SOTs). SETTING: Tertiary neurotology clinic. SUBJECTS: Men and women over age 18 with chronic vestibular disorders and chief complaints of unsteadiness, imbalance, and/or motion intolerance, and normal subjects. INTERVENTIONS: Pre- and post-rehabilitation assessment included CDP, vestibular disability, and activities of daily living questionnaires. Individualized rehabilitation plans were designed and implemented to address the subject's specific complaints and functional deficits. Supervised sessions were held at weekly intervals, and self-administered programs were devised for daily home use. MAIN OUTCOME MEASURES: CDP composite and SOT scores, number of falls on CDP, and self-assessment questionnaire results. RESULTS: Subjects who underwent rehabilitation (Group A) showed statistically significant improvements in SOTs, overall composite score, and reduction in falls compared with abnormal (Group B) control groups. Group A's performances after rehabilitation were not significantly different from those of normal subjects (Group C) in SOTs 3 through 6, and close to normal on SOTs 1 and 2. Subjects in Group A also reported statistically significant symptomatic improvement. CONCLUSIONS: Outcome measures of vestibular protocol physical therapy confirmed objective and subjective improvement in subjects with chronic peripheral vestibular disorders. These findings support results reported by other investigators.

NASA Discipline Neuroscience↗

ARM Research in the Equatorial Western Pacific: A Decade and Counting

The tropical western Pacific (TWP) is an important climatic region. Strong solar heating, warm sea surface temperatures, and the annual progression of the intertropical convergence zone (ITCZ) across this region generate abundant convective systems, which through their effects on the heat and water budgets have a profound impact on global climate and precipitation. In order to accurately evaluate tropical cloud systems in models, measurements of tropical clouds, the environment in which they reside, and their impact on the radiation and water budgets are needed. Because of the remote location, ground-based datasets of cloud, atmosphere, and radiation properties from the TWP region have come primarily from short-term field experiments. While providing extremely useful information on physical processes, these short-term datasets are limited in statistical and climatological information. To provide longterm measurements of the surface radiation budget in the tropics and the atmospheric properties that affect it, the Atmospheric Radiation Measurement program established a measurement site on Manus Island, Papua New Guinea, in 1996 and on the island republic of Nauru in late 1998. These sites provide unique datasets now available for more than 10 years on Manus and Nauru. This article presents examples of the scientific use of these datasets including characterization of cloud properties, analysis of cloud radiative forcing, model studies of tropical clouds and processes, and validation of satellite algorithms. New instrumentation recently installed at the Manus site will provide expanded opportunities for tropical atmospheric science.

tropical regions↗

K-Means Cluster Study for Radiofrequency Propagation Characterization

The objective of this study is to design a simple method for mining radio frequency (RF) propagation data. The study explored the characteristics of a large dataset of propagation experiments conducted over the span of years and using several ground stations around the world. Furthermore, this study developed simple predictive models that can be used for link characterization and overall propagation behavior description, without the need for physical measurements on-site. It is understood that such statistical learning has several drawbacks in terms of accuracy and precision. K-means clustering was used to characterize the data set in a way never explored before in an attempt to create useful tools that reduce cost, time and risk. K-means clustering was used to characterize the data set. Cosine distance was used as a method to determine the optimal number for clustering each feature. Dependence and independence analysis was performed to explore intra and inter-sensitivity between the presented features, with respect to each other and time. Several predicative models were generated and evaluated with respect to a test set to assess a measure of prediction accuracy and precision. A simple method for data analysis was developed and tested as the basis for further studies and future refinement to produce optimal performing models.

Cognitive↗

Modeling Spacecraft Safe Mode Events

Spacecraft enter a ‘safe mode’ to protect the vehicle when a potentially harmful anomaly occurs. This minimally functioning state isolates faults, establishes contact with Earth, and orients the vehicle into a power positive attitude until operators intervene. Though ‘safings’ are inherently unpredictable, mission teams build in time margin during operations to determine root causes and restore functionality. Planning and managing this margin is both critical and enabling on mission architectures dependent on near-continuous operability – such as a low-thrust electric propulsion mission. To better quantify the occurrences and severity of safe mode anomalies, the Jet Propulsion Laboratory (JPL) has assembled a database of safings from past and active missions. Currently nearly 240 records are captured from 21 beyond-Earth missions, stemming from a collaboration between teams at JPL, Ames Research Center, Goddard Space Flight Center, and the Johns Hopkins University Applied Physics Laboratory. This paper discusses the event database, explores a statistical approach in modeling the occurrences and severity of safing events, presents a simulation technique, and details recommendations and future work to benefit future concepts.

Nicholas, Austin↗

Detector noise statistics in the non-linear regime

The statistical behavior of an idealized linear detector in the presence of threshold and saturation levels is examined. It is assumed that the noise is governed by the statistical fluctuations in the number of photons emitted by the source during an exposure. Since physical detectors cannot have infinite dynamic range, our model illustrates that all devices have non-linear regimes, particularly at high count rates. The primary effect is a decrease in the statistical variance about the mean signal due to a portion of the expected noise distribution being removed via clipping. Higher order statistical moments are also examined, in particular, skewness and kurtosis. In principle, the expected distortion in the detector noise characteristics can be calibrated using flatfield observations with count rates matched to the observations. For this purpose, some basic statistical methods that utilize Fourier analysis techniques are described.

Shopbell, P. L.↗

Evaluation of the 29-km Eta Model: Objective Verification at Three Selected Stations - Part 1

This paper describes an objective verification of the National Centers for Environmental Prediction (NCEP) 29-km eta model from May 1996 through January 1998. The evaluation was designed to assess the model's surface and upper-air point forecast accuracy at three selected locations during separate warm (May - August) and cool (October - January) season periods. In order to enhance sample sizes available for statistical calculations, the objective verification includes two consecutive warm and cool season periods. Systematic model deficiencies comprise the larger portion of the total error in most of the surface forecast variables that were evaluated. The error characteristics for both surface and upper-air forecasts vary widely by parameter, season, and station location. At upper levels, a few characteristic biases are identified. Overall however, the upper-level errors are more nonsystematic in nature and could be explained partly by observational measurement uncertainty. With a few exceptions, the upper-air results also indicate that 24-h model error growth is not statistically significant. In February and August 1997, NCEP implemented upgrades to the eta model's physical parameterizations that were designed to change some of the model's error characteristics near the surface. The results shown in this paper indicate that these upgrades led to identifiable and statistically significant changes in forecast accuracy for selected surface parameters. While some of the changes were expected, others were not consistent with the intent of the model updates and further emphasize the need for ongoing sensitivity studies and localized statistical verification efforts. Objective verification of point forecasts is a stringent measure of model performance, but when used alone, is not enough to quantify the overall value that model guidance may add to the forecast process. Therefore, results from a subjective verification of the meso-eta model over the Florida peninsula are discussed in the companion paper by Manobianco and Nutter. Overall verification results presented here and in part two should establish a reasonable benchmark from which model users and developers may pursue the ongoing eta model verification strategies in the future.

Nutter, Paul A.↗

Generic results of the space physics community survey

This report summarizes the results of a survey of the members of the space physics research community conducted in 1990-1991 to ascertain demographic information on the respondents and information on their views on a number of facets of their space physics research. The survey was conducted by questionnaire and the information received was compiled in a database and analyzed statistically. The statistical results are presented for the respondent population as a whole and by four different respondent cross sections: individual disciplines of space physics, type of employers, age groups, and research techniques employed. Data from a brief corresponding survey of the graduate students of respondents are also included.

Sharma, Rikhi R.↗

The Strings of Eta Carina: The HST/STIS Spectra and [Ca II]

Long linear, filamentary ejecta, are found to move at very high velocity external to the Homunculus, the circumstellar hourglass-shaped ejecta surrounding Eta Carinae. The origin of the strings is a puzzle. As an example, the Weigelt Blobs have N at 10X solar and C, O at 0.01X solar abundance, along with He/H significantly enhanced. This abundance pattern is evidence for extreme CNO-processing. Similarly, the Strontium Filament has Ti/Ni at 100X solar, presumably due to the lack of oxygen to form Ti-oxide precipitates onto dust grains. We have obtained 2-D spectra with the HST/STIS of the Strontium Filament and a portion of a string. These deep spectral exposures, at moderate dispersion, span much of the near red spectral region from 5000 to 9000A. We have identified twelve emission lines in these spectra with proper velocities and spatial structure of this string and obtained line ratios for [Ca II] (7293/7325A) and [Fe Ill (7157/8619A) which are useful for determining physical conditions in this nebulosity. In an attempt to use the [Ca II] ratio to determine the physical parameters, and ultimately the abundances in the strings, we have constructed a statistical equilibrium model for Ca II , including radiative and collisional rates. These results incorporate our newly calculated atomic data for levels n = 3,4,5 and 6 configurations of Ca II. The aim is to compute the [Ca II] line ratios and use them as a diagnostic of the physical parameters. Using the [Fe II] ratio we find that for Te=10,000 K, the electron density is Ne approx.10(exp 6)/cu cm. We plan to use the [Ca II] ratio to confirm this result. Then, we will extend the use of this multilevel model Ca II atom to study the physical conditions of the Strontium filament where eight lines of Ca II, both allowed and forbidden, had been identified. With the physical conditions determined, we will be able to derive reliable estimates for the gas phase abundances in the strings.

Melendez, M. B.↗

Mapping of the Ronda peridotite massif (Spain) from AVIRIS spectro-imaging survey: A first attempt

In both AVIRIS and ISM data, through the use of mixing models, geological boundaries of the Ronda massif are identified with respect to the surrounding rocks. We can also yield first-order vegetation maps. ISM and AVIRIS instruments give consistent results. On the basis of endmember fraction images, it is then possible to discard areas highly vegetated or not belonging to the peridotite massif. Within the remaining part of the mosaic, spectro-mixing analysis reveals spectral variations in the peridotite massif between the well-exposed areas. Spatially organized units are depicted, related to differences in the relative depth of the absorption band at 1 micron, and it may be due to a different pyroxene content. At this stage, it is worth noting that, although mineralogical variations observed in the rocks are at a sub-pixel scale for the airborne analysis, we see an emerging spatial pattern in the distribution of spectral variations across the massif which might be prevailingly related to mineralogy. Although it is known from fieldwork that the Ronda peridotite massif exhibits mineralogical variations at local scale in the content of pyroxene, and at regional scale in different mineral facies, ranging from garnet-, to spinel- to plagioclase-lherzolites, no attempt has been done yet to produce a synoptic map relating the two scales of analysis. The present work is a first attempt to reach this objective, though a lot more work is still required. In particular, for the purpose of mineralogical interpretation, it is critical to relate the airborne observation to field work and laboratory spectra of Ronda rocks already obtained, with the use of image endmembers and associated reference endmembers. Also, the pretty rough linear mixing model used here is taken as a 'black-box' process which does not necessarily apply correctly to the physical situation at the sub-pixel level. One may think of using the ground-truth observations bearing on the sub-pixel statistical characteristics (texture, structural pattern, surface distribution and vegetation contribution (grass,..)) to produce a more advanced mixing model, physically appropriate to the geologic and environmental contexts.

Pinet, P. C.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in General Circulation Models (GCM'S). The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s). and passive scalars K(sub c). First, we took K(sub s) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derived using the latest model, the KPP model by the NCAR group. Presently, we are running the same Geophysical Fluid Dynamics Laboratory (GFDL) code relaxing the assumptive K(sub s) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R rho (= Beta delta S/delta z/alpha delta T/delta z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗

The Statistical Mechanics of Ideal MHD Turbulence

Turbulence is a universal, nonlinear phenomenon found in all energetic fluid and plasma motion. In particular. understanding magneto hydrodynamic (MHD) turbulence and incorporating its effects in the computation and prediction of the flow of ionized gases in space, for example, are great challenges that must be met if such computations and predictions are to be meaningful. Although a general solution to the "problem of turbulence" does not exist in closed form, numerical integrations allow us to explore the phase space of solutions for both ideal and dissipative flows. For homogeneous, incompressible turbulence, Fourier methods are appropriate, and phase space is defined by the Fourier coefficients of the physical fields. In the case of ideal MHD flows, a fairly robust statistical mechanics has been developed, in which the symmetry and ergodic properties of phase space is understood. A discussion of these properties will illuminate our principal discovery: Coherent structure and randomness co-exist in ideal MHD turbulence. For dissipative flows, as opposed to ideal flows, progress beyond the dimensional analysis of Kolmogorov has been difficult. Here, some possible future directions that draw on the ideal results will also be discussed. Our conclusion will be that while ideal turbulence is now well understood, real turbulence still presents great challenges.

Shebalin, John V.↗

Confronting Models with Data: The GEWEX Cloud Systems Study

The GEWEX Cloud System Study (GCSS; GEWEX is the Global Energy and Water Cycle Experiment) was organized to promote development of improved parameterizations of cloud systems for use in climate and numerical weather prediction models, with an emphasis on the climate applications. The strategy of GCSS is to use two distinct kinds of models to analyze and understand observations of the behavior of several different types of clouds systems. Cloud-system-resolving models (CSRMs) have high enough spatial and temporal resolutions to represent individual cloud elements, but cover a wide enough range of space and time scales to permit statistical analysis of simulated cloud systems. Results from CSRMs are compared with detailed observations, representing specific cases based on field experiments, and also with statistical composites obtained from satellite and meteorological analyses. Single-column models (SCMs) are the surgically extracted column physics of atmospheric general circulation models. SCMs are used to test cloud parameterizations in an un-coupled mode, by comparison with field data and statistical composites. In the original GCSS strategy, data is collected in various field programs and provided to the CSRM Community, which uses the data to "certify" the CSRMs as reliable tools for the simulation of particular cloud regimes, and then uses the CSRMs to develop parameterizations, which are provided to the GCM Community. We report here the results of a re-thinking of the scientific strategy of GCSS, which takes into account the practical issues that arise in confronting models with data. The main elements of the proposed new strategy are a more active role for the large-scale modeling community, and an explicit recognition of the importance of data integration.

Randall, David↗

Objective determination of image end-members in spectral mixture analysis

Spectral mixture analysis was shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. The essence of the first phase of the approach is to determine a set of image end-members that best account for the spectral variance in an image cube within a constrained, linear least squares mixing model. The selection of the image end-member is usually achieved using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. The problem of image end-member determination was approached objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.

Tompkins, Stefanie↗

Objective determination of image end-members in spectral mixture analysis of AVIRIS data

Spectral mixture analysis has been shown to be a powerful, multifaceted tool for analysis of multi- and hyper-spectral data. Applications of AVIRIS data have ranged from mapping soils and bedrock to ecosystem studies. During the first phase of the approach, a set of end-members are selected from an image cube (image end-members) that best account for its spectral variance within a constrained, linear least squares mixing model. These image end-members are usually selected using a priori knowledge and successive trial and error solutions to refine the total number and physical location of the end-members. However, in many situations a more objective method of determining these essential components is desired. We approach the problem of image end-member determination objectively by using the inherent variance of the data. Unlike purely statistical methods such as factor analysis, this approach derives solutions that conform to a physically realistic model.

Tompkins, Stefanie↗

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗