Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “randomization”

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,351 records · Page 75

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

Polarimetric backscattering from sea ice

Polarimetric backscattering from sea ice is presented. Theoretical models for backscattering are first presented for various ice types. Then, theoretical results are compared with experimental data for new thin ice, first-year ice, and multi-year ice. Sea ice is modeled as a layer medium containing random scatterers and rough interfaces. For multi-year sea ice with snow cover, the sea ice layer is modeled as an ice background with embedded spheroidal air bubbles and the snow layer as air with ice grains. The hummocky topography on multi-year ice is characterized by a Gaussian distribution which has an averaging effect on backscattering coefficients. First-year sea ice is described by an ice medium hosting ellipsoidal brine inclusions. These inclusions are oriented preferentially in the vertical direction due to the columnar structure of first-year sea ice. Azimuthally, the orientation of the brine pockets are random corresponding to the random c-axes in sea ice, unless the axes are oriented by sea currents. For thin ice in newly opened leads, it has been observed that there exists a thin brine layer with very high salinity on the top surface of the new ice. This brine layer is depicted as a medium with high permittivity which can significantly affect electromagnetic scattering signatures from the lower thin ice layer, with a higher fractional volume of brine inclusions due to higher salinity as compared to thick first-year sea ice. The rough medium interfaces are described as Gaussian rough surfaces characterized by root-mean-square heights and surface correlation lengths. The contribution from rough surface, calculated under the Kirchhoff approximation or small perturbation method, is assumed to be independent from volume scattering. The total loss including absorption and scattering losses in the scattering media is represented by the imaginary part of effective permittivity obtained from the strong fluctuation theory. The polarimetric scattering coefficients for different ice types are then derived under the distorted Born approximation.

Nghiem, S. V.↗

Interaction of the sonic boom with atmospheric turbulence

Theoretical research was carried out to study the effect of free-stream turbulence on sonic boom pressure fields. A new transonic small-disturbance model to analyze the interactions of random disturbances with a weak shock was developed. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. An alternative approach shows that the pressure field may be described by an equation that has an extended form of the classic nonlinear acoustics equation that describes the propagation of sound beams with narrow angular spectrum. The model shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed type elliptic-hyperbolic flows around the shock wave was also developed. Numerical calculations of shock wave interactions with various deterministic and random fluctuations will be presented in a future report.

Rusak, Zvi↗

System level mechanical testing of the Clementine spacecraft

This paper discusses the system level structural testing that was performed to qualify the Clementine Spacecraft for flight. These tests included spin balance, combined acoustic and axial random vibration, lateral random vibration, quasi-static loads, pyrotechnic shock, modal survey and on-orbit jitter simulation. Some innovative aspects of this effort were: the simultaneously combined acoustic and random vibration test; the mass loaded interface modal survey test; and the techniques used to assess how operating on board mechanisms and thrusters affect sensor vision.

Haughton, James↗

A probabilistic model of a porous heat exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Agrawal, O. P.↗

Probabilistic evaluation of fuselage-type composite structures

A methodology is developed to computationally simulate the uncertain behavior of composite structures. The uncertain behavior includes buckling loads, natural frequencies, displacements, stress/strain etc., which are the consequences of the random variation (scatter) of the primitive (independent random) variables in the constituent, ply, laminate and structural levels. This methodology is implemented in the IPACS (Integrated Probabilistic Assessment of Composite Structures) computer code. A fuselage-type composite structure is analyzed to demonstrate the code's capability. The probability distribution functions of the buckling loads, natural frequency, displacement, strain and stress are computed. The sensitivity of each primitive (independent random) variable to a given structural response is also identified from the analyses.

Shiao, Michael C.↗

A diffusion-limited aggregation model for the evolution of drainage networks

We propose a modified diffusion-limited aggregation (DLA) model for the evolution of fluvial drainage networks. Random walkers are introduced randomly on a grid, and each two-dimensional random walk proceeds until the walker finds a drainage network on which to accrete. This model for headward growth of drainage networks generates drainage patterns remarkably similar to actual drainages. The model also predicts statistical features which agree with actual networks, including the frequency-order (bifurcation) ratio (R(sub b) = 3.98) and the stream length-order (R(sub r) = 2.09). Using the definition of network fractal dimension D = log R(sub b)/log R(sub r), we find that our DLA model gives D = 1.87, near the observed range of D approximately equal to 1.80 - 1.85.

Masek, Jeffrey G.↗

Model studies of the beam-filling error for rain-rate retrieval with microwave radiometers

Low-frequency (less than 20 GHz) single-channel microwave retrievals of rain rate encounter the problem of beam-filling error. This error stems from the fact that the relationship between microwave brightness temperature and rain rate is nonlinear, coupled with the fact that the field of view is large or comparable to important scales of variability of the rain field. This means that one may not simply insert the area average of the brightness temperature into the formula for rain rate without incurring both bias and random error. The statistical heterogeneity of the rain-rate field in the footprint of the instrument is key to determining the nature of these errors. This paper makes use of a series of random rain-rate fields to study the size of the bias and random error associated with beam filling. A number of examples are analyzed in detail: the binomially distributed field, the gamma, the Gaussian, the mixed gamma, the lognormal, and the mixed lognormal ('mixed' here means there is a finite probability of no rain rate at a point of space-time). Of particular interest are the applicability of a simple error formula due to Chiu and collaborators and a formula that might hold in the large field of view limit. It is found that the simple formula holds for Gaussian rain-rate fields but begins to fail for highly skewed fields such as the mixed lognormal. While not conclusively demonstrated here, it is suggested that the notionof climatologically adjusting the retrievals to remove the beam-filling bias is a reasonable proposition.

Ha, Eunho↗

On the spatial distribution of magnetic fields on the solar surface

Recent measurements of solar surface magnetic fields suggest that the spatial distribution of these fields is fractal. In order to understand the physical basis for such geometric complexity, we study here the advection of magnetic flux tubes relatively simple random motions on the surface of a fluid and investigate the spatial statistics of the resulting surface field. While this study does not directly address the question of why solar surface fields have the observed spatial structure, it is designed to build our intuition about how surface flows lead to complex spatial structuring of magnetic fields. As part of our study, we discuss the various methods by which one can describe the spatial distribution of the surface magnetic flux and relate them mathematically; this turns out to be a crucial point of our work since, as we show, a number of previous analyses have misinterpreted the analysis procedures for determining fractal dimensions. Our principal result is the explicit demonstration that simple random flows lead to magnetic flux spatial distributions with a multifractal dimension spectrum. Furthermore, we demonstrate that this magnetic spatial structure is generic, i.e., is characteristic of a very large class of random flows.

Tao, L.↗

Estimation and classification by sigmoids based on mutual information

An estimate of the probability density function of a random vector is obtained by maximizing the mutual information between the input and the output of a feedforward network of sigmoidal units with respect to the input weights. Classification problems can be solved by selecting the class associated with the maximal estimated density. Newton's s method, applied to an estimated density, yields a recursive maximum likelihood estimator, consisting of a single internal layer of sigmoids, for a random variable or a random sequence. Applications to the diamond classification and to the prediction of a sun-spot process are demonstrated.

Baram, Yoram↗

Probabilistic Material Strength Degradation Model for Inconel 718 Components Subjected to High Temperature, High-Cycle and Low-Cycle Mechanical Fatigue, Creep and Thermal Fatigue Effects

The development of methodology for a probabilistic material strength degradation is described. The probabilistic model, in the form of a postulated randomized multifactor equation, provides for quantification of uncertainty in the lifetime material strength of aerospace propulsion system components subjected to a number of diverse random effects. This model is embodied in the computer program entitled PROMISS, which can include up to eighteen different effects. Presently, the model includes five effects that typically reduce lifetime strength: high temperature, high-cycle mechanical fatigue, low-cycle mechanical fatigue, creep and thermal fatigue. Results, in the form of cumulative distribution functions, illustrated the sensitivity of lifetime strength to any current value of an effect. In addition, verification studies comparing predictions of high-cycle mechanical fatigue and high temperature effects with experiments are presented. Results from this limited verification study strongly supported that material degradation can be represented by randomized multifactor interaction models.

Bast, Callie C.↗

Fluid Physics in a Fluctuating Acceleration Environment

We summarize several aspects of an ongoing investigation of the effects that stochastic residual accelerations (g-jitter) onboard spacecraft can have on experiments conducted in a microgravity environment. The residual acceleration field is modeled as a narrow band noise, characterized by three independent parameters: intensity (g(exp 2)), dominant angular frequency Omega, and characteristic correlation time tau. Realistic values for these parameters are obtained from an analysis of acceleration data corresponding to the SL-J mission, as recorded by the SAMS instruments. We then use the model to address the random motion of a solid particle suspended in an incompressible fluid subjected to such random accelerations. As an extension, the effect of jitter on coarsening of a solid-liquid mixture is briefly discussed, and corrections to diffusion controlled coarsening evaluated. We conclude that jitter will not be significant in the experiment 'Coarsening of solid-liquid mixtures' to be conducted in microgravity. Finally, modifications to the location of onset of instability in systems driven by a random force are discussed by extending the standard reduction to the center manifold to the stochastic case. Results pertaining to time-modulated oscillatory convection are briefly discussed.

Thomson, J. Ross↗

Simulation of Crack Propagation in Engine Rotating Components under Variable Amplitude Loading

The crack propagation life of tested specimens has been repeatedly shown to strongly depend on the loading history. Overloads and extended stress holds at temperature can either retard or accelerate the crack growth rate. Therefore, to accurately predict the crack propagation life of an actual component, it is essential to approximate the true loading history. In military rotorcraft engine applications, the loading profile (stress amplitudes, temperature, and number of excursions) can vary significantly depending on the type of mission flown. To accurately assess the durability of a fleet of engines, the crack propagation life distribution of a specific component should account for the variability in the missions performed (proportion of missions flown and sequence). In this report, analytical and experimental studies are described that calibrate/validate the crack propagation prediction capability ]or a disk alloy under variable amplitude loading. A crack closure based model was adopted to analytically predict the load interaction effects. Furthermore, a methodology has been developed to realistically simulate the actual mission mix loading on a fleet of engines over their lifetime. A sequence of missions is randomly selected and the number of repeats of each mission in the sequence is determined assuming a Poisson distributed random variable with a given mean occurrence rate. Multiple realizations of random mission histories are generated in this manner and are used to produce stress, temperature, and time points for fracture mechanics calculations. The result is a cumulative distribution of crack propagation lives for a given, life limiting, component location. This information can be used to determine a safe retirement life or inspection interval for the given location.

Bonacuse, P. J.↗

Probabilistic Evaluation of Fuselage-Type Composite Structures

A methodology is developed to simulate computationally the uncertain behavior of composite structures. The uncertain behavior includes buckling loads, natural frequencies, displacements, stress/strain, etc., which are the consequences of the random variation (scatter) of the primitive (independent random) variables in the constituent, plv. laminate and structural levels. This methodology is implemented in a computer code integrated probabilistic assessment of composite structures (IPACS). A fuselage-type composite structure is analyzed to demonstrate the code's capability . The probability distribution functions of the buckling loads, natural frequency, displacement, strain and stress are computed. The sensitivity of each primitive (independent random) variable to a given structural response is also identified from the analyses.

Shiao, Michael C.↗

Snow Crystal Orientation Effects on the Scattering of Passive Microwave Radiation

For this study, consideration is given to the role crystal orientation plays in scattering and absorbing microwave radiation. A discrete dipole scattering model is used to measure the passive microwave radiation, at two polarizations (horizontal and vertical), scattered by snow crystals oriented in random and non random positions, having various sizes (ranging between 1 micrometers to 10,000 micrometers in radius), and shapes (including spheroids, cylinders, hexagons). The model results demonstrate that for the crystal sizes typically found in a snowpack, crystal orientation is insignificant compared to crystal size in terms of scattering microwave energy in the 8,100 gm (37 GHz) region of the spectrum. Therefore, the assumption used in radiative transfer approaches, where snow crystals are modeled as randomly oriented spheres, is adequate to account for the transfer of microwave energy emanating from the ground and passing through a snowpack.

Foster, J. L.↗

Simulation of Crack Propagation in Engine Rotating Components Under Variable Amplitude Loading

The crack propagation life of tested specimens has been repeatedly shown to strongly depend on the loading history. Overloads and extended stress holds at temperature can either retard or accelerate the crack growth rate. Therefore, to accurately predict the crack propagation life of an actual component, it is essential to approximate the true loading history. In military rotorcraft engine applications, the loading profile (stress amplitudes, temperature, and number of excursions) can vary significantly depending on the type of mission flown. To accurately assess the durability of a fleet of engines, the crack propagation life distribution of a specific component should account for the variability in the missions performed (proportion of missions flown and sequence). In this report, analytical and experimental studies are described that calibrate/validate the crack propagation prediction capability for a disk alloy under variable amplitude loading. A crack closure based model was adopted to analytically predict the load interaction effects. Furthermore, a methodology has been developed to realistically simulate the actual mission mix loading on a fleet of engines over their lifetime. A sequence of missions is randomly selected and the number of repeats of each mission in the sequence is determined assuming a Poisson distributed random variable with a given mean occurrence rate. Multiple realizations of random mission histories are generated in this manner and are used to produce stress, temperature, and time points for fracture mechanics calculations. The result is a cumulative distribution of crack propagation lives for a given, life limiting, component location. This information can be used to determine a safe retirement life or inspection interval for the given location.

P J Bonacuse↗

Fluid Physics Under a Stochastic Acceleration Field

The research summarized in this report has involved a combined theoretical and computational study of fluid flow that results from the random acceleration environment present onboard space orbiters, also known as g-jitter. We have focused on a statistical description of the observed g-jitter, on the flows that such an acceleration field can induce in a number of experimental configurations of interest, and on extending previously developed methodology to boundary layer flows. Narrow band noise has been shown to describe many of the features of acceleration data collected during space missions. The scale of baroclinically induced flows when the driving acceleration is random is not given by the Rayleigh number. Spatially uniform g-jitter induces additional hydrodynamic forces among suspended particles in incompressible fluids. Stochastic modulation of the control parameter shifts the location of the onset of an oscillatory instability. Random vibration of solid boundaries leads to separation of boundary layers. Steady streaming ahead of a modulated solid-melt interface enhances solute transport, and modifies the stability boundaries of a planar front.

Vinals, Jorge↗

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗