Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational 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 649 records · Page 36

Theory and procedure for determining loads and motions in chine-immersed hydrodynamic impacts of prismatic bodies

A theoretical method is derived for the determination of the motions and loads during chine-immersed water landings of prismatic bodies. This method makes use of a variation of two-dimensional deflected water mass over the complete range of immersion, modified by a correction for three-dimensional flow. Equations are simplified through omission of the term proportional to the acceleration of the deflected mass for use in calculation of loads on hulls having moderate and heavy beam loading. The effects of water rise at the keel are included in these equations. In order to make a direct comparison of theory with experiment, a modification of the equations was made to include the effect of finite test-carriage mass. A simple method of computation which can be applied without reading the body of this report is presented as an appendix along with the required theoretical plots for determination of loads and motions in chine-immersed landings.

Schnitzer, Emanuel↗

Drawbar Pull (DP) Procedures for Off-Road Vehicle Testing

As NASA strives to explore the surface of the Moon and Mars, there is a continued need for improved tire and vehicle development. When tires or vehicles are being designed for off-road conditions where significant thrust generation is required, such as climbing out of craters on the Moon, it is important to use a standard test method for evaluating their tractive performance. The drawbar pull (DP) test is a way of measuring the net thrust generated by tires or a vehicle with respect to performance metrics such as travel reduction, sinkage, or power efficiency. DP testing may be done using a single tire on a traction rig, or with a set of tires on a vehicle; this report focuses on vehicle DP tests. Though vehicle DP tests have been used for decades, there are no standard procedures that apply to exploration vehicles. This report summarizes previous methods employed, shows the sensitivity of certain test parameters, and provides a body of knowledge for developing standard testing procedures. The focus of this work is on lunar applications, but these test methods can be applied to terrestrial and planetary conditions as well. Section 1.0 of this report discusses the utility of DP testing for off-road vehicle evaluation and the metrics used. Section 2.0 focuses on test-terrain preparation, using the example case of lunar terrain. There is a review of lunar terrain analogs implemented in the past and a discussion on the lunar terrain conditions created at the NASA Glenn Research Center, including methods of evaluating the terrain strength variation and consistency from test to test. Section 3.0 provides details of the vehicle test procedures. These consist of a review of past methods, a comprehensive study on the sensitivity of test parameters, and a summary of the procedures used for DP testing at Glenn.

wheel↗

Simultaneous Noncontact Precision Imaging of Microstructural and Thickness Variation in Dielectric Materials Using Terahertz Energy

This article describes a noncontact single-sided terahertz electromagnetic measurement and imaging method that simultaneously characterizes microstructural (egs. spatially-lateral density) and thickness variation in dielectric (insulating) materials. The method was demonstrated for two materials-Space Shuttle External Tank sprayed-on foam insulation and a silicon nitride ceramic. It is believed that this method can be used as an inspection method for current and future NASA thermal protection system and other dielectric material inspection applications, where microstructural and thickness variation require precision mapping. Scale-up to more complex shapes such as cylindrical structures and structures with beveled regions would appear to be feasible.

Roth, Don J.↗

A spherical harmonics method for asteroid pole determination

A new and general method is presented for determining pole orientations of unresolved solar system bodies that exhibit brightness variation with rotation. The accuracy of the method is limited only by the number of available apparitions, which imposes the degree where the Legendre expansion must be truncated. It is shown that power spectra of lightcurves provide enough information to extract pole coordinates from the data, independent of body shape and albedo variegation. The method has been tested using several computer-generated synthetic light-curves, and the accuracy is found to be about 5 deg in both the longitude and latitude of the pole, which is considerably better than that obtained from the standard triaxial ellipsoid model. Finally, the new spherical harmonics method has been applied to six well-observed asteroids: 16 Psyche, 22 Kalliope, 39 Laetitia, 44 Nysa, 45 Eugenia, and 349 Dembowska.

Lumme, K.↗

On the Numerical Formulation of Parametric Linear Fractional Transformation (LFT) Uncertainty Models for Multivariate Matrix Polynomial Problems

Robust control system analysis and design is based on an uncertainty description, called a linear fractional transformation (LFT), which separates the uncertain (or varying) part of the system from the nominal system. These models are also useful in the design of gain-scheduled control systems based on Linear Parameter Varying (LPV) methods. Low-order LFT models are difficult to form for problems involving nonlinear parameter variations. This paper presents a numerical computational method for constructing and LFT model for a given LPV model. The method is developed for multivariate polynomial problems, and uses simple matrix computations to obtain an exact low-order LFT representation of the given LPV system without the use of model reduction. Although the method is developed for multivariate polynomial problems, multivariate rational problems can also be solved using this method by reformulating the rational problem into a polynomial form.

Belcastro, Christine M.↗

High-Temperature Slow Crack Growth of Silicon Carbide Determined by Constant-Stress-Rate and Constant-Stress Testing

High-temperature slow-crack-growth behaviour of hot-pressed silicon carbide was determined using both constant-stress-rate ("dynamic fatigue") and constant-stress ("static fatigue") testing in flexure at 1300 C in air. Slow crack growth was found to be a governing mechanism associated with failure of the material. Four estimation methods such as the individual data, the Weibull median, the arithmetic mean and the median deviation methods were used to determine the slow crack growth parameters. The four estimation methods were in good agreement for the constant-stress-rate testing with a small variation in the slow-crack-growth parameter, n, ranging from 28 to 36. By contrast, the variation in n between the four estimation methods was significant in the constant-stress testing with a somewhat wide range of n= 16 to 32.

Choi, Sung H.↗

Evaluation of selected strapdown inertial instruments and pulse torque loops, volume 1

Design, operational and performance variations between ternary, binary and forced-binary pulse torque loops are presented. A fill-in binary loop which combines the constant power advantage of binary with the low sampling error of ternary is also discussed. The effects of different output-axis supports on the performance of a single-degree-of-freedom, floated gyroscope under a strapdown environment are illustrated. Three types of output-axis supports are discussed: pivot-dithered jewel, ball bearing and electromagnetic. A test evaluation on a Kearfott 2544 single-degree-of-freedom, strapdown gyroscope operating with a pulse torque loop, under constant rates and angular oscillatory inputs is described and the results presented. Contributions of the gyroscope's torque generator and the torque-to-balance electronics on scale factor variation with rate are illustrated for a SDF 18 IRIG Mod-B strapdown gyroscope operating with various pulse rebalance loops. Also discussed are methods of reducing this scale factor variation with rate by adjusting the tuning network which shunts the torque coil. A simplified analysis illustrating the principles of operation of the Teledyne two-degree-of-freedom, elastically-supported, tuned gyroscope and the results of a static and constant rate test evaluation of that instrument are presented.

Sinkiewicz, J. S.↗

The use of parabolic variations and the direct determination of stress intensity factors using the BIE method

Two advances in the numerical techniques of utilizing the BIE method are presented. The boundary unknowns are represented by parabolas over each interval which are integrated in closed form. These integrals are listed for easy use. For problems involving crack tip singularities, these singularities are included in the boundary integrals so that the stress intensity factor becomes just one more unknown in the set of boundary unknowns thus avoiding the uncertainties of plotting and extrapolating techniques. The method is applied to the problems of a notched beam in tension and bending, with excellent results.

Mendelson, A.↗

The use of parabolic variations and the direct determination of stress intensity factors using the BIE method

Two advances in the numerical techniques of utilizing the BIE method are presented. The boundary unknowns are represented by parabolas over each interval which are integrated in closed form. These integrals are listed for easy use. For problems involving crack tip singularities, these singularities are included in the boundary integrals so that the stress intensity factor becomes just one more unknown in the set of boundary unknowns thus avoiding the uncertainties of plotting and extrapolating techniques. The method is applied to the problems of a notched beam in tension and bending, with excellent results.

Mendelson, A.↗

Shape optimization of 3D continuum structures via force approximation techniques

The existing need to develop methods whereby the shape design efficiency can be improved through the use of high quality approximation methods is addressed. An efficient approximation method for stress constraints in 3D shape design problems is proposed based on expanding the nodal forces in Taylor series with respect to shape variations. The significance of this new method is shown through elementary beam theory calculations and via numerical computations using 3D solid finite elements. Numerical examples including the classical cantilever beam structure and realistic automotive parts like the engine connecting rod are designed for optimum shape using the proposed method. The numerical results obtained from these methods are compared with other published results, to assess the efficiency and the convergence rate of the proposed method.

Vanderplaats, Garret N.↗

Twenty years of interplanetary magnetic field variations with periods in the range of 10 days to 3 years

Twenty years of interplanetary magnetic field data collected primarily by the IMP-8 spacecraft near Earth has been analyzed by a dynamic periodogram method in search of significant periodicities in the range of 10 days to 3 years. The method has the advantage of detecting variations with time in the periodicities besides determining the power and phase of the dominant variations. It has been found that the well known periodicities near 1 year and 27 days are strongly modulated by the solar cycle. Both of these periodicities are only detected during solar minimum. During solar maximum. a number of unusual variations are observed. Special emphasis will be placed on the recently reported 1.3 year variation in solar wind parameters besides periods in the interplanetary magnetic field near 51, 73 and 154 days. Correlations with solar wind plasma and solar index variations will also be presented.

Szabo, A.↗

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗

Ultrasonic evaluation of the strength of unidirectional graphite-polyimide composites

An acoustic-ultrasonic method is described that was successful in ranking unidirectional graphite-polyimide composite specimens according to variations in interlaminar shear strength. Using this method, a quantity termed the stress wave factor was determined. It was found that this factor increases directly with interlaminar shear strength. The key variables in this investigation were composite density, fiber weight fraction, and void content. The stress wave factor and other ultrasonic factors that were studied were found to provide a powerful means for nondestructive evaluation of mechanical strength properties.

Vary, A.↗

Determination of Stress Intensity Factor Distributions for "Interface" Cracks in Incompressible, Dissimilar Materials

The present study was undertaken in order to develop test methods and procedures for measuring the variation of the stress intensity factor through the thickness in bimaterial specimens containing cracks within and parallel to the bond line using the frozen stress photoelastic method. Since stress freezing materials are incompressible above critical temperature, and since thick plates are to be employed which tend to produce a state of plane strain near the crack tip, the interface near tip fracture equations reduce to the classic form for homogeneous materials. Moreover, zero thickness interfaces do not exist when materials are bonded together. It was decided early on that it would be important to insure a uniform straight and accurate crack tip region through the thickness of the body to reduce scatter in the SIF distribution through the thickness. It was also observed that rubberlike materials which were desired to be modeled exhibited significant tip blunting prior to crack extension and that some blunting of the tip would provide a more realistic model. It should be noted that, in normal stress freezing photoelastic work, it is considered good practice to avoid utilizing data near bond lines in photoelastic models due to the bond line stresses which inevitably develop when two parts are bonded together. Thus, the present study involves certain exploratory aspects in deviating from standard practice in stress freezing work. With the above ideas in mind, several different test methods were investigated and are described in the following sections and appendices. The geometry selected for the program was a thick, edge cracked specimen containing a bond line.

Smith, C. W.↗

Quantifying Changes in the Arctic Shortwave Cloud Radiative Effects

The shortwave cloud radiative effect (SWCRE) is important on the Arctic surface radiation budget and the major source of inter-model spread in predictions of Arctic climate. To better understand the individual contributions of various radiative processes to changes in SWCRE, the paper presents the use of the extended APRP (Atmospheric Radiative Perturbation Potential) method. This involves adding the absorptivity for the upward beam and considering differences in reflectivity between upward and downward beams, as well as analyzing the cloud masking effect resulting from changes in surface albedo in more detail. Using data from the CMIP5 and CMIP6 climate models, the study decomposes the SWCRE over the Arctic surface and analyzes inter-model differences in quadrupled CO2 simulations. The study takes into account the fact that the response of SWCRE to Arctic warming is influenced by changes in surface albedo, cloud amount, and cloud microphysics. Results show that in the sunlight season, the reduction in surface albedo associated with sea ice loss is directly linked to strong negative SWCRE, which explains the considerable model discrepancy. Arctic clouds can hinder the positive surface albedo feedback by changing the albedo in two ways: (1) decreasing incoming shortwave radiation due to cloud reflection and (2) by decreasing the shortwave reaching the surface after being reflected by clouds. In addition, increased (decreased) cloud amount and cloud liquid water are shown to be less (more) incoming shortwave fluxes at the surface, but not dominating factors to the Arctic surface radiation budget and its inter-model variation. Overall, the extended APRP method offers a useful tool for analyzing the complex interactions between clouds and radiative process, reasonably decomposes the individual SWCRE responses at the Arctic surface, and emphasizes that considering not only the cloud amount or its properties, but also surface albedo change is critical for the prediction of SWCRE on the Arctic surface.

Arctic cloud↗

Latitudinal and altitudinal variation of size distribution of stratospheric aerosols inferred from SAGE aerosol extinction coefficient measurements at two wavelengths

A method of retrieving aerosol size distribution from the measured extinction of solar radiation at wavelengths of 0.45 microns and 1.0 microns has recently been proposed. This method is utilized to obtain latitudinal and altitudinal variations of size distributions of stratospheric aerosols from the Stratospheric Aerosol and Gas Experiment data for March 1979. Small particles are found in the lower stratosphere of the tropical region, and large particles are found at higher altitudes and latitudes in both hemispheres. Results of this study are consistent with the suggestion that the upper troposphere in tropical regions is a source of condensation nuclei in the stratosphere, and they become mature as they move to higher altitudes and latitude.

Yue, G. K.↗

Solid-propellant rocket motor internal ballistic performance variation analysis, phase 2

The Monte Carlo method was used to investigate thrust imbalance and its first time derivative throughtout the burning time of pairs of solid rocket motors firing in parallel. Results obtained compare favorably with Titan 3 C flight performance data. Statistical correlations of the thrust imbalance at various times with corresponding nominal trace slopes suggest several alternative methods of predicting thrust imbalance. The effect of circular-perforated grain deformation on internal ballistics is discussed, and a modified design analysis computer program which permits such an evaluation is presented. Comparisons with SRM firings indicate that grain deformation may account for a portion of the so-called scale factor on burning rate between large motors and strand burners or small ballistic test motors. Thermoelastic effects on burning rate are also investigated. Burning surface temperature is calculated by coupling the solid phase energy equation containing a strain rate term with a model of gas phase combustion zone using the Zeldovich-Novozhilov technique. Comparisons of solutions with and without the strain rate term indicate a small but possibly significant effect of the thermoelastic coupling.

Sforzini, R. H.↗

Three-dimensional modelling in magnetotelluric and magnetic variational sounding

The Galerkin finite-element method is used to obtain approximate solutions for the three-dimensional induction problem. A rectangular conductive prism is considered as an example, and solutions are obtained for linear and circularly polarized incident plane-wave fields. Magnetotelluric tensor impedances and magnetic transfer functions are computed. Polar diagrams of the tensor impedances and magnetic transfer functions along with their amplitude contour maps are presented. The dimensionality parameter, skew, is contoured at the surface of the earth. It is shown that the relative amplitudes and shapes of the additional and principal impedance polar diagrams can be used to determine the dimensionality of geoelectrical structures. Stations with skew values greater than 0.2 are significantly influenced by the three-dimensionality of the geoelectric structure. The amplitudes of the magnetic transfer function and the orientations of its polar diagrams exhibit large anomalies in the vicinity of the intersection of the lateral contacts.

Reddy, I. K.↗