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 937 records · Page 52

Polarimetric thermal emission from rough surfaces

Recent theoretical works have suggested the potential of passive polarimetry in the remote sensing of geophysical media. It was shown that the third Stokes parameter U of the thermal emission may become larger for azimuthally asymmetric fields of observation. In order to investigate the potential applicability of passive polarimetry to the remote sensing of ocean surface, measurements of the polarimetric thermal emission from a sinusoidal water surface and a numerical study of the polarimetric thermal emission from randomly rough ocean surfaces were performed. Measurements of sinusoidal water surface thermal emission were performed using a sinusoidal water surface which was created by placing a thin sheet of fiberglass with a sinusoidal profile in two dimensions extended infinitely in the third dimension onto a water surface. The theory of thermal emission from a 'two-layer' periodic surface is derived and the exact solution is performed using both the extended boundary condition method (EBC) and the method of moments (MOM). The theoretical predictions are found to be in good agreement with the experimental results once the effects of the radiometer antenna pattern are included and the contribution of background noise to the measurements is modeled. The experimental results show that the U parameter indicates the direction of periodicity of the water surface and can approach values of up to 30 K for the surface observed. Next, a numerical study of polarimetric thermal emission from randomly rough surfaces was performed. A Monte Carlo technique utilizing an exact method for calculating thermal emission was chosen for the study to avoid any of the limitations of the commonly used approximate methods in rough surface scattering. In this Monte Carlo technique, a set of finite rough surface profiles in two dimensions with desired statistics was generated and extended periodically. The polarimetric thermal emission from each surface of the set was then calculated using both the EBC and the MOM and the results were averaged. The surface statistics chosen were intended to model a wind perturbed ocean surface in the X to Ku band microwave region. The results indicate that the U parameter is sensitive to the azimuthal angle between the surface periodicity and the looking angle and to the rms height of the surface, and that the U parameter is fairly insensitive to variations in polar angle, permittivity, surface power law spectrum, and surface spectrum high frequency cutoff. These properties give further strength to the idea of using the U parameter to detect wind direction over the ocean.

Johnson, J. T.↗

Natural gas flow through critical nozzles

Empirical method for calculating both the mass flow rate and upstream volume flow rate through critical flow nozzles is determined. Method requires knowledge of the composition of natural gas, and of the upstream pressure and temperature.

Johnson, R. C.↗

The influence of line shape and band structure on temperatures in planetary atmospheres

Numerical experiments are performed to examine the effects of line shape and band structure on the radiative equilibrium temperature profile in planetary atmospheres. In order to accurately determine these effects, a method for calculating radiative terms is developed which avoids the usual approximations. It differs from the more commonly used methods in that it allows arbitrary dependence of the absorption coefficient on wave number, without requiring tedious line by line integration and without the constraints of band models. The present formulation is restricted to homogeneous atmospheres but the concept can be extended to the more general case. The numerical experiments reveal that the line shape and band structure of the absorbing gases have a large effect on temperatures in the higher layers of the atmosphere (corresponding to the stratosphere and mesosphere). The more nongrey the spectrum (that is, the higher the peaks and the deeper the troughs in the spectrum), the lower the temperature.

Arking, A.↗

Designing More-Efficient Spur Gears

Relatively simple method to calculate spur-gear system power loss for wide range of gear geometries and operating conditions developed. Method used to determine design requirements for efficient gearset. Effects of spur-gear size, pitch, ratio, pitch-line velocity and load on efficiency readily predictable with method. Analysis uses simple algebraic expressions to determine gear sliding, rolling, and windage losses and incorporates approximate ball-bearing power-loss expression. Predicted results show good agreement with published data.

Loewenthal, S. H.↗

Decoupled direct method for sensitivity analysis in combustion kinetics

An efficient, decoupled direct method for calculating the first order sensitivity coefficients of homogeneous, batch combustion kinetic rate equations is presented. In this method the ordinary differential equations for the sensitivity coefficients are solved separately from , but sequentially with, those describing the combustion chemistry. The ordinary differential equations for the thermochemical variables are solved using an efficient, implicit method (LSODE) that automatically selects the steplength and order for each solution step. The solution procedure for the sensitivity coefficients maintains accuracy and stability by using exactly the same steplengths and numerical approximations. The method computes sensitivity coefficients with respect to any combination of the initial values of the thermochemical variables and the three rate constant parameters for the chemical reactions. The method is illustrated by application to several simple problems and, where possible, comparisons are made with exact solutions and those obtained by other techniques.

Radhakrishnan, Krishnan↗

An equivalent domain integral method in the two-dimensional analysis of mixed mode crack problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies is presented. The details of the method and its implementation are presented for isoparametric elements. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented.

Raju, I. S.↗

Analysis and prediction of Multiple-Site Damage (MSD) fatigue crack growth

A technique was developed to calculate the stress intensity factor for multiple interacting cracks. The analysis was verified through comparison with accepted methods of calculating stress intensity factors. The technique was incorporated into a fatigue crack growth prediction model and used to predict the fatigue crack growth life for multiple-site damage (MSD). The analysis was verified through comparison with experiments conducted on uniaxially loaded flat panels with multiple cracks. Configuration with nearly equal and unequal crack distribution were examined. The fatigue crack growth predictions agreed within 20 percent of the experimental lives for all crack configurations considered.

Dawicke, D. S.↗

A theoretical investigation of the lateral oscillations of an airplane with free rudder with special reference to the effect of friction

Charts showing the variation in dynamic stability with the rudder hinge-moment characteristics are presented. A stabilizing rudder floating tendency combined with a high degree of aerodynamic balance is shown to lead to oscillations of increasing amplitude. This dynamic instability is increased by viscous-friction in the rudder control system. The presence of solid friction in the rudder control system will cause steady oscillations of constant amplitude if the floating angle of the rudder per unit angle of sideslip is stabilizing and greater than a certain critical value that depends on other airplane parameters, such as vertical-tail area and airplane moment of inertia about the vertical axis. The amplitude of the steady oscillation is proportional to the amount of friction and is generally quite small but increases as the condition of dynamic instability is approached. An approximate method of calculating the amplitudes of the steady oscillation is explained and is illustrated by a numerical example. A more accurate step-by-step calculation of the motion is also made and it is shown that the agreement with the approximate method is good.

Greenberg, Harry↗

SEADYN Analysis of a Tow Line for a High Altitude Towed Glider

The concept of using a system, consisting of a tow aircraft, glider and tow line, which would enable subsonic flight at altitudes above 24 km (78 kft) has previously been investigated. The preliminary results from these studies seem encouraging. Under certain conditions these studies indicate the concept is feasible. However, the previous studies did not accurately take into account the forces acting on the tow line. Therefore in order to investigate the concept further a more detailed analysis was needed. The code that was selected was the SEADYN cable dynamics computer program which was developed at the Naval Facilities Engineering Service Center. The program is a finite element based structural analysis code that was developed over a period of 10 years. The results have been validated by the Navy in both laboratory and at actual sea conditions. This code was used to simulate arbitrarily-configured cable structures subjected to excitations encountered in real-world operations. The Navy's interest was mainly for modeling underwater tow lines, however the code is also usable for tow lines in air when the change in fluid properties is taken into account. For underwater applications the fluid properties are basically constant over the length of the tow line. For the tow aircraft/glider application the change in fluid properties is considerable along the length of the tow line. Therefore the code had to be modified in order to take into account the variation in atmospheric properties that would be encountered in this application. This modification consisted of adding a variable density to the fluid based on the altitude of the node being calculated. This change in the way the code handled the fluid density had no effect on the method of calculation or any other factor related to the codes validation.

Colozza, Anthony J.↗

A New Study of Sea Spray Optical Properties from Multi-Sensor Spaceborne Observations

Retrievals of aerosol optical depth (AOD) from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) satellite sensor require the assumption of an extinction-to-backscatter ratio, also known as the lidar ratio. This paper evaluates a new method to calculate lidar ratio of sea spray aerosol using two independent sources: the AOD from Synergized Optical Depth of Aerosols (SODA) and the integrated attenuated backscatter from CALIOP. The method applied in this study allows particulate lidar ratio to be calculated for individual CALIOP retrievals of single aerosol layer columns over the ocean. Analyses are carried out using CALIOP level 2, 5km sea spray aerosol layer products and collocated SODA nighttime data from December 2007 to December 2009. The global mean lidar ratio for sea spray aerosols was found to be 26 sr, roughly 30 higher than the one prescribed by CALIOP. Data analysis also showed considerable spatiotemporal variability in calculated lidar ratio over different parts of the remote oceans. The calculated aerosol lidar ratios are shown to be inversely related to the mean ocean surface wind speed: increase in ocean surface wind speed (U10) from 0 to 15 ms-1 reduces the mean lidar ratios for sea spray particles from 32 sr (for 0 U10 4 ms-1) to 22 sr (for U10 15 ms-1). Such changes in the lidar ratio are expected to have a corresponding effect on sea spray AOD. The outcomes of this study are relevant for future improvements of the SODA and CALIOP operational product and could lead to more accurate retrievals of sea spray AOD.

aerosol optical depth↗

Noniterative grid generation using parabolic difference equations for fuselage-wing flow calculations

A fast method for generating three-dimensional grids for fuselage-wing transonic flow calculations using parabolic difference equations is described. No iterative scheme is used in the three-dimensional sense; grids are generated from one grid surface to the next starting from the fuselage surface. The computational procedure is similar to the iterative solution of the two-dimensional heat conduction equation. The proposed method is at least 10 times faster than the elliptic grid generation method and has much smaller memory requirements. Results are presented for a fuselage and wing of NACA-0012 section and thickness ratio of 10 percent. Although only H-grids are demonstrated, the present technique should be applicable to C-grids and O-grids in three dimensions.

Nakamura, S.↗

General theory of the transverse dielectric constant of III-V semiconducting compounds

A general model of the transverse dielectric constant of III-V compounds is developed using a hybrid method which combines the kp method with a nonlocal pseudopotential calculation. In this method the Brillouin zone is partitioned into three regions by expanding the energy bands and matrix elements about the F, X, and L symmetry points. The real and imaginary parts of the dielectric constant are calculated as a sum of the individual contributions of each region. By using this partition method, it is possible to get good insight into the dependence of the dielectric constant on the shape of the band structure.

Kahen, K. B.↗

Scattering from thin dielectric straps surrounding a perfectly conducting structure

A method to calculate the electromagnetic scattered fields from a dielectric strap wrapped around convex, conducting structure is presented. A moment method technique is used to find the current excited within the strap by the incident plane wave. Then, Uniform Geometrical Theory of Diffraction (UTD) is used to compute the fields scattered by the strap. Reasonable agreement was obtained between the computed and the measured results. The results found in this study are useful in evaluating straps as a target support structure for scattering measurements.

Al-Hekail, Zeyad↗

Incorporating finite element analysis into component life and reliability

A method for calculating a component's design survivability by incorporating finite element analysis and probabilistic material properties was developed. The method evaluates design parameters through direct comparisons of component survivability expressed in terms of Weibull parameters. The analysis was applied to a rotating disk with mounting bolt holes. The highest probability of failure occurred at, or near, the maximum shear stress region of the bolt holes. Distribution of failure as a function of Weibull slope affects the probability of survival. Where Weibull parameters are unknown for a rotating disk, it may be permissible to assume Weibull parameters, as well as the stress-life exponent, in order to determine the disk speed where the probability of survival is highest.

August, Richard↗

Comparison of two-dimensional and three-dimensional droplet trajectory calculations in the vicinity of finite wings

Computational predictions of ice accretion on flying aircraft most commonly rely on modeling in two dimensions (2D). These 2D methods treat an aircraft geometry either as wing-like with infinite span, or as an axisymmetric body. Recently, fully three dimensional (3D) methods have been introduced that model an aircrafts true 3D shape. Because 3D methods are more computationally expensive than 2D methods, 2D methods continue to be widely used. However, a 3D method allows us to investigate whether it is valid to continue applying 2D methods to a finite wing. The extent of disagreement between LEWICE, a 2D method, and LEWICE3D, a 3D method, in calculating local collection efficiencies at the leading edge of finite wings is investigated in this paper.

Mohler, Stanley R., Jr.↗

Backscattering by nonspherical hydrometeors as calculated by the coupled-dipole method - An application in radar meteorology

The severest test of a theory of scattering by particles is how well it calculates scattering in the backward direction. The coupled-dipole method can be used for accurately calculating backscattering at 94 GHz by hexagonal ice crystals. Backscattering by columns is markedly different from that by plates, which indicates that it might be possible to infer size and shape distributions of ice crystals using recently developed millimeter wave radar.

Dungey, Clifton E.↗

Cell volume control at a surface for three-dimensional grid generation packages

An alternate method of calculating the cell size for orthogonality control in the solution of Poisson's 3D space equations is presented. The method provides the capability to enforce a better initial guess for the grid distribution required for boundary layer resolution. This grid point distribution is accomplished by enforcing grid spacing from a grid block boundary where orthogonality is required. The actual grid spacing or cell size for that boundary is determined by the two or four adjacent boundaries in the grid block definition, which are two dimensional grids. These two dimensional grids are in turn defined by the user using insight into the flow field and boundary layer characteristics. The adjoining boundaries are extended using a multifunctional blending scheme, with user control of the blending and interpolating functions to be used. This grid generation procedure results in an enhanced computational fluid dynamics calculation by allowing a quicker resolution of the configuration's boundary layer and flow field and by limiting the number of grid re-adaptations. The cell size specification calculation was applied to a variety of configurations ranging from axisymmetric to complex three-dimensional configurations. Representative grids are shown for the Space Shuttle and the Langley Lifting Body (HL-20).

Alter, Stephen J.↗