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 415 records · Page 23

Finite-analytic numerical solution of heat transfer in two-dimensional cavity flow

Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.

Chen, C.-J.↗

Pair production in superstrong magnetic fields

The production of electron-positron pairs by single photons in magnetic fields 10 to the twelth power G was investigated in detail for photon energies near threshold as well as for the asymptotic limit of high photon energy. The exact attenuation coefficient, which is derived and then evaluated numerically, is strongly influenced by the discrete energy states of the electron and positron. Near threshold, it exhibits a sawtooth pattern as a function of photon energy, and its value is significantly below that predicted by the asymptotic expression for the attenuation coefficient. The energy distributions of the created pair are computed numerically near threshold and analytic expressions are derived in the asymptotic limit. These results indicate that as field strength and photon energy increase, it becomes increasingly probable for the pair to divide the photon energy unequally. This effect, as well as the threshold behavior of the attenuation coefficient, could have important consequences for pulsar models.

Daugherty, J. K.↗

Pair production in superstrong magnetic fields

The production of electron-positron pairs by single photons in magnetic fields of not less than about 10 to the 12th gauss has been investigated in detail for photon energy. The exact attenuation coefficient, which is derived and then evaluated numerically, is strongly influenced by the discrete energy states of the electron and positron. Near threshold, it exhibits a 'sawtooth' pattern as a function of photon energy and its value is significantly below that predicted by the asymptotic expression for the attenuation coefficient. The energy distributions of the created pair are computed numerically near threshold, and analytic expressions are derived in the asymptotic limit. These results indicate that as field strength and photon energy increase, it becomes increasingly probable for the pair to divide strength energy unequally. This effect, as well as the threshold behavior of the attenuation coefficient, could have important consequences to pulsar models.

Daugherty, J. K.↗

Reliability analysis of continuous fiber composite laminates

This paper describes two methods, the maximum distortion energy (MDE) and the principle of independent action (PIA), developed for the analysis of the reliability of a single continuous composite lamina. It is shown that, for the typical laminated plate structure, the individual lamina reliabilities can be combined in order to produce the upper and the lower bounds of reliability for the laminate, similar in nature to the bounds on properties produced from variational elastic methods. These limits were derived for both the interactive and the model failure considerations. Analytical expressions were also derived for the sensitivity of the reliability limits with respect to changes in the Weibull parameters and in loading conditions.

Thomas, David J.↗

Speckle interferometry of asteroids. I - 433 Eros

Analytical expressions are derived for the semimajor and semiminor axes and orientation angle of the ellipse projected by a triaxial asteroid, and the results are applied speckle-interferometry observations of the 433 Eros asteroid. The expressions were calculated as functions of the dimensions and pole of the body and of the asterocentric position of the earth and the sun. On the basis of the analytical expressions, the dimensions of 433 Eros are obtained. The light curve from December 18, 1981 is compared to the dimensions to obtain a geometric albedo of 0.156 (+ or - 0.010). A series of two-dimensional power spectra and autocorrelation functions for 433 Eros show that it is spinning in space.

Drummond, J. D.↗

Hyper Illumination of Exoplanets: Analytical and Numerical Approaches

This work describes the illumination of exoplanets whose orbits are close enough to their host star that the finite angular size of their host star causes hyper illumination, in which more than 50% of the planet receives light. Such exoplanets include the hot Jupiters KELT-9 b (64.5% illuminated) and Kepler-91 b (69.6% illuminated). We describe the geometry of three primary illumination zones: the fully illuminated zone, penumbral zone, and unilluminated zone. The integrals required to determine the incident radiation as a function of position from the substellar point on the exoplanet are explained and derived, and the analytical solution is presented within the fully illuminated zone. We find that the illumination predicted by our model is greater at the substellar point than the typical plane-parallel ray model used would suggest. In addition, it is greater within the region of the penumbral zone extending into the antistellar side of the exoplanet. Finally, we compare our model to that used in starry, an open-source software package used to create albedo maps. It appears that starry may be overestimating the illumination of closely orbiting exoplanets because the foreshortening of the area element of the host star is not included in its calculation.

Exoplanets↗

Cloud-top pressure retrieval with DSCOVR EPIC oxygen A- and B-band observations

An analytic transfer inverse model for Earth Polychromatic Imaging Camera (EPIC) observations is proposed to retrieve the cloud-top pressure (CTP) with the consideration of in-cloud photon penetration. In this model, an analytic equation was developed to represent the reflection at the top of the atmosphere from above cloud, in cloud, and below cloud. The coefficients of this analytic equation can be derived from a series of EPIC simulations under different atmospheric conditions using a nonlinear regression algorithm. With estimated cloud pressure thickness, the CTP can be retrieved from EPIC observation data by solving the analytic equation. To simulate the EPIC measurements, a program package using the double-k approach was developed. Compared to line-by-line calculation, this approach can calculate high-accuracy results with a 100-fold computation time reduction. During the retrieval processes, two kinds of retrieval results, i.e., baseline CTP and retrieved CTP, are provided. The baseline CTP is derived without considering in-cloud photon penetration, and the retrieved CTP is derived by solving the analytic equation, taking into consideration in-cloud and below-cloud interactions. The retrieved CTPs for the oxygen A and B bands are smaller than their related baseline CTP. At the same time, both baseline CTP and retrieved CTP at the oxygen B band are larger than those at the oxygen A band. Compared to the difference in baseline CTP between the B band and A band, the difference in retrieved CTP between these two bands is generally reduced. Out of around 10 000 cases, in retrieved CTP between the A and B bands we found an average bias of 93 mb with a standard deviation of 81 mb. The cloud layer top pressure from Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) measurements is used for validation. Under single-layer cloud situations, the retrieved CTPs for the oxygen A band agree well with the CTPs from CALIPSO, the mean difference of which within 5 mb in the case study. Under multiple-layer cloud situations, the CTPs derived from EPIC measurements may be larger than the CTPs of high-level thin clouds due to the effect of photon penetration.

DSCOVR↗

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗

Theoretical investigation of wave-vector-dependent analytical and numerical formulations of the interband impact-ionization transition rate for electrons in bulk silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E. Scholl, and M. Rudan (1993), and a more exact, numerical method of Y. Wang and K. F. Brennan (1994). Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formular, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with the completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to the overlap intergrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such the transition rate may be overestimated at those points for which 'near vertical,' small change in k, transitions occur.

Kolnik, Jan↗

Orbit-averaged behavior of magnetic control laws for momentum unloading

Analytical formulas are derived for orbit-averaged behavior of magnetic control laws for unloading the excess angular momentum of a spacecraft reaction wheel control system in the presence of secular environmental torques. The specific example of an axially symmetric spacecraft with an inertially fixed attitude for which the dominant environmental torque is the gravity-gradient torque is treated in detail, but extensions of the general approach to other inertially fixed and earth-pointing spacecraft are discussed. The analytical formulas are compared to detailed simulations performed for the Solar Maximum Mission spacecraft, and agreement to within 10% is found. The analytical formulas can be used in place of detailed simulations for preliminary studies, and can be used to find selected cases giving the most stringent tests of momentum unloading capability for which detailed simulations may be performed.

Camillo, P. J.↗

Design Considerations for Lightweight Space Radiators Based on Fabrication and Test Experience With a Carbon-Carbon Composite Prototype Heat Pipe

This report discusses the design implications for spacecraft radiators made possible by the successful fabrication and proof-of-concept testing of a graphite-fiber-carbon-matrix composite (i.e., carbon-carbon (C-C)) heat pipe. The prototype heat pipe, or space radiator element, consists of a C-C composite shell with integrally woven fins. It has a thin-walled furnace-brazed metallic (Nb-1%Zr) liner with end caps for containment of the potassium working fluid. A short extension of this liner, at increased wall thickness beyond the C-C shell, forms the heat pipe evaporator section which is in thermal contact with the radiator fluid that needs to be cooled. From geometric and thermal transport properties of the C-C composite heat pipe tested, a specific radiator mass of 1.45 kg/sq m can be derived. This is less than one-fourth the specific mass of present day satellite radiators. The report also discusses the advantage of segmented space radiator designs utilizing heat pipe elements, or segments, in their survivability to micrometeoroid damage. This survivability is further raised by the use of condenser sections with attached fins, which also improve the radiation heat transfer rate. Since the problem of heat radiation from a fin does not lend itself to a closed analytical solution, a derivation of the governing differential equation and boundary conditions is given in appendix A, along with solutions for rectangular and parabolic fin profile geometries obtained by use of a finite difference computer code written by the author.

Juhasz, Albert J.↗

Design Considerations for Lightweight Space Radiators Based on Fabrication and Test Experience with a Carbon-Carbon Composite Prototype Heat Pipe

This report discusses the design implications for spacecraft radiators made possible by the successful fabrication and Proof-of-concept testing of a graphite-fiber-carbon-matrix composite (i.e., carbon-carbon (C-C)) heat pipe. The proto-type heat pipe, or space radiator element, consists of a C-C composite shell with integrally woven fins. It has a thin-walled furnace-brazed metallic (Nb-1%Zr) liner with end caps for containment of the potassium working fluid. A short extension of this liner, at increased wall thickness beyond the C-C shell, forms the heat pipe evaporator section which is in thermal contact with the radiator fluid that needs to be cooled. From geometric and thermal transport properties of the C-C composite heat pipe tested, a specific radiator mass of 1.45 kg/m2 can be derived. This is less than one-fourth the specific mass of present day satellite radiators. The report also discusses the advantage of segmented space radiator designs utilizing heat pipe elements, or segments, in their survivability to micro-meteoroid damage. This survivability is further raised by the use of condenser sections with attached fins, which also improve the radiation heat transfer rate. Since the problem of heat radiation from a fin does not lend itself to a closed analytical solution, a derivation of the governing differential equation and boundary conditions is given in appendix A, along with solutions for rectangular and parabolic fin profile geometries obtained by use of a finite difference computer code written by the author.

Juhasz, Albert J.↗

The effect of wall interference upon the aerodynamic characteristics of an airfoil spanning a closed-throat circular wind tunnel

The results of a theoretical and experimental investigation of wall interference for an airfoil spanning a closed-throat circular wind tunnel are presented. Analytical equations are derived which relate the characteristics of an airfoil in the tunnel at subsonic speeds with the characteristics in free air. The analysis takes into consideration the effect of fluid compressibility and is based upon the assumption that the chord of the airfoil is small as compared with the diameter of the tunnel. The development is restricted to an untwisted, constant-chord airfoil spanning the middle of the tunnel. Brief theoretical consideration is also given to the problem of choking at high speeds. Results are then presented of tests to determine the low-speed characteristics of an NACA 4412 airfoil for two chord-diameter ratios. While, on the basis of these experiments, no appraisal is possible of the accuracy of the corrections at high speeds, the data indicate that at low Mach numbers the analytical results are valid, even for relatively large values of the chord-diameter ratio.

Vincenti, Walter G↗

Technique for lowering the noise figure in RF amplifiers

Low-level RF amplifiers are operated at noise figures at or near the theoretical minimum, while maintaining high-gain stability. Formulas for the minimum amplifier noise figure value are derived from the analytical noise-figure expression involving the parameters of the equivalent circuit representing the actual system.

Bellomo, A. F.↗

Pyrotechnic hazards classification and evaluation program. Phase 2, segment 3: Test plan for determining hazards associated with pyrotechnic manufacturing processes

A comprehensive test plan for determining the hazards associated with pyrotechnic manufacturing processes is presented. The rationale for each test is based on a systematic analysis of historical accounts of accidents and a detailed study of the characteristics of each manufacturing process. The most hazardous manufacturing operations have been determined to be pressing, mixing, reaming, and filling. The hazard potential of a given situation is evaluated in terms of the probabilities of initiation, communication, and transition to detonation (ICT). The characteristics which affect the ICT probabilities include the ignition mechanisms which are present either in normal or abnormal operation, the condition and properties of the pyrotechnic material, and the configuration of the processing equipment. Analytic expressions are derived which describe the physical conditions of the system, thus permitting a variety of processes to be evaluated in terms of a small number of experiments.

Source record↗

Ion acoustic waves in a multi-ion plasma.

An exact treatment of the multispecies ion acoustic dispersion relation is given for an argon/helium plasma. Phase velocity and damping are obtained as a function of ion-electron temperature ratio and relative densities of the two species. There are two important modes in the plasma, with quite different phase velocities, which are referred to as principal heavy ion mode and principal light ion mode. Which of these is dominant depends on the relative densities of the two components, but, in general, the light ion mode becomes important for surprisingly small light ion contamination. Approximate analytic expressions are derived from damping rates and phase velocities and their domains of validity are investigated. Relevance of the results for the investigation of collisionless shocks is discussed.

Fried, B. D.↗

Minimum weight passive insulation requirements for hypersonic cruise vehicles.

Analytical solutions are derived for two representative cases of the transient heat conduction equation to determine the minimum weight requirements for passive insulation systems of hypersonic cruise vehicles. The cases discussed are the wet wall case with the interior wall temperature held to that of the boiling point of the fuel throughout the flight, and the dry wall case where the heat transferred through the insulation is absorbed by the interior structure whose temperature is allowed to rise.

Ardema, M. D.↗

Transient ion neutralization by electrons.

The nonlinear initial-boundary-value problems describing the lateral neutralization of ion beams for the cases that (1) an auxiliary electric field accelerates the electrons into the ion space, and (2) the electrons are injected into the ion space at a prescribed current density are treated. Analytical solutions are derived which give the position and speed of the neutralization front as a function of time, and the temporal development of the electron density, velocity, and electric fields during the neutralization process.

Wilhelm, H. E.↗