Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “second moment method”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1923

General biplane theory

This report deals with the air forces on biplane cellule. The first part of the report deals with the two-dimensional problem neglecting viscosity. The variation of the section, chord, gap, stagger, and decalage are investigated, a great number of examples are calculated, and all numerical results are given in tables. For the biplane without stagger it is found that the loss of lift in consequence of the mutual influence of the two wing sections is only half as much if the lift is produced by the curvature of the section as it is when the lift is produced by the inclination of the chord to the direction of motion. The second part deals with the influence of the lateral dimensions. It is found that the loss of lift due to induction is almost unchanged, whether the biplane is staggered or not. In the third part conclusions from previous investigations are drawn, viscosity and experimental experience are brought in, and the method is simplified for practical application. Simple formulas give the drag, lift, and moment. In order to make use of the simple formulas more convenient, tables for the dynamical pressure, induced drag, and angle of attack are added so that practically no computation is needed for the application of the results.

Munk, Max M↗
NASA NTRS1973

Mars gravity derived from the long-period motion of Mariner 9

The spherical harmonic coefficients in the Mars gravity field are determined using a two-stage filtering process. In the first stage, a weighted least-squares orbit determination processor is applied to Mariner 9 radio data to estimate the mean Kepler elements of each orbit. A finite difference method is then used to generate orbital element rates. In the second stage, the gravity coefficients are determined using another weighted least-squares processor which fits the long period Lagrange perturbation equations to the estimated rates. The resulting gravity potential and associated equipotential surface is presented and error sources discussed. Gravity field results are applied to various physical properties of Mars, such as moments of inertia.-

Ferrari, A. J.↗
NASA NTRS2008

Comparing Local Starbursts to High-Redshift Galaxies: A Search for Lyman-Break Analogs

We compare the restframe far-ultraviolet (FUV) morphologies of 8 nearby interacting and starburst galaxies (Arp 269, M 82, Mrk 08, NGC 0520, NGC 1068, NGC 3079, NGC 3310, NGC 7673) with 54 galaxies at z approx.1.5 and 46 galaxies at z approx.4 in the Great Observatories Origins Deep Survey (GOODS) images taken with the Advanced Camera for Surveys onboard the Hubble Space Telescope. We calculate the Gini coefficient (G), the second order moment of 20% of the brightest pixels (M20), and the S ersic index (n). We find that 20% (11/54) of z approx.1.5 and 37% (17/46) of z approx.4 galaxies are bulge-like, using G and M20. We also find approx.70% of the z approx.1.5 and z approx.4 galaxies have exponential disks with n > 0.8. The 2D profile combined with the nonparametric methods provides more detail, concerning the nature of disturbed systems, such as merger and post-merger types. We also provide qualitative descriptions of each galaxy system and at each redshift. We conclude that Mrk 08, NGC 3079, and NGC 7673 have similar morphologies as the starburst FUV restframe galaxies and Lyman-break galaxies at z approx.1.5 and 4, and determine that they are Lyman-break analogs.

Petty, Sara M.↗
NASA NTRS2013

Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot

The optical extinction of a laser due to scattering of particles is a well-known phenomenon. In a laboratory environment, this physical principle is known as the Beer-Lambert law, and is often used to measure the concentration of scattering particles in a fluid or gas. This method has been experimentally shown to be a usable means to measure the dust density from a rocket plume interaction with the lunar surface. Using the same principles and experimental arrangement, this technique can be applied to hydrometeor size distributions, and for launch-pad operations, specifically as a passive hail detection and measurement system. Calibration of a hail monitoring system is a difficult process. In the past, it has required comparison to another means of measuring hydrometeor size and density. Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through hail (or dust in the rocket case) yields an estimate of the second moment of the particle cloud, and hydrometeor size distribution in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rain and hail fall make indirect measurements of the drop-size distribution. Instruments that scatter microwaves off of hydrometeors, such as the WSR-88D (Weather Surveillance Radar 88 Doppler), vertical wind profilers, and microwave disdrometers, measure the sixth moment of the drop size distribution (DSD). By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain and hail yield a measurement of the DSD's second moment by way of the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required. Depending on the intensity of the hail fall rate for moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. For hail fall only, the laser path may be shorter because of greater scattering due to the properties of hailstones versus raindrops. A photodetector may replace the camera in automated installations. Laser-based rain and hail measurement systems are available, but they are based on measuring the interruption of a thin laser beam, thus counting individual hydrometeors. These systems are true disdrometers since they also measure size and velocity. The method reported here is a simple method, requiring far less processing, but it is not a disdrometer.

Lane, John↗
NASA NTRS1975

Perform a gyro test of general relativity in a satellite and develop associated control technology

The progress accomplished in the Stanford Gyro Relativity program during the period November 1974 to October 1975 was described. Gyro developments were continued in the main laboratory dewar, concentrating on the operation of a three axis gyro readout and on improvements to the methods of canceling trapped fields in the rotor; these efforts culminated in the first successful observation of the London moment in the spinning gyro rotor in March 1975. Following a review meeting at that time, a new goal was formulated for the next 12 to 18 months, namely to operate a gyroscope in the new ultra-low field facility with readout resolution approaching 1 arc-second. The following other tasks were also completed: (1) sputtering work, (2) magnetometry, (3) construction and installation of the North Star simulator, (4) analysis of torques on the gyro, especially in inclined orbits, (5) equivalence principle accelerometer, and (6) analysis of a twin-satellite test of relativity.

Fairbank, W. M.↗
NASA NTRS1975

Aeroelastic stability of trimmed helicopter blades in forward flight

Equations for moderately large amplitude coupled flap-lag motion of a torsionally rigid hingeless elastic helicopter blade in forward flight are derived. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. By using Galerkin's method the spatial dependence of the problem is eliminated and the equations are linearized about a time dependent equilibrium position determined from the trimmed equilibrium position of the rotor in forward flight. In the first trim procedure the rotor is maintained at a fixed value of thrust coefficient with forward flight and horizontal and vertical force equilibrium is satisfied in addition to maintaining zero pitch and roll moments. The second trim procedure maintains only zero pitch and roll moment simulating conditions under which a rotor would be tested in the wind tunnel.

Friedmann, P.↗
DOE OSTI2025

Partnership Center for High-Fidelity Boundary Plasma Simulation (Final Report)

Within the Partnership Center for High-Fidelity Boundary Plasma Simulation (HBPS), work at UT-Austin was aimed at improved verification, validation, and uncertainty quantification (VVUQ) for edge plasma simulations and on performing gyrokinetics simulations of pedestal instabilities and turbulence in order to expand foundational understanding of pedestal transport. Regarding VVUQ, the accomplishments can be summarized as follows. First, it was shown that the Moment Preserving Constrained Resampling technique, when applied periodically in particle-in-cell simulations in the XGC code, can dramatically improve the accuracy of the simulation at essentially equivalent computational cost. Second, a technique for estimating model correlations, which are required to solve the model selection and sample allocation problem in multifidelity UQ techniques, without sampling the highest fidelity, most computationally expensive model, was developed and demonstrated. Third, previously developed methods for estimating statistical and discretization errors were applied to numerical methods relevant to edge plasma simulations, namely in particle-in-cell-based approaches, and shown to work. Finally, benchmark studies for comparing gyrokinetic codes were developed and performed, leading to reasonable agreement between four commonly used codes. Regarding physics studies, gyrokinetic simulations to investigate microtearing modes in the DIII-D pedestal were performed using the GENE code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Subsonic panel method for the efficient analysis of multiple geometry perturbations

An accurate and efficient method was developed for the aerodynamic analysis of a series of arbitrary small geometry perturbations to a given baseline configuration. The method is appropriate for wing-fuselage configurations in incompressible potential flow. Mathematical formulations are presented for three computer programs that are employed. The first program is a conventional surface panel method for calculating the baseline singularity distribution. The second program calculates a partial derivative matrix. Each element of the matrix is the rate of change of singularity strength at one point with respect to a surface coordinate of a different point. For each baseline configuration, the calculated quantities from the first two programs establish an input file for the third. The third program calculates the surface pressure distribution and forces and moments for a series of geometry perturbations.

Bristow, D. R.↗
NASA NTRS2004

Employing Sensitivity Derivatives for Robust Optimization under Uncertainty in CFD

A robust optimization is demonstrated on a two-dimensional inviscid airfoil problem in subsonic flow. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables), an approximate first-order statistical moment method is employed to represent the Computational Fluid Dynamics (CFD) code outputs as expected values with variances. These output quantities are used to form the objective function and constraints. The constraints are cast in probabilistic terms; that is, the probability that a constraint is satisfied is greater than or equal to some desired target probability. Gradient-based robust optimization of this stochastic problem is accomplished through use of both first and second-order sensitivity derivatives. For each robust optimization, the effect of increasing both input standard deviations and target probability of constraint satisfaction are demonstrated. This method provides a means for incorporating uncertainty when considering small deviations from input mean values.

Newman, Perry A.↗
NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗
NASA NTRS1981

Some applications of aerodynamic formulations to problems in aircraft dynamics

Two applications of mathematical modeling to aerodynamic problems are cussed. The first application is an investigation of the capacity of a nonlinear aerodynamic mathematical model to describe the aerodynamic reactions on an airfoil with a deflecting flap in transonic flow. Flow field computational methods are used to evaluate the nonlinear, unsteady aerodynamic data in terms of characteristic motions called for by the model. Histories of unconstrained motions of the flap are generated from the flap equations of motion, with the aerodynamic reactions specified by the mathematical model. In the second application wing rock is investigated. The most recent model accommodates experimental results wing rock by admitting the existence of aerodynamic hysteresis in the variation of the steady state rolling moment coefficient with roll angle is described.

Schiff, L. B.↗
NASA NTRS1961

Estimation of Rotary Stability Derivatives at Subsonic and Transonic Speeds

The first part of this paper pertains to the estimation of subsonic rotary stability derivatives of wings. The unsteady potential flow problem is solved by a superposition of steady flow solutions. Numerical results for the damping coefficients of triangular wings are presented as functions of aspect ratio and Mach number, and are compared with experimental results over the Mach number range 0 to 1. In the second part, experimental results are used. to point out a close correlation between the nonlinear variations with angle of attack of the static pitching-moment curve slope and the damping-in-pitch coefficient. The underlying basis for the correlation is found as a result of an analysis in which the indicial function concept and. the principle of super-position are adapted to apply to the nonlinear problem. The form of the result suggests a method of estimating nonlinear damping coefficients from results of static wind-tunnel measurements.

Tobak, Murray↗
NASA NTRS2001

Free-form Airfoil Shape Optimization Under Uncertainty Using Maximum Expected Value and Second-order Second-moment Strategies

Free-form shape optimization of airfoils poses unexpected difficulties. Practical experience has indicated that a deterministic optimization for discrete operating conditions can result in dramatically inferior performance when the actual operating conditions are different from the - somewhat arbitrary - design values used for the optimization. Extensions to multi-point optimization have proven unable to adequately remedy this problem of "localized optimization" near the sampled operating conditions. This paper presents an intrinsically statistical approach and demonstrates how the shortcomings of multi-point optimization with respect to "localized optimization" can be overcome. The practical examples also reveal how the relative likelihood of each of the operating conditions is automatically taken into consideration during the optimization process. This is a key advantage over the use of multipoint methods.

Huyse, Luc↗
NASA NTRS1990

Two-dimensional shape recognition using sparse distributed memory

Researchers propose a method for recognizing two-dimensional shapes (hand-drawn characters, for example) with an associative memory. The method consists of two stages: first, the image is preprocessed to extract tangents to the contour of the shape; second, the set of tangents is converted to a long bit string for recognition with sparse distributed memory (SDM). SDM provides a simple, massively parallel architecture for an associative memory. Long bit vectors (256 to 1000 bits, for example) serve as both data and addresses to the memory, and patterns are grouped or classified according to similarity in Hamming distance. At the moment, tangents are extracted in a simple manner by progressively blurring the image and then using a Canny-type edge detector (Canny, 1986) to find edges at each stage of blurring. This results in a grid of tangents. While the technique used for obtaining the tangents is at present rather ad hoc, researchers plan to adopt an existing framework for extracting edge orientation information over a variety of resolutions, such as suggested by Watson (1987, 1983), Marr and Hildreth (1980), or Canny (1986).

Kanerva, Pentti↗
NASA NTRS1992

Modeling methods for high-fidelity rotorcraft flight mechanics simulation

The cooperative effort being carried out under the agreements of the United States-Israel Memorandum of Understanding is discussed. Two different models of the AH-64 Apache Helicopter, which may differ in their approach to modeling the main rotor, are presented. The first model, the Blade Element Model for the Apache (BEMAP), was developed at Ames Research Center, and is the only model of the Apache to employ a direct blade element approach to calculating the coupled flap-lag motion of the blades and the rotor force and moment. The second model was developed at the Technion-Israel Institute of Technology and uses an harmonic approach to analyze the rotor. The approach allows two different levels of approximation, ranging from the 'first harmonic' (similar to a tip-path-plane model) to 'complete high harmonics' (comparable to a blade element approach). The development of the two models is outlined and the two are compared using available flight test data.

Mansur, M. Hossein↗
DOE OSTI2025

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

A unified perturbation expansion for surface scattering

Starting with the extinction theorem, a perturbation expansion which, to first and second orders, converges over a wider domain than the small perturbation expansion and the momentum transfer expansion is presented. It is shown that, in the appropriate limits, both of these theories, as well as the two-scale expansion, are recovered. There is no adjustable parameter, such as a spectral split, in the theory. This theory is applied to random rough surfaces and derive analytic expressions for the coherent field and the bistatic cross section. Finally, a numerical test of the theory against method of moments results for Gaussian random rough surfaces with a power law spectrum is given. These results show that the expansion is ramarkably accurate over a large range of surface heights and slopes for both horizontal and vertical polarization.

Rodriguez, Ernesto↗
DOE OSTI2026

Analytic Chromaticity Formulas for the Muon g -2 experiment at fermilab

The Muon g-2 Experiment (E989) at Fermilab measured the muon anomalous magnetic moment aμ with unprecedented precision of 127 ppb using a storage ring. The experiment pushed the limits in terms of accuracy and systematic errors. Achieving the 78-ppb total systematic uncertainty required precise beam dynamics modeling and corrections for effects on the measured muon precession frequency. Here, we derived the first analytic aberration formulas up to the second order for the muon g-2 storage ring’s combined-function electrostatic quadrupoles (superimposed magnetic dipole and electric quadrupole fields) using an order-by-order perturbation method. From these, we obtained the exact chromaticity formulas for three ring models of different granularity and validated them against numerical calculations using COSY INFINITY, achieving analytic-numerical agreement to ⁠$\mathscr{O}(10^{-10})$. This work resolved discrepancies between previous approximate derivations and provided essential beam dynamics results for Runs 4-6 analyses. We also calculated nonlinear chromaticities up to ninth order. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous muon g-2 experiment at Brookhaven National Laboratory, with the final result announced in June 2025.

aberrations↗