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 937 records · Page 52

Composite structure dynamics and space applicability

Guided wave modes in cladded or uncladded fiber-reinforced composite plates and tubes have been analyzed using a stiffness method in which the displacement variation through the thickness is approximated by polynomial interpolation functions. This allows for an arbitrary number of laminations and fiber orientations different from lamina to lamina. It is shown that dispersive behavior of guided modes depends significantly on the cladding, number of laminae, and interfaces between the adjacent laminae. A hybrid modeling technique is described in which an inner region containing cracks (or other defects) is discretized by finite elements, and the field in the exterior region is represented in terms of modes that are found using the stiffness method described above. It is found that the reflected and transmitted amplitudes of modes vary significantly with the size of a transverse or longitudinal (delamination) crack and frequency. We have also studied the impact response of a unidirectional fiber-reinforced plate. Received signals at the epicentral and other locations are shown. Strong longitudinal anisotropy of the graphite/epoxy plate causes the signal to be considerably different from that in an isotropic plate.

Datta, Subhendu K.↗

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗

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↗

Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle

Within this work, we establish a stochastic contact variational integrator and its discrete version via stochastic Herglotz variational principle for stochastic contact Hamiltonian systems. A general structure-preserving stochastic contact method is provided to seek the stochastic contact variational integrators. Numerical experiments are performed to verify the validity of this approach.

97 MATHEMATICS AND COMPUTING↗

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.↗

New data-driven approach to bridging power system protection gaps with deep learning

Protection is a critical function in power systems to avoid equipment damage, maintain personnel safety, and support system reliability. However, current protective relay technology cannot adequately protect equipment and personnel from effects of some events; these deficiencies are termed protection gaps. In this paper, a data-driven approach is proposed to complement traditional protection technology and distinguish fault conditions from transients caused by normal operations. A combined convolutional neural network and long short-term memory (CNN-LSTM) network is implemented to achieve data translation invariance and capture the temporal correlation of the time-series input data. As a result, the data-driven method can accurately detect system faults despite variation and noise in the input data. In addition, using the CNN-LSTM--based method avoids the complicated, manual feature extraction procedure required by many traditional data-driven methods. The effectiveness of the proposed approach is tested on two kinds of protection gaps: high-impedance faults and transformer inter-turn faults. Lastly, a transfer learning method is also proposed to address the common issue of data-driven methods for which real-world training data are scarce. Extensive study results demonstrate that the proposed approach can accurately bridge power system protection gaps.

42 ENGINEERING↗

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.↗

Introducing the SLICE Method for estimating pebble-bed reactor inventories at equilibrium operation with SCALE

This paper introduces the SCALE Leap-In method for Cores at Equilibrium (SLICE) for estimating pebble-bed reactor equilibrium core isotopic inventories using capabilities in the SCALE code system, requiring only a small computational cluster and a few days of computation. This method uses an iterative approach that relies on (1) a surrogate spectrum model that captures spatial and time-dependent spectral conditions, (2) a multi-pass model that captures the pebble’s evolving nuclide inventory as a function of location and time in the core, and (3) a full-core model that captures the core’s spatial neutron flux distribution. The SLICE approach is applied to a generic fluoride salt–cooled high-temperature reactor, demonstrating fuel inventory convergence through nuclide concentration inspection across iterations and comparisons for core realizations with varying discretizations. Results agree within ~5% with another state-of-the-art code, with differences attributed to input parameter or modeling assumption variations in the equilibrium generation methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Revisiting excitation gaps in the fractional quantum Hall effect

Recent systematic measurements of the quantum well width dependence of the excitation gaps of fractional quantum Hall states in high mobility samples open the possibility of a better quantitative understanding of this important issue. Here we present what we believe to be accurate theoretical gaps including the effects of finite width and Landau level (LL) mixing. While theory captures the width dependence, there still remains a deviation between the calculated and the measured gaps, presumably caused by disorder. It is customary to model the experimental gaps of the n/(2n ± 1) states as Δ n/(2n+1) = Ce 2 /[(2n ± 1)ϵl]-Γ, where ϵ is the dielectric constant of the background semiconductor and l is the magnetic length; the first term is interpreted as the cyclotron energy of composite fermions and Γ as a disorder-induced broadening of composite-fermion LLs. Fitting the gaps for various fractional quantum Hall states, we find that Γ can be nonzero even in the absence of disorder.

36 MATERIALS SCIENCE↗

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↗