SEARCH · Engineering Papers

Results for “gradient methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Renormalization group methods for the Reynolds stress transport equations

The Yakhot-Orszag renormalization group is used to analyze the pressure gradient-velocity correlation and return to isotropy terms in the Reynolds stress transport equations. The perturbation series for the relevant correlations, evaluated to lowest order in the epsilon-expansion of the Yakhot-Orszag theory, are infinite series in tensor product powers of the mean velocity gradient and its transpose. Formal lowest order Pade approximations to the sums of these series produce a rapid pressure strain model of the form proposed by Launder, Reece, and Rodi, and a return to isotropy model of the form proposed by Rotta. In both cases, the model constants are computed theoretically. The predicted Reynolds stress ratios in simple shear flows are evaluated and compared with experimental data. The possibility is discussed of deriving higher order nonlinear models by approximating the sums more accurately. The Yakhot-Orszag renormalization group provides a systematic procedure for deriving turbulence models. Typical applications have included theoretical derivation of the universal constants of isotropic turbulence theory, such as the Kolmogorov constant, and derivation of two equation models, again with theoretically computed constants and low Reynolds number forms of the equations. Recent work has applied this formalism to Reynolds stress modeling, previously in the form of a nonlinear eddy viscosity representation of the Reynolds stresses, which can be used to model the simplest normal stress effects. The present work attempts to apply the Yakhot-Orszag formalism to Reynolds stress transport modeling.

Rubinstein, R.↗

NASA NTRS1994

Measuring attitude with a gradiometer

This paper explores using a gravity gradiometer to measure the attitude of a satellite, given that the gravity field is accurately known. Since gradiometers actually measure a combination of the gradient and attitude rate and acceleration terms, the answer is far from obvious. The paper demonstrates that it can be done and at microradian accuracy. The technique employed is dynamic estimation, based on the momentum biased Euler equations. The satellite is assumed nominally planet pointed, and subject to control, gravity gradient, and partly radom drag torques. The attitude estimator is unusual. While the standard method of feeding back measurement residuals is used, the feedback gain matrix isn't derived from Kalman theory. instead, it's chosen to minimize a measure of the terminal covariance of the error in the estimate. This depends on the gain matrix and the power spectra of all the process and measurement noises. An integration is required over multiple solutions of Lyapunov equations.

Sonnabend, David↗

NASA NTRS1994

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗

NASA NTRS1995

Flight test validation of a frequency-based system identification method on an F-15 aircraft

A frequency-based performance identification approach was evaluated using flight data from the NASA F-15 Highly Integrated Digital Electronic Control aircraft. The approach used frequency separation to identify the effectiveness of multiple controls simultaneously as an alternative to independent control identification methods. Fourier transformations converted measured control and response data into frequency domain representations. Performance gradients were formed using multiterm frequency matching of control and response frequency domain models. An objective function was generated using these performance gradients. This function was formally optimized to produce a coordinated control trim set. This algorithm was applied to longitudinal acceleration and evaluated using two control effectors: nozzle throat area and inlet first ramp. Three criteria were investigated to validate the approach: simultaneous gradient identification, gradient frequency dependency, and repeatability. This report describes the flight test results. These data demonstrate that the approach can accurately identify performance gradients during simultaneous control excitation independent of excitation frequency.

Schkolnik, Gerard S.↗

NASA NTRS2007

Space Technology 5 Multi-point Measurements of Near-Earth Magnetic Fields: Initial Results

The Space Technology 5 (ST-5) mission successfully placed three micro-satellites in a 300 x 4500 km dawn-dusk orbit on 22 March 2006. Each spacecraft carried a boom-mounted vector fluxgate magnetometer that returned highly sensitive and accurate measurements of the geomagnetic field. These data allow, for the first time, the separation of temporal and spatial variations in field-aligned current (FAC) perturbations measured in low-Earth orbit on time scales of approximately 10 sec to 10 min. The constellation measurements are used to directly determine field-aligned current sheet motion, thickness and current density. In doing so, we demonstrate two multi-point methods for the inference of FAC current density that have not previously been possible in low-Earth orbit; 1) the "standard method," based upon s/c velocity, but corrected for FAC current sheet motion, and 2) the "gradiometer method" which uses simultaneous magnetic field measurements at two points with known separation. Future studies will apply these methods to the entire ST-5 data set and expand to include geomagnetic field gradient analyses as well as field-aligned and ionospheric currents.

Slavin, James A.↗

NASA NTRS2012

Asymmetric Uncertainty Expression for High Gradient Aerodynamics

When the physics of the flow around an aircraft changes very abruptly either in time or space (e.g., flow separation/reattachment, boundary layer transition, unsteadiness, shocks, etc), the measurements that are performed in a simulated environment like a wind tunnel test or a computational simulation will most likely incorrectly predict the exact location of where (or when) the change in physics happens. There are many reasons for this, includ- ing the error introduced by simulating a real system at a smaller scale and at non-ideal conditions, or the error due to turbulence models in a computational simulation. The un- certainty analysis principles that have been developed and are being implemented today do not fully account for uncertainty in the knowledge of the location of abrupt physics changes or sharp gradients, leading to a potentially underestimated uncertainty in those areas. To address this problem, a new asymmetric aerodynamic uncertainty expression containing an extra term to account for a phase-uncertainty, the magnitude of which is emphasized in the high-gradient aerodynamic regions is proposed in this paper. Additionally, based on previous work, a method for dispersing aerodynamic data within asymmetric uncer- tainty bounds in a more realistic way has been developed for use within Monte Carlo-type analyses.

Pinier, Jeremy T↗

NASA NTRS2002

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

NASA NTRSDate not supplied

Advanced Metal Separation (AMS) via Marangoni Effect-driven Separation

For practical, economical, and sustainable human missions to the Moon, processing of lunar regolith and metal oxides is necessary for localized production of oxygen (O2) and metals that can be repurposed as rocket propellant, life support O2, and in-space manufactured parts. The primary candidates for in-situ resource utilization (ISRU) regolith processing are carbothermic reductions and molten regolith electrolysis (MRE) which both produce oxygen and a reduced, deoxygenated metallic slag. A metal separation technique is still required to capture and purify metals from the deoxygenated slag for reuse in areas such as metallic part fabrication, structural or design frameworks, and metal catalyst production. The Marangoni Effect is best explained as mass transfer along an interface driven by surface tension gradients induced by temperature and concentration differences. Due to the temperature gradients, molten metal can be transported in this fashion to create thin films. This method can transport molten metal in the absence of gravity or mechanical motion and greatly increases the available surface area for a vacuum fractional separation process. In this study, a novel vacuum furnace induces the Marangoni Effect to fractionally decompose and separate valuable metals (e.g., aluminum) from thin films of molten metal via vacuum distillation. The proposed work will examine lunar regolith, metal-based recycling feedstocks, and metallic slags produced from the carbothermal reduction of regolith. Data and observations from this study could produce a simple yet effective method of refining lunar derived metals by taking advantage of the high vacuum environment present on the lunar surface.

ISRU↗

NASA NTRSDate not supplied

Design Under Uncertainty for Conceptual Aircraft Design Leveraging Analytical Gradients

The purpose of this paper is to extend previously demonstrated methodologies for design under uncertainty, leveraging analytical gradients to higher fidelity analysis for use in conceptual aircraft design. Previous work developed methods to generate analytical derivatives through polynomial chaos expansion, eliminating the need to estimate derivatives via complex step or finite difference. In this research, the authors build upon the methods to include physics-based aircraft design codes for aircraft design under uncertainty. This extends the previous work’s case study, which employed analytical aerodynamics and Breguet range estimations for wing design, to a higher fidelity level. In addition, this work extends previous work on interface development between the Uncertainty Quantification with Polynomial Chaos Expansion (UQPCE) software and Model-Based Systems Analysis and Engineering (MBSA&E) frameworks. This paper will discuss the development work necessary to perform multidisciplinary design under uncertainty as well as demonstrate the mechanics of interfacing UQPCE and conceptual aircraft design tools such as NASA’s Aviary code. In a case study, a conceptual aircraft design under uncertainty was conducted and compared against a traditional deterministic design. When given information about the uncertainty space from UQPCE, the optimizer was able to shape the output distribution and produce a more robust design

UQ↗

NASA NTRSDate not supplied

Design Under Uncertainty for Conceptual Aircraft Design Leveraging Analytical Gradients

UQ↗

NASA NTRS1998

Preston Probe Calibrations at High Reynolds Number

The overall goal of the research effort is to study the performance of two Preston probes designed by NASA Langley Research Center across an unprecedented range of Reynolds number (based on friction velocity and probe diameter), and perform an accurate calibration over the same Reynolds number range. Using the Superpipe facility in Princeton, two rounds of experiments were performed. In each round of experiments for each Reynolds number, the pressure gradient, static pressure from the Preston probes and the total pressure from the Preston probes were measured. In the first round, 3 Preston probes having outer diameters of 0.058 inches, 0.083 inches and 0.203 inches were tested over a large range of pipe Reynolds numbers. Two data reduction methods were employed: first, the static pressure measured on the Preston probe was used to calculate P (modified Preston probe configuration), and secondly, the static pressure measured at the reference pressure tap was used to calculate P (un-modified Preston probe configuration). For both methods, the static pressure was adjusted to correspond with the static pressure at the Preston probe tip using the pressure gradient. The measurements for Preston probes with diameters of 0.058 inches, and 0.083 inches respectively were performed in the test pipe before it was polished a second time. Therefore, the measurements at high pipe Reynolds numbers may have been affected by roughness. In the second round of experiments the 0.058 inches and 0.083 inches diameter, un-modified probes were tested after the pipe was polished and prepared to ensure that the surface was smooth. The average velocity was estimated by assuming that the connection between the centerline velocity and the average velocity was known, and by using a Pitot tube to measure the centerline velocity. A preliminary error estimate suggests that it is possible to introduce a 1% to 2% error in estimating the average velocity using this approach. The evidence on the errors attending the second data set is somewhat circumstantial, and the measurements have not been repeated using a better approach, it seems probable that the correlation given applies to un-modified Preston probes over the range 6.4 less than x* less than 11.3.

Smits, Alexander J.↗

NASA NTRS2011

Investigation of the Crust of the Pannonian Basin, Hungary Using Low-Altitude CHAMP Horizontal Gradient Magnetic Anomalies

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. It is some 600 by 500 km in area and centered on Hungary. This area was chosen since it has one of the thinnest continental crusts in Europe and is the region of complex tectonic structures. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The SWARM constellation, scheduled to be launched next year, will have two lower altitude satellites flying abreast, with a separation of between ca. 150 to 200 km. to record the horizontal magnetic gradient. Since the CHAMP satellite has been in orbit for eight years and has obtained an extensive range of data, both vertically and horizontally there is a large enough data base to compute the horizontal magnetic gradients over the Pannonian Basin region using these many CHAMP orbits. We recomputed a satellite magnetic anomaly map, using the spherical-cap method of Haines (1985), the technique of Alsdorf et al. (1994) and from spherical harmonic coefficients of MF6 (Maus et aI., 2008) employing the latest and lowest altitude CHAMP data. We then computed the horizontal magnetic anomaly gradients (Kis and Puszta, 2006) in order to determine how these component data will improve our interpretation and to preview what the SW ARM mission will reveal with reference to the horizontal gradient anomalies. The gradient amplitude of an 1000 km northeast-southwest profile through our horizontal component anomaly map varied from 0 to 0.025 nT/km with twin positive anomalies (0.025 and 0.023 nT/km) separated by a sharp anomaly negative at o nT/km. Horizontal gradient indicate major magnetization boundaries in the crust (Dole and Jordan, 1978 and Cordell and Grauch, 1985). Our gradient anomaly was modeled with a twodimensional body and the anomaly, of some 200 km, correlates with a 200 km area of crustal thinning in the southwestern Pannonian Basin.

Taylor, Patrick T.↗

NASA NTRS2003

Measuring Diffusion of Liquids by Common-Path Interferometry

A method of observing the interdiffusion of a pair of miscible liquids is based on the use of a common-path interferometer (CPI) to measure the spatially varying gradient of the index refraction in the interfacial region in which the interdiffusion takes place. Assuming that the indices of refraction of the two liquids are different and that the gradient of the index of refraction of the liquid is proportional to the gradient in the relative concentrations of either liquid, the diffusivity of the pair of liquids can be calculated from the temporal variation of the spatial variation of the index of refraction. This method yields robust measurements and does not require precise knowledge of the indices of refraction of the pure liquids. Moreover, the CPI instrumentation is compact and is optomechanically robust by virtue of its common- path design. The two liquids are placed in a transparent rectangular parallelepiped test cell. Initially, the interface between the liquids is a horizontal plane, above which lies pure liquid 2 (the less-dense liquid) and below which lies pure liquid 1 (the denser liquid). The subsequent interdiffusion of the liquids gives rise to a gradient of concentration and a corresponding gradient of the index of refraction in a mixing layer. For the purpose of observing the interdiffusion, the test cell is placed in the test section of the CPI, in which a collimated, polarized beam of light from a low-power laser is projected horizontally through a region that contains the mixing layer.

Rashidnia, Nasser↗

NASA NTRS1995

A General Method for Solving Systems of Non-Linear Equations

The method of steepest descent is modified so that accelerated convergence is achieved near a root. It is assumed that the function of interest can be approximated near a root by a quadratic form. An eigenvector of the quadratic form is found by evaluating the function and its gradient at an arbitrary point and another suitably selected point. The terminal point of the eigenvector is chosen to lie on the line segment joining the two points. The terminal point found lies on an axis of the quadratic form. The selection of a suitable step size at this point leads directly to the root in the direction of steepest descent in a single step. Newton's root finding method not infrequently diverges if the starting point is far from the root. However, the current method in these regions merely reverts to the method of steepest descent with an adaptive step size. The current method's performance should match that of the Levenberg-Marquardt root finding method since they both share the ability to converge from a starting point far from the root and both exhibit quadratic convergence near a root. The Levenberg-Marquardt method requires storage for coefficients of linear equations. The current method which does not require the solution of linear equations requires more time for additional function and gradient evaluations. The classic trade off of time for space separates the two methods.

Nachtsheim, Philip R.↗

NASA NTRS1985

A Method for Computing the Core Flow in Three-Dimensional Leading-Edge Vortices

A theory is presented for calculating the flow in the core of a separation-induced leading-edge vortex. The method is based on matching inner and outer representations of the vortex. The inner model of the vortex is based on the quasicylindrical Navier-Stokes equations; the flow is assumed to be steady, axially symmetric, and incompressible and in addition, gradients in the radial direction are assumed to be much larger then gradients in the axial direction. The outer model is based on the three-dimensional free-vortex-sheet theory, a higher-order panel method which solves the Prandtl-Glauert equation including nonlinear boundary conditions pertinent to the concentrated vorticity representation of the leading edge vortex. The resultant flow is evaluated a posteriori for evidence of incipient vortex breakdown and the critical helix angle concept, in conjunction with an adverse longitudinal pressure gradient, is found to correlate well with the occurrence of vortex breakdown at the trailing edge of delta, arrow, and diamond wings.

Luckring, J. M.↗

NASA NTRS2007

Some Scenarios for Transition on Turbomachinery Blading

Measurements on transition under different levels of adverse pressure gradient and free-stream turbulence level are described. This extensive series of investigations, which was predicated on intermittency measurement techniques, has resulted in correlations for transition length and turbulent spot formation rate. These correlations rae intended to be used in conjunction with boundary layer prediction methods and examples are given of such predictions. More effective predictions of the transition region, especially under conditions of variable pressure gradient, are dependent on a more comprehensive understanding of transition and spot behavior. It is expected that this will result in improved transition modeling.

Gostelow, J. P.↗

NASA NTRS1991

Description and evaluation of an interference assessment for a slotted-wall wind tunnel

A wind-tunnel interference assessment method applicable to test sections with discrete finite-length wall slots is described. The method is based on high order panel method technology and uses mixed boundary conditions to satisfy both the tunnel geometry and wall pressure distributions measured in the slotted-wall region. Both the test model and its sting support system are represented by distributed singularities. The method yields interference corrections to the model test data as well as surveys through the interference field at arbitrary locations. These results include the equivalent of tunnel Mach calibration, longitudinal pressure gradient, tunnel flow angularity, wall interference, and an inviscid form of sting interference. Alternative results which omit the direct contribution of the sting are also produced. The method was applied to the National Transonic Facility at NASA Langley Research Center for both tunnel calibration tests and tests of two models of subsonic transport configurations.

Kemp, William B., Jr.↗

NASA NTRS2001

Method Improvements in Thermal Analysis of Mach 10 Leading Edges

Several improvements have recently been made in the thermal analysis methods for leading edges of a hypersonic vehicle. The leading edges of this vehicle undergo exceptionally high heat loads that incorporate extreme spatial gradients as well as severe transients. Due to the varying flight conditions, complex geometry, and need for thermal loads at many points along the trajectory, full computational fluid dynamics (CFD) analysis of the aeroheating; loads is not feasible. Thus, engineering methods must be used to determine the aeroheating on the vehicle surfaces, and that must be utilized in the thermal analysis. Over the last year, the thermal analysis of a hypersonic vehicle has been enhanced in several ways. Two different engineering codes are used to predict aeroheating loads: one over the curve near the stagnation point, and the other on flat surfaces downstream of the leading edge. These two are matched together at the intersection point using a method that allows closer approximation of CFD results. User-developed FORTRAN, which is part of the thermal solver PATRAN Thermal, is used to accomplish this. The customizable FORTRAN code also allows use of many different time- and space-dependent factors, interpolation of the heat load in time and space, and inclusion of both highly swept and unswept grid structures. This FORTRAN is available to other PATRAN users who may want to accomplish a similar objective in analysis. Flux, rather than convective coefficient, is used to define heat loads, which allows more accurate analysis as well as better application of margins. Improvements have also been made in more efficient utilization of imported CAD geometry, by creating faces on solids to facilitate load application.

Amundsen, Ruth M.↗