Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta 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 109 records · Page 6

Cure kinetics of epoxy matrix resin by differential scanning calorimetry

A study was made on the cure kinetics of an epoxy neat-resin (Narmco 5208) using Differential Scanning Calorimetry (DSC). Two interrelated analytical methods were applied to dynamic DSC data for evaluating the kinetic parameters, such as activation energy, E, the order of reaction, n, and the total heat of polymerization (or crosslinking), delta H sub t. The first method was proposed by Ellerstein (1968), and uses a thorough differential-integral analysis of a single DSC curve to evaluate the kinetic parameters. The second method was proposed by Kissinger (1957), and uses multiple DSC curves obtained at various heating rates to evaluate E regardless of n. Kinetic analysis of Narmco 5208 epoxy resin showed that the reaction order, n, is substantially affected by the rate of heating; i.e., n is approximately 2 at slow scan rates but is reduced to 1.5 at higher scan rates. The activation energy, E, is not affected by the scan rate, and the average value of E is 25.6 + or - 1.8 kcal/mole.

Cizmecioglu, M.↗

Design considerations for large space electric power systems

As power levels of spacecraft rise to the 50 to 100 kW range, it becomes apparent that low voltage (28 V) dc power distribution and management systems will not operate efficiently at these higher power levels. The concept of transforming a solar array voltage at 150 V dc into a 1000 V ac distribution system operating at 20 kHz is examined. The transformation is accomplished with series-resonant inverter by using a rotary transformer to isolate the solar array from the spacecraft. The power can then be distributed in any desired method such as three phase delta to delta. The distribution voltage can be easily transformed to any desired load voltage and operating frequency. The reasons for the voltage limitations on the solar array due to plasma interactions and the many advantages of a high voltage, high frequency at distribution system are discussed.

Renz, D. D.↗

Integration of Libration Point Orbit Dynamics into a Universal 3-D Autonomous Formation Flying Algorithm

The autonomous formation flying control algorithm developed by the Goddard Space Flight Center (GSFC) for the New Millennium Program (NMP) Earth Observing-1 (EO-1) mission is investigated for applicability to libration point orbit formations. In the EO-1 formation-flying algorithm, control is accomplished via linearization about a reference transfer orbit with a state transition matrix (STM) computed from state inputs. The effect of libration point orbit dynamics on this algorithm architecture is explored via computation of STMs using the flight proven code, a monodromy matrix developed from a N-body model of a libration orbit, and a standard STM developed from the gravitational and coriolis effects as measured at the libration point. A comparison of formation flying Delta-Vs calculated from these methods is made to a standard linear quadratic regulator (LQR) method. The universal 3-D approach is optimal in the sense that it can be accommodated as an open-loop or closed-loop control using only state information.

Folta, David↗

Robust Trajectory Optimization for NRHO Rendezvous Using SPICE Kernel Relative Motion

In this paper, robust optimization is performed on trajectory correction maneuvers during the lunar lander return phase of an Artemis mission, treating the trajectory from one hour after low lunar orbit departure to arrival in the vicinity of the lunar Gateway as a relative motion problem. To enable rapid stochastic optimization techniques requiring many candidate trajectories, SPICE kernel relative motion as implemented by the Quadratic Interpolated State Transition (QIST) system is used as the underlying dynamics propagation. The optimization is performed with a genetic optimizer using linear covariance (LinCov) software in a simplified operational context, taking into account the availability of navigation sensors with varying measurement models, ranges, and accuracies. No numerical integration is used, since the relative motion around Gateway is fully characterized with the a priori computation of the QIST coefficients. Maneuver placements are computed to optimize the minimum 3σ delta-v of the trajectory, the position dispersion at a target point, and a convex combination of these two metrics. An order of magnitude runtime improvement is provided over legacy methods with less than 10% error introduced. All QIST results are shown to be in-family with legacy methods. The tradespace for optimal delta-v design is found to range from 77.0 to 93.9 m/s, while the range of optimal dispersion is between 1.4 and 11.7 km.

Relative Motion↗

Adaptive computations of flow around a delta wing with vortex breakdown

An adaptive unstructured mesh solution method for the three-dimensional Euler equations was used to simulate the flow around a sharp edged delta wing. Emphasis was on the breakdown of the leading edge vortex at high angle of attack. Large values of entropy, which indicate vortical regions of the flow, specified the region in which adaptation was performed. The aerodynamic normal force coefficients show excellent agreement with wind tunnel data measured by Jarrah, and demonstrate the importance of adaptation in obtaining an accurate solution. The pitching moment coefficient and the location of vortex breakdown are compared with experimental data measured by Hummel and Srinivasan, showing good agreement in cases in which vortex breakdown is located over the wing.

Modiano, David L.↗

A Globally Stable Lyapunov Pointing and Rate Controller for the Magnetospheric MultiScale Mission (MMS)

The Magnetospheric MultiScale Mission (MMS) is scheduled to launch in late 2014. Its primary goal is to discover the fundamental plasma physics processes of reconnection in the Earth's magnetosphere. Each of the four MMS spacecraft is spin-stabilized at a nominal rate of 3 RPM. Traditional spin-stabilized spacecraft have used a number of separate modes to control nutation, spin rate, and precession. To reduce the number of modes and simplify operations, the Delta-H control mode is designed to accomplish nutation control, spin rate control, and precession control simultaneously. A nonlinear design technique, Lyapunov's method, is used to design the Delta-H control mode. A global spin rate controller selected as the baseline controller for MMS, proved to be insufficient due to an ambiguity in the attitude. Lyapunov's design method was used to solve this ambiguity, resulting in a controller that meets the design goals. Simulation results show the advantage of the pointing and rate controller for maneuvers larger than 90 deg and provide insight into the performance of this controller.

Shah, Neerav↗

Theoretical and experimental longitudinal aerodynamic characteristics of an aspect ratio 0.25 sharp-edge delta wing at subsonic, supersonic, and hypersonic speeds

The suction analogy concept of Polhamus for predicting vortex lift in conjunction with an appropriate potential-flow solution is called the present method. This method is applied to an aspect ratio 0.25 sharp-edge delta wing from a Mach number of 0.143 to 10.4 in free air and at 0.074 in ground effect, and also to an aspect ratio 0.35 triangular crosssectional body at a Mach number of 6.9. The models had subsonic leading edges at the test Mach numbers. Vortex-flow effects could be neither confirmed nor denied to exist at high speeds because of the lack of flow visualization above a Mach number of 0.143. The data, however, could be better predicted by including a vortex-flow effect, although not always to the extent predicted from the present method because of the presence of actual and hypothesized unmodeled flow situations. The method of Nenni and Tung (NASA CR-1860) tended to confirm the existence of vortex flow at hypersonic speeds. The hypersonic-tangent-cone method predicted best the delta-wing results over the test angle-of-attack range and hypersonic Mach number range and did equally as well as the present method for the triangular body.

Fox, C. H., Jr.↗

Streamlined Convergence Acceleration for CFD Codes

Enigma, a simplified interface to the PETSc library, is shown to enable the rapid solution of discrete partial differential equations. Two CFD codes, LAURA and HyperSolve, use Enigma to compute steady solutions of the Navier-Stokes equations. Using PETSc, Enigma is shown to provide a Jacobian-Free Newton-Krylov method (JFNK), globalized with pseudotransient continuation, that improves efficiency over the point-implicit relaxation method traditionally used by LAURA. It is shown that iterative error has a large impact on surface heat transfer predicted by LAURA on an axisymmetric sphere-cone geometry. Also, the convergence rate of HyperSolve simulating subsonic flow over a delta wing geometry with the JFNK method is shown to be more efficient than employing a defect correction method as the nonlinear solver.

Thompson, Kyle B.↗

Accurate simulation of vortical flows

Even with recent supercomputers having a large memory, Navier-Stokes simulations for vortical flows do not provide satisfactory results because of the lack of grid resolution to accurately simulate the strength of separation vortices. To overcome this problem, a zonal method is proposed to increase the number of grid points locally. Interface scheme which is critical for an efficient and stable zonal method is based on the Fortified Navier-Stokes concept. Application to both two-dimensional conical and three-dimensional delta wing problems indicates this simple zonal method can improve the accuracy of vortical flow simulations.

Fujii, Kozo↗

Long Term Measurement of the Vapor Pressure of Gold in the Au-C System

Incorporating the {Au(s,l) + graphite} reference in component activity measurements made with the multiple effusion-cell vapor source mass spectrometry (multicell KEMS) technique provides a fixed temperature defining ITS-90 (T(sub mp)(Au) = 1337.33K) and a systematic method to check accuracy. Over a 2 year period delta H sub(298)Au was determined by the 2nd and 3rd law methods in 25 separate experiments and were in the ranges 362.2 plus or minus 3.3 kJmol(sup -1) and 367.8 plus or minus 1.1 kJmol(sup -1), respectively. This 5 kJmol-1 discrepancy is transferred directly to the measured activities. This is unacceptable and the source of this discrepancy needs to be understood and corrected. Accepting the 2nd law value increases p(Au) by about 50 percent, brings the 2nd and 3rd law values into agreement and removes the T dependence in the 3rd law values. While compelling, there is no way to independently determine instrument sensitivities, S(sub Au), with T in a single experiment with KEMS. This lack of capability is stopping a deeper understanding of this problem. In addition, the Au-C phase diagram suggests a eutectic invariant reaction: L-Au(4.7at%C) = FCC-Au(0.08at%C) + C(graphite) at T(sub e) approximately 1323K. This high C concentration in Au(l) must reduce p(Au) in equilibrium with {Au(s,l) + graphite} and raises some critical questions about the Gibbs free energy functions of Au(s,l) and the Au fixed point (T(sub mp)(Au) = 1337.33K) which is always measured in graphite.

Copland, Evan H.↗

Using the radiative temperature difference at 3.7 and 11 microns to track dust outbreaks

The radiative temperature difference between 3.7 and 11 microns (Delta T) is investigated as a possible method for tracking dust outbreaks. Theoretical calculations indicate that the technique would be most sensitive to dust loading during the day when the Delta Ts are enhanced by reflection of 3.7-micron solar radiation. The feasibility of tracking dust outbreaks is demonstrated by comparing satellite observations with surface estimates of visibility. Theoretical calculations also demonstrate the potential of inferring the dust layer optical depth from these spectral measurements.

Ackerman, Steven A.↗

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" (DR) method is introduced for analyzing lightning field changes (DELTAEs) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an Overdetermined Fixed Matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all 8 parameters (x, y, z, Q, x', y', z', Q') are found, yet a numerical search of only 4 parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused DeltaEs derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the Lightning Detection And Ranging (LDAR) system at the KSC-ER, and from least squares error estimation theory, and the method is shown to be a useful "stand-alone" charge retrieval tool. Since more than one charge distribution describes a finite set of DELTAEs (i.e., solutions are non-unique), and since there can exist appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence, an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions, thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning in several Florida storms.

Koshak, W. J.↗

Numerical and experimental determination of secondary separation on delta wings in subsonic flow

A quasi-three-dimensional boundary layer method for conical inviscid flows is presented. The model was developed to characterize the flow fields over slender delta wings, particularly the shed vortex sheet. A steady, incompressible, laminar boundary layer is assumed and a solution is obtained with the Smith (1966) finite difference code. Predictions are compared with oil flow patterns and surface pressure measurements for 74 and 80 deg delta wings at 0-20 deg angles of attack. The model is capable of accurately predicting the secondary separation lines. Inviscid velocities derived from panel methods lead, however, to inaccurate pressure distributions.

Dejarnette, F. R.↗

An Approximate Method of Calculating the Deformations of Wings Having Swept, M or W, A, and Swept-Tip Plan Forms

An approximate method of calculating the deformations of wings of uniform thickness having swept, M or W, Delta, and swept-tip plan forms is presented. The method employs an adjustment to the elementary beam theory to account for the effect of the triangular root portion of a swept wing on the deformation of the outboard section of the wing. To demonstrate the general applicability of the method, the modified elementary theory is applied to the more complex M or W, Delta, and swept-tip plan forms as well as to swept plan forms. For the purpose of calculating angles of attack, it is shown that the unmodified elementary beam theory applied to that part of the wing outboard of the root triangle produces satisfactory results. However, for calculating deflections it is necessary to include the effects of the root-triangle deformation.

Zender, George W.↗

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" ("DR") method is introduced for analyzing lightning field changes ((Delta)Es) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an overdetermined fixed matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all eight parameters Ix. y, z, Q, x', y', z', Q') are found, yet a numerical search of only four parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused (Delta)Es derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and U.S. Air Force Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the lightning detection and ranging (LDAR) system at the KSC-ER. and from least squares error estimation theory, and the method is shown to be a useful "stand alone" charge retrieval tool. Since more than one charge distribution describes a finite set of (Delta)Es (i.e., solutions are nonunique), and since there can be appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence. an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions. thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning m several Florida storms.

Koshak, W. J.↗

Three dimensional steady and unsteady asymmetric flow past wings of arbitrary planforms

The nonlinear discrete vortex method is extended to treat the problem of asymmetric flows past a wing with leading edge separation, including steady and unsteady flows. The problem is formulated in terms of a body fixed frame of reference and the nonlinear-discrete vortex method is modified accordingly. Although the method is general, only examples of flows past delta wings are presented due to the availability of experimental data as well as approximate theories. Comparison of results with experimental results for a delta wing undergoing a steady rolling motion at zero angle of attack demonstrate the superiority of the present method over existing approximate theories in obtaining highly accurate loads. Numerical results for yawed wings at large angles of attack are also presented. In all cases, total load coefficients, pressure distributions, and shapes of the free vortex sheets are shown.

Kandil, O. A.↗

Lateral aerodynamics of delta wings with leading edge separation

An unsteady vortex lattice method is presented for the calculation of the aerodynamic forces acting on lifting surfaces undergoing complex three dimensional motion. For the present case the nonsymmetric motion of a slender delta wing was considered and the resulting lateral characteristics were calculated. The flow separation line was specified along the wing leading edge and the emanating vortex sheet shape and rollup was then calculated. Numerical results are presented for the combined high angle of attack and side slip condition and for the wing constant roll and coning motions.

Katz, J.↗

A method for computing instantaneous downtrack error from TDRS Doppler residuals

A simple method has been developed which calculates the instantaneous downtrack error in an orbiting vehicle's state vector from TDRS (Tracking and Data Relay Satellite) Doppler observation residuals. The downtrack error is expressed as a delta-time correction to the timetag of the vehicle's state vector. The method requires only a few easily obtained parameters and is easy to implement. This paper details the derivation of computing a delta-time correction from TDRS Doppler observation residuals and evaluates the method's assumptions and limitations. The results presented show that typically the method computes delta-time corrections that are less than five percent off from the true timetag error for Doppler observation residuals within twenty minutes of closest approach to the TDRS.

Haas, Franklin M., Jr.↗