Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Directional 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 55 records · Page 3

A variable-metric algorithm employing linear and quadratic penalties

A variable-metric algorithm is described that uses both linear and quadratic penalty terms for handling nonlinear constraints. Quadratic penalty coefficients are adjusted in a process which maintains a positive-definite matrix of second partial derivatives of the function without generating the large positive eigenvalues which cause zigzagging and slow convergence. The schemes suggested use inferred second-order properties not only in terms of the variable metric of the Davidson-Fletcher-Powell algorithm (or its relatives) but by estimating of second directional derivatives by fitting cubics to various functions along search directions.

Kelley, H. J.

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.

Flight Investigation of Prescribed Simultaneous Independent Surface Excitations for Real-Time Parameter Identification

Near real-time stability and control derivative extraction is required to support flight demonstration of Intelligent Flight Control System (IFCS) concepts being developed by NASA, academia, and industry. Traditionally, flight maneuvers would be designed and flown to obtain stability and control derivative estimates using a postflight analysis technique. The goal of the IFCS concept is to be able to modify the control laws in real time for an aircraft that has been damaged in flight. In some IFCS implementations, real-time parameter identification (PID) of the stability and control derivatives of the damaged aircraft is necessary for successfully reconfiguring the control system. This report investigates the usefulness of Prescribed Simultaneous Independent Surface Excitations (PreSISE) to provide data for rapidly obtaining estimates of the stability and control derivatives. Flight test data were analyzed using both equation-error and output-error PID techniques. The equation-error PID technique is known as Fourier Transform Regression (FTR) and is a frequency-domain real-time implementation. Selected results were compared with a time-domain output-error technique. The real-time equation-error technique combined with the PreSISE maneuvers provided excellent derivative estimation in the longitudinal axis. However, the PreSISE maneuvers as presently defined were not adequate for accurate estimation of the lateral-directional derivatives.

Moes, Timothy R.

Pareto-optimal target definition for multi-axis random vibration testing

In random vibration testing with multiple control channels, existing control laws require specification of a complete spectral density matrix at each control frequency. Spectral density matrices include autospectral densities on the diagonal and cross-spectral densities on the off-diagonal. In practice, the off-diagonal terms are often unknown, and recent vibration testing research has focused on fixing the diagonal and specifying the off-diagonal to minimize the required control energy, subject to a constraint that the target matrix is positive semidefinite. This paper shows that, even with a fixed diagonal, off-diagonal terms strongly affect control residuals. This overlooked effect occurs in both square and rectangular systems. By jointly considering input energy and control residuals, open-loop inputs are derived directly from the diagonal without specifying the off-diagonal terms. Vibration targets that can be used in closed-loop control are then derived using the optimal inputs, with positive semidefinite constraints applied during the derivation. The result is a set of Pareto-optimal control solutions. For each solution in the set, any other possible solution produces greater control error, greater input energy, or both. A balanced solution is selected automatically, though others can be chosen based on test needs. Simulations and experiments show that the proposed method outperforms state-of-the-art energy-minimizing approaches, achieving significant reductions in both control error and input energy.

Autospectral density

Design for prevention of acoustic fatigue

It is pointed out that new noise prediction methods and acoustic life estimation methods have matured to the point where they can be combined into a unified engineering procedure. "Life derivatives" can be extracted from parametric charts to furnish design data for preventing acoustic fatigue. The acoustic fatigue life is shown to be sensitive to the damping ratio through the use of life derivatives. The localized nature of propeller noise can be quantified with an "isodecibel" contour diagram.Even though the peak sound pressure level may be high, the directional derivatives show the noise decay rates with distance. Acoustic fatigue design is discussed from the overall design methodology and is shown to be similar to other structural design problems. While nonlinearities present a formidable design engineering problem, they are manageable by proven semi-empirical techniques. For new design problems, it is imperative to determine whether the data base completely spans the design variables.

Smith, H. W.

Sensitivity of optimum solutions to problem parameters

Derivation of the sensitivity equations that yield the sensitivity derivatives directly, which avoids the costly and inaccurate perturb-and-reoptimize approach, is discussed and solvability of the equations is examined. The equations apply to optimum solutions obtained by direct search methods as well as those generated by procedures of the sequential unconstrained minimization technique class. Applications are discussed for the use of the sensitivity derivatives in extrapolation of the optimal objective function and design variable values for incremented parameters, optimization with multiple objectives, and decomposition of large optimization problems.

Sobieszczanski-Sobieski, J.

Expansion-based passive ranging

A new technique of passive ranging which is based on utilizing the image-plane expansion experienced by every object as its distance from the sensor decreases is described. This technique belongs in the feature/object-based family. The motion and shape of a small window, assumed to be fully contained inside the boundaries of some object, is approximated by an affine transformation. The parameters of the transformation matrix are derived by initially comparing successive images, and progressively increasing the image time separation so as to achieve much larger triangulation baseline than currently possible. Depth is directly derived from the expansion part of the transformation. To a first approximation, image-plane expansion is independent of image-plane location with respect to the focus of expansion (FOE) and of platform maneuvers. Thus, an expansion-based method has the potential of providing a reliable range in the difficult image area around the FOE. In areas far from the FOE the shift parameters of the affine transformation can provide more accurate depth information than the expansion alone, and can thus be used similarly to the way they were used in conjunction with the Inertial Navigation Unit (INU) and Kalman filtering. However, the performance of a shift-based algorithm, when the shifts are derived from the affine transformation, would be much improved compared to current algorithms because the shifts - as well as the other parameters - can be obtained between widely separated images. Thus, the main advantage of this new approach is that, allowing the tracked window to expand and rotate, in addition to moving laterally, enables one to correlate images over a very long time span which, in turn, translates into a large spatial baseline - resulting in a proportionately higher depth accuracy.

Barniv, Yair

Expansion-based passive ranging

This paper describes a new technique of passive ranging which is based on utilizing the image-plane expansion experienced by every object as its distance from the sensor decreases. This technique belongs in the feature/object-based family. The motion and shape of a small window, assumed to be fully contained inside the boundaries of some object, is approximated by an affine transformation. The parameters of the transformation matrix are derived by initially comparing successive images, and progressively increasing the image time separation so as to achieve much larger triangulation baseline than currently possible. Depth is directly derived from the expansion part of the transformation. To a first approximation, image-plane expansion is independent of image-plane location with respect to the focus of expansion (FOE) and of platform maneuvers. Thus, an expansion-based method has the potential of providing a reliable range in the difficult image area around the FOE. In areas far from the FOE the shift parameters of the affine transformation can provide more accurate depth information than the expansion alone, and can thus be used similarly to the way they have been used in conjunction with the Inertial Navigation Unit (INU) and Kalman filtering. However, the performance of a shift-based algorithm, when the shifts are derived from the affine transformation, would be much improved compared to current algorithms because the shifts--as well as the other parameters--can be obtained between widely separated images. Thus, the main advantage of this new approach is that, allowing the tracked window to expand and rotate, in addition to moving laterally, enables one to correlate images over a very long time span which, in turn, translates into a large spatial baseline resulting in a proportionately higher depth accuracy.

Barniv, Yair

Chlorophyll Fluorescence Emissions of Vegetation Canopies From High Resolution Field Reflectance Spectra

A two-year experiment was performed on corn (Zea mays L.) crops under nitrogen (N) fertilization regimes to examine the use of hyperspectral canopy reflectance information for estimating chlorophyll fluorescence (ChlF) and vegetation production. Fluorescence of foliage in the laboratory has proven more rigorous than reflectance for correlation to plant physiology. Especially useful are emissions produced from two stable red and far-red chlorophyll ChlF peaks centered at 685V10 nm and 735V5 nm. Methods have been developed elsewhere to extract steady state solar induced fluorescence (SF) from apparent reflectance of vegetation canopies/landscapes using the Fraunhofer Line Depth (FLD) principal. Our study utilized these methods in conjunction with field-acquired high spectral resolution canopy reflectance spectra obtained in 2004 and 2005 over corn crops, as part of an ongoing multi-year experiment at the USDA/Agriculture Research Service in Beltsville, MD. A spectroradiometer (ASD-FR Fieldspec Pro, Analytical Spectral Devices, Inc., Boulder, CO) was used to measure canopy radiances 1 m above plant canopies with a 22deg field of view and a 0deg nadir view zenith angle. Canopy and plant measurements were made at the R3 grain fill reproductive stage on 3-4 replicate N application plots provided seasonal inputs of 280, 140, 70, and 28 kg N/ha. Leaf level measurements were also made which included ChlF, photosynthesis, and leaf constituents (photosynthetic pigment, carbon (C), and N contents). Crop yields were determined at harvest. SIF intensities for ChlF were derived directly from canopy reflectance spectra in specific narrowband regions associated with atmospheric oxygen absorption features centered at 688 and 760 nm. The red/far-red S F ratio derived from these field reflectance spectra successfully discriminated foliar pigment levels (e.g., total chlorophyll, Chl) associated with N application rates in both corn crops. This canopy-level spectral ratio was also positively correlated to the foliar C/N ratio (r = 0.89, n = go), as was a leaf-level steady state fluorescence ratio (Fs/Chl, r = 0.92). The latter ratio was inversely correlated with crop grain yield (Kg 1 ha) (r = 0.9). This study has relevance to future passive satellite remote sensing approaches to monitoring C dynamics from space.

Middleton, E. M.

Reconstructing Quasar Spectra and Measuring the Lyα Forest with SpenderQ

Quasar spectra carry the imprint of foreground intergalactic medium (IGM) through absorption features. In particular, absorption caused by neutral hydrogen gas, the "Lyα forest," is a key spectroscopic tracer for cosmological analyses used to measure cosmic expansion and test physics beyond the standard model. Despite their importance, current methods for measuring Lyα absorption cannot directly derive the intrinsic quasar continuum and make strong assumptions on its shape, thus distorting the measured Lyα clustering. We present SpenderQ , a ML-based approach for directly reconstructing the intrinsic quasar spectra and measuring the Lyα forest from observations. SpenderQ uses the Spender spectrum autoencoder to learn a compact and redshift-invariant latent encoding of quasar spectra, combined with an iterative procedure to identify and mask absorption regions. To demonstrate its performance, we apply SpenderQ to 400,000 synthetic quasar spectra created to validate the Dark Energy Spectroscopic Instrument Year 1 Lyα cosmological analyses. SpenderQ accurately reconstructs the true intrinsic quasar spectra, including the broad Lyβ, Lyα, SiIV, CIV, and CIII emission lines. Redward of Lyα, SpenderQ provides percent-level reconstructions of the true quasar spectra. Blueward of Lyα, SpenderQ reconstructs the true spectra to < 5%. SpenderQ reproduces the shapes of individual quasar spectra more robustly than the current state-of-the-art. We, thus, expect it will significantly reduce biases in Lyα clustering measurements and enable studies of quasars and their physical properties. SpenderQ also provides informative latent variable encodings that can be used to, e.g., classify quasars with Broad Absorption Lines. Overall, SpenderQ provides a new data-driven approach for unbiased Lyα forest measurements in cosmological, quasar, and IGM studies.

Hahn, ChangHoon [Arizona U., Astron. Dept. - Stewa

Journal bearing impedance descriptions for rotordynamic applications

The paper deals with the development of analytic descriptions for plain circumferentially-symmetric fluid journal bearings, which are suitable for use in rotor dynamic analysis. The bearing impedance vector is introduced, which defines the bearing reaction force components as a function of the bearing motion. Impedances are derived directly for the Ocvirk (short) and Sommerfeld (long) bearings, and the relationships between the impedance vector and the more familiar mobility vector are developed and used to derive analytic impedance for finite-length bearings. The static correctness of the finite-length cavitating impedance is verified. Analytic stiffness and damping coefficient definitions are derived in terms of an impedance vector for small motion around an equilibrium position and demonstrated for the finite-length cavitating impedance. Nonlinear transient rotordynamic simulations are presented for the short pi and 2-pi impedances and the finite-length cavitating impedance. It is shown that finite-length impedance yields more accurate results for substantially less computer time than the short-bearing numerical-pressure-integration approach.

Childs, D.

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Several Different Wings Using a Discrete Vortex Method

Calculated numerical values for some aerodynamic terms and stability Derivatives for several different wings in unseparated inviscid incompressible flow were made using a discrete vortex method involving a limited number of horseshoe vortices. Both longitudinal and lateral-directional derivatives were calculated for steady conditions as well as for sinusoidal oscillatory motions. Variables included the number of vortices used and the rotation axis/moment center chordwise location. Frequencies considered were limited to the range of interest to vehicle dynamic stability (kb <.24 ). Comparisons of some calculated numerical results with experimental wind-tunnel measurements were in reasonable agreement in the low angle-of-attack range considering the differences existing between the mathematical representation and experimental wind-tunnel models tested. Of particular interest was the presence of induced drag for the oscillatory condition.

Riley, Donald R.

Calculation of the flow field including boundary layer effects for supersonic mixed compression inlets at angles of attack

The flow field in supersonic mixed compression aircraft inlets at angle of attack is calculated. A zonal modeling technique is employed to obtain the solution which divides the flow field into different computational regions. The computational regions consist of a supersonic core flow, boundary layer flows adjacent to both the forebody/centerbody and cowl contours, and flow in the shock wave boundary layer interaction regions. The zonal modeling analysis is described and some computational results are presented. The governing equations for the supersonic core flow form a hyperbolic system of partial differential equations. The equations for the characteristic surfaces and the compatibility equations applicable along these surfaces are derived. The characteristic surfaces are the stream surfaces, which are surfaces composed of streamlines, and the wave surfaces, which are surfaces tangent to a Mach conoid. The compatibility equations are expressed as directional derivatives along streamlines and bicharacteristics, which are the lines of tangency between a wave surface and a Mach conoid.

Vadyak, J.

Io - Geochemistry of sulfur

The evidence from Voyager imaging, Earth-based spectral reflectivity studies, and thermal emission measurements combine to suggest an extremely fresh, volcanically recycled sulfur-rich crust for Io, with very shallow large-scale melting. Two present styles of volcanism are possible, depending on the thickness of local deposits of sulfur: shallow liquid sulfur magma generation with quiescent flooding, and high-temperature volcanism with violet eruption of a sulfur-iron magma driven by SO2. Evolutionary considerations preclude direct derivation of Io's lithosphere from any metal-bearing chondritic source material. Metal-free C3V- or C2M-type parent material of either primary or secondary origin is the most plausible direct antecedent of the present sulfur-rich crust. Sulfates are almost certainly important constituents of the mantle, and can participate in the recycling of reduced, dense sulfide species to prevent total extraction of sulfur into the core.

Lewis, J. S.

Effect of front and rear incident proton irradiation on silicon solar cells

Four solar cell types of current manufacture were irradiated through the front and rear surfaces with protons in the energy range between 1 and 10 MeV. The solar cell parameters varied for this study were cell thickness and back surface field (BSF) vs. no BSF. Some cells were irradiated at normal incidence and an equal number were irradiated with simulated isotropic fluences. The solar cell electrical characteristics were measured under simulated AM0 illumination after each fluence. Using the normal incidence data, proton damage coefficients were computed for all four types of cells for both normal and omnidirectional radiation fields. These were found to compare well with the omnidirectional damage coefficients derived directly from the rear-incidence radiation data. Similarly, the rear-incidence omnidirectional radiation data were used to compute appropriate damage coefficients. A method for calculating the effect of a spectrum of energies is derived from these calculations. It is suitable for calculating the degradation of cells in space when they have minimal rear-surface shielding.

Anspaugh, Bruce

Finite difference operators on unstructured triangular meshes

A new form of the Laplace operator is derived which is shown to be second-order accurate when calculated at mesh triangle nodes if the cells corresponding to the nodes have a high degree of symmetry. It is demonstrated that the Laplace operator directly derived from Stokes' integral theorem is only zeroth-order accurate on these same cells. The construction of gradient and Laplace operators on unstructured triangular grids are discussed from the variational point of view. Two discrete representations of the Laplacian are compared to each other, and conclusions are drawn regarding their accuracy on a pointwise basis. Finally, geometric relations valid on arbitrary polygons are given for completeness in an appendix.

Erlebacher, G.

Shape sensitivity analysis of flutter response of a laminated wing

A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates' modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities, but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other. The sensitivities of the eigenvalues were used to predict the flutter speed, frequency, and reduced frequency. These approximations were found to be in good agreement with those obtained using a complete reanalysis.

Bergen, Fred D.

Shape sensitivity analysis of flutter response of a laminated wing

A method is presented for calculating the shape sensitivity of a wing aeroelastic response with respect to changes in geometric shape. Yates' modified strip method is used in conjunction with Giles' equivalent plate analysis to predict the flutter speed, frequency, and reduced frequency of the wing. Three methods are used to calculate the sensitivity of the eigenvalue. The first method is purely a finite difference calculation of the eigenvalue derivative directly from the solution of the flutter problem corresponding to the two different values of the shape parameters. The second method uses an analytic expression for the eigenvalue sensitivities of a general complex matrix, where the derivatives of the aerodynamic, mass, and stiffness matrices are computed using a finite difference approximation. The third method also uses an analytic expression for the eigenvalue sensitivities, but the aerodynamic matrix is computed analytically. All three methods are found to be in good agreement with each other.

Kapania, Rakesh K.