Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation 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 1,135 records · Page 63

Integrated aerodynamic-structural-control wing design

The aerodynamic-structural-control design of a forward-swept composite wing for a high subsonic transport aircraft is considered. The structural analysis is based on a finite-element method. The aerodynamic calculations are based on a vortex-lattice method, and the control calculations are based on an output feedback control. The wing is designed for minimum weight subject to structural, performance/aerodynamic and control constraints. Efficient methods are used to calculate the control-deflection and control-effectiveness sensitivities which appear as second-order derivatives in the control constraint equations. To suppress the aeroelastic divergence of the forward-swept wing, and to reduce the gross weight of the design aircraft, two separate cases are studied: (1) combined application of aeroelastic tailoring and active controls; and (2) aeroelastic tailoring alone. The results of this study indicated that, for this particular example, aeroelastic tailoring is sufficient for suppressing the aeroelastic divergence, and the use of active controls was not necessary.

Rais-Rohani, M.↗

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.↗

Relativistic core–valence-separated equation-of-motion coupled-cluster singles and doubles method: Efficient implementation and benchmark calculations

An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast and sensitive measurements of sub-3 nm particles using Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM)

New particle formation (NPF) is the atmospheric process whereby gas molecules react and nucleate to form detectable particles. NPF has a strong impact on Earth's radiative balance as it produces roughly half of global cloud condensation nuclei. However, the time resolution and sensitivity of current instrumentation are inadequate in measuring the size distribution of sub-3 nm particles, the particles critical for understanding NPF. Here we present the Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM), a method to measure the concentrations of freshly nucleated particles. The CPC FARM consists of five CPCs operating in parallel, each configured to operate at different detectable particle sizes between 1–3 nm. This study explores two methods to calculate the size distribution from the differential measurements across the CPC channels. The performance of both inversion methods was tested against the size distribution measured by a pair of stepping particle mobility sizers (SMPSs) during an ambient air sampling study in Pittsburgh, PA. Observational results indicate that the CPC FARM is more accurate with higher time resolution and sensitivity in the sub-3 nm range compared to the SMPS.

Cheng, Darren [Carnegie Mellon Univ., Pittsburgh, ↗

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↗