Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial differential equation”

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.

1,139 records · Page 64

A new approach to radiative transfer theory using Jones's vectors. I Fundamentals

Radiative transfer of partially polarized radiation in an anisotropically scattering, inhomogeneous atmosphere containing an arbitrary polydispersion of particles is described using Jones's (1941) amplitude vectors and matrices. This approach exploits the close analogy between the quantum-mechanical states of spin 1/2 systems and the polarization states of electromagnetic radiation described by Jones's vector. The complete equivalence between the transport equation for Jones's vectors and the classical radiative-transfer equation for Stokes's intensity vectors is demonstrated in two independent ways after deriving the transport equations for the polarization coherency matrices and for the quaternions corresponding to the Jones's vectors. A compact-operator formulation of the theory is provided and used to derive the necessary equations for both a local and a global description of the transport of Jones's vectors. The integro-differential equations for the amplitude reflection and transmission matrices are derived and related to the usual corresponding equations.

Fymat, A. L.↗

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗

Comparing Consider-Covariance Analysis with Sigma-Point Consider Filter and Linear-Theory Consider Filter Formulations

Recent literature in applied estimation theory reflects growing interest in the sigma-point (also called unscented ) formulation for optimal sequential state estimation, often describing performance comparisons with extended Kalman filters as applied to specific dynamical problems [c.f. 1, 2, 3]. Favorable attributes of sigma-point filters are described as including a lower expected error for nonlinear even non-differentiable dynamical systems, and a straightforward formulation not requiring derivation or implementation of any partial derivative Jacobian matrices. These attributes are particularly attractive, e.g. in terms of enabling simplified code architecture and streamlined testing, in the formulation of estimators for nonlinear spaceflight mechanics systems, such as filter software onboard deep-space robotic spacecraft. As presented in [4], the Sigma-Point Consider Filter (SPCF) algorithm extends the sigma-point filter algorithm to the problem of consider covariance analysis. Considering parameters in a dynamical system, while estimating its state, provides an upper bound on the estimated state covariance, which is viewed as a conservative approach to designing estimators for problems of general guidance, navigation and control. This is because, whether a parameter in the system model is observable or not, error in the knowledge of the value of a non-estimated parameter will increase the actual uncertainty of the estimated state of the system beyond the level formally indicated by the covariance of an estimator that neglects errors or uncertainty in that parameter. The equations for SPCF covariance evolution are obtained in a fashion similar to the derivation approach taken with standard (i.e. linearized or extended) consider parameterized Kalman filters (c.f. [5]). While in [4] the SPCF and linear-theory consider filter (LTCF) were applied to an illustrative linear dynamics/linear measurement problem, in the present work examines the SPCF as applied to nonlinear sequential consider covariance analysis, i.e. in the presence of nonlinear dynamics and nonlinear measurements. A simple SPCF for orbit determination, exemplifying an algorithm hosted in the guidance, navigation and control (GN&C) computer processor of a hypothetical robotic spacecraft, was implemented, and compared with an identically-parameterized (standard) extended, consider-parameterized Kalman filter. The onboard filtering scenario examined is a hypothetical spacecraft orbit about a small natural body with imperfectly-known mass. The formulations, relative complexities, and performances of the filters are compared and discussed.

Lisano, Michael E.↗

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗

Small Particle Response to Fluid Motion Using Tethered Particles to Simulate Microgravity

This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Through an innovative application of fractional differential equations, two members of our team produced the first analytical solution of a fundamental equation of motion, which had only been solved numerically or by approximation before. The general solution predicts that the usually neglected history term becomes important in particle response to a sinusoidal fluid movement when the characteristic viscous time is in the same order as the fluid oscillation period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We have developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the virtual apparatus for on-earth microgravity experimentation. We can quantify precisely the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. We are examining the response of particles up to 2 mm radius to fluid oscillation at frequencies up to 80 Hz with amplitudes up to 200 microns. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle to within microns during recording, and techniques for achieving this are one of the project challenges. Focused microscopic images and diffraction patterns are used. To make the experiment more versatile, the spaceflight system will record holograms both on film and electronically. A cross correlation procedure enables sub pixel accuracies for electronic recordings, partially accommodating the lower spatial resolution of CCDs. The electronic holograms can be down linked providing real time data. Results of the ground experiments, the flight experiment design, and data analysis procedures are reported.

Trolinger, James↗