Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Randomized methods”

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 703 records · Page 39

Test particle calculations of pick-up ions in the vicinity of Comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked-up by the solar wind. The cometary ion pick-up process for Comet Giacobini-Zinner is studied using a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density. The cometary ion distribution function was calculated at several locations upstream of the bow shock. Several types of magnetic fluctuations were included. Circularly and linearly polarized monochromatic waves result in less pitch-angle scattering of pick-up ions than turbulent fluctuations which have about the same amplitude. Some anisotropy remains in the distribution function just upstream of the shock, even for the turbulent fluctuations.

Cravens, T. E.↗

Theoretical calculations of ion acceleration in the vicinity of comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked up by the solar wind. The cometary ion pickup process for comet Giacobini-Zinner is studied in two ways: (1) with a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.

Mckenzie, M. L.↗

Soil Sampling Techniques For Alabama Grain Fields

Characterizing the spatial variability of nutrients facilitates precision soil sampling. Questions exist regarding the best technique for directed soil sampling based on a priori knowledge of soil and crop patterns. The objective of this study was to evaluate zone delineation techniques for Alabama grain fields to determine which method best minimized the soil test variability. Site one (25.8 ha) and site three (20.0 ha) were located in the Tennessee Valley region, and site two (24.2 ha) was located in the Coastal Plain region of Alabama. Tennessee Valley soils ranged from well drained Rhodic and Typic Paleudults to somewhat poorly drained Aquic Paleudults and Fluventic Dystrudepts. Coastal Plain s o i l s ranged from coarse-loamy Rhodic Kandiudults to loamy Arenic Kandiudults. Soils were sampled by grid soil sampling methods (grid sizes of 0.40 ha and 1 ha) consisting of: 1) twenty composited cores collected randomly throughout each grid (grid-cell sampling) and, 2) six composited cores collected randomly from a -3x3 m area at the center of each grid (grid-point sampling). Zones were established from 1) an Order 1 Soil Survey, 2) corn (Zea mays L.) yield maps, and 3) airborne remote sensing images. All soil properties were moderately to strongly spatially dependent as per semivariogram analyses. Differences in grid-point and grid-cell soil test values suggested grid-point sampling does not accurately represent grid values. Zones created by soil survey, yield data, and remote sensing images displayed lower coefficient of variations (8CV) for soil test values than overall field values, suggesting these techniques group soil test variability. However, few differences were observed between the three zone delineation techniques. Results suggest directed sampling using zone delineation techniques outlined in this paper would result in more efficient soil sampling for these Alabama grain fields.

Thompson, A. N.↗

Estimation of Sonic Fatigue by Reduced-Order Finite Element Based Analyses

A computationally efficient, reduced-order method is presented for prediction of sonic fatigue of structures exhibiting geometrically nonlinear response. A procedure to determine the nonlinear modal stiffness using commercial finite element codes allows the coupled nonlinear equations of motion in physical degrees of freedom to be transformed to a smaller coupled system of equations in modal coordinates. The nonlinear modal system is first solved using a computationally light equivalent linearization solution to determine if the structure responds to the applied loading in a nonlinear fashion. If so, a higher fidelity numerical simulation in modal coordinates is undertaken to more accurately determine the nonlinear response. Comparisons of displacement and stress response obtained from the reduced-order analyses are made with results obtained from numerical simulation in physical degrees-of-freedom. Fatigue life predictions from nonlinear modal and physical simulations are made using the rainflow cycle counting method in a linear cumulative damage analysis. Results computed for a simple beam structure under a random acoustic loading demonstrate the effectiveness of the approach and compare favorably with results obtained from the solution in physical degrees-of-freedom.

Rizzi, Stephen A.↗

VAPEPS Management Center

A VibroAcoustic Payload Environment Prediction System (VAPEPS) Management Center is established and being maintained at JPL. VAPEPS is a computer program that utilizes statistical energy analysis (SEA) theory and SEA extrapolation methods to predict the vibroacoustic environments of complex systems. The program is commonly utilized to establish random vibration and acoustic test requirements for aerospace components exposed to launch environments. It can also be used in the payload or launch vehicle design process for control or mitigation of the launch or on-orbit vibration and acoustic environments, as well as to calculate acoustically induced stresses in structures such as solar panels and reflectors.

VAPEPS↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Carbon Nanotube – Polysiloxane Nanocomposites for Ablative Thermal Protection Systems

Ultra-high temperature resins (UHTR) such as polysiloxane possess excellent characteristics for ablation including low thermal conductivity, lightweight composition, high thermal resistance, and the ability to form a durable char layer. To enhance its mechanical strength, the aerospace industry utilizes fiber reinforcement, such as carbon fiber reinforced polysiloxane (CF-polysiloxane), due to its high strength-to-weight ratio. Carbon nanotubes (CNTs) have advantages over microfibers including higher specific strength, lower density, and customizable thermal conductivity. Their nano-scale dimensions allow for a more homogenous material compared to traditional microfibers. Conventional methods of incorporating CNTs into resins often have structural defects such as voids, CNT agglomerations, and randomly dispersed CNTs, resulting in marginal improvements in properties. A novel approach known as bulk nanocomposite laminating (BNL) has been introduced involving CNT densification, resin polymer infiltration, and pressurized curing to achieve a high packing density exceeding 35% by volume. Characterization and testing to simulate re-entry conditions is performed with a sample of CNT-polysiloxane nanocomposite cured onto a CF-polysiloxane backing. These samples are subjected to ablation by an oxyacetylene flame while collecting temperature data to evaluate thermal protection properties.

Palak B Patel↗

Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. I: Approach Based on the Far-Field Foldy Equations

In this paper, we revisit, with further enhancements and clarifications, the self-consistent first-principles approach developed previously for deriving the vector radiative transfer theory for a discrete random medium with a sparse concentration of particles. We specifically consider the case of a plane-parallel particulate layer embedded in an otherwise homogeneous unbounded medium. The solution method is based on the far-field Foldy equations, an order-of-scattering expansion for the total field derived under the Twersky approximation, the computation of the coherent field by assuming that the positions of the particles are uncorrelated, and the ladder approximation for the coherency dyadic. The latter yields an integral equation for the diffuse specific coherency dyadic, defined through an angular spectrum repre- sentation for the coherency dyadic, which in turn, gives the vector radiative transfer equation for the diffuse specific intensity column vector. We analyze specifically the computation of the coherent field for inhomogeneous particulate media and multiple species of particles, the continuous extension of the far-field representation to the near field, the Foldy approximation, and the Foldy integral equation for the coherent field. Finally, we discuss the transition from the vector to the scalar radiative transfer equation.

Doicu, Adrian↗

Optical devices and methods employing nanoparticles, microcavities, and semicontinuous metal films

An optical sensing enhancing material (and corresponding method of making) comprising: a medium, the medium comprising a plurality of aggregated nanoparticles comprising fractals; and a microcavity, wherein the medium is located in a vicinity of the microcavity. Also an optical sensor and sensing method comprising: providing a doped medium, the medium comprising a plurality of aggregated nanoparticles comprising fractals, with the material; locating the doped medium in the vicinity of a microcavity; exciting the doped medium with a light source; and detecting light reflected from the doped medium. Also an optical sensing enhancing material comprising a medium, the medium comprising a semicontinuous metal film of randomly distributed metal particles and their clusters at approximately their percolation threshold. The medium preferably additionally comprises a microcavity/microresonator. Also devices and methods employing such material.

Armstrong, Robert L.↗

Validation of a Deterministic Vibroacoustic Response Prediction Model

This report documents the recently completed effort involving validation of a deterministic theory for the random vibration problem of predicting the response of launch pad structures in the low-frequency range (0 to 50 hertz). Use of the Statistical Energy Analysis (SEA) methods is not suitable in this range. Measurements of launch-induced acoustic loads and subsequent structural response were made on a cantilever beam structure placed in close proximity (200 feet) to the launch pad. Innovative ways of characterizing random, nonstationary, non-Gaussian acoustics are used for the development of a structure's excitation model. Extremely good correlation was obtained between analytically computed responses and those measured on the cantilever beam. Additional tests are recommended to bound the problem to account for variations in launch trajectory and inclination.

Caimi, Raoul E.↗

Layered Electrode Architectures Enabled by Holey Graphene Hosts

Battery electrodes are typically prepared using a slurry-based method, where high boiling point organic solvents are typically used with a polymer binder to disperse active electrode materials and conductive carbon additives followed by casting into current collector sheets. In this process, the use of organic solvents are environmentally hazardous, energy intensive, and time-consuming, while polymeric binders add weight and can undergo parasitic reactions during battery cycling. Here we discuss the use of a unique conductive carbon nanomaterial, namely holey graphene, as a multifunctional host that can be processed via a facile, solvent-free, binder-free, cold pressing method to prepare electrodes with unique architectures. This approach can accommodate a broad range of battery chemistries, such as lithium oxygen/carbon dioxide (Li-O2/CO2), lithium sulfur/selenium (Li-S/Se), and various lithium ion chemistries. The successful use of dry, cold compression allows the convenient modification of electrode mass loading with unique architectures containing active materials randomly mixed or layered with single or multiple layers of one or more active materials. In the conventional slurry-based electrode fabrication method, it is difficult, if not impossible, to achieve such layered architectures. In particular, the layered electrode architectures from the dry, cold press approach are easy to fabricate without significantly sacrificing device performance, and are amenable to continuous roll-to-roll dry processing. In addition, these architectures offer unprecedented perspectives on the battery material loading vs. utilization efficiency, as well as enabling unique characterization opportunities to advance mechanistic understanding of various battery chemistries. Potential applications of this unique electrode fabrication approach for solid-state batteries will also be discussed.

holey graphene, batteries, electrode architecture,↗

Effects of Dose Error and Sample Size on Sonic Boom Dose-Response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

Meteorological Measurements from Satellite Platforms

The stability of stochastic satellites and the stability and control of flexible satellites were investigated. The effects of random environmental torques and noises in the moments of inertia of spinning and three-axes stabilized satellites were first compared analytically by four methods and by analog simulations. Among the analytical methods, it was shown that the Fokker-Planck formulation yields predictions which most coincide with the simulation results. It was then shown that the required stability criterion of a satellite is quite different from that obtained by a deterministic approach, under the assumption that the environmental and control torques experienced by the satellite are random. Finally, it was demonstrated that, by monitoring the deformations of the flexible elements of a satellite, the effectiveness of the satellite control system can be increased considerably.

Suomi, V. E.↗

Assessment of Bulk Absorber Properties for Multi-Layer Perforates in Porous Honeycomb Liners

CONTINUING progress in materials technology provides potential for improved acoustic liners for attenuating broadband fan noise emissions from aircraft engine nacelles. Conventional liners (local-reacting perforate-over-honeycomb structures) provide significant narrow-band attenuation, but limited attenuation over wide bandwidths. Two approaches for increasing attenuation bandwidth are to (1) replace the honeycomb structure with bulk material, or (2) cascade multiple layers of perforate/honeycomb structures. Usage of the first approach is limited because of mechanical and maintenance reasons, while multi-layer liners are limited to about three layers because of their additional mechanical complexity, depth and weight. The current research concerns a novel approach reported by the University of Cincinnati, in which a single-layer conventional liner is converted into an extended-reaction, broadband absorber by making the honeycomb core structure porous. This modified single-layer liner requires no increase in depth and weight, and minimal increase in mechanical complexity. Langley has initiated research to identify potential benefits of liner structures with porous cell walls. This research has two complementary goals: (1) develop and validate experimental techniques for treating multi-layer perforates (representative of the internal cells of a liner with porous cell walls) as 1-D bulk materials, and (2) develop analytical approaches to validate this bulk material assumption. If successful, the resultant model can then be used to design optimized porous honeycomb liners. The feasibility of treating an N-layer perforate system (N porous plates separated by uniform air gaps) as a one-dimensional bulk absorber is assessed using the Two-Thickness Method (TTM), which is commonly used to educe bulk material intrinsic acoustic parameters. Tests are conducted with discrete tone and random noise sources, over an SPL range sufficient to determine the nonlinearity of the test specimens, for test specimens consisting of 5, 10 and 15% porous plates. Measured impedances for two liner thicknesses (e.g., 12 and 24 layers) are used as input to the TTM to determine the characteristic impedance and propagation constant that characterize these liners as bulk absorbers. These parameters are then used to calculate the predicted impedance of liners with different thicknesses (e.g., 36 layers), and a comparison of predicted and measured impedances for these other thicknesses is used to determine the efficacy of this approach. Finally, an independent method is used to educe the propagation constant for a single representative sample, and excellent comparison between the results for this method and those for the TTM provides increased confidence in the results achieved with the TTM. In general, the results demonstrate these multi-layer perforates can be acceptably treated as bulk absorbers.

Jones, Michael G.↗

Extending unbiased stereology of brain ultrastructure to three-dimensional volumes

OBJECTIVE: Analysis of brain ultrastructure is needed to reveal how neurons communicate with one another via synapses and how disease processes alter this communication. In the past, such analyses have usually been based on single or paired sections obtained by electron microscopy. Reconstruction from multiple serial sections provides a much needed, richer representation of the three-dimensional organization of the brain. This paper introduces a new reconstruction system and new methods for analyzing in three dimensions the location and ultrastructure of neuronal components, such as synapses, which are distributed non-randomly throughout the brain. DESIGN AND MEASUREMENTS: Volumes are reconstructed by defining transformations that align the entire area of adjacent sections. Whole-field alignment requires rotation, translation, skew, scaling, and second-order nonlinear deformations. Such transformations are implemented by a linear combination of bivariate polynomials. Computer software for generating transformations based on user input is described. Stereological techniques for assessing structural distributions in reconstructed volumes are the unbiased bricking, disector, unbiased ratio, and per-length counting techniques. A new general method, the fractional counter, is also described. This unbiased technique relies on the counting of fractions of objects contained in a test volume. A volume of brain tissue from stratum radiatum of hippocampal area CA1 is reconstructed and analyzed for synaptic density to demonstrate and compare the techniques. RESULTS AND CONCLUSIONS: Reconstruction makes practicable volume-oriented analysis of ultrastructure using such techniques as the unbiased bricking and fractional counter methods. These analysis methods are less sensitive to the section-to-section variations in counts and section thickness, factors that contribute to the inaccuracy of other stereological methods. In addition, volume reconstruction facilitates visualization and modeling of structures and analysis of three-dimensional relationships such as synaptic connectivity.

NASA Discipline Neuroscience↗

A Monte Carlo study of Weibull reliability analysis for space shuttle main engine components

The incorporation of a number of additional capabilities into an existing Weibull analysis computer program and the results of Monte Carlo computer simulation study to evaluate the usefulness of the Weibull methods using samples with a very small number of failures and extensive censoring are discussed. Since the censoring mechanism inherent in the Space Shuttle Main Engine (SSME) data is hard to analyze, it was decided to use a random censoring model, generating censoring times from a uniform probability distribution. Some of the statistical techniques and computer programs that are used in the SSME Weibull analysis are described. The methods documented in were supplemented by adding computer calculations of approximate (using iteractive methods) confidence intervals for several parameters of interest. These calculations are based on a likelihood ratio statistic which is asymptotically a chisquared statistic with one degree of freedom. The assumptions built into the computer simulations are described. The simulation program and the techniques used in it are described there also. Simulation results are tabulated for various combinations of Weibull shape parameters and the numbers of failures in the samples.

Abernethy, K.↗

System and method for character recognition

A character recognition system is disclosed in which each character in a retina, defining a scanning raster, is scanned with random lines uniformly distributed over the retina. For each type of character to be recognized the system stores a probability density function (PDF) of the random line intersection lengths and/or a PDF of the random line number of intersections. As an unknown character is scanned, the random line intersection lengths and/or the random line number of intersections are accumulated and based on a comparison with the prestored PDFs a classification of the unknown character is performed.

Hong, J. P.↗

NASTRAN computer system level 12.1

Program uses finite element displacement method for solving linear response of large, three-dimensional structures subject to static, dynamic, thermal, and random loadings. Program adapts to computers of different manufacture, permits up-dating and extention, allows interchange of output and input information between users, and is extensively documented.

Butler, T. G.↗