Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic derivatives”

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 271 records · Page 15

Maneuver Planning for Conjunction Risk Mitigation with Ground-track Control Requirements

The planning of conjunction Risk Mitigation Maneuvers (RMM) in the presence of ground-track control requirements is analyzed. Past RMM planning efforts on the Aqua, Aura, and Terra spacecraft have demonstrated that only small maneuvers are available when ground-track control requirements are maintained. Assuming small maneuvers, analytical expressions for the effect of a given maneuver on conjunction geometry are derived. The analytical expressions are used to generate a large trade space for initial RMM design. This trade space represents a significant improvement in initial maneuver planning over existing methods that employ high fidelity maneuver models and propagation.

McKinley, David↗

Vibration characteristics of 1/8-scale dynamic models of the space-shuttle solid-rocket boosters

Vibration tests and analyses of six 1/8 scale models of the space shuttle solid rocket boosters are reported. Natural vibration frequencies and mode shapes were obtained for these aluminum shell models having internal solid fuel configurations corresponding to launch, midburn (maximum dynamic pressure), and near endburn (burnout) flight conditions. Test results for longitudinal, torsional, bending, and shell vibration frequencies are compared with analytical predictions derived from thin shell theory and from finite element plate and beam theory. The lowest analytical longitudinal, torsional, bending, and shell vibration frequencies were within + or - 10 percent of experimental values. The effects of damping and asymmetric end skirts on natural vibration frequency were also considered. The analytical frequencies of an idealized full scale space shuttle solid rocket boosted structure are computed with and without internal pressure and are compared with the 1/8 scale model results.

Leadbetter, S. A.↗

Partially coherent scattering in stellar chromospheres. II - The first-order escape probability method. III - A second-order escape probability method

Approximate analytic expressions are derived for resonance-line wing diagnostics, accounting for frequency redistribution effects, for homogeneous slabs, and slabs with a constant Planck function gradient. Resonance-line emission profiles from a simplified conceptual standpoint are described in order to elucidate the basic physical parameters of the line-forming layers prior to the performance of detailed numerical calculations. An approximate analytic expression is derived for the dependence on stellar surface gravity of the location of the Ca II and Mg II resonance-line profile peaks. An approximate radiative transfer equation using generalized second-order escape probabilities, applicable even in the presence of nearly coherent scattering in the damping wings of resonance lines, is derived. Approximate analytic solutions that can be applied in special regimes and achieve good agreement with accurate numerical results are found.

Gayley, K. G.↗

Nonlinear feedback control for high alpha flight

Analytical aerodynamic models are derived from a high alpha 6 DOF wind tunnel model. One detail model requires some interpolation between nonlinear functions of alpha. One analytical model requires no interpolation and as such is a completely continuous model. Flight path optimization is conducted on the basic maneuvers: half-loop, 90 degree pitch-up, and level turn. The optimal control analysis uses the derived analytical model in the equations of motion and is based on both moment and force equations. The maximum principle solution for the half-loop is poststall trajectory performing the half-loop in 13.6 seconds. The agility induced by thrust vectoring capability provided a minimum effect on reducing the maneuver time. By means of thrust vectoring control the 90 degrees pitch-up maneuver can be executed in a small place over a short time interval. The agility capability of thrust vectoring is quite beneficial for pitch-up maneuvers. The level turn results are based currently on only outer layer solutions of singular perturbation. Poststall solutions provide high turn rates but generate higher losses of energy than that of classical sustained solutions.

Stalford, Harold↗

Design sensitivity derivatives for isoparametric elements by analytical and semi-analytical approaches

This paper presents an analytical approach for incorporating design sensitivity calculations directly into the finite element analysis. The formulation depends on the implicit differentiation approach and requires few additional calculations to obtain the design sensitivity derivatives. In order to evaluate this approach, it is compared with the semi-analytical approach which is based on commonly used finite difference formulations. Both approaches are implemented to calculate the design sensitivities for continuum and structural isoparametric elements. To demonstrate the accuracy and robustness of the developed analytical approach compared to the semi-analytical approach, some test cases using different structural and continuum element types are presented.

Zumwalt, Kenneth W.↗

Approximate solutions of radiative transfer in dusty nebulae. 2: Hydrogen and helium

The ionization structure of hydrogen and helium in dusty nebulae is discussed. The general equations for pure hydrogen nebulae are modified to include helium. An approximate analytic solution is presented for Stroemgren radii of hydrogen and helium in the absence of dust. It is shown that these results, with simple modifications, are also applicable to dusty nebulae where the effective absorption cross section of dust grains varies slowly with frequency in the 1000 to 200A range. No analytic solutions are possible if this cross section varies rapidly with frequency. In this case, however, simple coupled differential equations are derived. Approximate analytic expressions are presented for evaluation of the variation of the fraction of ionizing radiation absorbed by dust and the ratio of the volume emission measures of He II to H II regions with the spectrum of the ionizing source, helium abundance, and absorption properties of dust. The effects of dust on the He III zone are discussed. Results are restricted to spherically symmetric nebulae, but nonuniform gas and dust distributions and clumpiness can be taken into account.

Dana, R. A.↗

A theory for creep crack growth

The paper derives an analytical expression for creep crack growth rate based on a model in which the higher and concentrated stresses ahead of the crack tip enhance creep deformation, thus progressively causing stress rupture and continuous crack advance. The equation derived for creep crack growth rate is expressed in terms of a geometrical factor in the stress intensity expression and in terms of the crack tip radius whose lower limit estimate is the fracture mechanics crack-tip opening displacement. The functional features of the derived creep crack growth rate equation can be compared with those of a pertinent empirical equation. The analytical expression provides some guidelines for alloy and microstructure design.

Purushothaman, S.↗

Analytic expressions for perturbations and partial derivatives of range and range rate of a spacecraft with respect to the coefficient of the second harmonic

Closed-form analytic expressions for the time variations of instantaneous orbital parameters and of the topocentric range and range rate of a spacecraft moving in the gravitational field of an oblate large body are derived using a first-order variation of parameters technique. In addition, the closed-form analytic expressions for the partial derivatives of the topocentric range and range rate are obtained, with respect to the coefficient of the second harmonic of the potential of the central body (J sub 2). The results are applied to the motion of a point-mass spacecraft moving in the orbit around the equatorially elliptic, oblate sun, with J sub 2 approximately equal to .000027.

Georgevic, R. M.↗

Emission, absorption and polarization of gyrosynchrotron radiation of mildly relativistic particles

Approximate analytic expressions are presented for the emissivity and absorption coefficient of synchrotron radiation of mildly relativistic particles with an arbitrary energy spectrum and pitch angle distribution. From these, an expression for the degree of polarization is derived. The analytic results are compared with numerical results for both thermal and non-thermal (power law) distributions of particles.

Petrosian, V.↗

Convergent-tapered annular seals - Analysis and testing for rotordynamic coefficients

A combined analytical-computational method has been developed to calculate the pressure field and dynamic coefficients for tapered high-pressure annular seals. Completely developed turbulent flow is assumed in both the circumferential and axial directions, according to Hirs' (1973) bulk-flow turbulent-lubrication equations. In a numerical experiment with the method, a short bearing approximation is used to derive an analytical expression for the first-order (dynamic) pressure gradient. This expression is integrated numerically to define dynamic coefficients of the seal. Numerical results of the integration are compared with previous results for straight and tapered seals. It is found that direct stiffness and leakage coefficients increase in the present seal, while the remaining coefficients decrease. An optimal taper angle is shown to exist with respect to: (1) the direct stiffness and (2) the ratio of direct stiffness to leakage. On the basis of the optimal taper angle calculations, stiffness increases on the order of 40 to 50 percent are predicted. Experimental results for three taper angles are presented, and are found to be in good agreement with the theoretical predictions.

Childs, D. W.↗

Effect of blunt trailing edge on rotor broadband noise

The production of high-frequency broadband noise by turbulent vortex shedding from rotor blades with blunt trailing edges is investigated analytically. The derivation of the governing equations, analogous to that of Kim and George (1982) for boundary-layer/trailing-edge noise, is explained, and numerical results are compared with the experimental data of Hubbard et al. (1981) and Lowson et al. (1972) in graphs. It is shown that vortex-shedding noise is a significant component of blunt-trailing-edge rotor broadband noise and that the analytical method employed gives reasonable predictions. The need for a better empirical expression for the normalized spectrum and for more measurements of surface pressure fluctuations near blunt trailing edges is indicated.

Chou, S.-T.↗

Shape sensitivity analysis of wing static aeroelastic characteristics

A method is presented to calculate analytically the sensitivity derivatives of wing static aeroelastic characteristics with respect to wing shape parameters. The wing aerodynamic response under fixed total load is predicted with Weissinger's L-method; its structural response is obtained with Giles' equivalent plate method. The characteristics of interest include the spanwise distribution of lift, trim angle of attack, rolling and pitching moments, wind induced drag, as well as the divergence dynamic pressure. The shape parameters considered are the wing area, aspect ratio, taper ratio, sweep angle, and tip twist angle. Results of sensitivity studies indicate that: (1) approximations based on analytical sensitivity derivatives can be used over wide ranges of variations of the shape parameters considered, and (2) the analytical calculation of sensitivity derivatives is significantly less expensive than the conventional finite-difference alternative.

Barthelemy, Jean-Francois M.↗

Analytical Approach Validation for the Spin-Stabilized Satellite Attitude

An analytical approach for spin-stabilized spacecraft attitude prediction is presented for the influence of the residual magnetic torques and the satellite in an elliptical orbit. Assuming a quadripole model for the Earth s magnetic field, an analytical averaging method is applied to obtain the mean residual torque in every orbital period. The orbit mean anomaly is used to compute the average components of residual torque in the spacecraft body frame reference system. The theory is developed for time variations in the orbital elements, giving rise to many curvature integrals. It is observed that the residual magnetic torque does not have component along the spin axis. The inclusion of this torque on the rotational motion differential equations of a spin stabilized spacecraft yields conditions to derive an analytical solution. The solution shows that the residual torque does not affect the spin velocity magnitude, contributing only for the precession and the drift of the spin axis of the spacecraft. The theory developed has been applied to the Brazilian s spin stabilized satellites, which are quite appropriated for verification and comparison of the theory with the data generated and processed by the Satellite Control Center of Brazil National Research Institute. The results show the period that the analytical solution can be used to the attitude propagation, within the dispersion range of the attitude determination system performance of Satellite Control Center of Brazil National Research Institute.

Zanardi, Maria Cecilia F. P. S.↗

Shape sensitivity analysis of wing static aeroelastic characteristics

A method is presented to calculate analytically the sensitivity derivatives of wing static aeroelastic characteristics with respect to wing shape parameters. The wing aerodynamic response under fixed total load is predicted using Weissinger's L-method; its structural response is obtained with Giles' equivalent plate method. The characteristics of interest in this study include the spanwise distribution of lift, trim angle of attack, rolling and pitching moments, wing induced drag as well as the divergence dynamic pressure. The shape parameters considered are the wing area, aspect ratio, taper ratio, sweep angle and tip twist angle. Results of sensitivity studies indicate that (1) approximations based on analytical sensitivity derivatives can be used for wide ranges of variations of the shape parameters considered and (2) the analytical calculation of sensitivity derivatives is significantly less expensive than the conventional finite difference alternative.

Barthelemy, Jean-Francois M.↗

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

Long-term evolution of crustal neutron star magnetic fields

We have derived an analytic solution to the asymptotic behavior of dipolar magnetic fields that are generated in the crusts of neutron stars. We show that if the conductivity is due to impurity scattering, as expected for late stages of evolution, the surface field strength at the magnetic pole declines with the power law B(sub p) approximately = (t/t(sub 0))(exp -2/3). The results are shown to be qualitatively consistent with detailed numerical calculations. These latter results are consistent with some recent analyses of pulsar statistics and the magnetic fields of several binary pulsars with white dwarf companions whose ages have been determined. The dependence of the surface magnetic field on spin period of the pulsar is derived.

Urpin, V. A.↗

Monte Carlo simulation of particle acceleration at astrophysical shocks

A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.

Campbell, Roy K.↗