Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational statistical mechanics”

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 73 records · Page 4

Reliability analysis of composite structures

A probabilistic static stress analysis methodology has been developed to estimate the reliability of a composite structure. Closed form stress analysis methods are the primary analytical tools used in this methodology. These structural mechanics methods are used to identify independent variables whose variations significantly affect the performance of the structure. Once these variables are identified, scatter in their values is evaluated and statistically characterized. The scatter in applied loads and the structural parameters are then fitted to appropriate probabilistic distribution functions. Numerical integration techniques are applied to compute the structural reliability. The predicted reliability accounts for scatter due to variability in material strength, applied load, fabrication and assembly processes. The influence of structural geometry and mode of failure are also considerations in the evaluation. Example problems are given to illustrate various levels of analytical complexity.

Kan, Han-Pin↗

Simulation Analysis of Computer-Controlled pressurization for Mixture Ratio Control

A procedural code (C++) simulation was developed to investigate potentials for mixture ratio control of pressure-fed spacecraft rocket propulsion systems by measuring propellant flows, tank liquid quantities, or both, and using feedback from these measurements to adjust propellant tank pressures to set the correct operating mixture ratio for minimum propellant residuals. The pressurization system eliminated mechanical regulators in favor of a computer-controlled, servo- driven throttling valve. We found that a quasi-steady state simulation (pressure and flow transients in the pressurization systems resulting from changes in flow control valve position are ignored) is adequate for this purpose. Monte-Carlo methods are used to obtain simulated statistics on propellant depletion. Mixture ratio control algorithms based on proportional-integral-differential (PID) controller methods were developed. These algorithms actually set target tank pressures; the tank pressures are controlled by another PID controller. Simulation indicates this approach can provide reductions in residual propellants.

Alexander, Leslie A.↗

Finite element methods in probabilistic mechanics

Probabilistic methods, synthesizing the power of finite element methods with second-order perturbation techniques, are formulated for linear and nonlinear problems. Random material, geometric properties and loads can be incorporated in these methods, in terms of their fundamental statistics. By construction, these methods are applicable when the scale of randomness is not too large and when the probabilistic density functions have decaying tails. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. Applications showing the effects of combined random fields and cyclic loading/stress reversal are studied and compared with Monte Carlo simulation results.

Liu, Wing Kam↗

At least some errors are randomly generated (Freud was wrong)

An experiment was carried out to expose something about human error generating mechanisms. In the context of the experiment, an error was made when a subject pressed the wrong key on a computer keyboard or pressed no key at all in the time allotted. These might be considered, respectively, errors of substitution and errors of omission. Each of seven subjects saw a sequence of three digital numbers, made an easily learned binary judgement about each, and was to press the appropriate one of two keys. Each session consisted of 1,000 presentations of randomly permuted, fixed numbers broken into 10 blocks of 100. One of two keys should have been pressed within one second of the onset of each stimulus. These data were subjected to statistical analyses in order to probe the nature of the error generating mechanisms. Goodness of fit tests for a Poisson distribution for the number of errors per 50 trial interval and for an exponential distribution of the length of the intervals between errors were carried out. There is evidence for an endogenous mechanism that may best be described as a random error generator. Furthermore, an item analysis of the number of errors produced per stimulus suggests the existence of a second mechanism operating on task driven factors producing exogenous errors. Some errors, at least, are the result of constant probability generating mechanisms with error rate idiosyncratically determined for each subject.

Sellen, A. J.↗

The use of statistical characteristics of reducer vibrations as diagnostic symptoms

The results of a statistical analysis of the vibrations of the experimental RS-1 reducer stand, with a spiral-gear transmission, operating on a closed circuit, are presented. The analysis was carried out on the Minsk-2 and Minsk-32 digital computers, with two-channel analog-digital converter, built in the Institute of the Science of Mechanics. Two-dimensional distribution patterns, conditional dispersions and dispersion ratios were calculated. The octave-band-filtered first harmonics of the tooth frequency f sub z of the vibrations at two different measurement points were considered as the components of the vibration process to be analyzed. The regression lines, corresponding to different values of the loading torque, are presented. Since it was not the gear drive parameters which were determined by diagnostic methods, but the characteristics most sensitive to change in state of the object of the investigation, the loading torque, which is the simplest and most accessible for measurement, was chosen as the condition parameter.

Balitskiy, F. Y.↗

Empirical Noise Modeling of Internally Mixed Exhaust Systems

There appear to be no non-proprietary methods to predict the noise of internally mixed exhaust systems, and no guidance for how to adapt known jet noise models for these configurations. This paper surveys literature and historical databases acquired at NASA Glenn's Aero-Acoustic Propulsion Lab to give such guidance. The core premise is that an exhaust system with well-designed mixer produces noise that is to first approximation the same as a fully-mixed jet flow. Refinement of what is meant by a “fully mixed jet” can lead to more accurate prediction of the main jet noise. Additional noise is often generated within the nozzle, typically at high frequencies, whose source mechanism(s) are not obvious. However, a noise prediction method can be established that captures some aspects of the excess noise and provide an estimate of the total jet noise. Explorations of source mechanisms associated with the internal mixer have led to a new noise model which includes the effects of having an external plug nozzle, a feature desirable for near-term supersonic aircraft. Statistical analysis of historical data is provided to estimate the uncertainty in using this method given the variations found that cannot be directly computed without detailed mixer geometry.

Noise Prediction↗

Empirical Noise Modeling of Internally Mixed Exhaust Systems

There appears to be no non-proprietary methods to predict the noise of internally mixed exhaust systems, and no guidance for how to adapt known jet noise models for these configurations. This paper surveys literature and historical databases acquired at NASA Glenn's Aero-Acoustic Propulsion Lab to give such guidance. The core premise is that an exhaust system with well-designed mixer produces noise that is to first approximation the same as a fully-mixed jet flow. Refinement of what is meant by a “fully mixed jet” can lead to more accurate prediction of the main jet noise. Additional noise is often generated within the nozzle, typically at high frequencies, whose source mechanism(s) are not obvious. However, a noise prediction method can be established that captures some aspects of the excess noise and provide an estimate of the total jet noise. Explorations of source mechanisms associated with the internal mixer have led to a new noise model which includes the effects of having an external plug nozzle, a feature desirable for near-term supersonic aircraft. Statistical analysis of historical data is provided to estimate the uncertainty in using this method given the variations found that cannot be directly computed without detailed mixer geometry.

Noise Prediction↗

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗

Signal Detection Techniques for Diagnostic Monitoring of Space Shuttle Main Engine Turbomachinery

An investigation to develop, implement, and evaluate signal analysis techniques for the detection and classification of incipient mechanical failures in turbomachinery is reviewed. A brief description of the Space Shuttle Main Engine (SSME) test/measurement program is presented. Signal analysis techniques available to describe dynamic measurement characteristics are reviewed. Time domain and spectral methods are described, and statistical classification in terms of moments is discussed. Several of these waveform analysis techniques have been implemented on a computer and applied to dynamc signals. A laboratory evaluation of the methods with respect to signal detection capability is described. A unique coherence function (the hyper-coherence) was developed through the course of this investigation, which appears promising as a diagnostic tool. This technique and several other non-linear methods of signal analysis are presented and illustrated by application. Software for application of these techniques has been installed on the signal processing system at the NASA/MSFC Systems Dynamics Laboratory.

Coffin, Thomas↗

Statistical mechanics of light elements at high pressure. VI - Liquid-state calculations with Thomas-Fermi-Dirac theory

A model free energy is developed for hydrogen-helium mixtures based on solid-state Thomas-Fermi-Dirac calculations at pressures relevant to the interiors of giant planets. Using a model potential similar to that for a two-component plasma, effective charges for the nuclei (which are in general smaller than the actual charges because of screening effects) are parameterized, being constrained by calculations at a number of densities, compositions, and lattice structures. These model potentials are then used to compute the equilibrium properties of H-He fluids using a charged hard-sphere model. The results find critical temperatures of about 0 K, 500 K, and 1500 K, for pressures of 10, 100, and 1000 Mbar, respectively. These phase separation temperatures are considerably lower (approximately 6,000-10,000 K) than those found from calculations using free electron perturbation theory, and suggest that H-He solutions should be stable against phase separation in the metallic zones of Jupiter and Saturn.

Macfarlane, J. J.↗

Ablation effects in oxygen-lead fragmentation at 2.1 GeV/nucleon

The mechanism of particle evaporation was used to examine ablation effects in the fragmentation of 2.1 GeV/nucleon oxygen nuclei by lead targets. Following the initial abrasion process, the excited projectile prefragment is assumed to statistically decay in a manner analogous to that of a compound nucleus. The decay probabilities for the various particle emission channels are calculated by using the EVAP-4 Monte Carlo computer program. The input excitation energy spectrum for the prefragment is estimated from the geometric ""clean cut'' abrasion-ablation model. Isotope production cross sections are calculated and compared with experimental data and with the predictions from the standard geometric abrasion-ablation fragmentation model.

Townsend, L. W.↗

Error Representation in Time For Compressible Flow Calculations

Time plays an essential role in most real world fluid mechanics problems, e.g. turbulence, combustion, acoustic noise, moving geometries, blast waves, etc. Time dependent calculations now dominate the computational landscape at the various NASA Research Centers but the accuracy of these computations is often not well understood. In this presentation, we investigate error representation (and error control) for time-periodic problems as a prelude to the investigation of feasibility of error control for stationary statistics and space-time averages. o These statistics and averages (e.g. time-averaged lift and drag forces) are often the output quantities sought by engineers. o For systems such as the Navier-Stokes equations, pointwise error estimates deteriorate rapidly which increasing Reynolds number while statistics and averages may remain well behaved.

Barth, Timothy J.↗

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗

Calculation of the thermodynamic properties of a mixture of gases as a function of temperature and pressure

The evaluation of the thermodynamic properties of a gas mixture can be performed using a generalized correlation which makes use of the second virial coefficient. This coefficient is based on statistical mechanics and is a function of temperature and composition, but not of pressure. The method provides results accurate to within 3 percent for gases which are nonpolar or only slightly polar. When applied to highly polar gases, errors of 5 to 10 percent may result. For gases which associate, even larger errors are possible. The sequences of calculations can be routinely programmed for a digital computer. The thermodynamic properties of a mixture of neon, argon and ethane were calculated by such a program. The result will be used for the design of the gas replenishment system for the Energetic Gamma Ray Experiment Telescope.

Colon, G.↗

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗

Space Radiation Effects on Human Cells: Modeling DNA Breakage, DNA Damage Foci Distribution, Chromosomal Aberrations and Tissue Effects

Future long-tem space travel will face challenges from radiation concerns as the space environment poses health risk to humans in space from radiations with high biological efficiency and adverse post-flight long-term effects. Solar particles events may dramatically affect the crew performance, while Galactic Cosmic Rays will induce a chronic exposure to high-linear-energy-transfer (LET) particles. These types of radiation, not present on the ground level, can increase the probability of a fatal cancer later in astronaut life. No feasible shielding is possible from radiation in space, especially for the heavy ion component, as suggested solutions will require a dramatic increase in the mass of the mission. Our research group focuses on fundamental research and strategic analysis leading to better shielding design and to better understanding of the biological mechanisms of radiation damage. We present our recent effort to model DNA damage and tissue damage using computational models based on the physics of heavy ion radiation, DNA structure and DNA damage and repair in human cells. Our particular area of expertise include the clustered DNA damage from high-LET radiation, the visualization of DSBs (DNA double strand breaks) via DNA damage foci, image analysis and the statistics of the foci for different experimental situations, chromosomal aberration formation through DSB misrepair, the kinetics of DSB repair leading to a model-derived spectrum of chromosomal aberrations, and, finally, the simulation of human tissue and the pattern of apoptotic cell damage. This compendium of theoretical and experimental data sheds light on the complex nature of radiation interacting with human DNA, cells and tissues, which can lead to mutagenesis and carcinogenesis later in human life after the space mission.

Ponomarev, A. L.↗

The Transition from Thick to Thin Plate Wake Physics: Whither Vortex Shedding?

The near and very near wake of a flat plate with a circular trailing edge is investigated with data from direct numerical simulations. Computations were performed for six different combinations of the Reynolds numbers based on plate thickness (D) and boundary layer momentum thickness upstream of the trailing edge (theta). Unlike the case of the cylinder, these Reynolds numbers are independent parameters for the flat plate. The separating boundary layers are turbulent in all the cases investigated. One objective of the study is to understand the changes in the wake vortex shedding process as the plate thickness is reduced (increasing theta/D). The value of D varies by a factor of 16 and that of theta by approximately 5 in the computations. Vortex shedding is vigorous in the low theta/D cases with a substantial decrease in shedding intensity in the large theta/D cases. Other shedding characteristics are also significantly altered with increasing theta/D. A visualization of the shedding process in the different cases is provided and discussed. The basic shedding mechanism is explored in depth. The effect of changing theta/D on the time-averaged, near-wake velocity statistics is also discussed. A functional relationship between the shedding frequency and the Reynolds numbers mentioned above is obtained.

CFD↗