Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CALCULATION 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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Estimation of effective hydrologic properties of soils from observations of vegetation density

An existing one-dimensional model of the annual water balance is reviewed. Slight improvements are made in the method of calculating the bare soil component of evaporation, and in the way surface retention is handled. A natural selection hypothesis, which specifies the equilibrium vegetation density for a given, water limited, climate-soil system, is verified through comparisons with observed data and is employed in the annual water balance of watersheds in Clinton, Ma., and Santa Paula, Ca., to estimate effective areal average soil properties. Comparison of CDF's of annual basin yield derived using these soil properties with observed CDF's provides excellent verification of the soil-selection procedure. This method of parameterization of the land surface should be useful with present global circulation models, enabling them to account for both the non-linearity in the relationship between soil moisture flux and soil moisture concentration, and the variability of soil properties from place to place over the Earth's surface.

Tellers, T. E.↗
NASA NTRS1983

Computing the projected area of a spacecraft solar array with matrices

Matrix methods for computing the projected area of a solar array as a function of the rotational position of a spacecraft and of array position on a spacecraft are presented. Formulas are derived which provide the optimum solar array pitch, cant and tracking angles for a given spacecraft configuration and orbit. These formulas are general and applicable to many spacecraft. Formulas are also provided for determining the energy output from an array for a given orbit and the resultant energy available to spacecraft loads. Results are simply obtained and realizable with a hand calculator. The methods above can be extended to the case of computer analysis of solar array shadowing. How this can be done is outlined and results are presented from a spacecraft study program.

E. M. Gaddy↗
NASA NTRS2006

Maximum von Mises Stress in the Loading Environment of Mass Acceleration Curve

Method for calculating stress due to acceleration loading: 1) Part has been designed by FEA and hand calculation in one critical loading direction judged by the analyst; 2) Maximum stress can be due to loading in another direction; 3) Analysis procedure to be presented determines: a) The maximum Mises stress at any point; and b) The direction of maximum loading associated with the "stress". Concept of Mass Acceleration Curves (MAC): 1) Developed by JPL to perform preliminary structural sizing (i.e. Mariners, Voyager, Galileo, Pathfinder, MER,...MSL); 2) Acceleration of physical masses are bounded by a curve; 3) G-levels of vibro-acoustic and transient environments; 4) Convergent process before the couple loads cycle; and 5) Semi-empirical method to effectively bound the loads, not a simulation of the actual response.

von Mises stress↗
NASA NTRS2012

Causal Correlation Functions and Fourier Transforms: Application in Calculating Pressure Induced Shifts

By adopting a concept from signal processing, instead of starting from the correlation functions which are even, one considers the causal correlation functions whose Fourier transforms become complex. Their real and imaginary parts multiplied by 2 are the Fourier transforms of the original correlations and the subsequent Hilbert transforms, respectively. Thus, by taking this step one can complete the two previously needed transforms. However, to obviate performing the Cauchy principal integrations required in the Hilbert transforms is the greatest advantage. Meanwhile, because the causal correlations are well-bounded within the time domain and band limited in the frequency domain, one can replace their Fourier transforms by the discrete Fourier transforms and the latter can be carried out with the FFT algorithm. This replacement is justified by sampling theory because the Fourier transforms can be derived from the discrete Fourier transforms with the Nyquis rate without any distortions. We apply this method in calculating pressure induced shifts of H2O lines and obtain more reliable values. By comparing the calculated shifts with those in HITRAN 2008 and by screening both of them with the pair identity and the smooth variation rules, one can conclude many of shift values in HITRAN are not correct.

correlation function↗
NASA NTRS2017

Kepler Planet Detection Metrics: Per-Target Detection Contours for Data Release 25

A necessary input to planet occurrence calculations is an accurate model for the pipeline completeness (Burke et al., 2015). This document describes the use of the Kepler planet occurrence rate products in order to calculate a per-target detection contour for the measured Data Release 25 (DR25) pipeline performance. A per-target detection contour measures for a given combination of orbital period, Porb, and planet radius, Rp, what fraction of transit signals are recoverable by the Kepler pipeline (Twicken et al., 2016; Jenkins et al., 2017). The steps for calculating a detection contour follow the procedure outlined in Burke et al. (2015), but have been updated to provide improved accuracy enabled by the substantially larger database of transit injection and recovery tests that were performed on the final version (i.e., SOC 9.3) of the Kepler pipeline (Christiansen, 2017; Burke Catanzarite, 2017a). In the following sections, we describe the main inputs to the per-target detection contour and provide a worked example of the python software released with this document (Kepler Planet Occurrence Rate Tools KeplerPORTs)1 that illustrates the generation of a detection contour in practice. As background material for this document and its nomenclature, we recommend the reader be familiar with the previous method of calculating a detection contour (Section 2 of Burke et al.,2015), input parameters relevant for describing the data quantity and quality of Kepler targets (Burke Catanzarite, 2017b), and the extensive new transit injection and recovery tests of the Kepler pipeline (Christiansen et al., 2016; Burke Catanzarite, 2017a; Christiansen, 2017).

Planet Detection Metrics↗
NASA NTRS2000

Effect of Negative Ions on the Conductivity of the Titan Atmosphere

In an earlier paper, Borucki et al (1987) calculated the electrical conductivity and electrical charge on aerosols in Titan's atmosphere due to the ionization by galactic cosmic rays and electron precipitation from Saturn's magnetosphere. The lower atmosphere was predicted to be substantially more conducting than the atmospheres of Earth and Venus because of the high concentration of free electrons. The prediction of a high conductivity is based on the lack of electrophillic species which form negative ions with low mobility and which reduce the number of free electrons. At that time, no molecular species capable of forming negative ions in concentrations sufficient to perturb the atmospheric conductivity were identified. Recently, E. Bakes and her colleagues have been investigating the formation of nitrogenous macromolecules using quantum mechanical methods. Their calculations indicate that the molecules will be highly electrophillic and are likely to be present in the atmosphere at mixing ratios of order 10(exp -7). This mixing ratio is sufficiently large that a substantial reduction in the conductivity is expected at altitudes below 100 km. Revision of the atmospheric model to accommodate the presence of negative ions and to increase the fidelity of the modeling will be described.

Borucki, W. J.↗
NASA NTRS1989

Estimation of velocity perturbations in satellite fragmentation events

The magnitude, variance and directionality of the velocity perturbations of the fragments of a satellite can shed valuable information regarding the nature and intensity of the fragmentation. Up until now, the only method used to calculate the three orthogonal components of the velocity change consisted of inverting the process of evaluating the changes in the orbital elements of the fragments due to velocity perturbing forces. But the traditional method failed in five different cases: (1) when the parent satellite's orbit was circular; (2) and (3) when the true anomaly of the parent was either 0 deg or 180 deg; and (4) and (5) when the argument of latitude of the parent was 90 deg or 270 deg. Described here is a new method of calculating the velocity perturbations which is free from the shortcomings of the traditional method and could be used in all occasions, provided the fragmentation data and the orbital elements data are consistent with one another.

Tan, Arjun↗
NASA NTRS1972

A theoretical and experimental investigation of the hypersonic turbulent boundary layer subject to normal and longitudinal pressure gradients and cross flow along a windward plane of symmetry.

A theoretical and experimental investigation for the hypersonic turbulent boundary layer undergoing both normal and longitudinal pressure gradients and cross flow along the plane of symmetry on the windward side of a body of revolution is presented. The purpose of this investigation is to provide a method for calculating boundary layers that occur on the centerline of symmetry of an inlet. The three-dimensional compressible integral equations are written along the symmetry plane and integrated numerically, and the effect of the cross flow on the behavior of the boundary layer is also studied. The experiments were performed at Mach number 6 and Reynolds number of 44,000,000/ft for the axisymmetric case and for the case with cross flow. The velocity profiles and the boundary layer parameters agree relatively well with theoretical values with and without cross flow effects. The presence of the cross flow alters significantly the value of the profiles. Total temperature-velocity correlations deviate from the linear Crocco distribution. Heat transfer rates agree with theoretical values and values obtained by the flat plate reference enthalpy method. Cross flow influence on the heat transfer rates is of minor relevance.

Zakkay, V.↗
NASA NTRS2005

Simplified Solar Modulation Model of Inner Trapped Belt Proton Flux As a Function of Atmospheric Density

No simple algorithm seems to exist for calculating proton fluxes and lifetimes in the Earth's inner, trapped radiation belt throughout the solar cycle. Most models of the inner trapped belt in use depend upon AP8 which only describes the radiation environment at solar maximum and solar minimum in Cycle 20. One exception is NOAAPRO which incorporates flight data from the TIROS/NOAA polar orbiting spacecraft. The present study discloses yet another, simple formulation for approximating proton fluxes at any time in a given solar cycle, in particular between solar maximum and solar minimum. It is derived from AP8 using a regression algorithm technique from nuclear physics. From flux and its time integral fluence, one can then approximate dose rate and its time integral dose. It has already been published in this journal that the absorbed dose rate, D, in the trapped belts exhibits a power law relationship, D = A(rho)(sup -n), where A is a constant, rho is the atmospheric density, and the index n is weakly dependent upon shielding. However, that method does not work for flux and fluence. Instead, we extend this idea by showing that the power law approximation for flux J is actually bivariant in energy E as well as density rho. The resulting relation is J(E,rho)approx.(sum of)A(E(sup n))rho(sup -n), with A itself a power law in E. This provides another method for calculating approximate proton flux and lifetime at any time in the solar cycle. These in turn can be used to predict the associated dose and dose rate.

Wilson, Thomas L.↗
NASA NTRS1993

Preconditioned conjugate-gradient methods for low-speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations is integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the Lower-Upper Successive Symmetric Over-Relaxation iterative scheme is more efficient than a preconditioner based on Incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional Line Gauss-Seidel Relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS2009

Optimization of Layer Densities for Spacecraft Multilayered Insulation Systems

Numerous tests of various multilayer insulation systems have indicated that there are optimal densities for these systems. However, the only method of calculating this optimal density was by a complex physics based algorithm developed by McIntosh. In the 1970's much data were collected on the performance of these insulation systems with many different variables analyzed. All formulas generated included number of layers and layer density as geometric variables in solving for the heat flux, none of them was in a differentiable form for a single geometric variable. It was recently discovered that by converting the equations from heat flux to thermal conductivity using Fourier's Law, the equations became functions of layer density, temperatures, and material properties only. The thickness and number of layers of the blanket were merged into a layer density. These equations were then differentiated with respect to layer density. By setting the first derivative equal to zero, and solving for the layer density, the critical layer density was determined. Taking a second derivative showed that the critical layer density is a minimum in the function and thus the optimum density for minimal heat leak, this is confirmed by plotting the original function. This method was checked and validated using test data from the Multipurpose Hydrogen Testbed which was designed using McIntosh's algorithm.

Johnson, W. L.↗
NASA NTRS2012

Wavefront Sensing Analysis of Grazing Incidence Optical Systems

Wavefront sensing is a process by which optical system errors are deduced from the aberrations in the image of an ideal source. The method has been used successfully in near-normal incidence, but not for grazing incidence systems. This innovation highlights the ability to examine out-of-focus images from grazing incidence telescopes (typically operating in the x-ray wavelengths, but integrated using optical wavelengths) and determine the lower-order deformations. This is important because as a metrology tool, this method would allow the integration of high angular resolution optics without the use of normal incidence interferometry, which requires direct access to the front surface of each mirror. Measuring the surface figure of mirror segments in a highly nested x-ray telescope mirror assembly is difficult due to the tight packing of elements and blockage of all but the innermost elements to normal incidence light. While this can be done on an individual basis in a metrology mount, once the element is installed and permanently bonded into the assembly, it is impossible to verify the figure of each element and ensure that the necessary imaging quality will be maintained. By examining on-axis images of an ideal point source, one can gauge the low-order figure errors of individual elements, even when integrated into an assembly. This technique is known as wavefront sensing (WFS). By shining collimated light down the optical axis of the telescope and looking at out-of-focus images, the blur due to low-order figure errors of individual elements can be seen, and the figure error necessary to produce that blur can be calculated. The method avoids the problem of requiring normal incidence access to the surface of each mirror segment. Mirror figure errors span a wide range of spatial frequencies, from the lowest-order bending to the highest order micro-roughness. While all of these can be measured in normal incidence, only the lowest-order contributors can be determined through this WFS technique.

Rohrbach, Scott↗
NASA NTRS2011

Packing Optimization of Sorbent Bed Containing Dissimilar and Irregular Shaped Media

The Fire Cartridge is a packed bed air filter with two different and separate layers of media designed to provide respiratory protection from combustion products after a fire event on the International Space Station (ISS). The first layer of media is a carbon monoxide catalyst and the second layer of media is universal carbon. During development of Fire Cartridge prototypes, the two media beds were noticed to have shifted inside the cartridge. The movement of media within the cartridge can cause mixing of the bed layers, air voids, and channeling, which could cause preferential air flow and allow contaminants to pass through without removal. An optimally packed bed mitigates these risks and ensures effective removal of contaminants from the air. In order to optimally pack each layer, vertical, horizontal, and orbital agitations were investigated and a packed bulk density was calculated for each method. Packed bulk density must be calculated for each media type to accommodate variations in particle size, shape, and density. Additionally, the optimal vibration parameters must be re-evaluated for each batch of media due to variations in particle size distribution between batches. For this application it was determined that orbital vibrations achieve an optimal pack density and the two media layers can be packed by the same method. Another finding was media with a larger size distribution of particles achieve an optimal bed pack easier than media with a smaller size distribution of particles.

Holland, Nathan↗
NASA NTRS1999

Sensitivity of Lumped Constraints Using the Adjoint Method

Adjoint sensitivity calculation of stress, buckling and displacement constraints may be much less expensive than direct sensitivity calculation when the number of load cases is large. Adjoint stress and displacement sensitivities are available in the literature. Expressions for local buckling sensitivity of isotropic plate elements are derived in this study. Computational efficiency of the adjoint method is sensitive to the number of constraints and, therefore, the method benefits from constraint lumping. A continuum version of the Kreisselmeier-Steinhauser (KS) function is chosen to lump constraints. The adjoint and direct methods are compared for three examples: a truss structure, a simple HSCT wing model, and a large HSCT model. These sensitivity derivatives are then used in optimization.

Akgun, Mehmet A.↗
NASA NTRS1990

A study of the limitations of linear theory methods as applied to sonic boom calculations

Current sonic boom minimization theories have been reviewed to emphasize the capabilities and flexibilities of the methods. Flexibility is important because it is necessary for the designer to meet optimized area constraints while reducing the impact on vehicle aerodynamic performance. Preliminary comparisons of sonic booms predicted for two Mach 3 concepts illustrate the benefits of shaping. Finally, for very simple bodies of revolution, sonic boom predictions were made using two methods - a modified linear theory method and a nonlinear method - for signature shapes which were both farfield N-waves and midfield waves. Preliminary analysis on these simple bodies verified that current modified linear theory prediction methods become inadequate for predicting midfield signatures for Mach numbers above 3. The importance of impulse is sonic boom disturbance and the importance of three-dimensional effects which could not be simulated with the bodies of revolution will determine the validity of current modified linear theory methods in predicting midfield signatures at lower Mach numbers.

Darden, Christine M.↗
NASA NTRS2016

Covariance Manipulation for Conjunction Assessment

Use of probability of collision (Pc) has brought sophistication to CA. Made possible by JSpOC precision catalogue because provides covariance. Has essentially replaced miss distance as basic CA parameter. Embrace of Pc has elevated methods to 'manipulate' covariance to enable/improve CA calculations. Two such methods to be examined here; compensation for absent or unreliable covariances through 'Maximum Pc' calculation constructs, projection (not propagation) of epoch covariances forward in time to try to enable better risk assessments. Two questions to be answered about each; situations to which such approaches are properly applicable, amount of utility that such methods offer.

Probability↗
NASA NTRS1963

RADIATIVE ENERGY TRANSFER ON ENTRY INTO MARS AND VENUSQUARTERLY REPORT NO. 4

The emissivity of carbon dioxide in the 2.7 4 µ region has been measured at temperatures from 2000°K to 3500°K, for optical densities less than 5 atm-cm. The high temperatures were obtained by passing a shock wave through the test gas. A total gas pressure of one atm was maintained by varying the carbon dioxide concentration in an inert gas. Emissivities obtained from absorption and emission measurements are in agreement. results are compared to an approximate method for calculating emissivity, The which considers the molecular population concentrated in a single band whose center is temperature dependent. A comparison with previous results suggests that this single band model will overestimate the peak and underestimate the width of the observed values. The observed and calculated emissivities are in agreement if the intensities of the (n 1 +1,n 2 , n 3 + 1 → n 1 , n 2 , n 3 ) and (n 1 ,n 2 +2, n 3 + 1 → n 1 , n 2 , n 3 ) transitions are assumed proportional to (n 1 +1) (n 3 +1) and (n 2 +1)(n 3 +1)/2. The use of this single band model which was originally intended to provide a rough prediction of emissivity, indicates that the use of the appropriate terms of the electric dipole expansion is satisfactory for the (n 1 +1,n 2 , n 3 + 1 → n 1 , n 2 , n 3 ) transition, but is not for (n 1 ,n 2 +2, n 3 + 1 → n 1 , n 2 , n 3 ) transition.

radiative heat transfer↗