A lyapunov method for the estimation of statistical averages
Liapunov method for estimation of statistical averages - application to control theory and theory of random vibration
CONTROL SYSTEM↗
SEARCH · Search NASA
Results for “method of averaging”
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.
Minimum-fuel, power-limited transfers between coplanar elliptical orbits
Two techniques of computing coplanar, minimum-fuel, power-limited transfers are presented, based on approximate solutions obtained by the averaging method. For the first technique, the average solution provides estimates of the initial adjoint variables, and the second provides approximations of the optimal controls in feedback form. The accuracy of the techniques for computing coplanar, minimum-fuel, power-limited transfers is evaluated for a number of initial and final orbits.
Mease, Kenneth D.↗
Improved uranium isotopic ratio determinations for the liquid sampling-atmospheric pressure glow discharge orbitrap mass spectrometer by use of moving average processing
The liquid sampling-atmospheric pressure glow discharge (LS-APGD) is a versatile combined atomic and molecular (CAM) ionization source capable of ionizing elemental species, small polar compounds, low-polarity polycyclic aromatic hydrocarbons, and proteins. While the LS-APGD has been proven capable of determining the 235U/238U isotope ratio in enriched and natural uranium, recent efforts have strived to reduce the analysis time of individual measurements by employing higher-order data processing techniques, specifically, moving average methods. Additionally, the use of a moving average (of data windows of various widths) improves the precision of the ratio measurements and reduces the number of scans needed to be collected to generate high-precision results. Additionally, use of the moving average minimizes the quantity of sample that needs to be analyzed while reducing the measurement time. These reductions have allowed isotope ratio (IR) determinations to be made using injections instead of direct infusion of the analyte. For example, employing a 40-point window width in the moving average to direct infusion data reduced the percent relative standard deviation (% RSD) of the 235U/238U values from 3.4% to 0.5%. For the injection of a 100 μL aliquot, the % RSD was reduced from 5.3% to 0.7%.
74 ATOMIC AND MOLECULAR PHYSICS↗
A comparison of the dielectric function for p-point water models
We compute the longitudinal dielectric function ɛl(k) for p-point water models, p = 3 (SPC/E, TIP3P), p = 4 (TIP4P-EW, OPC), and p = 5 (TIP5P-E) from the charge–charge fluctuation function SZZ(k), obtained from two methods: average of the charge in k-space (method:sum) and from the pair distribution function in Fourier space (method:integral). The latter requires a continuation to small k-values, which we thoroughly discuss. We conclude with a detailed comparison of the longitudinal dielectric function for the different p-point models and its main characteristics: small k and large k limits, position of the poles, and its (complex) zeros.
Shevy, Noah [Ames National Laboratory, Department ↗
Energy-dependent Orbital Modulation of X-rays and Constraints on Emission of the Jet in Cyg X-3
We study orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and RXTE. Using the wealth of the presently available data and an improved averaging method, we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. We find that above ∼5 keV, the modulation depth decreases with the increasing energy, which is consistent with the modulation being caused by both bound-free absorption and Compton scattering in the stellar wind of the donor, with minima corresponding to the highest optical depth, which occurs around the superior conjunction. We find a decrease of the depth below ∼3 keV, which appears to be due to re-emission of the absorbed continuum by the wind in soft X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet in Cyg X-3, which emits γ-rays detected at energies > 0.1 GeV in soft spectral states, is found to be minor up to ∼100 keV. This implies the presence of a rather sharp low-energy break in the jet MeV-range spectrum.We also calculate phase-resolved RXTE X-ray spectra, and show the difference between the spectra corresponding to phases around the superior and inferior conjunctions can indeed be accounted for by a combined effect of bound-free absorption in an ionized medium and Compton scattering.
Energy-dependent Orbital↗
Overview of recent turbulence studies across multiple confinement modes at the ASDEX Upgrade tokamak using the Correlation Electron Cyclotron Emission diagnostic
This work presents an overview of recent and ongoing experimental measurements of core and edge turbulence across multiple confinement regimes using the Correlation Electron Cyclotron Emission (CECE) diagnostic at the ASDEX Upgrade (AUG) tokamak. A common goal among these investigations is to identify how the properties of the turbulent electron temperature fluctuations measured by CECE influence and regulate the unique transport characteristics of each confinement regime, including L-mode, I-mode, ELMy H-mode, and ELM-free H-mode. Optics and signal processing methods to aid in the analysis and interpretation of experimental turbulence results are also presented. These methods, and particularly the down-sampling and ensemble averaging method, are relevant to a wide variety of fusion and non-fusion applications.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of a Rotating Aerodynamic Probe on the Flow Field of a Compressor Rotor
An investigation of distortions of the rotor exit flow field caused by an aerodynamic probe mounted in the rotor is described in this paper. A rotor total pressure Kiel probe, mounted on the rotor hub and extending up to the mid-span radius of a rotor blade channel, generates a wake that forms additional flow blockage. Three types of high-response aerodynamic probes were used to investigate the distorted flow field behind the rotor. These probes were: a split-fiber thermo-anemometric probe to measure velocity and flow direction, a total pressure probe, and a disk probe for in-flow static pressure measurement. The signals acquired from these high-response probes were reduced using an ensemble averaging method based on a once per rotor revolution signal. The rotor ensemble averages were combined to construct contour plots for each rotor channel of the rotor tested. In order to quantify the rotor probe effects, the contour plots for each individual rotor blade passage were averaged into a single value. The distribution of these average values along the rotor circumference is a measure of changes in the rotor exit flow field due to the presence of a probe in the rotor. These distributions were generated for axial flow velocity and for static pressure.
Lepicovsky, Jan↗
Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems
An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.
54 ENVIRONMENTAL SCIENCES↗
Micromechanical models for textile structural composites
The objective is to develop micromechanical models for predicting the stiffness and strength properties of textile composite materials. Two models are presented to predict the homogeneous elastic constants and coefficients of thermal expansion of a textile composite. The first model is based on rigorous finite element analysis of the textile composite unit-cell. Periodic boundary conditions are enforced between opposite faces of the unit-cell to simulate deformations accurately. The second model implements the selective averaging method (SAM), which is based on a judicious combination of stiffness and compliance averaging. For thin textile composites, both models can predict the plate stiffness coefficients and plate thermal coefficients. The finite element procedure is extended to compute the thermal residual microstresses, and to estimate the initial failure envelope for textile composites.
Marrey, Ramesh V.↗
A note on obtaining an approximate value for the period of a periodic solution of x-double prime = -V prime /x/
The first integral of x-double prime = - V prime (x) yields an integral for the period of a periodic solution, if such exists. In general, this integral cannot be evaluated. By means of an approximate solution along with the minimization of a mean-square error, one can obtain an approximate value for the period in terms of the amplitude of the motion. The calculated period agrees very well with the period obtained by means of numerical integration for the case of orbit-orbit resonance involving the motion of two satellites of a planet. The same method is applied to obtain an alternative derivation of the first Krylov-Bogoliuboff averaging method in nonlinear mechanics.
Lass, H.↗
Micromechanical models for textile composites
Numerical and analytical micromechanical models are presented to predict the thermoelastic behavior of a textile composite. In the numerical model, the unit-cell is discretized with finite elements, and periodic boundary conditions are imposed between opposite faces of the unit-cell. For a thin textile composite, stress gradients effects through the thickness are demonstrated. The consequent difference in the stiffness and strength behavior of thick and thin composites are discussed. The numerical model is implemented to predict 3-D thermo-elastic constants for a thick textile composite, and the plate thermo-mechanical properties for a thin textile composite. The numerical model is extended to compute the thermal residual microstresses due to processing to predict the composite failure envelopes. An analytical model - Selective Averaging Method (SAM) - is proposed, which is based on a judicious combination of stiffness and compliance averaging to estimate the 3-D elastic constants. Both the models are tested and verified for several examples by comparing the stiffness properties with elasticity solutions and available results.
Sankar, Bhavani V.↗
MIRO Continuum Calibration for Asteroid Mode
MIRO (Microwave Instrument for the Rosetta Orbiter) is a lightweight, uncooled, dual-frequency heterodyne radiometer. The MIRO encountered asteroid Steins in 2008, and during the flyby, MIRO used the Asteroid Mode to measure the emission spectrum of Steins. The Asteroid Mode is one of the seven modes of the MIRO operation, and is designed to increase the length of time that a spectral line is in the MIRO pass-band during a flyby of an object. This software is used to calibrate the continuum measurement of Steins emission power during the asteroid flyby. The MIRO raw measurement data need to be calibrated in order to obtain physically meaningful data. This software calibrates the MIRO raw measurements in digital units to the brightness temperature in Kelvin. The software uses two calibration sequences that are included in the Asteroid Mode. One sequence is at the beginning of the mode, and the other at the end. The first six frames contain the measurement of a cold calibration target, while the last six frames measure a warm calibration target. The targets have known temperatures and are used to provide reference power and gain, which can be used to convert MIRO measurements into brightness temperature. The software was developed to calibrate MIRO continuum measurements from Asteroid Mode. The software determines the relationship between the raw digital unit measured by MIRO and the equivalent brightness temperature by analyzing data from calibration frames. The found relationship is applied to non-calibration frames, which are the measurements of an object of interest such as asteroids and other planetary objects that MIRO encounters during its operation. This software characterizes the gain fluctuations statistically and determines which method to estimate gain between calibration frames. For example, if the fluctuation is lower than a statistically significant level, the averaging method is used to estimate the gain between the calibration frames. If the fluctuation is found to be statistically significant, a linear interpolation of gain and reference power is used to estimate the gain between the calibration frames.
Lee, Seungwon↗
Analytical Approach Validation for the Spin-Stabilized Satellite Attitude
An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques and the satellite in an elliptical orbit. Assuming a quadripole model for the Earth s magnetic field, an analytical averaging method is applied to obtain the mean residual torque in every orbital period. The orbit mean anomaly is used to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that the residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. The theory developed has been applied to the Brazilian s spin stabilized satellites, which are quite appropriated for verification and comparison of the theory with the data generated and processed by the Satellite Control Center of Brazil National Research Institute. The results show the period that the analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of Brazil National Research Institute.
Zanardi, Maria Cecilia F. P. S.↗
Clumping index estimation with 30°-tilted cameras in row crops: Evaluation of methods and segment size effects
The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Furthermore, applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.
30°-tilted camera↗
Poisson equations of rotational motion for a rigid triaxial body with application to a tumbling artificial satellite
A mathematical model is developed for studying the effects of gravity gradient torque on the attitude stability of a tumbling triaxial rigid satellite. Poisson equations are used to investigate the rotation of the satellite (which is in elliptical orbit about an attracting point mass) about its center of mass. An averaging method is employed to obtain an intermediate set of differential equations for the nonresonant, secular behavior of the osculating elements which describe the rotational motions of the satellite, and the averaged equations are then integrated to obtain long-term secular solutions for the osculating elements.
Liu, J. J. F.↗
Current Trends in Modeling Research for Turbulent Aerodynamic Flows
The engineering tools of choice for the computation of practical engineering flows have begun to migrate from those based on the traditional Reynolds-averaged Navier-Stokes approach to methodologies capable, in theory if not in practice, of accurately predicting some instantaneous scales of motion in the flow. The migration has largely been driven by both the success of Reynolds-averaged methods over a wide variety of flows as well as the inherent limitations of the method itself. Practitioners, emboldened by their ability to predict a wide-variety of statistically steady, equilibrium turbulent flows, have now turned their attention to flow control and non-equilibrium flows, that is, separation control. This review gives some current priorities in traditional Reynolds-averaged modeling research as well as some methodologies being applied to a new class of turbulent flow control problems.
Gatski, Thomas B.↗
Adaptive Inverse Control for Rotorcraft Vibration Reduction
This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.
Jacklin, Stephen A.↗
Tracking blobs in the turbulent edge plasma of a tokamak fusion device
Abstract The analysis of turbulence in plasmas is fundamental in fusion research. Despite extensive progress in theoretical modeling in the past 15 years, we still lack a complete and consistent understanding of turbulence in magnetic confinement devices, such as tokamaks. Experimental studies are challenging due to the diverse processes that drive the high-speed dynamics of turbulent phenomena. This work presents a novel application of motion tracking to identify and track turbulent filaments in fusion plasmas, called blobs, in a high-frequency video obtained from Gas Puff Imaging diagnostics. We compare four baseline methods (RAFT, Mask R-CNN, GMA, and Flow Walk) trained on synthetic data and then test on synthetic and real-world data obtained from plasmas in the Tokamak à Configuration Variable (TCV). The blob regime identified from an analysis of blob trajectories agrees with state-of-the-art conditional averaging methods for each of the baseline methods employed, giving confidence in the accuracy of these techniques. By making a dataset and benchmark publicly available, we aim to lower the entry barrier to tokamak plasma research, thereby greatly broadening the community of scientists and engineers who might apply their talents to this endeavor.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗