Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random structure generation”

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 127 records · Page 7

Step-by-Step Simulation of Radiation Chemistry Using Green Functions for Diffusion-Influenced Reactions

Radiolytic species are formed approximately 1 ps after the passage of ionizing radiation through matter. After their formation, they diffuse and chemically react with other radiolytic species and neighboring biological molecules, leading to various oxidative damage. Therefore, the simulation of radiation chemistry is of considerable importance to understand how radiolytic species damage biological molecules [1]. The step-by-step simulation of chemical reactions is difficult, because the radiolytic species are distributed non-homogeneously in the medium. Consequently, computational approaches based on Green functions for diffusion-influenced reactions should be used [2]. Recently, Green functions for more complex type of reactions have been published [3-4]. We have developed exact random variate generators of these Green functions [5], which will allow us to use them in radiation chemistry codes. Moreover, simulating chemistry using the Green functions is which is computationally very demanding, because the probabilities of reactions between each pair of particles should be evaluated at each timestep [2]. This kind of problem is well adapted for General Purpose Graphic Processing Units (GPGPU), which can handle a large number of similar calculations simultaneously. These new developments will allow us to include more complex reactions in chemistry codes, and to improve the calculation time. This code should be of importance to link radiation track structure simulations and DNA damage models.

Plante, Ianik↗

Statistical characterization of life drivers for a probabilistic design analysis

This paper discusses the issue of statistical characterization of life drivers for a probabilistic design analysis (PDA) approach to support the conventional deterministic structural design methods that are currently used. The probabilistic approach takes into consideration the modeling inadequacies and uncertainties in many design variables such as loads, environments, and material properties. The importance of the distributional assumption is motivated by illustrating an example where the results differ substantially due to the distribution selected. Different types of distributions are discussed and techniques for estimating the parameters are given. Given this information, procedures are outlined for selecting the appropriate distribution based on the particular type of variable (i.e., dimensional, performance) as well as the information that is available (i.e., test data, engineering analysis). Finally, techniques are given for generating random numbers from these selected distributions within the PDA process.

Fox, Eric P.↗

Planetary image conversion task

The Planetary Image Conversion Task group processed 12,500 magnetic tapes containing raw imaging data from JPL planetary missions and produced an image data base in consistent format on 1200 fully packed 6250-bpi tapes. The output tapes will remain at JPL. A copy of the entire tape set was delivered to US Geological Survey, Flagstaff, Ariz. A secondary task converted computer datalogs, which had been stored in project specific MARK IV File Management System data types and structures, to flat-file, text format that is processable on any modern computer system. The conversion processing took place at JPL's Image Processing Laboratory on an IBM 370-158 with existing software modified slightly to meet the needs of the conversion task. More than 99% of the original digital image data was successfully recovered by the conversion task. However, processing data tapes recorded before 1975 was destructive. This discovery is of critical importance to facilities responsible for maintaining digital archives since normal periodic random sampling techniques would be unlikely to detect this phenomenon, and entire data sets could be wiped out in the act of generating seemingly positive sampling results. Reccomended follow-on activities are also included.

Martin, M. D.↗

Investigation of Near Shannon Limit Coding Schemes

Turbo codes can deliver performance that is very close to the Shannon limit. This report investigates algorithms for convolutional turbo codes and block turbo codes. Both coding schemes can achieve performance near Shannon limit. The performance of the schemes is obtained using computer simulations. There are three sections in this report. First section is the introduction. The fundamental knowledge about coding, block coding and convolutional coding is discussed. In the second section, the basic concepts of convolutional turbo codes are introduced and the performance of turbo codes, especially high rate turbo codes, is provided from the simulation results. After introducing all the parameters that help turbo codes achieve such a good performance, it is concluded that output weight distribution should be the main consideration in designing turbo codes. Based on the output weight distribution, the performance bounds for turbo codes are given. Then, the relationships between the output weight distribution and the factors like generator polynomial, interleaver and puncturing pattern are examined. The criterion for the best selection of system components is provided. The puncturing pattern algorithm is discussed in detail. Different puncturing patterns are compared for each high rate. For most of the high rate codes, the puncturing pattern does not show any significant effect on the code performance if pseudo - random interleaver is used in the system. For some special rate codes with poor performance, an alternative puncturing algorithm is designed which restores their performance close to the Shannon limit. Finally, in section three, for iterative decoding of block codes, the method of building trellis for block codes, the structure of the iterative decoding system and the calculation of extrinsic values are discussed.

Kwatra, S. C.↗

Investigation of spectral analysis techniques for randomly sampled velocimetry data

It is well known that velocimetry (LV) generates individual realization velocity data that are randomly or unevenly sampled in time. Spectral analysis of such data to obtain the turbulence spectra, and hence turbulence scales information, requires special techniques. The 'slotting' technique of Mayo et al, also described by Roberts and Ajmani, and the 'Direct Transform' method of Gaster and Roberts are well known in the LV community. The slotting technique is faster than the direct transform method in computation. There are practical limitations, however, as to how a high frequency and accurate estimate can be made for a given mean sampling rate. These high frequency estimates are important in obtaining the microscale information of turbulence structure. It was found from previous studies that reliable spectral estimates can be made up to about the mean sampling frequency (mean data rate) or less. If the data were evenly samples, the frequency range would be half the sampling frequency (i.e. up to Nyquist frequency); otherwise, aliasing problem would occur. The mean data rate and the sample size (total number of points) basically limit the frequency range. Also, there are large variabilities or errors associated with the high frequency estimates from randomly sampled signals. Roberts and Ajmani proposed certain pre-filtering techniques to reduce these variabilities, but at the cost of low frequency estimates. The prefiltering acts as a high-pass filter. Further, Shapiro and Silverman showed theoretically that, for Poisson sampled signals, it is possible to obtain alias-free spectral estimates far beyond the mean sampling frequency. But the question is, how far? During his tenure under 1993 NASA-ASEE Summer Faculty Fellowship Program, the author investigated from his studies on the spectral analysis techniques for randomly sampled signals that the spectral estimates can be enhanced or improved up to about 4-5 times the mean sampling frequency by using a suitable prefiltering technique. But, this increased bandwidth comes at the cost of the lower frequency estimates. The studies further showed that large data sets of the order of 100,000 points, or more, high data rates, and Poisson sampling are very crucial for obtaining reliable spectral estimates from randomly sampled data, such as LV data. Some of the results of the current study are presented.

Sree, Dave↗

How thick are lunar mare basalts

It is argued that De Hon's estimates of the thickness of lunar mare basalts, made by analyzing 'ghost' craters on mare surfaces, were inflated as the result of the crater morphometric data of Pike (1977) to reconstruct rim heights of degraded craters. Crater rim heights of 82 randomly selected highland craters of various states of degradation were determined, and median rim height was compared to that of corresponding fresh impact structures. Results indicate that the thickness estimates of De Hon may be reduced by a factor of 2, and that the total volume of mare basalt produced throughout lunar history could be as little as 1-2 million cubic kilometers. A survey of geochemical and petrographic evidence indicates that lateral transport of regolith components over distances of much greater than 10 km is relatively inefficient; it is suggested that vertical mixing of a highland substrate underlying the basaltic fill may have had a primordial role in generating the observed mare width distributions and high concentrations of exotic components in intrabasin regoliths.

Hoerz, F.↗

Microscopic and Metallurgical Aspects of the Space Shuttle Columbia Accident Investigation and Reconstruction

The Space Shuttle Columbia was descending for a landing at the Kennedy Space Center (KSC) on February 1, 2003. Approximately 20 minutes prior to touchdown, the Columbia began disintegrating over the western United States; the majority of debris eventually impacted in eastern Texas and western Louisiana. A monumental effort eventually recovered approximately 84,000 pieces of debris, approximately 38% of the Orbiter's original dry weight. The debris was transported to KSC, where the items were catalogued and evaluated. Critical areas of interest, such as the left and right leading edge surfaces and the underside of the ship, were placed upon a grid to aid in the reconstruction. Items of interest included metallic structures, reinforced carbon-carbon composites, and ceramic heat insulation tiles. Many of the leading edge elements had re-solidified metallic deposits spattered on them. These deposits became known as slag and were one of the main focuses of the investigation. In order to help determine the sequence of events inside the left wing during the accident, the slag's composition, layering order, and directionality of deposition were studied. A myriad of analytical tests were performed in an attempt to ascertain the compositional and depositional characteristics of selected slag deposits, including the ordering of deposited layers within each individual slag deposit harvested. Initially, Scanning Electron Microscopy and Energy Dispersive X-Ray Spectroscopy (SEM/EDX) were performed to quickly characterize the overall composition of individual slag deposits: SEM utilizes a narrowlyfocused high-energy electron beam impinging upon a specimen. The incident beam excites and liberates lower energy secondary electrons, which are detected and analyzed, providing a visual representation of the sample's surface topography. EDX also relies on an incident electron beam, except an EDX unit measures X-ray energies generated by the impinging beam. Each element generates a unique X-ray signature; the EDX detector measures these discreet energies. EDX actually penetrates approximately 2 microns into the bulk of the sample. However, random examination of various portions of slag, coupled with the semiquantitative nature of the SEM/EDX analysis, did not yield convincingly pertinent data. Therefore, X-ray dot mapping was conducted, which provided more understandable data, both in terms of slag layering and composition. An X-ray dot map is generated by performing numerous EDX scans for individual elements, then compiling the scans in a visual representation. Eventually, specimens consisting of not only the slag, but of the adjacent RCC substrate as well were cross-sectioned. X-Ray dot mapping of the materialographicallymounted and -polished cross- sections provided a visual representation of both the layering sequence and compositional characteristics of the slag. Contemporaneously, Electron Spectroscopy for Chemical Analysis/X-Ray Photoelectron Spectroscopy (ESCA/XPS) and powdered X-Ray Diffraction (XRD) were performed to further characterize the deposits and to attempt to identify what, if any, compounds were present. The ESCA/XPS analysis allowed the analyst to "sputter" into the sample with an electron gun, aiding in the identification of the layering sequence. XPS uses photons, rather than electrons, which impinge upon the surface of the sample. XPS measures the electrons emitted from within the first 5 nm of the sample's surface. The XRD measures the scatter angles of incident X-rays; the angle and intensity of scatter depend upon the crystalline structure of the pulverized sample. XRD is considered a qualitative rather than quantitative technique. ESCA/XPS revealed that the final layer to deposit was predominantly carbonaceous. XRD was successful in identifying specific compounds, such as Al 2O3, Al and/or Al3 21SiO47, mullite (3(Al2)O3 -SiO2), and nickel-aluminides. Eventually, Electron MicroProbe Analysis (EMPA) was conducted on the marialographically-prepared cross- sections of selected slag deposits. Microprobe combines SEM and Wavelength Dispersive X-Ray Spectroscopy (WDS), and, like EDX, uses a narrowly-focused high-energy electron beam impinging upon a specimen to elicit, in the case of EPMA, characteristic X-rays with specific wavelengths. This quantitative, analytical tool proved the most useful in determining depositional layering and composition of the slag deposits. This information was utilized in verifying the location of the breach in the left leading edge of the wing of the Columbia.

McDaniels, Steven J.↗

Numerical Modeling of Turbulent Combustion

The work in numerical modeling is focused on the use of the random vortex method to treat turbulent flow fields associated with combustion while flame fronts are considered as interfaces between reactants and products, propagating with the flow and at the same time advancing in the direction normal to themselves at a prescribed burning speed. The latter is associated with the generation of specific volume (the flame front acting, in effect, as the locus of volumetric sources) to account for the expansion of the flow field due to the exothermicity of the combustion process. The model was applied to the flow in a channel equipped with a rearward facing step. The results obtained revealed the mechanism of the formation of large scale turbulent structure in the wake of the step, while it showed the flame to stabilize on the outer edges of these eddies.

Ghoneim, A. F.↗

Allowing for Horizontally Heterogeneous Clouds and Generalized Overlap in an Atmospheric GCM

While fully accounting for 3D effects in Global Climate Models (GCMs) appears not realistic at the present time for a variety of reasons such as computational cost and unavailability of 3D cloud structure in the models, incorporation in radiation schemes of subgrid cloud variability described by one-point statistics is now considered feasible and is being actively pursued. This development has gained momentum once it was demonstrated that CPU-intensive spectrally explicit Independent Column Approximation (lCA) can be substituted by stochastic Monte Carlo ICA (McICA) calculations where spectral integration is accomplished in a manner that produces relatively benign random noise. The McICA approach has been implemented in Goddard's GEOS-5 atmospheric GCM as part of the implementation of the RRTMG radiation package. GEOS-5 with McICA and RRTMG can handle horizontally variable clouds which can be set via a cloud generator to arbitrarily overlap within the full spectrum of maximum and random both in terms of cloud fraction and layer condensate distributions. In our presentation we will show radiative and other impacts of the combined horizontal and vertical cloud variability on multi-year simulations of an otherwise untuned GEOS-5 with fixed SSTs. Introducing cloud horizontal heterogeneity without changing the mean amounts of condensate reduces reflected solar and increases thermal radiation to space, but disproportionate changes may increase the radiative imbalance at TOA. The net radiation at TOA can be modulated by allowing the parameters of the generalized overlap and heterogeneity scheme to vary, a dependence whose behavior we will discuss. The sensitivity of the cloud radiative forcing to the parameters of cloud horizontal heterogeneity and comparisons of CERES-derived forcing will be shown.

Lee, D.↗

Vibration Modal Characterization of a Stirling Convertor via Base-Shake Excitation

The U.S. Department of Energy (DOE), Lockheed Martin (LM), Stirling Technology Company (STC), and NASA John H. Glenn Research Center (GRC) are currently developing a high-efficiency Stirling convertor for use in a Stirling Radioisotope Generator (SRG). NASA and DOE have identified the SRG for potential use as an advanced power system for future NASA Space Science missions, providing spacecraft onboard electric power for deep space missions and power for unmanned Mars rovers. Low-level, baseshake sine vibration tests were conducted on the Stirling Technology Demonstration Convertor (TDC), at NASA GRC's Structural Dynamics Laboratory, in February 2001, as part of the development of this Stirling technology. The purpose of these tests was to provide a better understanding of the TDC's internal dynamic response to external vibratory base excitations. The knowledge obtained can therein be used to help explain the success that the TDC enjoyed in its previous random vibration qualification tests (December 1999). This explanation focuses on the TDC s internal dynamic characteristics in the 50 to 250 Hz frequency range, which corresponds to the maximum input levels of its qualification random vibration test specification. The internal dynamic structural characteristics of the TDC have now been measured in two separate tests under different motoring and dynamic loading conditions: (1) with the convertor being electrically motored, under a vibratory base-shake excitation load, and (2) with the convertor turned off, and its alternator internals undergoing dynamic excitation via hammer impact loading. This paper addresses the test setup, procedure and results of the base-shake vibration testing conducted on the motored TDC, and will compare these results with those results obtained from the dynamic impact tests (May 2001) on the nonmotored TDC.

Suarez, Vicente J.↗

On the phase-locked polarization variations in Cygnus X-1

The paper presents data on the accumulated V-band orbital phase curves of polarization in Cygnus X-1 which disagree with the analysis of Simmons, Aspin, and Brown (1980). No curves which conform to the canonical model were observed, and a noncanonical structure is indicated by F tests on third, fourth, and fifth harmonics, which have a probability of 0.85 as a group. The fifth harmonic has a probability of 0.99; in contrast, the fundamental harmonic has a probability of only 0.6, which is below the detection threshold. Therefore, major conclusions of Simmons et al based on a nonsymmetrical canonical model assuming that fundamental harmonics are present, but harmonics with n greater than 2 are absent, are unacceptable. Analysis of the author's data by Simmons et al is criticized on the grounds that he did not take into account the cumulative aspects of the data and did not recognize the need to average out the random noise generated by Cygnus X-1 over many orbital cycles.

Kemp, J. C.↗

Turbulent Compressible Convection with Rotation: Mean Flows and Differential Rotation - 2

The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Brummell, Nicholas H.↗

Some Interesting Applications of Probabilistic Techiques in Structural Dynamic Analysis of Rocket Engines

Numerical and Analytical methods developed to determine damage accumulation in specific engine components when speed variation included. Dither Life Ratio shown to be well over factor of 2 for specific example. Steady‐State assumption shown to be accurate for most turbopump cases, allowing rapid calculation of DLR. If hot‐fire speed data unknown, Monte Carlo method developed that uses speed statistics for similar engines. Application of techniques allow analyst to reduce both uncertainty and excess conservatism. High values of DLR could allow previously unacceptable part to pass HCF criteria without redesign. Given benefit and ease of implementation, recommend that any finite life turbomachine component analysis adopt these techniques. Probability Values calculated, compared, and evaluated for several industry‐proposed methods for combining random and harmonic loads. Two new excel macros written to calculate combined load for any specific probability level. Closed form Curve fits generated for widely used 3(sigma) and 2(sigma) probability levels. For design of lightweight aerospace components, obtaining accurate, reproducible, statistically meaningful answer critical.

Brown, Andrew M.↗

Interstellar turbulence, random density variations, and scintillation measurements

The presence of random electron variations suggests that the ionized interstellar medium is turbulent. In the interstellar plasma the presence of power spectra of such variations extending to spatial scales much less than a Coulomb mean free path, Lambda sub c, is required by analyses of measurements of scintillation and angular broadening of pulsar radio signals. The existence of corresponding variations in magnetic field strength could efficiently scatter cosmic rays and thus constrain cosmic-ray propagation. Unfortunately both the origin of the electron density variations and mechanisms by which these variations couple to fluctuations in magnetic field strength are unknown. It is conjectured that the small-scale density variations are generated by the convective distortion of initially large-scale isobaric entropy structures in the turbulent interstellar plasma. An investigation of the spectra of turbulent entropy structures, velocity, and magnetic fields at small spatial scales is made. The modifier small is employed to characterize length scales much less than the dimension, L, containing the bulk of the turbulent energy.

Higdon, J. C.↗

Understanding the Scalability of Bayesian Network Inference Using Clique Tree Growth Curves

One of the main approaches to performing computation in Bayesian networks (BNs) is clique tree clustering and propagation. The clique tree approach consists of propagation in a clique tree compiled from a Bayesian network, and while it was introduced in the 1980s, there is still a lack of understanding of how clique tree computation time depends on variations in BN size and structure. In this article, we improve this understanding by developing an approach to characterizing clique tree growth as a function of parameters that can be computed in polynomial time from BNs, specifically: (i) the ratio of the number of a BN s non-root nodes to the number of root nodes, and (ii) the expected number of moral edges in their moral graphs. Analytically, we partition the set of cliques in a clique tree into different sets, and introduce a growth curve for the total size of each set. For the special case of bipartite BNs, there are two sets and two growth curves, a mixed clique growth curve and a root clique growth curve. In experiments, where random bipartite BNs generated using the BPART algorithm are studied, we systematically increase the out-degree of the root nodes in bipartite Bayesian networks, by increasing the number of leaf nodes. Surprisingly, root clique growth is well-approximated by Gompertz growth curves, an S-shaped family of curves that has previously been used to describe growth processes in biology, medicine, and neuroscience. We believe that this research improves the understanding of the scaling behavior of clique tree clustering for a certain class of Bayesian networks; presents an aid for trade-off studies of clique tree clustering using growth curves; and ultimately provides a foundation for benchmarking and developing improved BN inference and machine learning algorithms.

Mengshoel, Ole J.↗

Microwave backscattering and emission model for grass canopies

A two-layer model is developed that treats the grass canopy as a collection of randomly oriented elliptical dielectric discs over a layer of thatch with underlying soil surface. The distorted Born approximation in conjunction with the Peake formulation is used to calculate the backscattering coefficient and the emissivity from the canopy. Two particular features of this model which are unique for grass canopies are the variation of the canopy structure and the presence of the thatch layer. The basic parameters in the model such as the size and orientation of grass blades, dielectric constant of soil and vegetation, and thickness and water content of the thatch layer have been obtained from ground truth data. To interpret the available experimental observations of grasslands, numerical results from both passive and active models at L-band (1.4 GHz) are generated and various scattering and emission properties of the grass canopies are discussed.

Saatchi, Sasan S.↗

Incorporation of Half-Cycle Theory Into Ko Aging Theory for Aerostructural Flight-Life Predictions

The half-cycle crack growth theory was incorporated into the Ko closed-form aging theory to improve accuracy in the predictions of operational flight life of failure-critical aerostructural components. A new crack growth computer program was written for reading the maximum and minimum loads of each half-cycle from the random loading spectra for crack growth calculations and generation of in-flight crack growth curves. The unified theories were then applied to calculate the number of flights (operational life) permitted for B-52B pylon hooks and Pegasus adapter pylon hooks to carry the Hyper-X launching vehicle that air launches the X-43 Hyper-X research vehicle. A crack growth curve for each hook was generated for visual observation of the crack growth behavior during the entire air-launching or captive flight. It was found that taxiing and the takeoff run induced a major portion of the total crack growth per flight. The operational life theory presented can be applied to estimate the service life of any failure-critical structural components.

Ko, William L.↗

Development and application of an atmospheric turbulence model for use in flight simulators in flight simulators

The influence of simulated turbulence on aircraft handling qualities was investigated. Pilot opinion of the handling qualities of a light general aviation aircraft were evaluated in a motion-base simulator using a simulated turbulence environment. A realistic representation of turbulence disturbances is described in terms of rms intensity and scale length and their random variations with time. The time histories generated by the proposed turbulence models showed characteristics which appear to be more similar to real turbulence than the frequently-used Gaussian turbulence model. In addition, the proposed turbulence models can flexibly accommodate changes in atmospheric conditions and be easily implemented in flight simulator studies. Six turbulence time histories, including the conventional Gaussian model, were used in an IFR-tracking task. The realism of each of the turbulence models and the handling qualities of the simulated airplane were evaluated. Analysis of pilot opinions shows that at approximately the same rms intensities of turbulence, the handling quality ratings transit from the satisfactory level, for the simple Gaussian model, to an unacceptable level for more realistic and compositely structured turbulence models.

Jacobson, I. D.↗