Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “random number generators”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS2011

Origins of Protein Functions in Cells

In modern organisms proteins perform a majority of cellular functions, such as chemical catalysis, energy transduction and transport of material across cell walls. Although great strides have been made towards understanding protein evolution, a meaningful extrapolation from contemporary proteins to their earliest ancestors is virtually impossible. In an alternative approach, the origin of water-soluble proteins was probed through the synthesis and in vitro evolution of very large libraries of random amino acid sequences. In combination with computer modeling and simulations, these experiments allow us to address a number of fundamental questions about the origins of proteins. Can functionality emerge from random sequences of proteins? How did the initial repertoire of functional proteins diversify to facilitate new functions? Did this diversification proceed primarily through drawing novel functionalities from random sequences or through evolution of already existing proto-enzymes? Did protein evolution start from a pool of proteins defined by a frozen accident and other collections of proteins could start a different evolutionary pathway? Although we do not have definitive answers to these questions yet, important clues have been uncovered. In one example (Keefe and Szostak, 2001), novel ATP binding proteins were identified that appear to be unrelated in both sequence and structure to any known ATP binding proteins. One of these proteins was subsequently redesigned computationally to bind GTP through introducing several mutations that introduce targeted structural changes to the protein, improve its binding to guanine and prevent water from accessing the active center. This study facilitates further investigations of individual evolutionary steps that lead to a change of function in primordial proteins. In a second study (Seelig and Szostak, 2007), novel enzymes were generated that can join two pieces of RNA in a reaction for which no natural enzymes are known. Recently it was found that, as in the previous case, the proteins have a structure unknown among modern enzymes. In this case, in vitro evolution started from a small, non-enzymatic protein. A similar selection process initiated from a library of random polypeptides is in progress. These results not only allow for estimating the occurrence of function in random protein assemblies but also provide evidence for the possibility of alternative protein worlds. Extant proteins might simply represent a frozen accident in the world of possible proteins. Alternative collections of proteins, even with similar functions, could originate alternative evolutionary paths.

Seelig, Burchard↗
NASA NTRS2000

A Parallel Genetic Algorithm for Automated Electronic Circuit Design

Parallelized versions of genetic algorithms (GAs) are popular primarily for three reasons: the GA is an inherently parallel algorithm, typical GA applications are very compute intensive, and powerful computing platforms, especially Beowulf-style computing clusters, are becoming more affordable and easier to implement. In addition, the low communication bandwidth required allows the use of inexpensive networking hardware such as standard office ethernet. In this paper we describe a parallel GA and its use in automated high-level circuit design. Genetic algorithms are a type of trial-and-error search technique that are guided by principles of Darwinian evolution. Just as the genetic material of two living organisms can intermix to produce offspring that are better adapted to their environment, GAs expose genetic material, frequently strings of 1s and Os, to the forces of artificial evolution: selection, mutation, recombination, etc. GAs start with a pool of randomly-generated candidate solutions which are then tested and scored with respect to their utility. Solutions are then bred by probabilistically selecting high quality parents and recombining their genetic representations to produce offspring solutions. Offspring are typically subjected to a small amount of random mutation. After a pool of offspring is produced, this process iterates until a satisfactory solution is found or an iteration limit is reached. Genetic algorithms have been applied to a wide variety of problems in many fields, including chemistry, biology, and many engineering disciplines. There are many styles of parallelism used in implementing parallel GAs. One such method is called the master-slave or processor farm approach. In this technique, slave nodes are used solely to compute fitness evaluations (the most time consuming part). The master processor collects fitness scores from the nodes and performs the genetic operators (selection, reproduction, variation, etc.). Because of dependency issues in the GA, it is possible to have idle processors. However, as long as the load at each processing node is similar, the processors are kept busy nearly all of the time. In applying GAs to circuit design, a suitable genetic representation 'is that of a circuit-construction program. We discuss one such circuit-construction programming language and show how evolution can generate useful analog circuit designs. This language has the desirable property that virtually all sets of combinations of primitives result in valid circuit graphs. Our system allows circuit size (number of devices), circuit topology, and device values to be evolved. Using a parallel genetic algorithm and circuit simulation software, we present experimental results as applied to three analog filter and two amplifier design tasks. For example, a figure shows an 85 dB amplifier design evolved by our system, and another figure shows the performance of that circuit (gain and frequency response). In all tasks, our system is able to generate circuits that achieve the target specifications.

Long, Jason D.↗
NASA NTRS2007

Automated Knowledge Discovery From Simulators

A computational method, SimLearn, has been devised to facilitate efficient knowledge discovery from simulators. Simulators are complex computer programs used in science and engineering to model diverse phenomena such as fluid flow, gravitational interactions, coupled mechanical systems, and nuclear, chemical, and biological processes. SimLearn uses active-learning techniques to efficiently address the "landscape characterization problem." In particular, SimLearn tries to determine which regions in "input space" lead to a given output from the simulator, where "input space" refers to an abstraction of all the variables going into the simulator, e.g., initial conditions, parameters, and interaction equations. Landscape characterization can be viewed as an attempt to invert the forward mapping of the simulator and recover the inputs that produce a particular output. Given that a single simulation run can take days or weeks to complete even on a large computing cluster, SimLearn attempts to reduce costs by reducing the number of simulations needed to effect discoveries. Unlike conventional data-mining methods that are applied to static predefined datasets, SimLearn involves an iterative process in which a most informative dataset is constructed dynamically by using the simulator as an oracle. On each iteration, the algorithm models the knowledge it has gained through previous simulation trials and then chooses which simulation trials to run next. Running these trials through the simulator produces new data in the form of input-output pairs. The overall process is embodied in an algorithm that combines support vector machines (SVMs) with active learning. SVMs use learning from examples (the examples are the input-output pairs generated by running the simulator) and a principle called maximum margin to derive predictors that generalize well to new inputs. In SimLearn, the SVM plays the role of modeling the knowledge that has been gained through previous simulation trials. Active learning is used to determine which new input points would be most informative if their output were known. The selected input points are run through the simulator to generate new information that can be used to refine the SVM. The process is then repeated. SimLearn carefully balances exploration (semi-randomly searching around the input space) versus exploitation (using the current state of knowledge to conduct a tightly focused search). During each iteration, SimLearn uses not one, but an ensemble of SVMs. Each SVM in the ensemble is characterized by different hyper-parameters that control various aspects of the learned predictor - for example, whether the predictor is constrained to be very smooth (nearby points in input space lead to similar output predictions) or whether the predictor is allowed to be "bumpy." The various SVMs will have different preferences about which input points they would like to run through the simulator next. SimLearn includes a formal mechanism for balancing the ensemble SVM preferences so that a single choice can be made for the next set of trials.

Burl, Michael↗
NASA NTRS1993

New perspectives for high accuracy SLR with second generation geodesic satellites

This paper reports on the accuracy limitations imposed by geodesic satellite signatures, and on the potential for achieving millimetric performances by means of alternative satellite concepts and an optimized 2-color system tradeoff. Long distance laser ranging, when performed between a ground (emitter/receiver) station and a distant geodesic satellite, is now reputed to enable short arc trajectory determinations to be achieved with an accuracy of 1 to 2 centimeters. This state-of-the-art accuracy is limited principally by the uncertainties inherent to single-color atmospheric path length correction. Motivated by the study of phenomena such as postglacial rebound, and the detailed analysis of small-scale volcanic and strain deformations, the drive towards millimetric accuracies will inevitably be felt. With the advent of short pulse (less than 50 ps) dual wavelength ranging, combined with adequate detection equipment (such as a fast-scanning streak camera or ultra-fast solid-state detectors) the atmospheric uncertainty could potentially be reduced to the level of a few millimeters, thus, exposing other less significant error contributions, of which by far the most significant will then be the morphology of the retroreflector satellites themselves. Existing geodesic satellites are simply dense spheres, several 10's of cm in diameter, encrusted with a large number (426 in the case of LAGEOS) of small cube-corner reflectors. A single incident pulse, thus, results in a significant number of randomly phased, quasi-simultaneous return pulses. These combine coherently at the receiver to produce a convolved interference waveform which cannot, on a shot to shot basis, be accurately and unambiguously correlated to the satellite center of mass. This paper proposes alternative geodesic satellite concepts, based on the use of a very small number of cube-corner retroreflectors, in which the above difficulties are eliminated while ensuring, for a given emitted pulse, the return of a single clean pulse with an adequate cross-section.

Lund, Glenn↗
NASA NTRSDate not supplied

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗
NASA NTRS1989

PRIME: A bottom-up approach to probabilistic rule development

PRIME is a system to be used by an intelligent machine to allow it to operate in an abstract but uncertain (or stochastic) environment. It maintains a model of the effects of the machine's actions in the form of a rule base, which is induced from experience. This bottom-up approach to rule development allows the model to adapt to changes in the environment. Each rule consists of a condition under which the rule is active, an action, the effect of the action on the environment, and an estimate of the probability of this effect occurring. The effect probabilities are used to model the uncertainty in the environment, permitting multiple possible effects for a single action under a particular set of conditions. The objective of the intelligent machine is to satisfy user-specified goals with maximum probability of success. PRIME fulfills this requirement in two ways: it continuously updates the rule base with the most recent information, to ensure the validity of the model; and it generates plans which have the maximum probability of achieving the goals, based on the probability estimates in the rule base. PRIME is composed of three main processes: exploration, generalization, and planning. In exploration, the machine executes various randomly chosen actions, observes the effects on the environment, and updates the rule base accordingly. This process is used to develop the rule base in simulation, as well as to supplement the current knowledge during normal operation. Generalization is the procedure used to induce general rules from experience, which is encoded in the form of specific rules. These general rules extend the machine's knowledge to situations which have not been encountered yet, thereby increasing the capability of the machine to plan effectively. Planning is the process of constructing an optimal sequence of actions to satisfy a goal, using the rule base to predict the effects of these actions and to determine the probability of success of the plan. The rule representation and many other data structures were specifically chosen to maximize the efficiency of these processes. A simulated environment was designed to test the performance of PRIME. The results of experimentation were largely negative. The main problem was that the domain coverage of the rules was inadequate for the number of rules stored in the rule base, due to redundancies in general rules and numerous rules covering ineffective actions. It was determined that a more efficient generalization, and some form of goal-directed exploration, are necessary in order to solve most of the current deficiencies in PRIME.

Miller, Scott A.↗
NASA NTRS1995

Computation of vertical profiles of longwave radiative cooling over the equatorial Pacific

An important quantity whose magnitude has not been throughly examined is the vertical distribution of heating in the Tropics. The details of the vertical distribution of heating have a significant impact on a number of phenomena, including the 30-60 day oscillation, sometimes known as the intraseasonal oscillation. Prior attempts to establish the structure of the heating relied on limited field data or assimilated data, coupled with climatological radiative heating parameters. The availability of high quality global-scale datasets has made it possible to make more accurate calculations than were possible a few years ago. An important component of the apparent heat budget is the longwave radiative cooling, which in this paper is found by using the ECMWF/WCRP/TOGA Archive 2 and ISCCP C1 datasets, together with a well-established parameterization scheme. A method is developed that can be used to estimate the vertical structure of cloud amounts based on top-of-atmosphere cloud observations, and the results are used with a wide-band longwave parameterization to produce longwave cooling rates over the tropical Pacific Ocean. Outgoing longwave radiation is calculated and compared the ERBE results. The calculated values are generally higher than those from ERBE, though the spatial distributions are similar. Some significant problems exist with the ECMWF upper-tropospheric water vapor amounts, which could imply uncertainties of 0.5 C/day in the calculated cooling rates. This is comparable to the differences associated with the minimum or random overlap assumptions used to generate cloud profiles.

Ramsey, Perry G.↗
NASA NTRS1999

Transduction-like gene transfer in the methanogen Methanococcus voltae

Strain PS of Methanococcus voltae (a methanogenic, anaerobic archaebacterium) was shown to generate spontaneously 4.4-kbp chromosomal DNA fragments that are fully protected from DNase and that, upon contact with a cell, transform it genetically. This activity, here called VTA (voltae transfer agent), affects all markers tested: three different auxotrophies (histidine, purine, and cobalamin) and resistance to BES (2-bromoethanesulfonate, an inhibitor of methanogenesis). VTA was most effectively prepared by culture filtration. This process disrupted a fraction of the M. voltae cells (which have only an S-layer covering their cytoplasmic membrane). VTA was rapidly inactivated upon storage. VTA particles were present in cultures at concentrations of approximately two per cell. Gene transfer activity varied from a minimum of 2 x 10(-5) (BES resistance) to a maximum of 10(-3) (histidine independence) per donor cell. Very little VTA was found free in culture supernatants. The phenomenon is functionally similar to generalized transduction, but there is no evidence, for the time being, of intrinsically viral (i.e., containing a complete viral genome) particles. Consideration of VTA DNA size makes the existence of such viral particles unlikely. If they exist, they must be relatively few in number;perhaps they differ from VTA particles in size and other properties and thus escaped detection. Digestion of VTA DNA with the AluI restriction enzyme suggests that it is a random sample of the bacterial DNA, except for a 0.9-kbp sequence which is amplified relative to the rest of the bacterial chromosome. A VTA-sized DNA fraction was demonstrated in a few other isolates of M. voltae.

Methanococcus/genetics/growth & development/metabo↗
NASA NTRS2004

Effects of Elevated Free-Stream Turbulence on Active Control of a Separation Bubble

The Effects of elevated free-stream turbulence (FST) on the natural and periodically excited separation bubbles were studied experimentally, due to the relevance of this flow to low-pressure turbine blades at low Reynolds numbers. A bubble was formed at the leading edge of a flat plate and the FST level was altered by placing a grid across the flow at different locations upstream of the plate. The mixing across the separated shear-layer, forming the free boundary of the bubble, increased due to the elevated FST and due to nominally two-dimensional periodic excitation, both flattening and shortening the bubble. Periodic excitation at frequencies that were at least an order of magnitude lower than those associated with the initial shear-layer instability, were very effective at low FST, because the amplitudes of the excitation frequency and its harmonic were amplified over the bubble. High frequency excitation (F+ 3, based on the length of the baseline low FST bubble) had a major effect close to the separation location, while farther downstream the excited fluctuations rapidly decayed in the reattachment region. Low frequency excitation, that generated waves comparable to the length of the unperturbed bubble (F+ 1) were less effective and their magnitude decayed at a slower rate downstream of reattachment. An increase in the level of the FST reduced the net effect of the periodic excitation on the mixing enhancement and subsequent reattachment process, probably due to a destructive interference between the nominally 2D excitation and the random (in space and time) FST, reducing the spanwise coherence and therefore the effectiveness of the current control strategy. However, even at the reduced effectiveness of 2D periodic excitation at elevated FST, it accelerated the reattachment process and the recovery rate of the reattached boundary layer, enhancing the boundary layer resistance to repeat separation and reducing its momentum loss further downstream.

Ashpis, D.↗
NASA NTRS2006

Earth Satellite Population Instability: Underscoring the Need for Debris Mitigation

A recent study by NASA indicates that the implementation of international orbital debris mitigation measures alone will not prevent a significant increase in the artificial Earth satellite population, beginning in the second half of this century. Whereas the focus of the aerospace community for the past 25 years has been on the curtailment of the generation of long-lived orbital debris, active remediation of the current orbital debris population should now be reconsidered to help preserve near-Earth space for future generations. In particular, we show in this paper that even if launch operations were to cease today, the population of space debris would continue to grow. Further, proposed remediation techniques do not appear to offer a viable solution. We therefore recommend that, while the aerospace community maintains the current debris-limiting mission regulations and postmission disposal procedures, future emphasis should be placed on finding new remediation technologies for solving this growing problem. Since the launch of Sputnik 1, space activities have created an orbital debris environment that poses increasing impact risks to existing space systems, including human space flight and robotic missions (1, 2). Currently, more than 9,000 Earth orbiting man-made objects (including many breakup fragments), with a combined mass exceeding 5 million kilograms, are tracked by the US Space Surveillance Network and maintained in the US satellite catalog (3-5). Three accidental collisions between cataloged satellites during the period from late 1991 to early 2005 have already been documented (6), although fortunately none resulted in the creation of large, trackable debris clouds. Several studies conducted during 1991-2001 demonstrated, with assumed future launch rates, the unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects (7-13). In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new satellites (i.e., debris) due to collisions exceeds the loss of objects due to orbital decay. NASA s evolutionary satellite population model LEGEND (LEO-to-GEO Environment Debris model), developed by the Orbital Debris Program Office at the NASA Lyndon B. Johnson Space Center, is a high fidelity three-dimensional physical model that is capable of simulating the historical satellite environment, as well as the evolution of future debris populations (14, 15). The subject study assumed no rocket bodies and spacecraft were launched after December 2004, and no future disposal maneuvers were allowed for existing spacecraft, few of which currently have such a capability. The rate of satellite explosions would naturally decrease to zero within a few decades as the current satellite population ages. The LEGEND future projection adopts a Monte Carlo approach to simulate future on-orbit explosions and collisions. Within a given projection time step, once the explosion probability is estimated for an intact object, a random number is drawn and compared with the probability to determine if an explosion would occur. A similar procedure is applied to collisions for each pair of target and projectile involved within the same time step. Due to the nature of the Monte Carlo process, multiple projection runs must be performed and analyzed before one can draw reliable and meaningful conclusions from the outcome. A total of fifty, 200-year future projection Monte Carlo simulations were executed and evaluated (16).

Liou, Jer-chyi↗
NASA NTRS2011

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engines is to decompose a CFD-generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non-harmonic excitation sources that become present in complex flows. A substantial effort has been made to account for this denser spatial Fourier content in frequency response analysis (described in another paper by the author), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements. Six loading cases were generated by varying a baseline harmonic excitation in different ways based upon cold-flow testing from Heritage Fuel Air Turbine Test. It was hypothesized that the randomness and other variation from the standard harmonic excitation would reduce the blade structural response, but the results showed little reduction. The second study was of a realistic model of a bladed-disk excited by the same CFD used in the J2X engine program. It was hypothesized that enforcing periodicity in the CFD (inherent in the frequency response technique) would overestimate the response. The results instead showed that the transient analysis results were up to 10% higher for "clean" nodal diameter excitations and six times larger for "messy" excitations, where substantial Fourier content around the main harmonic exists. Because the bulk of resonance problems are due to the "clean" excitations, a 10% underprediction is not necessarily a problem, especially since the average response in the transient is similar to the frequency response result, and so in a realistic finite life calculation, the life would be same. However, in the rare cases when the "messy" excitations harmonics are identified as the source of potential resonance concerns, this research does indicate that frequency response analysis is inadequate for accurate characterization of blade structural capability.

Brown, Andrew M.↗
NASA NTRS2003

Non-Genomic Origins of Proteins and Metabolism

It is proposed that evolution of inanimate matter to cells endowed with a nucleic acid- based coding of genetic information was preceded by an evolutionary phase, in which peptides not coded by nucleic acids were able to self-organize into networks capable of evolution towards increasing metabolic complexity. Recent findings that truly different, simple peptides (Keefe and Szostak, 2001) can perform the same function (such as ATP binding) provide experimental support for this mechanism of early protobiological evolution. The central concept underlying this mechanism is that the reproduction of cellular functions alone was sufficient for self-maintenance of protocells, and that self- replication of macromolecules was not required at this stage of evolution. The precise transfer of information between successive generations of the earliest protocells was unnecessary and, possibly, undesirable. The key requirement in the initial stage of protocellular evolution was an ability to rapidly explore a large number of protein sequences in order to discover a set of molecules capable of supporting self- maintenance and growth of protocells. Undoubtedly, the essential protocellular functions were carried out by molecules not nearly as efficient or as specific as contemporary proteins. Many, potentially unrelated sequences could have performed each of these functions at an evolutionarily acceptable level. As evolution progressed, however proteins must have performed their functions with increasing efficiency and specificity. This, in turn, put additional constraints on protein sequences and the fraction of proteins capable of performing their functions at the required level decreased. At some point, the likelihood of generating a sufficiently efficient set of proteins through a non-coded synthesis was so small that further evolution was not possible without storing information about the sequences of these proteins. Beyond this point, further evolution required coupling between proteins and informational polymers that is characteristic to all known forms of life. The emergence of such coupling must be postulated in any scenario of the origin of life, no matter whether it starts with RNA or proteins. To examine the evolutionary potential of non-genomic systems, a simple, computationally tractable model, which is still capable of capturing the essential features of the real system, has been studied computationally. Both constructive and destructive processes have been introduced into the model in a stochastic manner. Instead of assuming random reaction sets, only a suite of protobiologically plausible reactions has been considered. Peptides have been explicitly considered as protoenzymes and their catalytic efficiencies have been assigned on the basis of biochemical principles and experimental estimates. Simulations have been carried out using a novel approach (The Next Reaction Method) that is appropriate even for very low concentrations of reactants. Studies have focused on global autocatalytic processes and their diversity.

Pohorille, Andrew↗
NASA NTRS1999

An Analysis of Gravity-Field Estimation Based on Intersatellite Dual-1-Way Biased Ranging

The GRACE (Gravity Recovery And Climate Experiment) mission is designed to make global, highly accurate measurements of the Earth's gravity field with high spatial resolution. Ancillary GPS occultation measurements are also to be carried out for atmospheric monitoring. In the dual-1-way biased ranging of this mission, the range between two satellites separated by 100 to 200 km in nearly polar, coplanar, circular orbits, is measured to very high precision, to within an additive constant, through the exchange of K- and Ka-band sinusoidal signals. Such biased ranging data, along with GPS L-band range and phase data, can be processed and fit over successive multiday intervals to obtain accurate estimates of the Earth's gravity field. This report approximately models and analyzes this process, from the generation of the RF signals at the two satellites through the extraction of the geopotential. The steps include generation of the transmitted signals, processing the received signals to extract high-rate baseband phase, carrying out a dual-1-way combination of baseband phase to extract high-rate biased range for each band, combining K- and Ka-band ranges to correct for the ionosphere effect, and processing the resulting high-rate biased range values to extract three types of reduced-rate observables: biased range, range rate and range acceleration. The version of dual-1-way biased ranging developed by this report improves upon previous versions in a number of ways: highly accurate satellite-timetag corrections derived from concurrent GPS data, better baseband phase extraction using highly digital processing, highly accurate USO-rate calibration derived from concurrent GPS data, an improved method for extracting high-rate biased range from baseband phase, improved filtering for extracting reduced- rate observables from high-rate biased range, and parallel extraction of three observable types.

SINUSOIDAL SIGNALS↗
NASA NTRS1996

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗
NASA NTRS2016

Recurrence Rate and Magma Effusion Rate for the Latest Volcanism on Arsia Mons, Mars

Magmatism and volcanism have evolved the Martian lithosphere, surface, and climate throughout the history of Mars. Constraining the rates of magma generation and timing of volcanism on the surface clarifies the ways in which magma and volcanic activity have shaped these Martian systems. The ages of lava flows on other planets are often estimated using impact crater counts, assuming that the number and size-distribution of impact craters per unit area reflect the time the lava flow has been on the surface and exposed to potential impacts. Here we show that impact crater age model uncertainty is reduced by adding stratigraphic information observed at locations where neighboring lavas abut each other, and demonstrate the significance of this reduction in age uncertainty for understanding the history of a volcanic field comprising 29 vents in the 110-kilometer-diameter caldera of Arsia Mons, Mars. Each vent within this caldera produced lava flows several to tens of kilometers in length; these vents are likely among the youngest on Mars, since no impact craters in their lava flows are larger than 1 kilometer in diameter. First, we modeled the age of each vent with impact crater counts performed on their corresponding lava flows and found very large age uncertainties for the ages of individual vents, often spanning the estimated age for the entire volcanic field. The age model derived from impact crater counts alone is broad and unimodal, with estimated peak activity in the field around 130Ma (megaannum, 1 million years). Next we applied our volcano event age model (VEAM), which uses a directed graph of stratigraphic relationships and random sampling of the impact crater age determinations to create alternative age models. Monte Carlo simulation was used to create 10,000 possible vent age sets. The recurrence rate of volcanism is calculated for each possible age set, and these rates are combined to calculate the median recurrence rate of all simulations. Applying this approach to the 29 volcanic vents, volcanism likely began around 200-300Ma then first peaked around 150Ma, with an average production rate of 0.4 vents per Myr (million years). The recurrence rate estimated including stratigraphic data is distinctly bimodal, with a second, lower peak in activity around 100Ma. Volcanism then waned until the final vents were produced 10-90Ma. Based on this model, volume flux is also bimodal, reached a peak rate of 1-8 cubic kilometers per million years by 150Ma and remained above half this rate until about 90Ma, after which the volume flux diminished greatly. The onset of effusive volcanism from 200-150Ma might be due to a transition of volcanic style away from explosive volcanism that emplaced tephra on the western flank of Arsia Mons, while the waning of volcanism after the 150Ma peak might represent a larger-scale diminishing of volcanic activity at Arsia Mons related to the emplacement of flank apron lavas.

Richardson, Jacob A.↗
NASA NTRS2003

Large-Eddy Simulation Code Developed for Propulsion Applications

A large-eddy simulation (LES) code was developed at the NASA Glenn Research Center to provide more accurate and detailed computational analyses of propulsion flow fields. The accuracy of current computational fluid dynamics (CFD) methods is limited primarily by their inability to properly account for the turbulent motion present in virtually all propulsion flows. Because the efficiency and performance of a propulsion system are highly dependent on the details of this turbulent motion, it is critical for CFD to accurately model it. The LES code promises to give new CFD simulations an advantage over older methods by directly computing the large turbulent eddies, to correctly predict their effect on a propulsion system. Turbulent motion is a random, unsteady process whose behavior is difficult to predict through computer simulations. Current methods are based on Reynolds-Averaged Navier- Stokes (RANS) analyses that rely on models to represent the effect of turbulence within a flow field. The quality of the results depends on the quality of the model and its applicability to the type of flow field being studied. LES promises to be more accurate because it drastically reduces the amount of modeling necessary. It is the logical step toward improving turbulent flow predictions. In LES, the large-scale dominant turbulent motion is computed directly, leaving only the less significant small turbulent scales to be modeled. As part of the prediction, the LES method generates detailed information on the turbulence itself, providing important information for other applications, such as aeroacoustics. The LES code developed at Glenn for propulsion flow fields is being used to both analyze propulsion system components and test improved LES algorithms (subgrid-scale models, filters, and numerical schemes). The code solves the compressible Favre-filtered Navier- Stokes equations using an explicit fourth-order accurate numerical scheme, it incorporates a compressible form of Smagorinsky s model for the subgrid-scale turbulence, and it uses generalized curvilinear coordinates to allow analysis of a wide range of geometries. The code runs in parallel on shared memory multiprocessor computers and is written in Fortran 90 with dynamic memory allocation. A sample result for a Mach-1.4 round jet is presented in the figure. Instantaneous Mach number contours in several cross-planes downstream of the nozzle exit are shown, illustrating how an LES captures the large unsteady three-dimensional turbulent structures present in the jet.

DeBonis, James R.↗
NASA NTRS2019

Using Deep Learning to Automate Inference of Meteoroid Pre-Entry Properties

Properly assessing the asteroid threat depends on the knowledge of asteroid pre-entry parameters, such as size, velocity, mass, density, and strength. Although a vast number of possible bodies to study exist, such characterization of asteroid populations is currently limited by substantial costs associated with space rendezvous missions and rare meteorite findings. As asteroids fragment, ablate, and decelerate in the atmosphere, they emit light detectable by ground-based and space-borne instruments. Earth’s atmosphere, thus, becomes an accessible laboratory that enables impactor risk assessments by facilitating inference of the pre-entry parameters. These asteroid pre-entry conditions are typically deduced by modeling the entry and breakup physics that best reproduce the observed light or energy deposition curve. However, this process requires extensive manual trial-and-error of uncertain modeling parameters. Automating meteor modeling and inference would improve property distributions used in risk assessments and enable population characterization as more light curves become more readily available through the presence of space assets and ground-based camera networks. We previously developed a genetic algorithm to automate meteor modeling by using the fragment-cloud model (FCM) to search for the values of the FCM input parameters (e.g., diameter) that generate energy deposition profiles that match the observed one. Now, we apply deep learning to infer asteroid diameter, velocity, and density from observed energy deposition curves. We trained and tested our neural network models with synthetic energy deposition curves modeled using the FCM rubble pile implementation. We present an application of a 1D convolutional neural network and compare its performance to other attempted regressors and machine learning techniques, such as a fully connected neural network and Random Forest regression, to demonstrate its capabilities. We validate our model weights and approach using the Chelyabinsk, Tagish Lake, Benešov, Košice, and Lost City meteors.

Tarano, Ana Maria↗