Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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 469 records · Page 26

Minimal modeling of the extratropical general circulation

The ability of low-order, two-layer models to reproduce basic features of the mid-latitude general circulation is investigated. Changes in model behavior with increased spectral resolution are examined in detail. Qualitatively correct time-mean heat and momentum balances are achieved in a beta-plane channel model which includes the first and third meridional modes. This minimal resolution also reproduces qualitatively realistic surface and upper-level winds and mean meridional circulations. Higher meridional resolution does not result in substantial changes in the latitudinal structure of the circulation. A qualitatively correct kinetic energy spectrum is produced when the resolution is high enough to include several linearly stable modes. A model with three zonal waves and the first three meridional modes has a reasonable energy spectrum and energy conversion cycle, while also satisfying heat and momentum budget requirements. This truncation reproduces the basic mechanisms and zonal circulation features that are obtained at higher resolution. The model performance improves gradually with higher resolution and is smoothly dependent on changes in external parameters.

O'Brien, Enda↗

Atmospheric drag model calibrations for spacecraft lifetime prediction

Although solar activity prediction uncertainty normally dominates decay prediction error budget for near-Earth spacecraft, the effect of drag force modeling errors for given levels of solar activity needs to be considered. Two atmospheric density models, the modified Harris-Priester model and the Jacchia-Roberts model, to reproduce the decay histories of the Solar Mesosphere Explorer (SME) and Solar Maximum Mission (SMM) spacecraft in the 490- to 540-kilometer altitude range were analyzed. Historical solar activity data were used in the input to the density computations. For each spacecraft and atmospheric model, a drag scaling adjustment factor was determined for a high-solar-activity year, such that the observed annual decay in the mean semimajor axis was reproduced by an averaged variation-of-parameters (VOP) orbit propagation. The SME (SMM) calibration was performed using calendar year 1983 (1982). The resulting calibration factors differ by 20 to 40 percent from the predictions of the prelaunch ballistic coefficients. The orbit propagations for each spacecraft were extended to the middle of 1988 using the calibrated drag models. For the Jaccia-Roberts density model, the observed decay in the mean semimajor axis of SME (SMM) over the 4.5-year (5.5-year) predictive period was reproduced to within 1.5 (4.4) percent. The corresponding figure for the Harris-Priester model was 8.6 (20.6) percent. Detailed results and conclusions regarding the importance of accurate drag force modeling for lifetime predictions are presented.

Binebrink, A. L.↗

Soft X-ray emission from electron-beam-heated solar flares

Using time-dependent numerical simulations and Solar Maximum Mission observations of a solar flare on 1985 January 23, a study is conducted of the ability of an electron-beam-heating model to reproduce the rise phase of a flare as observed in soft X-ray lines of Ca XIX. The electron beam is parameterized by a peak flux, a low-energy cutoff, and a spectral index, and has a time dependence similar to the observed hard X-ray burst. For a spectral index of 6, only models with a low-energy cutoff of 20 keV reproduce the observed peak emission in the Ca XIX line complex. All models with a low-energy cutoff of 15 keV produce too much emission, while all models with a 25-keV cutoff too little emission. None of the models reproduces the temporal behavior of the soft X-ray emission. The electron-beam-heated component is theorized to only represent a small fraction of the energy released in the impulsive phase of this flare.

Mariska, John T.↗

A global potential energy surface for ArH2

We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.

Schwenke, David W.↗

Non-LTE radiating acoustic shocks and Ca II K2V bright points

We present, for the first time, a self-consistent solution of the time-dependent 1D equations of non-LTE radiation hydrodynamics in solar chromospheric conditions. The vertical propagation of sinusoidal acoustic waves with periods of 30, 180, and 300 s is calculated. We find that departures from LTE and ionization recombination determine the temperature profiles of the shocks that develop. In LTE almost all the thermal energy goes into ionization, so the temperature rise is very small. In non-LTE, the finite transition rates delay the ionization to behind the shock front. The compression thus goes into thermal energy at the shock front leading to a high temperature amplitude. Further behind the shock front, the delayed ionization removes energy from the thermal pool, which reduces the temperature, producing a temperature spike. The 180 s waves reproduce the observed temporal changes in the calcium K line profiles quite well. The observed wing brightening pattern, the violet/red peak asymmetry and the observed line center behavior are all well reproduced. The short-period waves and the 5 minute period waves fail especially in reproducing the observed behavior of the wings.

Carlsson, Mats↗

A global potential energy surface for ArH2

We describe a simple analytic representation of the ArH2 potential energy surface which well reproduces the results of extensive ab initio electronic structure calculations. The analytic representation smoothly interpolates between the dissociated H2 and strong bonding limits. In the fitting process, emphasis is made on accurately reproducing regions of the potential expected to be important for high temperature (ca. 3000 K) collision processes. Overall, the anisotropy and H2 bond length dependence of the analytic representation well reproduce the input data.

Schwenke, David W.↗

Material Data Representation of Hysteresis Loops for Hastelloy X Using Artificial Neural Networks

The artificial neural network (ANN) model proposed by Rumelhart, Hinton, and Williams is applied to develop a functional approximation of material data in the form of hysteresis loops from a nickel-base superalloy, Hastelloy X. Several different ANN configurations are used to model hysteresis loops at different cycles for this alloy. The ANN models were successful in reproducing the hysteresis loops used for its training. However, because of sharp bends at the two ends of hysteresis loops, a drift occurs at the corners of the loops where loading changes to unloading and vice versa (the sharp bends occurred when the stress-strain curves were reproduced by adding stress increments to the preceding values of the stresses). Therefore, it is possible only to reproduce half of the loading path. The generalization capability of the network was tested by using additional data for two other hysteresis loops at different cycles. The results were in good agreement. Also, the use of ANN led to a data compression ratio of approximately 22:1.

Alam, Javed↗

A coronal magnetic field model with horizontal volume and sheet currents

When globally mapping the observed photospheric magnetic field into the corona, the interaction of the solar wind and magnetic field has been treated either by imposing source surface boundary conditions that tacitly require volume currents outside the source surface or by limiting the interaction to thin current sheets between oppositely directed field regions. Yet observations and numerical Magnetohydrodynamic (MHD) calculations suggest the presence of non-force-free volume currents throughout the corona as well as thin current sheets in the neighborhoods of the interfaces between closed and open field lines or between oppositely directed open field lines surrounding coronal helmet-streamer structures. This work presents a model including both horizontal volume currents and streamer sheet currents. The present model builds on the magnetostatic equilibria developed by Bogdan and Low and the current-sheet modeling technique developed by Schatten. The calculation uses synoptic charts of the line-of-sight component of the photospheric magnetic field measured at the Wilcox Solar Observatory. Comparison of an MHD model with the calculated model results for the case of a dipole field and comparison of eclipse observations with calculations for CR 1647 (near solar minimum) show that this horizontal current-current-sheet model reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces polar plumes and axes of corona streamers better than the source-surface model and reproduces coro nal helmet structures better than the current-sheet model.

Zhao, Xuepu↗

A two-dimensional model study of the QBO signal in SAGE II NO2 and O3

Calculations of the quasi biennial oscillation (QBO) signal in Stratospheric Aerosol and Gas Experiment (SAGE) II O3 and NO2 data between 1984 and 1991 are presented and have been investigated by using a two-dimensional model. The isentropic 2D model is a fully interactive radiative-dynamical-chemical model in which the eddy fluxes of chemical species are calculated in a consistent manner. The QBO in the model has been forced by relaxing the equatorial zonal wind toward the observations at Singapore allowing the comparison of the model with observations from specific years. The model reproduces the observed vertical structure of the equatorial ozone anomaly with the well-known transition from dynamical to photochemical control at around 28km. The model also reproduces the observed vertical structure of the SAGE II observed NO2 anomaly. The model studies have shown that it is the QBO modulation of NO2 which the main cause of QBO signal in O3 above 30 km. The model also reproduces the observed latitudinal structure of the QBO signals in O3 and NO2. Due to the differing horizontal distribution of O3 and NO(y) the ozone signal shows a distinct phase change in the subtropics whereas the NO2 anomaly gives a broader signal.

Chipperfield, M. P.↗

Mirages and the nature of Pluto's atmosphere

We present model occultation lightcurves demonstrating that a strong thermal inversion layer at the base of Pluto's stratosphere can reproduce the minimum flux measured by the Kuiper Airborne Observatory (KAO) during the 1988 occultation of a star by Pluto. The inversion layer also forms the occultation equivalent of a mirage at a radius of 1198 km, which is capable of hiding tropospheres of significant depth. Pluto's surface lies below 1198 km, its radius depending on the depth of the troposphere. We begin by computing plausible temperature structures for Pluto's lower atmosphere, constrained by a calculation of the temperature of the atmosphere near the surface. We then trace rays from the occulted star through the model atmosphere, computing the resultant bending of the ray. Model light curves are obtained by summing the contribution of individual rays within the shadow of Pluto on Earth. We find that we can reproduce the KAO lightcurve using model atmospheres with a temperature inversion and no haze. We have explored models with tropospheres as deep as 40 km (implying a Pluto radius of 1158 km) that reproduce the suite of occultation data. Deeper tropospheres can be fitted to the data, but the mutual event radius of 1150 km probably provides a lower bound. If Pluto has a shallow or nonexistent troposphere, its density is consistent with formation in the solar nebula with modest water loss due to impact ejection. If the troposhere is relatively deep, implying a smaller radius and larger density, significant amounts of water loss are required.

Stansberry, J. A.↗

Development of an Antimicrobial Susceptibility Testing Method Suitable for Performing During Space Flight

Very little is known regarding the affects of the microgravity environment of space flight upon the action of antimicrobial agents on bacterial pathogens. This study was undertaken to develop a simple method for conducting antibacterial susceptibility tests during a Space Shuttle mission. Specially prepared susceptibility test research cards (bioMerieux Vitek, Hazelwood, MO) were designed to include 6-11 serial two-fold dilutions of 14 antimicrobial agents, including penicillins, cephalosporins, a Beta-lactamase inhibitor, vancomycin, erythromycin, tetracycline, gentamicin, ciprofloxacin, and trimethoprim/sulfamethoxazole. Minimal inhibitory concentrations (MICS) of the drugs were determined by visual reading of color endpoints in the Vitek research cards made possible by incorporation of a colorimetric growth indicator (alamarBlue(Trademark), Accumed International, Westlake, OH). This study has demonstrated reproducible susceptibility results when testing isolates of Staphylococcus aurezis, Group A Streptococcus, Enterococcusfaecalis, Escherichia coli (beta-lactamase positive and negative strains), Klebsiella pneumoniae, Enterobacter cloacae, and Pseudomoiias aeruginosa. In some instances, the MICs were comparable to those determined using a standard broth microdilution method, while in some cases the unique test media and format yielded slightly different values, that were themselves reproducible. The proposed in-flight experiment will include inoculation of the Vitek cards on the ground prior to launch of the Space Shuttle, storage of inoculated cards at refrigeration temperature aboard the Space Shuttle until experiment initiation, then incubation of the cards for 18-48 h prior to visual interpretation of MICs by the mission's astronauts. Ground-based studies have shown reproducible MICs following storage of inoculated cards for 7 days at 4-8 C to accommodate the mission's time schedule and the astronauts' activities. For comparison, ground-based control (normal gravity) MIC values will be generated by simultaneous inoculation and incubation of a second set of test cards in a laboratory at the launch site. This procedure can provide a safe and compact experiment that should yield new information on the affects of microgravity on the biological activities of various classes of antibiotics.

Jorgensen, James H.↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

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↗

Comparison of TOPEX/Poseidon Sea Level and Linear Model Results forced by Various Wind Products for the Tropical Pacific

The goal of this paper is to compare TOPEX/Posaidon (T/P) sea level with sea level results from linear ocean model experiments forced by several different wind products for the tropical Pacific. During the period of this study (October 1992 - October 1995), available wind products include satellite winds from the ERS-1 scatterometer product of [HALP 97] and the passive microwave analysis of SSMI winds produced using the variational analysis method (VAM) of [ATLA 91]. In addition, atmospheric GCM winds from the NCEP reanalysis [KALN 96], ECMWF analysis [ECMW94], and the Goddard EOS-1 (GEOS-1) reanalysis experiment [SCHU 93] are available for comparison. The observed ship wind analysis of FSU [STRI 92] is also included in this study. The linear model of [CANE 84] is used as a transfer function to test the quality of each of these wind products for the tropical Pacific. The various wind products are judged by comparing the wind-forced model sea level results against the T/P sea level anomalies. Correlation and RMS difference maps show how well each wind product does in reproducing the T/P sea level signal. These results are summarized in a table showing area average correlations and RMS differences. The large-scale low-frequency temporal signal is reproduced by all of the wind products, However, significant differences exist in both amplitude and phase on regional scales. In general, the model results forced by satellite winds do a better job reproducing the T/P signal (i.e. have a higher average correlation and lower RMS difference) than the results forced by atmospheric model winds.

Hackert, Eric C.↗

Modeling the Diurnal Tides in the MLT Region with the Doppler Spread Parameterization of Gravity Waves

Special Session: SA03 The mesosphere/lower thermosphere region: Structure, dynamics, composition, and emission. Ground based and satellite observations in the upper mesosphere and lower thermosphere (MLT) reveal large seasonal variations in the horizontal wind fields of the diurnal and semidiurnal tides. To provide an understanding of the observations, we discuss results obtained with our Numerical Spectral Model (NMS) that incorporates the gravity wave Doppler Spread Parameterization (DSP) of Hines. Our model reproduces many of the salient features observed, and we discuss numerical experiments that delineate the important processes involved. Gravity wave momentum deposition and the seasonal variations in the tidal excitation contribute primarily to produce the large equinoctial amplitude maxima in the diurnal tide. Gravity wave induced variations in eddy viscosity, not accounted for in the model, have been shown by Akmaev to be important too. For the semidiurnal tide, with amplitude maximum observed during winter solstice, these processes also contribute, but filtering by the mean zonal circulation is more important. A deficiency of our model is that it cannot reproduce the observed seasonal variations in the phase of the semidiurnal tide, and numerical experiments are being carried out to diagnose the cause and to alleviate this problem. The dynamical components of the upper mesosphere are tightly coupled through non-linear processes and wave filtering, and this may constrain the model and require it to reproduce in detail the observed phenomenology.

Mayr, H. G.↗

Who or What? Self-Replication and Function-Reproduction in the Origin of Life

In this presentation, we will present results on the fundamental properties of two classes of replicating systems: autocatalytic replicators that reproduce exact copies of a template molecule, and function reproducers that maintain a set of essential functions without replicating the identities of the functional moieties. We will describe the stability and behavior in-the-large of autocatalytic replicators. Most importantly, we have found no sharp distinction between an autocatalytic and a non-autocatalytic domain. We will also present a new derivation of von Kiedrowski's square-root rate law. Function - reproducers are proposed as an important component of protocells and we will present theoretical results on a simple model system that incorporates known peptide biophysics. For a wide range of parameters, we have shown that this type of system can improve its overall performance, even in the absence of any method for information storage. This type of system improvement is defined to be non-genomic evolution.

New, Michael H.↗