Search NASA⌕ Search

SEARCH · Search NASA

Results for “partial derivative”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Thrust law effects on the longitudinal stability of hypersonic cruise

In the present analytical treatment of the cruising flight stability of an aerospacecraft in nearly circular orbit, on the basis of a thrust law which arbitrarily depends on altitude, speed, and angle of attack, attention is given to thrust law effects on the translational (height and phugoid) and rotational (angle of attack) modes. The partial derivatives of the propulsive forces in conjunction with the aerodynamic forces, with respect to speed and altitude, are noted to exert a major influence on the stability of translational dynamics; the partial derivative of the component of propulsive and aerodynamic forces which is perpendicular to vehicle velocity (with respect to the angle of attack) is a primary determinant of the damping of angle-of-attack oscillations, while the partial derivative of the sum about the center of mass of aerodynamic and propulsive pitching moments, with respect to the angle of attack, determines the corresponding period.

Markopoulos, Nikos↗
NASA NTRS2007

Influences of Local Sea-Surface Temperatures and Large-scale Dynamics on Monthly Precipitation Inferred from Two 10-year GCM-Simulations

Two parallel sets of 10-year long: January 1, 1982 to December 31, 1991, simulations were made with the finite volume General Circulation Model (fvGCM) in which the model integrations were forced with prescribed sea-surface temperature fields (SSTs) available as two separate SST-datasets. One dataset contained naturally varying monthly SSTs for the chosen period, and the oth& had the 12-monthly mean SSTs for the same period. Plots of evaporation, precipitation, and atmosphere-column moisture convergence, binned by l C SST intervals show that except for the tropics, the precipitation is more strongly constrained by large-scale dynamics as opposed to local SST. Binning data by SST naturally provided an ensemble average of data contributed from disparate locations with same SST; such averages could be expected to mitigate all location related influences. However, the plots revealed: i) evaporation, vertical velocity, and precipitation are very robust and remarkably similar for each of the two simulations and even for the data from 1987-ENSO-year simulation; ii) while the evaporation increased monotonically with SST up to about 27 C, the precipitation did not; iii) precipitation correlated much better with the column vertical velocity as opposed to SST suggesting that the influence of dynamical circulation including non-local SSTs is stronger than local-SSTs. The precipitation fields were doubly binned with respect to SST and boundary-layer mass and/or moisture convergence. The analysis discerned the rate of change of precipitation with local SST as a sum of partial derivative of precipitation with local SST plus partial derivative of precipitation with boundary layer moisture convergence multiplied by the rate of change of boundary-layer moisture convergence with SST (see Eqn. 3 of Section 4.5). This analysis is mathematically rigorous as well as provides a quantitative measure of the influence of local SST on the local precipitation. The results were recast to examine the dependence of local rainfall on local SSTs; it was discernible only in the tropics. Our methodology can be used for computing relationship between any forcing function and its effect(s) on a chosen field.

Sud, Y. C.↗
NASA NTRS2007

Scope of Gradient and Genetic Algorithms in Multivariable Function Optimization

Global optimization of a multivariable function - constrained by bounds specified on each variable and also unconstrained - is an important problem with several real world applications. Deterministic methods such as the gradient algorithms as well as the randomized methods such as the genetic algorithms may be employed to solve these problems. In fact, there are optimization problems where a genetic algorithm/an evolutionary approach is preferable at least from the quality (accuracy) of the results point of view. From cost (complexity) point of view, both gradient and genetic approaches are usually polynomial-time; there are no serious differences in this regard, i.e., the computational complexity point of view. However, for certain types of problems, such as those with unacceptably erroneous numerical partial derivatives and those with physically amplified analytical partial derivatives whose numerical evaluation involves undesirable errors and/or is messy, a genetic (stochastic) approach should be a better choice. We have presented here the pros and cons of both the approaches so that the concerned reader/user can decide which approach is most suited for the problem at hand. Also for the function which is known in a tabular form, instead of an analytical form, as is often the case in an experimental environment, we attempt to provide an insight into the approaches focusing our attention toward accuracy. Such an insight will help one to decide which method, out of several available methods, should be employed to obtain the best (least error) output. *

Shaykhian, Gholam Ali↗
NASA NTRS1994

The 630 nm MIG and the vertical neutral wind in the low latitude nighttime thermosphere

It is shown that large negative divergences (gradients) in the horizontal neutral wind in the equatorial thermosphere can support downward neutral winds in excess of 20 m/s. With attention to the meridional and vertical winds only, the pressure tendency equation is used to derive the expression U(sub z0) approximately equals (Partial derivative U(sub y)/Partial derivative y)H for the vertical wind U(sub z0) at the reference altitude for the pressure tendency equation; H is the atmospheric density scale height, and (Partial derivative U(sub y)/Partial derivative y) is the meridional wind gradient. The velocity gradient associated with the Meridional Intensity Gradient (MIG) of the O((sup 1)D) emission (630 nm) at low latitudes is used to estimate the vertical neutral wind in the MIG region. Velocity gradients derived from MIG data are about 0.5 (m/s)/km) or more, indicating that the MIG region may contain downward neutral winds in excess of 20 m/s. Though direct measurements of the vertical wind are scarce, Fabry-Perot interferometer data of the equatorial F-region above Natal, Brazil, showed downward winds of 30 m/s occurring during a strong meridional wind convergence in 1982. In-situ measurements with the WATS instrument on the DE-2 satellite also show large vertical neutral winds in the equatorial region.

Herrero, F. A.↗
NASA NTRS1981

Influence of velocity gradients on measurements of velocity and streamwise vorticity with hot-wire X-array probes

An analysis and measurement of the effects of the streamwise velocity gradients partial derivative of U with respect to y and partial derivative of U with respect to z, on the velocity components, U, v, and w, and the streamwise vorticity component, omega sub x measured in turbulent flow with a pair of orthogonal hot-wire X arrays, is presented. It is shown that these gradients, which can have the same order of magnitude instantaneously as the mean shear stress at the wall, cause extremely large errors in the measured instantaneous cross-stream velocity and streamwise vorticity components.

Vukoslavcevic, P.↗
NASA NTRS2017

Mapping of the Resistance of a Superconducting Transition Edge Sensor as a Function of Temperature, Current, and Applied Magnetic Field

We have measured the resistance R (T, I, B(sub ext) of a superconducting transition edge sensor over the entire transition region on a fine scale, producing a four-dimensional map of the resistance surface. The dimensionless temperature and current sensitivities (alpha equivalence partial derivative log R/partial derivative log T|(sub I) and beta equivalence partial derivative log R/partial derivative log I|(sub T) of the TES resistance have been determined at each point. alpha and beta are closely related to the sensor performance, but show a great deal of complex, large amplitude fine structure over large portions of the surface that is sensitive to the applied magnetic field. We discuss the relation of this structure to the presence of Josephson weak link fringes.

Zhang, Shou↗
NASA NTRS1989

Kinematic sensitivity of robot manipulators

Kinematic sensitivity vectors and matrices for open-loop, n degrees-of-freedom manipulators are derived. First-order sensitivity vectors are defined as partial derivatives of the manipulator's position and orientation with respect to its geometrical parameters. The four-parameter kinematic model is considered, as well as the five-parameter model in case of nominally parallel joint axes. Sensitivity vectors are expressed in terms of coordinate axes of manipulator frames. Second-order sensitivity vectors, the partial derivatives of first-order sensitivity vectors, are also considered. It is shown that second-order sensitivity vectors can be expressed as vector products of the first-order sensitivity vectors.

Vuskovic, Marko I.↗
NASA NTRS1994

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗
NASA NTRS1985

Navigation systems

The elements of the measurement and communications network comprising the global deep space navigation system (DSN) for NASA missions are described. Among the measurement systems discussed are: VLBI, two-way Doppler and range measurements, and optical measurements carried out on board the spacecraft. Processing of navigation measurement is carried out using two modules: an N-body numerical integration of the trajectory (and state transition partial derivatives) based on pre-guessed initial conditions; and partial derivatives of simulated observables corresponding to each actual observation. Calculations of velocity correction parameters is performed by precise modelling of all physical phenomena influencing the observational measurements, including: planetary motions; tracking station locations, gravity field structure, and transmission media effects. Some of the contributions to earth-relative orbit estimate errors for the Doppler/range system on board Voyager are discussed in detail. A line drawing of the DSN navigation system is provided.

Jordan, J. F.↗
NASA NTRS1989

A direct procedure for interpolation on a structured curvilinear two-dimensional grid

A direct procedure is presented for locally bicubic interpolation on a structured, curvilinear, two-dimensional grid. The physical (Cartesian) space is transformed to a computational space in which the grid is uniform and rectangular by a generalized curvilinear coordinate transformation. Required partial derivative information is obtained by finite differences in the computational space. The partial derivatives in physical space are determined by repeated application of the chain rule for partial differentiation. A bilinear transformation is used to analytically transform the individual quadrilateral cells in physical space into unit squares. The interpolation is performed within each unit square using a piecewise bicubic spline.

Zingg, David W.↗
NASA NTRS1972

Closed-form solutions for atomspheric flight with applications to shuttle guidance

Closed-form solutions for the motion of a rocket-powered vehicle during atmospheric ascent and closed-form solutions for unpowered atmospheric reentry are developed. These closed-form solutions are then used to develop a simplified guidance scheme and to develop a variation-of-parameters integration of more accurate equations of motion with the closed-form solutions as base solutions. The variation-of-parameters integration of the more accurate equations of motion also allows the transition partial derivative matrices associated with these equations to be easily developed. Then the partial derivative transition matrices are used to develop a guidance scheme based on the more accurate equations of motion instead of the less accurate closed-form solutions.

Ingram, H. L.↗
NASA NTRS2000

Detached and Floating-Zone Growth of Semiconductor Crystals on the ISS

Understanding the mechanism of detached Bridgman growth and establishing the growth of large scale germanium-silicon crystals by the float-zone technique are the key points of the project "RDGS - Reduction of Defects in Germanium-Silicon". The contact angle of the melt and the growth angle of the crystal are essential parameters which allow a controlled use of detached growth. The contact angle was determined for a variety of different substrates and melt compositions; pBN showed the highest value for pure germanium as well as for germanium-rich GeSi melts. The growth angle of Ge(sub 0.95) Si(sub 0.05) was measured to be 8.5-10.5 degrees which concurs with the values of pure germanium and silicon, respectively. The temperature dependence and the concentration dependence of the surface tension were determined for concentrations up to 10at% silicon (partial derivative (gamma)/partial derivative T=-0.08 (raised dot) 10(exp -3)N/m (raised dot) K, partial derivative (gamma)/partial derivative (C)=2.2 (raised dot) 10(exp -3)N/m (raised dot) at%). Using these values, the critical Marangoni number indicating the transition to time-dependent thermocapillary flow will be exceeded for the growth of large scale float-zone crystals onboard the ISS. Therefore, suitable tools for flow control are required.

Dold, P.↗
NASA NTRS1986

Geometerial description for a proposed aeroassist flight experiment vehicle

One geometry currently under consideration for the Aeroassist Flight Experiment (AFE) vehicle is composed of several segments of simple general conics: an ellipsoidal nose tangent to an elliptical cone and a base skirt with the base plane raked relative to the body axis. An analytic representation for the body coordinates and first and second partial derivatives of this configuration has been developed. Equations are given which define the body radius and partial derivatives for a prescribed axial and circumferential position on the vehicle. The results for a sample case are tabulated and presented graphically.

Cheatwood, F. M.↗
NASA NTRS1994

Motion models in attitude estimation

Attitude estimator use observations from different times to reduce the effects of noise. If the vehicle is rotating, the attitude at one time needs to be propagated to that at another time. If the vehicle measures its angular velocity, attitude propagating entails integrating a rotational kinematics equation only. If a measured angular velocity is not available, torques can be computed and an additional rotational dynamics equation integrated to give the angular velocity. Initial conditions for either of these integrations come from the estimation process. Sometimes additional quantities, such as gyro and torque parameters, are also solved for. Although the partial derivatives of attitude with respect to initial attitude and gyro parameters are well known, the corresponding partial derivatives with respect to initial angular velocity and torque parameters are less familiar. They can be derived and computed numerically in a way that is analogous to that used for the initial attitude and gyro parameters. Previous papers have demonstrated the feasibility of using dynamics models for attitude estimation but have not provided details of how each angular velocity and torque parameters can be estimated. This tutorial paper provides some of that detail, notably how to compute the state transition matrix when closed form expressions are not available. It also attempts to put dynamics estimation in perspective by showing the progression from constant to gyro-propagated to dynamics-propagated attitude motion models. Readers not already familiar with attitude estimation will find this paper an introduction to the subject, and attitude specialists may appreciate the collection of heretofore scattered results brought together in a single place.

Chu, D.↗
DOE OSTI2026

Approximation of refrigerant thermophysical properties using neural networks to speed up transient thermofluid simulations

Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.

Ma, Jiacheng↗
NASA NTRS1989

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗
NASA NTRS1972

Quadrature errors in the partical derivatives required for the direct recovery of gravity anomalies from satellite observations

The equations of motion of a geodetic satellite in the earth's gravitational field expressed by gravity anomalies require the evaluation, amongst others, of the partial derivatives of the disturbing force with respect to individual gravity anomalies. Data are discussed on how anomaly blocks should be subdivided so that the partial derivatives may be numerically evaluated for each subdivision, and then finally meaned to give the value representative of the whole blocks, with accuracies better than 2 to 3 percent for all blocks. The number of subdivisions is large for the blocks nearest to the satellite subpoint and decreases away from it. The actual values of this spherical distance and the actual subdivision of the mean gravity anomaly blocks was determined numerically for 184 15 deg x 15 deg equal area blocks. Satellite heights above the earth of 400 km, 800 km and 1600 km were considered. The computer times for the suggested scheme were compared with alternative solutions.

Hajela, D. P.↗