Search NASA⌕ Search

SEARCH · Search NASA

Results for “dependable computing”

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 451 records · Page 25

On the solution of problems involving impact type loading

The response of beam and plate members to pulse and impact loading is investigated, and numerical methods are analyzed and compared. The equations of motion, the finite difference and finite element methods, the time integration/Runge-Kutta techniques, and material modeling are discussed in detail. It is found that both the finite difference and the finite element methods could be accurately employed to discretize the spatial variation in the displacements. Central differences or a fourth order Runge-Kutta algorithm could be used for the time integration. The total energy of the system would give the stability and accuracy of the solution. Results showed that the finite element method provided better efficiency in obtaining accurate solution than the finite difference method when a scalar processor is used. However, the finite difference method was more efficient on the vector processor. Therefore, the most efficient method of solution depends on the type of computer present for the analysis.

Moyer, E. T., Jr.↗

Joint nonlinearity effects in the design of a flexible truss structure control system

Nonlinear effects are introduced in the dynamics of large space truss structures by the connecting joints which are designed with rather important tolerances to facilitate the assembly of the structures in space. The purpose was to develop means to investigate the nonlinear dynamics of the structures, particularly the limit cycles that might occur when active control is applied to the structures. An analytical method was sought and derived to predict the occurrence of limit cycles and to determine their stability. This method is mainly based on the quasi-linearization of every joint using describing functions. This approach was proven successful when simple dynamical systems were tested. Its applicability to larger systems depends on the amount of computations it requires, and estimates of the computational task tend to indicate that the number of individual sources of nonlinearity should be limited. Alternate analytical approaches, which do not account for every single nonlinearity, or the simulation of a simplified model of the dynamical system should, therefore, be investigated to determine a more effective way to predict limit cycles in large dynamical systems with an important number of distributed nonlinearities.

Mercadal, Mathieu↗

Theoretical study of the electron temperature in the high-latitude ionosphere for solar maximum and winter conditions

The T(e) variation in the high-latitude ionosphere at altitudes between 120 and 800 km has been modeled for solar maximum, winter solstice, and strong magnetic activity conditions. The calculated electron temperatures are consistent with the plasma densities and ion temperatures computed from a time-dependent ionospheric model. Heating rates for both solar EUV and auroral precipitation were included. In general, the predicted UT variation of the electron temperature that results from the displacement between the magnetic and geographic poles is only a few hundred degrees. However, in sunlit trough regions, T(e) hot spots develop, and these hot spots show a marked UT variation, by as much as 2500 K. The dominant parameter controlling the T(e) variation above 200 km is the magnetospheric heat flux into the ionosphere, which is essentially unknown. For realistic values of the magnetospheric heat flux, the maximum electron temperature ranges from 5000 to 10,000 K at 800 km. A magnetospheric heat flux is particularly effective in enhancing trough electron temperatures. In general, the electron heat flux at high altitudes is uniquely related to the electron temperature and gradient, except on auroral field lines where thermoelectric heat flow is important.

Schunk, R. W.↗

A computer model of global thermospheric winds and temperatures

Output data from the NCAR Thermospheric GCM and a vector-spherical-harmonic (VSH) representation of the wind field are used in constructing a computer model of time-dependent global horizontal vector neutral wind and temperature fields at altitude 130-300 km. The formulation of the VSH model is explained in detail, and some typical results obtained with a preliminary version (applicable to December solstice at solar maximum) are presented graphically. Good agreement with DE-2 satellite measurements is demonstrated.

Killeen, T. L.↗

General method for assessing the geographically correlated error for an orbiter with applications to Topex

An estimation based computational algorithm for evaluating geographically dependent orbit determination errors is presented. A derivation is included of the generalized geographically correlated error analysis algorithm and the algorithm is applied to a Topex simulation. Then a general numerical technique which provides the capability to analyze any type of mismodeled parameter with the geopotential as a special case is described. The numerical algorithm is compared with the analytic method of Rosborough (1986) and it is shown that both methods generate consistent geographically dependent orbit error profiles when the ephemeris error is dominated by the mismodeled geopotential.

Wolff, P. J.↗

Topology of three-dimensional, variable density flows

This paper is concerned with the interpretation of unsteady, variable-density flow fields. The topology of the flow is determined by finding critical points and identifying the character of local solution trajectories. The time evolution of the flow is studied by following the paths of the critical points in the three-dimensional space of invariants of the local deformations tensor. The methodology can be applied to any smooth vector field and its associated gradient tensor including the vorticity and pressure gradient fields. This approach provides a framework for describing the geometry of complex flow patterns. Concisely summarizing that geometry in the space of invariants of the local gradient tensor may be a useful way of gaining insight into time-dependent processes described by large computational data bases. Applications to the descriptions of a flickering diffusion flame and a compressible wake are discussed.

Cantwell, Brian↗

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗

Asynchronous transfer mode link performance over ground networks

The results of an experiment to determine the feasibility of using asynchronous transfer mode (ATM) technology to support advanced spacecraft missions that require high-rate ground communications and, in particular, full-motion video are reported. Potential nodes in such a ground network include Deep Space Network (DSN) antenna stations, the Jet Propulsion Laboratory, and a set of national and international end users. The experiment simulated a lunar microrover, lunar lander, the DSN ground communications system, and distributed science users. The users were equipped with video-capable workstations. A key feature was an optical fiber link between two high-performance workstations equipped with ATM interfaces. Video was also transmitted through JPL's institutional network to a user 8 km from the experiment. Variations in video depending on the networks and computers were observed, the results are reported.

Chow, E. T.↗

Finite difference time domain analysis of chirped dielectric gratings

The finite difference time domain (FDTD) method for solving Maxwell's time-dependent curl equations is accurate, computationally efficient, and straight-forward to implement. Since both time and space derivatives are employed, the propagation of an electromagnetic wave can be treated as an initial-value problem. Second-order central-difference approximations are applied to the space and time derivatives of the electric and magnetic fields providing a discretization of the fields in a volume of space, for a period of time. The solution to this system of equations is stepped through time, thus, simulating the propagation of the incident wave. If the simulation is continued until a steady-state is reached, an appropriate far-field transformation can be applied to the time-domain scattered fields to obtain reflected and transmitted powers. From this information diffraction efficiencies can also be determined. In analyzing the chirped structure, a mesh is applied only to the area immediately around the grating. The size of the mesh is then proportional to the electric size of the grating. Doing this, however, imposes an artificial boundary around the area of interest. An absorbing boundary condition must be applied along the artificial boundary so that the outgoing waves are absorbed as if the boundary were absent. Many such boundary conditions have been developed that give near-perfect absorption. In this analysis, the Mur absorbing boundary conditions are employed. Several grating structures were analyzed using the FDTD method.

Hochmuth, Diane H.↗

Non-local sub-characteristic zones of influence in unsteady interactive boundary-layers

The properties of incompressible, unsteady, interactive, boundary layers are examined for a model hypersonic boundary layer and internal flow past humps or, equivalently, external flow past short-scaled humps. Using a linear high frequency analysis, it is shown that the domains of dependence within the viscous sublayer may be a strong function of position within the sublayer and may be strongly influenced by the pressure displacement interaction, or the prescribed displacement condition. Detailed calculations are presented for the hypersonic boundary layer. This effect is found to carry over directly to the fully viscous problem as well as the nonlinear problem. In the fully viscous problem, the non-local character of the domains of dependence manifests itself in the sub-characteristics. Potential implications of the domain of dependence structure on finite difference computations of unsteady boundary layers are briefly discussed.

Rothmayer, A. P.↗

Modeling coronal evolution

Helmet streamers on the sun were observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. A two dimensional (azimuthally symmetric) helmet streamer configuration was developed by computing solutions of the time dependent magnetohydrodynamic (MHD) equations, for a specified magnetic flux distribution on the sun. The helmet streamer is not symmetric about the equator. The evolution of the configuration, when differential rotation is applied, was investigated. It was found that after many rotations the configuration does not reach a steady state, but disrupts recurrently with the ejection of a plasmoid. These results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated.

Linker, Jon A.↗

On the decay of outbursts in dwarf novae nad X-ray novae

We perform computations using a time-dependent model for the accretion disk limit-cycle mechanism to examine the decay of the optical light following the peak of a dwarf nova outburst. We present the results of a parameter study of the physical input variables which affect the decay rate. In the model, the decay is brought about by a cooling transition front which begins at large radii in the disk and moves inward. The nature of the decay is strongly influenced by the radial dependence of the accretion disk viscosity parameter alpha. To obtain exponential decays for typical dwarf nova parameters, we require alpha proportional to r(exp epsilon(sub 0)), where epsilon(sub 0) approximately = 0.3-0.4. The exact value of epsilon(sub 0) which produces exponential decays depends on factors such as the mass of the accreting star and the inner radius of the accretion disk. Therefore, the observed ubiquity of exponential decays in two different types of systems (dwarf novae and X-ray novae) leads us to believe that alpha is an unnatural scaling for the viscosity. The physics of the cooling transition front must be self-regulating in that the timescale (-parital derivative of lnSigma(r)/partial derivative +)(exp -1) (where Sigma is the surface density) for mass extraction across the front remains constant. This may be consistent with a scaling alpha proportional to (h/r)(exp n), where h is the local disk semi-thickness and n approximately 1-2. As regards the speed of the cooling front, we find v(sub F)(r) proportional to r(exp p), where p approximately 3 at large radii, with an abrupt transition to p approximately 0 at some smaller radius. The r(exp 3) dependence is much steeper than has been found by previous workers and appears to result from the strong variation of specific heat within the cooling front when the front resides at a large radius in the disk. The outflow of disk material across the cooling front causes a significant departure of dln T(sub dff0/dln r from the standard value of -0.75 (expected from steady state accretion) within about 0.2 dex in radius of the break associated with the cooling front -- T(sub eff) aproximately 10(exp 3.9) K (r/10(exp 10 cm)) (exp -0.1). These effects should be observable with eclipse mapping. Finally, it appears that the relatively slow decay rate for the optical flux in the 1975 outburst of A0620-00 can be accounted for if the primary is a approximately 10 Solar mass black hole.

Cannizzo, John K.↗

The high-state/low-state transition in V794 Aquilae

We present a V-magnitude light curve of the cataclysmic binary V794 Quilae covering an 800 day time span. The system shows variations of ups to 3 mag. In particular, there are two dips of approximately 1 mag which last for about 50 days, and are followed by a rapid recovery to the original brightness, and a third dip of approximately 3 mag lasting about 100 days. These fluctuations are thought to be caused by the response of the accretion disk in the system to the cessation of mass transfer from the mass-losing secondary star. We present computations using a time-dependent accretion disk code to delineate a range of allowed behavior for th accretion disk. To model the observed light curve, we require that the model parameters alpha(sub cold) and alpha(sub hot), which characterize the degree of coupling of the viscous dissipation to the orbital shear in the low and high states of the accretion disk, be smaller than they have been inferred to be in the dwarf novae, and that the ratio alpha(sub hot)/alpha(sub cold) also be smaller. The fact that alpha seems to vary with the rate of mass transfer may provide support for the model of Vishniac and Diamond, in which the impact of the mass transfer stream onto the outer edge of the accretion disk excites internal waves which transport angular momentum outward and provide the viscous dissipation.

Honeycutt, R. Kent↗

Carbon Monoxide Formation in SN 1987A

The formation of carbon monoxide in the supernova SN 1987A at 200 days is investigated. Rather than the one-zone that have previous been employed, this work utilizes a radially dependent radiative transport model to compute the temperature and ionization structure and the destruction of CO and CO+ by ultraviolet radiation. The CO profile is computed assuming chemical equilibrium. Two models for the density and abundance profiles are examined: model 10H (unmixed) and model 10HMM (mixed) of Pinto & Woosley. Sensitivity to adopted rates and temperature as it might be affected by CO cooling is examined. The models give a total mass of 10-4 M_sun at 200 days for the unmixed model, which is comparable with the estimated observed abundance by Syromilio et al. but nearly 2 orders of magnitude less than the estimated observed and that computed in thermal-chemical models by Liu & Dalgarno. The effect of different model assumptions and results concerning ionization structure and radiative transfer are discussed. We confirm that CO+ is not expected to be produced in significant amounts and that the amount of CO is sensitive to the degree of the mixing of the composition of various elements present in the ejecta.

Gearhart, Rob A.↗

Development of a Dynamically Configurable, Object-Oriented Framework for Distributed, Multi-modal Computational Aerospace Systems Simulation: Second Year Progress Report

Mesh generation has long been recognized as a bottleneck in the CFD process. While much research on automating the volume mesh generation process have been relatively successful,these methods rely on appropriate initial surface triangulation to work properly. Surface discretization has been one of the least automated steps in computational simulation due to its dependence on implicitly defined CAD surfaces and curves. Differences in CAD peometry engines manifest themselves in discrepancies in their interpretation of the same entities. This lack of "good" geometry causes significant problems for mesh generators, requiring users to "repair" the CAD geometry before mesh generation. The problem is exacerbated when CAD geometry is translated to other forms (e.g., IGES )which do not include important topological and construction information in addition to entity geometry. One technique to avoid these problems is to access the CAD geometry directly from the mesh generating software, rather than through files. By accessing the geometry model (not a discretized version) in its native environment, t h s a proach avoids translation to a format which can deplete the model of topological information. Our approach to enable models developed in the Denali software environment to directly access CAD geometry and functions is through an Application Programming Interface (API) known as CAPRI. CAPRI provides a layer of indirection through which CAD-specific data may be accessed by an application program using CAD-system neutral C and FORTRAN language function calls. CAPRI supports a general set of CAD operations such as truth testing, geometry construction and entity queries.

Afjeh, Abdollah A.↗

Kinetic dissection of individual steps in the poly(C)-directed oligoguanylate synthesis from guanosine 5'-monophosphate 2-methylimidazolide

A kinetic study of oligoguanylate synthesis on a polycytidylate template, poly(C), as a function of the concentration of the activated monomer, guanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpG, is reported. Reactions were run with 0.005-0.045 M 2-MeImpG in the presence of 0.05 M poly(C) at 23 degrees C. The kinetic results are consistent with a reaction scheme (eq 1) that consists of a series of consecutive steps, each step representing the addition of one molecule of 2-MeImpG to the growing oligomer. This scheme allows the calculation of second-order rate constants for every step by analyzing the time-dependent growth of each oligomer. Computer simulations of the course of reaction based on the determined rate constants and eq 1 are in excellent agreement with the product distributions seen in the HPLC profiles. In accord with an earlier study (Fakhrai, H.; Inoue, T.; Orgel, L. E. Tetrahedron 1984, 40, 39), rate constants, ki, for the formation of the tetramer and longer oligomers up to the 16-mer were found to be independent of length and somewhat higher than k3 (formation of trimer), which in turn is much higher than k2 (formation of dimer). The ki (i > or = 4), k3, and k2 values are not true second-order rate constants but vary with monomer concentration. Mechanistic models for the dimerization (Scheme I) and elongation reactions (Scheme II) are proposed that are consistent with our results. These models take into account that the monomer associates with the template in a cooperative manner. Our kinetic analysis allowed the determination of rate constants for the elementary processes of covalent bond formation between two monomers (dimerization) and between an oligomer and a monomer (elongation) on the template. A major conclusion from our study is that bond formation between two monomer units or between a primer and a monomer is assisted by the presence of additional next-neighbor monomer units. This is consistent with recent findings with hairpin oligonucleotides (Wu, T.; Orgel, L. E. J. Am. Chem. Soc. 1992, 114, 317). Our study is the first of its kind that shows the feasibility of a thorough kinetic analysis of a template-directed oligomerization and provides a detailed mechanistic model of these reactions.

Non-NASA Center↗

Thermophysical Properties of Selected II-VI Semiconducting Melts

Thermophysical properties are essential for the accurate predication of the crystal growth process by computational modeling. Currently, the temperature dependent thermophysical property data for the II-VI semiconductor melts are scarce. This paper reports the results of the temperature dependence of melt density, viscosity and electrical conductivity of selected II-VI compounds, including HgTe, HgCdTe and HgZnTe. The melt density was measured using a pycnometric method, and the viscosity and electrical conductivity were measured by a transient torque method. The results were compared with and showed good agreement with the existing data in the literature.

Li, C.↗

Modeling of Radiative Transfer in Protostellar Disks

This program implements a spectral line, radiative transfer tool for interpreting Spitzer Space Telescope observations by matching them with models of protostellar disks for improved understanding of planet and star formation. The Spitzer Space Telescope detects gas phase molecules in the infrared spectra of protostellar disks, with spectral lines carrying information on the chemical composition of the material from which planets form. Input to the software includes chemical models developed at JPL. The products are synthetic images and spectra for comparison with Spitzer measurements. Radiative transfer in a protostellar disk is primarily affected by absorption and emission processes in the dust and in molecular gases such as H2, CO, and HCO. The magnitude of the optical absorption and emission is determined by the population of the electronic, vibrational, and rotational energy levels. The population of the molecular level is in turn determined by the intensity of the radiation field. Therefore, the intensity of the radiation field and the population of the molecular levels are inter-dependent quantities. To meet the computational challenges of solving for the coupled radiation field and electronic level populations in disks having wide ranges of optical depths and spatial scales, the tool runs in parallel on the JPL Dell Cluster supercomputer with C++ and Fortran compiler with a Message Passing Interface. Because this software has been developed on a distributed computing platform, the modeling of systems previously beyond the reach of available computational resources is possible.

VonAllmen, Paul↗