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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Nebular properties from far-infrared spectrosopy

We describe a semiempirical methodology-based on measurements of far-infrared (FIR) lines-that yields information on electron densities in regions where various ionic species exist, effective temperatures (T(sub eff)) for stars ionizing H II regions, and gas-phase heavy element abundances. Although this capability has long been available via optical data, the special features of FIR lines-relative insensitivity to extinction and electron temperature variations-extend the analysis ability. Several line ratios serve as diagnostics of electron density, N(sub e), probing different ionization conditions and different density regimes. The more N(sub e)-diagnostic observations made, the more reliable will be the deciphering of the actual variation in density throughout a nebula. A method to estimate T(sub eff) from the FIR (N III)/(N II) line ratio requires that the nebula be ionization bounded and that substantially all of the flux from the revevant lines be observed. However, to estimate T(sub eff) by a second method that uses the ratio of FIR (S III)/(O III) lines, an ionization-bounded nebula is a sufficient, but not necessary, condition. These restrictions are unnecessary for estimating densities and heavy element abundances. We show that a fairly general determination of metallicity, via the S/H ratio, may be made for H II regions with observations of just two lines-(S III) 19 micron and a hydrogen recombination line (or appropriate substitute). These techniques are applied to recent FIR data for the G333.6-0.2 H II region, including application to the recently measured (N II) 122 and 205 micron lines.

Rubin, Robert H.↗
NASA NTRS1995

Spectroscopic study of the seasonal variation of carbon monoxide vertical distribution above Kitt Peak

A method to retrieve elements of the carbon monoxide (CO) vertical distribution from ground-based high-resolution infrared solar spectra has been developed. The method is based on the fact that the total column amount retrieved by nonlinear least squares spectral fitting techniques depends on the shape of the assumed a priori profile and this dependence is a function of the absorption line intensity and the lower state energy of the transition. Four CO lines between 2057 and 2159/cm have been selected and the method has been tested on synthetic spectra. The CO total column content and average concentrations in two atmospheric layers (surface to 400 mbar and 400 mbar to the top of the atmosphere) can be retrieved with precisions of about 1% and less than 10%, respectively. Solar spectra recorded at Kitt Peak from 1982 to 1993 have been analyzed. The CO total column and the average concentration in the two layers show an asymmetrical seasonal cycle with extreme values of (1.1-2.1) x 10(exp 18) molecules cm(exp -2), (50-80) parts per billion by volume (ppbv) in the top layer and (80-160) ppbv in the bottom layer, and precisions of 1, 3, and 6%, respectively; a spring maximum and late summer minimum are observed.

Pougatchev, N. S.↗
NASA NTRS1994

Accuracy assessment of the large-scale dynamic ocean topography from TOPEX/POSEIDON altimetry

The quality of TOPEX/POSEIDON determinations of the global scale dynamic ocean topography have been assessed by determining mean topography solutions for successive 10-day repeat cycles and by examining the temporal changes in the sea surface topography to identify known features. The assessment is based on the analysis of TOPEX altimeter data cycles 1 through 36. Important errors in the tide model used to correct the altimeter data have been identified. The errors were reduced significantly by use of a new tide model derived with the TOPEX/POSEIDON measurements. Maps of the global 1-year mean topography, produced using four of the most accurate of the marine geoid, show that the largest error in the dynamic ocean topography show expected features, such as the known annual hemispherical sea surface rise and fall and the seasonal variability due to monsoon influence in the Indian Ocean. Changes in the sequence of 10-day topography maps show the development and propagation of an equatorial Kelvin wave in the Pacific beginning in December 1992 with a propagation velocity of approximately 3 m/s. The observations are consistent with observed changes in the equatorial trade winds, and with tide gauge and other in situ observations of the strengthening of the El Nino. Comparison of TOPEX-determine sea surface height at points near oceanic tide gauges shows agreement at the 4 cm root-mean-square (RMS) level over the tropical Pacific. The results show that the TOPEX altimeter data set can be used to map the ocean surface with a temporal resolution of 10 days and an accuracy which is insonsistent with traditional in situ methods for the determination of sea level variations.

Tapley, B. D.↗
NASA NTRS1995

Nonperiodic variations in astrophysical systems: Investigating frequency evolution

We present a method related to the wavelet transform, the Gabor transform, for investigating astronomical time series containing nonconsistent frequencies. Instances in which such data sets may arise include variable star light curves, numerical studies of the gravitational three-body problem, X-ray binaries, and signals from more exotic objects such as planets around pulsars and mass infall from accretion structures onto compact objects. As an illustration of its power, we apply the technique to a numerical data set of a gravitational three-body interaction and to photometry of the rapidly oscillating peculiar A star HD 60435. In the three-body example, the method provides an insightful shorthand that allows for the determination of episodes where the system behaves as two nearly Keplerian orbits. For HD 60435, the power in the main frequency exhibits unusual evolution over the duration of the observation.

Boyd, Patricia T.↗
NASA NTRS2004

Determining Cloud Parameters with the Curve-Of-Growth: Application Eta Car

We have investigated the NUV part of the Eta Car spectrum, using data with high spatial and high spectral resolving power obtained with the HST/STIS under the Treasury Program. The NUV spectrum of Eta Car Shows a great contribution of absorption features from neutral and singly ionized elements along the line-of-sight. A large number of velocity systems have been observed. The two most prominent, with Doppler shifts corresponding to -146 and -513 km/s respectively, are shown to be useful for investigations of the gaseous environments responsible for the absorption. The -146 and the -513 km/s velocity systems display different characteristics regarding the ionization state and spectral line width, which suggest that they originate at different distances from the central object. We have investigated the absorption structures before the spectroscopic minimum, occurring during the summer of 2003, with a standard curve-of-growth. We have independently derived the column density and the b-value for the Fe II (-146 km/s) and Ti II (-513 km/s) velocity systems. The excitation temperature has been determined for the -146 km/s velocity system using the photo-ionization code \textsc(cloudy). The -146 km/s velocity structure shows noticeable variation over the spectroscopic minimum. The sudden appearance and disappearance of Ti II and V II are astonishing. We have made an attempt to analyze these variations with the curve-of-growth method and will present preliminary results.

Vieira, G. L.↗
NASA NTRS2011

Dynamic Shape Reconstruction of Three-Dimensional Frame Structures Using the Inverse Finite Element Method

A robust and efficient computational method for reconstructing the three-dimensional displacement field of truss, beam, and frame structures, using measured surface-strain data, is presented. Known as shape sensing , this inverse problem has important implications for real-time actuation and control of smart structures, and for monitoring of structural integrity. The present formulation, based on the inverse Finite Element Method (iFEM), uses a least-squares variational principle involving strain measures of Timoshenko theory for stretching, torsion, bending, and transverse shear. Two inverse-frame finite elements are derived using interdependent interpolations whose interior degrees-of-freedom are condensed out at the element level. In addition, relationships between the order of kinematic-element interpolations and the number of required strain gauges are established. As an example problem, a thin-walled, circular cross-section cantilevered beam subjected to harmonic excitations in the presence of structural damping is modeled using iFEM; where, to simulate strain-gauge values and to provide reference displacements, a high-fidelity MSC/NASTRAN shell finite element model is used. Examples of low and high-frequency dynamic motion are analyzed and the solution accuracy examined with respect to various levels of discretization and the number of strain gauges.

Gherlone, Marco↗
NASA NTRS2015

Composite Characterization Using Laser Doppler Vibrometry and Multi-Frequency Wavenumber Analysis

NASA has recognized the need for better characterization of composite materials to support advances in aeronautics and the next generation of space exploration vehicles. An area of related research is the evaluation of impact induced delaminations. Presented is a non-contact method of measuring the ply depth of impact delamination damage in a composite through use of a Scanning Laser Doppler Vibrometer (SLDV), multi-frequency wavenumber analysis, and a wavenumber-ply correlation algorithm. A single acquisition of a chirp excited lamb wavefield in an impacted composite is post-processed into a numerous single frequency excitation wavefields through a deconvolution process. A spatially windowed wavenumber analysis then extracts local wavenumbers from the wavefield, which are then correlated to theoretical dispersion curves for ply depth determination. SLDV based methods to characterize as-manufactured composite variation using wavefield analysis will also be discussed.

Juarez, Peter↗
NASA NTRS2017

Asymmetric Expansion of the Youngest Galactic Supernova Remnant G1.9+0.3

The youngest Galactic supernova remnant (SNR) G1.9+0.3, produced by a (probable) SN Ia that exploded approximately 1900 CE, is strongly asymmetric at radio wavelengths, much brighter in the north, but bilaterally symmetric in X-rays. We present the results of X-ray expansion measurements that illuminate the origin of the radio asymmetry. We confirm the mean expansion rate (2011-2015) of 0.58% per yr, but large spatial variations are present. Using the nonparametric 'Demons' method, we measure the velocity field throughout the entire SNR, finding that motions vary by a factor of 5, from 0.''09 to 0.''44 per yr. The slowest shocks are at the outer boundary of the bright northern radio rim, with velocities v(sub s) as low as 3600 km per sec (for an assumed distance of 8.5 kpc), much less than v(sub s) = 12,000-13,000 km per sec along the X-ray-bright major axis. Such strong deceleration of the northern blast wave most likely arises from the collision of SN ejecta with a much denser than average ambient medium there. This asymmetric ambient medium naturally explains the radio asymmetry. In several locations, significant morphological changes and strongly nonradial motions are apparent. The spatially integrated X-ray flux continues to increase with time. Based on Chandra observations spanning 8.3 yr, we measure its increase at 1.3% +/- 0.8% per yr. The SN ejecta are likely colliding with the asymmetric circumstellar medium ejected by the SN progenitor prior to its explosion.

ISM: individual objects (G1.9+0.3) ISM: individual↗
NASA NTRS2020

Remote monitoring of agricultural systems using NDVI time series and machine learning methods: a tool for an adaptive agricultural policy

This study aims to provide accurate information about changes in agricultural systems (AS) using phenological metrics derived from the NDVI time series. Use of such information could help land managers optimize land use choices and monitor the status of agricultural lands, under a variety of environmental and socioeconomic conditions. For this purpose, the Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI data were used to derive phenological metrics over the Oum Er-Rbia basin (central Morocco). Random forest (RF), support vector machine (SVM), and K-nearest neighbor (KNN) classifiers were explored and compared on their ability to classify AS classes over the study area. Four main AS classes have been considered: (1) irrigated annual crop (IAC), (2) irrigated perennial crop (IPC), (3) rainfed area (RA), and (4) fallow (FA). By comparing the accuracy of the three classifiers, the RF method showed the best performance with an overall accuracy of 0.97 and kappa coefficient of 0.96.The RF method was then chosen to examine time variations in AS over a 16-year period (2000–2016). The AS main variations were detected and evaluated for the four AS classes. These variations have been found to be linked well with other indicators of local agricultural land management, as well as the historical agricultural drought changes over the study area. Overall, the results present a tool for decision makers to improve agricultural management and provide a different perspective in understanding the spatiotemporal dynamics of agricultural systems.

Youssef Lebrini↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of A Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
NASA NTRSDate not supplied

A Physics-Based Work-Energy Formulation for Real-Time Trajectory Guidance of a Lunar Lander

Throughout the years, many researchers have calculated and optimized trajectory solutions for lunar landing systems by employing sophisticated mathematical methods, that include: Hamilton’s Principle of Variation, Pontryagin’s maximum principle, and well known convex-optimization techniques among others. Many of these approaches typically require expensive computational resources to achieve convergence in the solution. In an effort to reduce complexity and the computational load required to obtain real-time guidance commands, a simple physics-based work-energy approach has been formulated. This approach is based on the dissipation of the mechanical energy of the vehicle to its final desired energy state required to achieve a safe landing. The rocket engine(s) employed during landing (among other maneuvers) dissipates mechanical energy by both doing work against the velocity vector of the vehicle (thus defining the trajectory path), and by jettisoning mass. Therefore, by solving the energy dissipation problem at every step of the maneuver, a much simpler formulation that naturally and quickly attains convergence is obtained. This formulation is not limited to approach, landing, and divert maneuvers, but in principle it can be employed during de-orbiting, braking burn, ascent, as well as orbit insertion.

Guidance↗
DOE OSTI2021

Comparison of temperature adaptive calibration methods for laser induced fluorescence based fuel-in-oil instrument

Fuel dilution of engine oil (or fuel-in-oil, FiO) is an important issue as multiple and late-cycle fuel injection, integral to many combustion efficiency and emissions improvements (e.g., downsized boosted gasoline engines and catalyst thermal management) increases FiO rate. In addition to causing general wear and corrosion in engine due to decreased oil viscosity and pH buffering, FiO is also believed to cause destructive low speed pre-ignition (or super knock) in boosted SI engines. To understand the effects of engine operating conditions on the FiO rate, an optical diagnostic capable of measuring transient FiO on minute timescales was recently developed and demonstrated on a modified GM Ecotech engine system (Neupane et al., Applied Spectroscopy 2021). The measurement is based on adding a dye to the fuel and monitoring for its presence in oil via laser-induced fluorescence (LIF). The measured LIF signal is related to FiO concentration via pre-determined calibration factors using a multivariate classical least square (CLS) method.Since fluorescence quantum yield is a function of temperature, measured LIF intensity not only depends on FiO concentration but also oil temperature. To expand the applicability of the FiO diagnostic to transient oil-temperature conditions (e.g., cold start in practical engines), this study develops a method to account for oil-temperature variations. The effect of oil temperature (20°C - 95°C) on the LIF spectra of eight FiO samples ranging from ~0.8-15% was investigated. LIF intensities of the FiO samples decreased linearly with increasing temperature; the reductions being more significant at dye peaks. We develop a new calibration model (T-adaptive CLS) incorporating the temperature (T) effects on LIF intensity that enables simultaneous calculation of FiO and oil temperature. The improved FiO diagnostic with T-adaptive calibration is more robust, and applicable to varying oil-temperature conditions. For example, when strategies such as multiple/late fuel injections are applied to overcome cold-start instability due to use of low vaporization bio-based fuels such as ethanol, FiO rate is expected to be very high; the improved T-adaptive FiO diagnostic is hence relevant for engine- and fuel-system calibration and optimization. The diagnostic could also provide validation data for flow-field and spray interaction CFD models, further broadening the diagnostic’s utility for advancing engine technology and efficiency.

Neupane, Sneha↗
DOE OSTI2020

The German Site Selection Process - Implementation and Challenges - 20528

Federal Company for Radioactive Waste Disposal (Bundesgesellschaft fuer Endlagerung mbH; hereinafter referred to as BGE) is a state-owned company established in September 2016 following the 'Act on the Organisational Restructuring in the Field of Radioactive Waste Storage'. 'Repository Site Selection Act' (hereinafter referred to as StandAG) defines the site selection procedure, aiming to ensure the best possible safety for storing high-level radioactive waste in Germany for the duration of one million years. The selection procedure is a participatory, transparent, learning and self-questioning process based on scientific expertise. Possible repository sites are sought all over Germany, without excluding or predetermining any potential site at the beginning. The procedure consists of three phases. After each phase, the legislative decides on the regions for further consideration in the next phase of the site selection procedure. The current first phase (Phase I) is divided into two steps. Objective of Step 1 is to determine subareas of favorable geological conditions based on available data. Step 2 of Phase I will lead to identification of specific regions that are considered for further surface-based geoscientific and geophysical exploration. Potential host rocks are rock salt, clay rock and crystalline rock. As part of the site selection procedure, an overview of the status of the planning and implementation of R and D activities is given. Topics of particular interest are, for example, the adaptation of the maximum temperature of 100 deg. C at the container surface, the impact of long-term climate variations as well as minimally invasive exploration methods. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Automated Waterbox Inspection for Nuclear Power Plants Using Computer Vision - Based Change Detection

Nuclear power plant waterboxes require regular inspection for leaks, missing components, and structural damage during maintenance outages. Traditional manual inspection is time-consuming and poses safety risks from confined space entry. We developed an automated computer vision system for drone-based waterbox inspection in partnership with Florida Light and Power. Our approach uses feature detection and matching to identify critical changes between baseline and current inspection images, automatically flagging additions (leaks/debris), removals (missing plugs), and translations (displaced components) while compensating for drone movement and environmental variations. We systematically evaluated six feature matching methods, from classical approaches (SIFT+BF) to state-of-the-art neural networks (SuperPoint+SuperGlue), using both standard benchmarks (HPatches) and waterbox-specific validation with real-world augmentations. SuperPoint+SuperGlue achieved superior performance with 7.82 pixels RMSE and 100% success rate—2.8x better accuracy than our baseline. While the pre-trained model has commercial licensing restrictions for nuclear deployment, our findings validate this architecture for custom training. We implemented a real-time GUI demonstrating the SIFT+BF approach for immediate deployment, processing drone feeds at 30 FPS with color-coded change visualization. Future work includes training a custom SuperPoint+SuperGlue model on waterbox data and integrating Vision-Language Models for automated reporting and maintenance guidance.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Live 129I-129Xe dating

A technique of cosmic ray exposure age dating using cosmic ray produced I-129 and Xe-129 components is discussed. The live I-129 - Xe-129 method provides an ideal monitor for cosmic ray flux variations on the 10(7)y - 10(8)y time-scale. It is based on low-energy neutron reactions on Te, and these data, when coupled to those from other methods, may facilitate the detection of complex exposure histories.

Marti, K.↗
NASA NTRS1986

On the flight dynamics of aeroelastic vehicles

The nonlinear equations of motion for an elasic airplane are developed from first principles. Lagrange's Equation and the Principle of Virtual Work are used to generate the equations of motion and strip theory is then employed to obtain closed form integral expressions for the generalized forces. The inertial coupling is minimized by appropriate choice of the body reference axes and by making use of free vibration modes of the body. In addition, particular attention is paid to the simplifying assumptions used during the development of the equations of motion. A unique aspect of this modeling process is that since the generalized aerodynamic forces are determined from closed form, analytic expressions, this method can be used to gain insight into the effects of parameter variations not easily obtained from numerical models. A numerical example is also presented in which the modeling method is applied to a representative elastic aircraft. The model is used to address the effects of aerodynamic coupling which occurs between the rigid body degrees of freedom and the elastic degrees of freedom. Finally, model simplification is addressed and two methods are evaluated. The resulting frequency responses are compared.

Schmidt, D. K.↗
NASA NTRS1985

Variation in efficiency of parallel algorithms

The present study has the objective to investigate some iterative parallel-processor linear equation solving algorithms with respect to efficiency for analyses of typical linear engineering systems. Attention is given to a set of n linear equations, Ku = p, where K = an n x n positive definite, sparsely populated, symmetric matrix, u = an n x 1 vector of unknown responses, and p = an n x 1 vector of prescribed constants. This study is concerned with a hybrid method in which iteration is used to solve the problem, while a direct method is used on the local processor level. Variations in the efficiency of parallel algorithms are explored. Measures of the efficiency are based on computer experiments regarding the algorithms. For all the algorithms, the wall clock time is found to decrease as the number of processors increases.

Hayashi, A.↗