Search NASA⌕ Search

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.

At least 307 records · Page 17

Fast two-stream method for computing diurnal-mean actinic flux in vertically inhomogeneous atmospheres

Based on a derivation of the two-stream daytime-mean equations of radiative flux transfer, a method for computing the daytime-mean actinic fluxes in the absorbing and scattering vertically inhomogeneous atmosphere is suggested. The method applies direct daytime integration of the particular solutions of the two-stream approximations or the source functions. It is valid for any duration of period of averaging. The merit of the method is that the multiple scattering computation is carried out only once for the whole averaging period. It can be implemented with a number of widely used two-stream approximations. The method agrees with the results obtained with 200-point multiple scattering calculations. The method was also tested in runs with a 1-km cloud layer with optical depth of 10, as well as with aerosol background. Comparison of the results obtained for a cloud subdivided into 20 layers with those obtained for a one-layer cloud with the same optical parameters showed that direct integration of particular solutions possesses an 'analytical' accuracy. In the case of the source function interpolation, the actinic fluxes calculated above the one-layer and 20-layer clouds agreed within 1%-1.5%, while below the cloud they may differ up to 5% (in the worst case). The ways of enhancing the accuracy (in a 'two-stream sense') and computational efficiency of the method are discussed.

Filyushkin, V. V.↗

Status of Turbulence Modeling for Hypersonic Propulsion Flowpaths

This report provides an assessment of current turbulent flow calculation methods for hypersonic propulsion flowpaths, particularly the scramjet engine. Emphasis is placed on Reynolds-averaged Navier-Stokes (RANS) methods, but some discussion of newer meth- ods such as Large Eddy Simulation (LES) is also provided. The report is organized by considering technical issues throughout the scramjet-powered vehicle flowpath including laminar-to-turbulent boundary layer transition, shock wave / turbulent boundary layer interactions, scalar transport modeling (specifically the significance of turbulent Prandtl and Schmidt numbers) and compressible mixing. Unit problems are primarily used to conduct the assessment. In the combustor, results from calculations of a direct connect supersonic combustion experiment are also used to address the effects of turbulence model selection and in particular settings for the turbulent Prandtl and Schmidt numbers. It is concluded that RANS turbulence modeling shortfalls are still a major limitation to the accuracy of hypersonic propulsion simulations, whether considering individual components or an overall system. Newer methods such as LES-based techniques may be promising, but are not yet at a maturity to be used routinely by the hypersonic propulsion community. The need for fundamental experiments to provide data for turbulence model development and validation is discussed.

Georgiadis, Nicholas J.↗

Newmark-Type Pseudo-Three-Dimensional Back-Analysis of Co-Seismic Landslides in Egkremnoi, Lefkada, Greece

A pseudo-three-dimensional (pseudo-3D) methodology to back-analyze coseismic landslides was developed and applied to 68mapped landslides, which occurred over approximately 1 km2area in Egkremnoi, Greece, during the 2015 Mw6.5 Lefkada earthquake. The methodology is based on a one-dimensional (1D) Newmark-type sliding block model to assess instability and a spatial projection in 3D topography to derive landslide geometry. The strength parameters for modeled landslides that best match the landslide location, area, and volume were derived through an iterative scheme that optimizes the match using predefined matching criteria. The range of different-sized landslides produced shear strength estimates from ~10 to 300 kPa and led to the derivation of a regionally averaged strength envelope characterized by a cohesion of 6 kPa and a friction angle of 53°for the highly fractured limestones that are encountered in this area. Compared to previous full 3Dslope stability analyses in this area, the friction angle using this methodology was found to be generally consistent, but the cohesion was lower. The presented methodology can provide a computationally efficient method to estimate the average shear strength of a geologic unit over large areas, especially where extensive field and laboratory tests on the materials are unavailable or difficult to conduct.

Weibing Gong↗

Evaluation of the ERBE scene identification algo-rithm

The sensitivity of the radiation budget parameters and the scene selection process to different factors is evaluated. The use of ERB-7 CLE models provides instantaneous albedo and longwave flux values which are in essential agreement in the 70 deg cutoff case with those from the SAB method. The increase in albedo with increasing satellite zenith angle cutoff is not apparent at the target area level for different surface types. GOES models seem to show an increase in instantaneous albedo at the global level with satellite zenith angle of the same nature as the ERB-7 models investigated in a previous report and are, therefore, probably not an improvement. The use of recently derived NCLE models did not make a noticeable change in the budget parameters; but the cloud classification did show hemispherical differences and caused a day-night readjustment of cloudiness amounts. Preliminary results for an additional month (December 1979) indicate good agreement between the SAB and MLE methods. Additional work is required to establish the agreement for all seasons. Comparison of derived cloud amounts with other data sets such as THIR/TOMS indicates good zonal agreement. Sampling adequacy investigations at different temporal averaging intervals with the SAB method, indicate large uncertainties in small time average cases. For a 6-day averaging period, sampling is very poor with only 23% of the globe contributing toward the global mean albedo. Effect of modifications to the MLE procedure seem to have little effect on the derived budget quantities in an averaged sense. Significant differences in the flux values due to differences in scene selection are apparent in individual target areas studies. A modification of the MLE procedure to mimic the perpendicular bisector algorithm indicated no effect on the gross radiation budget quantities.

Vemury, S. K.↗

Detached Eddy Simulation Results for a Space Launch System Configuration at Liftoff Conditions and Comparison with Experiment

Computational simulations for a Space Launch System configuration at liftoff conditions for incidence angles from 0 to 90 degrees were conducted in order to generate integrated force and moment data and longitudinal lineloads. While the integrated force and moment coefficients can be obtained from wind tunnel testing, computational analyses are indispensable in obtaining the extensive amount of surface information required to generate proper lineloads. However, beyond an incidence angle of about 15 degrees, the effects of massive flow separation on the leeward pressure field is not well captured with state of the art Reynolds Averaged Navier-Stokes methods, necessitating the employment of a Detached Eddy Simulation method. Results from these simulations are compared to the liftoff force and moment database and surface pressure data derived from a test in the NASA Langley 14- by 22-Foot Subsonic Wind Tunnel.

Computational fluid dynamics↗

An investigation of gear mesh failure prediction techniques

A study was performed in which several gear failure prediction methods were investigated and applied to experimental data from a gear fatigue test apparatus. The primary objective was to provide a baseline understanding of the prediction methods and to evaluate their diagnostic capabilities. The methods investigated use the signal average in both the time and frequency domain to detect gear failure. Data from eleven gear fatigue tests were recorded at periodic time intervals as the gears were run from initiation to failure. Four major failure modes, consisting of heavy wear, tooth breakage, single pits, and distributed pitting were observed among the failed gears. Results show that the prediction methods were able to detect only those gear failures which involved heavy wear or distributed pitting. None of the methods could predict fatigue cracks, which resulted in tooth breakage, or single pits. It is suspected that the fatigue cracks were not detected because of limitations in data acquisition rather than in methodology. Additionally, the frequency response between the gear shaft and the transducer was found to significantly affect the vibration signal. The specific frequencies affected were filtered out of the signal average prior to application of the methods.

Zakrajsek, James J.↗

Dynamical Evolution of Meteoroid Streams, Developments Over the Last 30 Years

As soon as reliable methods for observationally determining the heliocentric orbits of meteoroids and hence the mean orbit of a meteoroid stream in the 1950s and 60s, astronomers strived to investigate the evolution of the orbit under the effects of gravitational perturbations from the planets. At first, the limitations in the capabilities of computers, both in terms of speed and memory, placed severe restrictions on what was possible to do. As a consequence, secular perturbation methods, where the perturbations are averaged over one orbit became the norm. The most popular of these is the Halphen- Goryachev method which was used extensively until the early 1980s. The main disadvantage of these methods lies in the fact that close encounter can be missed, however they remain useful for performing very long-term integrations. Direct integration methods determine the effects of the perturbing forces at many points on an orbit. This give a better picture of the orbital evolution of an individual meteoroid, but many meteoroids have to be integrated in order to obtain a realistic picture of the evolution of a meteoroid stream. The notion of generating a family of hypothetical meteoroids to represent a stream and directly integrate the motion of each was probably first used by Williams Murray & Hughes (1979), to investigate the Quadrantids. Because of computing limitations, only 10 test meteoroids were used. Only two years later, Hughes et. al. (1981) had increased the number of particles 20-fold to 200 while after a further year, Fox Williams and Hughes used 500 000 test meteoroids to model the Geminid stream. With such a number of meteoroids it was possible for the first time to produce a realistic cross-section of the stream on the ecliptic. From that point on there has been a continued increase in the number of meteoroids, the length of time over which integration is carried out and the frequency with which results can be plotted so that it is now possible to produce moving images of the stream. As a consequence, over recent years, emphasis has moved to considering stream formation and the role fragmentation plays in this.

Williams, I. P.↗

Probe shapes that measure time-averaged streamwise momentum and cross-stream turbulence intensity

A method and apparatus for directly measuring the time-averaged streamwise momentum in a turbulent stream use a probe which has total head response which varies as the cosine-squared of the angle of incidence. The probe has a nose with a slight indentation on its front face for providing the desired response. The method of making the probe incorporates unique design features. Another probe may be positioned in a side-by-side relationship to the first probe to provide a direct measurement of the total pressure. The difference between the two pressures yields the sum of the squares of the cross-stream components of the turbulence level.

Rossow, Vernon J.↗

Mapping Bone Mineral Density Obtained by Quantitative Computed Tomography to Bone Volume Fraction

Methods for relating or mapping estimates of volumetric Bone Mineral Density (vBMD) obtained by Quantitative Computed Tomography to Bone Volume Fraction (BVF) are outlined mathematically. The methods are based on definitions of bone properties, cited experimental studies and regression relations derived from them for trabecular bone in the proximal femur. Using an experimental range of values in the intertrochanteric region obtained from male and female human subjects, age 18 to 49, the BVF values calculated from four different methods were compared to the experimental average and numerical range. The BVF values computed from the conversion method used data from two sources. One source provided pre bed rest vBMD values in the intertrochanteric region from 24 bed rest subject who participated in a 70 day study. Another source contained preflight vBMD values from 18 astronauts who spent 4 to 6 months on the ISS. To aid the use of a mapping from BMD to BVF, the discussion includes how to formulate them for purpose of computational modeling. An application of the conversions would be used to aid in modeling of time varying changes in vBMD as it relates to changes in BVF via bone remodeling and/or modeling.

Pennline, James A.↗

Exploiting OMI NO2 Satellite Observations to Infer Fossil-fuel CO2 Emissions from U.S. Megacities

Fossil-fuel CO2 emissions and their trends in eight U.S. megacities during 2006–2017 are inferred by combining satellite-derived NOX emissions with bottom-up city-specific NOX-to-CO2 emission ratios. A statistical model is fit to a collection NO2 plumes observed from the Ozone Monitoring Instrument (OMI), and is used to calculate top-down NOX emissions. Decreases in OMI-derived NOX emissions are observed across the eight cities from 2006 to 2017 (−17% in Miami to −58% in Los Angeles), and are generally consistent with long-term trends of bottom-up inventories (−25% in Miami to −49% in Los Angeles), but there are some interannual discrepancies. City-specific NOX-to-CO2 emission ratios, used to calculate inferred CO2, are estimated through annual bottom-up inventories of NOX and CO2 emissions disaggregated to 1 × 1 km2 resolution. Over the study period, NOX-to-CO2 emission ratios have decreased by ~40% nationwide (−24% to −51% for our studied cities), which is attributed to a faster reduction in NOX when compared to CO2 due to policy regulations and fuel type shifts. Combining top-down NOX emissions and bottom-up NOX-to-CO2 emission ratios, annual fossil-fuel CO2 emissions are derived. Inferred OMI-based top-down CO2 emissions trends vary between +7% in Dallas to −31% in Phoenix. For 2017, we report annual fossil-fuel CO2 emissions to be: Los Angeles 113 ± 49 Tg/yr; New York City 144 ± 62 Tg/yr; and Chicago 55 ± 24 Tg/yr. A study in the Los Angeles area, using independent methods, reported a 2013–2016 average CO2 emissions rate of 104 Tg/yr and 120 Tg/yr, which suggests that the CO2 emissions from our method are in good agreement with other studies' top-down estimates. We anticipate future remote sensing instruments – with better spatial and temporal resolution – will better constrain the NOX-to-CO2 ratio and reduce the uncertainty in our method.

Ozone Monitoring Instrument (OMI)↗

The energy separation between the classical and nonclassical isomers of protonated acetylene - An extensive study in one- and n-particle space saturation

The energy separation between the classical and nonclassical forms of protonated acetylene has been reinvestigated in light of the recent experimentally deduced lower bound to this value of 6.0 kcal/mol. The objective of the present study is to use state-of-the-art ab initio quantum mechanical methods to establish this energy difference to within chemical accuracy (i.e., about 1 kcal/mol). The one-particle basis sets include up to g-type functions and the electron correlation methods include single and double excitation coupled-cluster (CCSD), the CCSD(T) extension, multireference configuration interaction, and the averaged coupled-pair functional methods. A correction for zero-point vibrational energies has also been included, yielding a best estimate for the energy difference between the classical and nonclassical forms of 3.7 + or - 1.3 kcal/mol.

Lindh, Roland↗

Investigation of dispersion-relation-preserving scheme and spectral analysis methods for acoustic waves

Important characteristics of the aeroacoustic wave propagation are mostly encoded in their dispersion relations. Hence, a computational aeroacoustic (CAA) algorithm, which reasonably preserves these relations, was investigated. It was derived using an optimization procedure to ensure, that the numerical derivatives preserved the wave number and angular frequency of the differential terms in the linearized, 2-D Euler equations. Then, simulations were performed to validate the scheme and a compatible set of discretized boundary conditions. The computational results were found to agree favorably with the exact solutions. The boundary conditions were transparent to the outgoing waves, except when the disturbance source was close to a boundary. The time-domain data generated by such CAA solutions were often intractable until their spectra was analyzed. Therefore, the relative merits of three different methods were included in the study. For simple, periodic waves, the periodogram method produced better estimates of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window was more effective when used with the weighted-overlapped-segment-averaging and Blackman-Tukey methods gave better results than the periodogram method. Finally, it was demonstrated that the representation of time domain-data was significantly dependent on the particular spectral analysis method employed.

Vanel, Florence O.↗

High-latitude truncation errors of box-type primitive equation models

The 'box-type' finite-difference method includes a weighted average of the pressure gradient with weights proportional to the surface of the grid walls. It is shown that this averaging introduces first-order truncation errors near the poles. An example is shown in which the relative error is of zero order and the scheme produces large distortions in the solution at high latitudes.

Kalnay-Rivas, E.↗

Lunar gravity - A long-term Keplerian rate method

Nearside tracking data from Apollo subsatellites and Lunar Orbiter 5 are used to construct a model of the farside lunar gravity field. An indirect method of calculating the average effect of all the gravity anomalies which affect spacecraft orbits is reduced to a linear estimation problem by eliminating the integration of spacecraft dynamics. Doppler frequency shifts between the lunar satellite and earth-based tracking station are used to measure the line of sight velocity component between the orbiter and the tracker; gravity determination is then performed in a weighted least squares estimator which employs the orbital rates as data and formulates the normal equations using the associated Keplerian state. Simulations performed demonstrate that farside gravity features can be economically determined by use of gravity estimates obtained from long-period models of nearside orbiting tracks.

Ferrari, A. J.↗

Long range orbital error estimation for applications satellites

A method of optimum orbital averaging was employed to study the long range accuracy potential of polar orbiting applications satellites. This approach involved the determination of the boundary conditions of one set of differential equations of motion by adjusting the initial conditions in a least square sense with the use of data generated by another set of differential equations of motion.

Bonavito, N. L.↗

The mechanics of vortex pairing in a circular jet under controlled excitation

The vortex pairing phenomenon in the near field circular jet under controlled excitation is monitored by the method of phase-averaging. Hot-wire measurements were made in a 7.62 cm diameter air at Re(D) = 32,000, excited at St(D) = 0.85 with an excitation amplitude of 3%. The trajectories of the vortex centers and their longitudinal and transverse convection velocities during the pairing program are measured. Spatial distributions of the average velocity, vorticity, coherent Reynolds numbers, streamlines and pseudo-streamfunction are presented at a particular phase of the pairing. These represent field measurements at one particular time and do not involve the use of Taylor hypothesis.

Zaman, K. B. M. Q.↗

A quantitative method for evaluating alternatives

When faced with choosing between alternatives, people tend to use a number of criteria (often subjective, rather than objective) to decide which is the best alternative for them given their unique situation. The subjectivity inherent in the decision-making process can be reduced by the definition and use of a quantitative method for evaluating alternatives. This type of method can help decision makers achieve degree of uniformity and completeness in the evaluation process, as well as an increased sensitivity to the factors involved. Additional side-effects are better documentation and visibility of the rationale behind the resulting decisions. General guidelines for defining a quantitative method are presented and a particular method (called 'hierarchical weighted average') is defined and applied to the evaluation of design alternatives for a hypothetical computer system capability.

Forthofer, M. J.↗

Limits to the burster repetition rate as deduced from the 2nd catalog of the interplanetary network

The burster repetition rate is an important parameter in many gamma ray burst models. The localizations of the interplanetary network, which have a relatively small combined surface area, may be used to estimate the average repetition rate. The method consists of (1) estimating the number of random overlaps between error boxes expected in the catalog and comparing this number to that actually observed; (2) modeling the response of the detectors in the network, so that the probability of detecting a burst can be estimated; and (3) simulating the arrival of bursts at the network assuming that burster repetition is governed by a Poisson process. The application of this method for many different burster luminosity functions shows that (1) the lower limit to the burster repetition rate depends strongly upon the assumed luminosity function; (2) the best lower limit to the repetition period obtainable from the data of the network is about 100 months; and (3) that a luminosity function for all bursters similar to that of the 1979 Mar 5 burster is inconsistent with the data.

Atteia, J. L.↗