Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stochastic modelling”

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 1,081 records · Page 60

Particle-laden swirling free jets - Measurements and predictions

A theoretical and experimental investigation of single-phase and particle-laden weakly swirling jets was conducted. The jets were injected vertically downward from a 19 mm diameter tube with swirl numbers ranging from 0 to 0.33. The particle-laden jets had a single loading ratio (0.2) with particles having an SMD of 39 micrometers. Mean and fluctuating properties of both phases were measured using nonintrusive laser based methods while particle mass flux was measured using an isokinetic sampling probe. The continuous phase was analyzed using both a baseline k-epsilon turbulence model and an extended version with modifications based on the flux Richardson number to account for effects of streamline curvature. To highlight effects of interphase transport rates and particle/turbulence interactions, effects of the particles were analyzed as follows: (1) locally homogeneous flow (LHF) analysis, where interphase transport rates are assumed to be infinitely fast; (2) deterministic separated flow (DSF) analysis, where finite interphase transport rates are considered but particle/turbulence interactions are ignored; and (3) stochastic separated flow (SSF) analysis, where both effects are considered using random-walk computations.

Bulzan, D. L.↗

Precipitating electron energy flux and auroral zone conductances - An empirical model

Data from the low energy electron (LEE) experiments on the Atmosphere Explorer C and D satellites have been used to determine the average global distribution of the energy flux of precipitating auroral electrons and their average energy for different levels of geomagnetic activity. Measurements from the Atmosphere Explorer unified abstract file (15-s resolution) have been binned according to invariant latitude (in the range 50-90 deg), magnetic local time, and geomagnetic activity as measured by the Kp and auroral electrojet (AE) indices, separately. Bin-averaged values of precipitating energy flux and average energy have been calculated, and a smoothing algorithm used to reduce stochastic variations in the raw data. The results indicate that, for the parameters studied, the AE inces does a superior job of ordering the data with regard to geomagnetic activity. The global distribution of the auroral enhancement porition of the Pedersen and Hall conductances were inferred from the data by means of an empirical fit to detailed energy deposition calculations.

Spiro, R. W.↗

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

Design and evaluation of a Stochastic Optimal Feed-forward and Feedback Technology (SOFFT) flight control architecture

This paper describes the design and evaluation of a stochastic optimal feed-forward and feedback technology (SOFFT) control architecture with emphasis on the feed-forward controller design. The SOFFT approach allows the designer to independently design the feed-forward and feedback controllers to meet separate objectives and then integrate the two controllers. The feed-forward controller has been integrated with an existing high-angle-of-attack (high-alpha) feedback controller. The feed-forward controller includes a variable command model with parameters selected to satisfy level 1 flying qualities with a high-alpha adjustment to achieve desired agility guidelines, a nonlinear interpolation approach that scales entire matrices for approximation of the plant model, and equations for calculating feed-forward gains developed for perfect plant-model tracking. The SOFFT design was applied to a nonlinear batch simulation model of an F/A-18 aircraft modified for thrust vectoring. Simulation results show that agility guidelines are met and that the SOFFT controller filters undesired pilot-induced frequencies more effectively during a tracking task than a flight controller that has the same feedback control law but does not have the SOFFT feed-forward control.

Ostroff, Aaron J.↗

Towards a Probabilistic Criterion for Preliminary Shell Design

Thin-walled stiffened or unstiffened, metallic or composite shells are widely used structural elements in aeronautical and space applications. Buckling strength, which is a major concern in all these areas, is affected by the uncertainties in the definition of loads, material properties, geometric variables, engineering models and the accuracy of the analysis tools used in the design phase. The NASA design criteria monographs from the late sixties account for these design uncertainties by the use of a lump sum safety factor or the so-called "knockdown" factor gamma, which usually results in an overly conservative design. In the present paper a new, reliability based, probabilistic design procedure for buckling critical imperfect isotropic shells is proposed. It essentially consists of a stochastic approach based on a new, improved "knockdown" factor lambda(sub a) that is not as conservative as the traditional one. It is felt that by quantifying and understanding the "problem uncertainties" such as initial imperfections and their influence on the design variables, one can develop a better engineered, better designed and safer system.

Arbocz, Johann↗

The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates

Cumulus clouds can become tilted or elongated in the presence of wind shear. Nevertheless, most studies of the interaction of cumulus clouds and radiation have assumed these clouds to be isotropic. This paper describes an investigation of the effect of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and atmospheric heating rate profiles. A stochastic field generation algorithm was used to produce twenty three-dimensional liquid water content fields based on the statistical properties of cloud scenes from a large eddy simulation. Progressively greater degrees of x-z plane tilting and horizontal stretching were imposed on each of these scenes, so that an ensemble of scenes was produced for each level of distortion. The resulting scenes were used as input to a three-dimensional Monte Carlo radiative transfer model. Domain-average transmission, reflection, and absorption of broadband solar radiation were computed for each scene along with the average heating rate profile. Both tilt and horizontal stretching were found to significantly affect calculated fluxes, with the amount and sign of flux differences depending strongly on sun position relative to cloud distortion geometry. The mechanisms by which anisotropy interacts with solar fluxes were investigated by comparisons to independent pixel approximation and tilted independent pixel approximation computations for the same scenes. Cumulus anisotropy was found to most strongly impact solar radiative transfer by changing the effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface perpendicular to the direction of the incident solar beam.

Hinkelman, Laura M.↗

The NANOGrav Nine-Year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr(exp -1) for a GWB from supermassive black hole binaries of A(sub gw) less than 1.5 x 10(exp -15). We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, OMEGA(sub gw) (f) h squared less than 4.2 x 10(exp -10). Our limit on the cosmic string GWB, OMEGA(sub gw) (f) h squared less than 2.2 x 10(exp -10), translates to a conservative limit on the cosmic string tension with G mu less than 3.3 x 10(exp -8), a factor of four better than the joint Planck and high-l‚ cosmic microwave background data from other experiments.

Arzoumanian, Z.↗

Modulating terrestrial impacts from Oort cloud comets by the adiabatically changing galactic tides

Time modulation of the flux of new Jupiter-dominated Oort cloud comets is the subject of interest here. The major perturbation of these comets during the present epoch is due to the tidal field of the relatively smooth distribution of matter in the galactic disk. A secondary source of the near-parabolic comet flux are stars penetrating the inner Oort cloud and providing impulses that create brief comet showers. Substantial stellar-induced showers occur approximately every 100 m.y. Less frequent (but stronger) impulses due to giant molecular clouds can also perturb comets from the inner cloud. These occur on timescales of approximately equal to 500 m.y. In contrast to these infrequent stochastic shower phenomena is the continuously varying tidal-induced flux due to the galaxy. As the Sun orbits the galactic center it undergoes quasiharmonic motion about the galactic midplane, which is superimposed on the small eccentricity, near-Keplerian motion in the plane having epicycle period approximately equal to 150 m.y. In the process the galactic tidal field on the Sun/cloud system will vary causing a modulation of the observable Oort cloud flux. We have created a model of the galactic matter distribution as it affects the solar motion over a time interval ranging from 300 m.y. in the past to 100 m.y. into the future. As constraints on the disk's compact dark matter component we require consistency with the following: (1) the observed galactic rotation curve, (2) today's flux distribution of new comets, (3) the studies of K-giant distributions, and (4) the periodicity found in the terrestrial cratering record. The adiabatically varying galactic tidal torque is then determined and used to predict the time dependence of the flux. We find that a model in which approximately half the disk matter is compact is consistent with these constraints. Under such circumstances the peak-to-trough flux variation will be approx. equal to 5:1 with a full width of 9 m.y. This variability will manifest in the terrestrial cratering record and is consistent with the observed cratering periodicity, if over half of the impacts on Earth are caused by comets or asteroids that originate in the outer Oort cloud.

Matese, J. J.↗

Determination of Nonlinear Stiffness Coefficients for Finite Element Models with Application to the Random Vibration Problem

In this paper, a method for obtaining nonlinear stiffness coefficients in modal coordinates for geometrically nonlinear finite-element models is developed. The method requires application of a finite-element program with a geometrically non- linear static capability. The MSC/NASTRAN code is employed for this purpose. The equations of motion of a MDOF system are formulated in modal coordinates. A set of linear eigenvectors is used to approximate the solution of the nonlinear problem. The random vibration problem of the MDOF nonlinear system is then considered. The solutions obtained by application of two different versions of a stochastic linearization technique are compared with linear and exact (analytical) solutions in terms of root-mean-square (RMS) displacements and strains for a beam structure.

Muravyov, Alexander A.↗

Probabilistic evaluation of SSME structural components

The application is described of Composite Load Spectra (CLS) and Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) family of computer codes to the probabilistic structural analysis of four Space Shuttle Main Engine (SSME) space propulsion system components. These components are subjected to environments that are influenced by many random variables. The applications consider a wide breadth of uncertainties encountered in practice, while simultaneously covering a wide area of structural mechanics. This has been done consistent with the primary design requirement for each component. The probabilistic application studies are discussed using finite element models that have been typically used in the past in deterministic analysis studies.

Rajagopal, K. R.↗

The Integrated Medical Model: A Risk Assessment and Decision Support Tool for Space Flight Medical Systems

The Integrated Medical Model (IMM) is a decision support tool that is useful to mission planners and medical system designers in assessing risks and designing medical systems for space flight missions. The IMM provides an evidence based approach for optimizing medical resources and minimizing risks within space flight operational constraints. The mathematical relationships among mission and crew profiles, medical condition incidence data, in-flight medical resources, potential crew functional impairments, and clinical end-states are established to determine probable mission outcomes. Stochastic computational methods are used to forecast probability distributions of crew health and medical resource utilization, as well as estimates of medical evacuation and loss of crew life. The IMM has been used in support of the International Space Station (ISS) medical kit redesign, the medical component of the ISS Probabilistic Risk Assessment, and the development of the Constellation Medical Conditions List. The IMM also will be used to refine medical requirements for the Constellation program. The IMM outputs for ISS and Constellation design reference missions will be presented to demonstrate the potential of the IMM in assessing risks, planning missions, and designing medical systems. The implementation of the IMM verification and validation plan will be reviewed. Additional planned capabilities of the IMM, including optimization techniques and the inclusion of a mission timeline, will be discussed. Given the space flight constraints of mass, volume, and crew medical training, the IMM is a valuable risk assessment and decision support tool for medical system design and mission planning.

Kerstman, Eric↗

Probabilistic Path Planning of Montgolfier Balloons in Strong, Uncertain Wind Fields

Lighter-than-air vehicles such as hot-air balloons have been proposed for exploring Saturn s moon Titan, as well as other bodies with significant atmospheres. For these vehicles to navigate effectively, it is critical to incorporate the effects of surrounding wind fields, especially as these winds will likely be strong relative to the control authority of the vehicle. Predictive models of these wind fields are available, and previous research has considered problems of planning paths subject to these predicted forces. However, such previous work has considered the wind fields as known a priori, whereas in practical applications, the actual wind vector field is not known exactly and may deviate significantly from the wind velocities estimated by the model. A probabilistic 3D path-planning algorithm was developed for balloons to use uncertain wind models to generate time-efficient paths. The nominal goal of the algorithm is to determine what altitude and what horizontal actuation, if any is available on the vehicle, to use to reach a particular goal location in the least expected time, utilizing advantageous winds. The solution also enables one to quickly evaluate the expected time-to-goal from any other location and to avoid regions of large uncertainty. This method is designed for balloons in wind fields but may be generalized for any buoyant vehicle operating in a vector field. To prepare the planning problem, the uncertainty in the wind field is modeled. Then, the problem of reaching a particular goal location is formulated as a Markov decision process (MDP) using a discretized space approach. Solving the MDP provides a policy of what actuation option (how much buoyancy change and, if applicable, horizontal actuation) should be selected at any given location to minimize the expected time-to-goal. The results provide expected time-to-goal values from any given location on the globe in addition to the action policy. This stochastic approach can also provide insights not accessible by deterministic methods; for example, one can evaluate variability and risk associated with different scenarios, rather than only viewing the expected outcome.

Wolf, Michael↗

Flare gamma ray continuum emission from neutral pion decay

We investigate, in detail, the production of solar flare gamma ray emission above 100 MeV via the interaction of high energy protons with the ambient solar atmosphere. We restrict our considerations to the broadband gamma ray spectrum resulting from the decay of neutral pions produced in p-H reactions. Thick-target calculations are performed to determine the photon fluences. However, proton transport is not considered. Inferences about the form of the proton spectrum at 10-100 MeV have already been drawn from de-excitation gamma ray lines. Our aim is to constrain the proton spectrum at higher energies. Thus, the injected proton spectrum is assumed to have the form of a Bessel Function, characteristics of stochastic energy at higher energies. The detailed shape of the gamma ray spectra around 100 MeV is found to have a strong dependence on the spectral index of the power law and on the turnover energy (from Bessel function to power law). As would be expected, the harder the photon spectrum the wider the 100 MeV feature. The photon spectra are to be compared with observations and used to place limits upon the number of particles accelerated and to constrain acceleration models.

Alexander, David↗

An adaptive tracking observer for failure-detection systems

The design problem of adaptive observers applied to linear, constant and variable parameters, multi-input, multi-output systems, is considered. It is shown that, in order to keep the observer's (or Kalman filter) false-alarm rate (FAR) under a certain specified value, it is necessary to have an acceptable proper matching between the observer (or KF) model and the system parameters. An adaptive observer algorithm is introduced in order to maintain desired system-observer model matching, despite initial mismatching and/or system parameter variations. Only a properly designed adaptive observer is able to detect abrupt changes in the system (actuator, sensor failures, etc.) with adequate reliability and FAR. Conditions for convergence for the adaptive process were obtained, leading to a simple adaptive law (algorithm) with the possibility of an a priori choice of fixed adaptive gains. Simulation results show good tracking performance with small observer output errors and accurate and fast parameter identification, in both deterministic and stochastic cases.

Sidar, M.↗

Navigation between the planets

Recent advances in spacecraft tracking, chronometry, ephemerides, and orbit and trajectory determinations are reviewed. Improvements in timekeeping are reviewed, as well as precision distance and range measurements; orbit determinations, trajectory-correction maneuvers, flight path optimization, and information provided by rotation of the tracking station with the earth's surface. Doppler and tropospheric wave propagation effects are dealt with. Nongravitational perturbations (solar radiation pressure, release of gases from the spacecraft, stochastic unmodeled accelerations and sequential estimation to cope with them), the effect of the target planet's gravitational field upon close approach, and navigation problems in the outer reaches of the solar system (TV data telemetered back for inertial navigation) are covered. By-products of the research include: refined data on the mass of planets, on planetary mass distributions, planet configurations, on physical properties of the atmospheres and ionospheres of planets, and opportunities for refined tests of gravitation and relativity theories and models.

Melbourne, W. G.↗

Chaotic appearance of the AE index

Results are reported from a stochastic analysis of a 5-day time series of the geomagnetic AE index during an active period. The original data, the power spectrum, and the autocorrelation function are shown, and the steps in the analysis are described in detail. It is found that the autocorrelation time scale is about 50 min, giving a correlation dimension (for the construction of a time series of m-dimensional vectors) of 2.4. This result is consistent with either colored-noise or deterministic-chaos magnetosphere models, indicating the need for further investigation.

Shan, L.-H.↗

RITRACKS: A Software for Simulation of Stochastic Radiation Track Structure, Micro and Nanodosimetry, Radiation Chemistry and DNA Damage for Heavy Ions

The code RITRACKS (Relativistic Ion Tracks) has been developed over the last few years at the NASA Johnson Space Center to simulate the effects of ionizing radiations at the microscopic scale, to understand the effects of space radiation at the biological level. The fundamental part of this code is the stochastic simulation of radiation track structure of heavy ions, an important component of space radiations. The code can calculate many relevant quantities such as the radial dose, voxel dose, and may also be used to calculate the dose in spherical and cylindrical targets of various sizes. Recently, we have incorporated DNA structure and damage simulations at the molecular scale in RITRACKS. The direct effect of radiations is simulated by introducing a slight modification of the existing particle transport algorithms, using the Binary‐Encounter‐Bethe model of ionization cross sections for each molecular orbitals of DNA. The simulation of radiation chemistry is done by a step‐by‐step diffusion‐reaction program based on the Green's functions of the diffusion equation]. This approach is also used to simulate the indirect effect of ionizing radiation on DNA. The software can be installed independently on PC and tablets using the Windows operating system and does not require any coding from the user. It includes a Graphic User Interface (GUI) and a 3D OpenGL visualization interface. The calculations are executed simultaneously (in parallel) on multiple CPUs. The main features of the software will be presented.

Plante, I↗

Fundamental and harmonic radiation in type III solar radio bursts

Type III solar radio bursts are investigated by modeling the propagation of the electron beam and the generation and subsequent propagation of waves to the observer. Predictions from this model are compared in detail with particle, Langmuir wave, and radio data from the International Sun Earth Explorer-3 (ISSE-3) spacecraft and with other observations to clarify the roles of fundamental and harmonic emission in type III radio bursts. Langmuir waves are seen only after the arrival of the beam, in accord with the standard theory. These waves persist after a positive beam slope is last resolved, implying that sporadic positive slopes persist for some time, unresolved but in accord with the predictions of stochastic growth theory. Local electromagnetic emission sets in only after Langmuir waves are seen, in accord with the standard theory, which relies on nonlinear processes involving Langmuir waves. In the events investigated here, fundamental radiation appears to dominate early in the event, followed and/or accompanied by harmonic radiation after the peak, with a long-lived tail of multiply scattered fundamental or harmonic emission extending long afterwards. These results are largely independent of, but generally consistent with, the conclusions of earlier works.

Robinson, P. A.↗