Engineering Papers⌕ Search

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.

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

At least 1,117 records · Page 62

An experimental study of the turbulent boundary layer on a transport wing in subsonic and transonic flow

The upper surface boundary layer on a transport wing model was extensively surveyed with miniature yaw probes at a subsonic and a transonic cruise condition. Additional data were obtained at a second transonic test condition, for which a separated region was present at mid-semispan, aft of mid-chord. Significant variation in flow direction with distance from the surface was observed near the trailing edge except at the wing root and tip. The data collected at the transonic cruise condition show boundary layer growth associated with shock wave/boundary layer interaction, followed by recovery of the boundary layer downstream of the shock. Measurements of fluctuating surface pressure and wingtip acceleration were also obtained. The influence of flow field unsteadiness on the boundary layer data is discussed. Comparisons among the data and predictions from a variety of computational methods are presented. The computed predictions are in reasonable agreement with the experimental data in the outboard regions where 3-D effects are moderate and adverse pressure gradients are mild. In the more highly loaded mid-span region near the trailing edge, displacement thickness growth was significantly underpredicted, except when unrealistically severe adverse pressure gradients associated with inviscid calculations were used to perform boundary layer calculations.

Spaid, Frank W.↗

AXAF SIM focus mechanism study

The design requirements and initial design concept for the AXAF-I Science Instrument Module (SIM) were reviewed at Ball on September 29, 1993. The concept design SIM focus mechanism utilizes a planetary gearset, with redundant motors, to drive a large ring (called 'main housing bearing') via a spur gearset. This large drive ring actuates three tangent bar links (called 'push rods'), which in turn actuate three levers (called 'pin levers'). Each of the three pin levers rotates an 'eccentric pin,' which in turn moves the base of a bipod flexure in both the radial (normal to optical axis) and axial (focus along optical axis) directions. Three bipod flexures are employed, equally spaced at 120 degrees apart, the base of each being translated in the two directions as described above. A focus adjustment is made by rotating the drive ring, which drives the push rods and therefore the pin levers, which in turn rotate the eccentric pins, finally imparting the two motions to the base of each of the bipod flexures. The axial translation (focus adjustment) of the focused structure is the sum of the direct axial motion plus axial motion which comes from uniformly squeezing the three bipod bases radially inward. SAO documented the following concerns regarding the focus mechanism in memo WAP-FY94-001, dated October 7, 1993: (1) The focus adjustment depends, in large part, on the structural properties (stiffnesses and end fixities) of the bipod flexures, push rods, pin levers and eccentric pins. If these properties are not matched very well, then lateral translations as well as unwanted rotations of the focussed structure will accompany focus motion. In addition, the stackup of linkage tolerances and any nonuniform wear in the linkages will result in the same unwanted motions. Thermal gradients will also affect these motions. At the review Ball did not present supporting analyses to support their choice of this design concept. (2) The proposed 'primary' method of measuring focus is by counting motor steps. The 'backup' method is by a pot mounted on the drive ring. Neither method provides for a direct measurement of the quantity desired (focus position). This is of concern because of the long and indirect relationship between focus and the sensed quantity (drive ring rotation). There are three sinusoidal relationships and structural stiffness in the path, and the resulting calibration is likely to be highly nonlinear. These methods would require an accurate ground calibration. (3) Ground calibration (and verification) of focus vs. drive position must be done in 1-g on the ground. This calibration will be complicated by both the structural characteristics of the bipods and the fact that the CG of the translating portion of the SIM is not on the optical axis (thereby causing unwated rotations and changing the focus position vs. motor step and pot readout relationships). The SIM translating weight could be offloaded, but the calibration then becomes sensitive to any errors in offloading (both magnitude and direction). There are concerns as to whether a calibration to the required accuracy can be accomplished on the ground. (4) The choice of a potentiometer as the focus position sensor is questionable in terms of reliability for a five year mission. The results of SAO's study of items 1, 2 and 3 described above are presented in this report.

Tananbaum, H. D.↗

Effects of solar radiation pressure torque on the rotational motion of an artificial satellite

The motion of an artificial satellite about its center of mass is studied considering torques due to the gravity gradient and direct solar radiation pressure. A model for direct solar radiation torque is derived for a circular cylindrical satellite. An analytical solution is obtained by the method of variation of the parameters. This solution shows that the angular variables have secular variation but that the modulus of the rotational angular momentum, the projection of rotational angular momentum on the z axis of the moment of inertia and inertial axis z, suffer only periodic variations. Considering a hypothetical artificial satellite, a numerical application is demonstrated.

Zanardi, Maria Cecilia F. P. S.↗

Automatic extraction of planetary image features

A method for the extraction of Lunar data and/or planetary features is provided. The feature extraction method can include one or more image processing techniques, including, but not limited to, a watershed segmentation and/or the generalized Hough Transform. According to some embodiments, the feature extraction method can include extracting features, such as, small rocks. According to some embodiments, small rocks can be extracted by applying a watershed segmentation algorithm to the Canny gradient. According to some embodiments, applying a watershed segmentation algorithm to the Canny gradient can allow regions that appear as close contours in the gradient to be segmented.

LeMoigne-Stewart, Jacqueline J.↗

Verifying Implementation of the Dryden Turbulence Model and MIL-F-8785 Gust Gradient

Turbulence modeling in human-in-the-loop simulation is important to assessing aircraft handling qualities and pilot performance and to provide additional realism for pilot training. In the simulation community, the Dryden turbulence spectra is a popular choice for modeling the linear turbulent gusts because its rational form is efficiently reproduced by passing white noise through linear filters. The MIL-F-8785 gust gradients similarly use additional linear filters to model the gradient of the turbulent gust over the wing, and it represents the gust gradients as perturbations to the air-relative rotational rates. The Cockpit Motion Facility at NASA Langley Research Center (LaRC) models continuous random turbulence using the Dryden one-dimensional spectra and MIL-F-8785 gust gradient. The facility recently reviewed and updated its verification of these models as part of an initiative to improve motion cueing under turbulence. This exercise introduced improved methods for verifying the turbulence models and led to rediscovery of model assumptions that informed improvements to implementation.

Madden, Michael M.↗

Advancement of the General Aviation Synthesis Program Using Python to Enable Optimization-Based Hybrid-Propulsion Aircraft Design

In support of the Electrified Powertrain Flight Demonstrator and Advanced Air Transport Technologies projects at NASA, a new tool has been developed at NASA’s Ames and Glenn Research Centers to enable coupled engine and airframe optimization and analysis. The new tool combines the engineering-level analysis methods and empirical models of the FORTRAN General Aviation Synthesis Program (GASP) with the Python-based OpenMDAO framework to provide a modular framework for efficient gradient-based optimization with the aim of incorporating new subsystem models for unconventional configurations. The tool has been verified against GASP analyses of several aircraft models and mission formulations. Preliminary efforts have been made to integrate pyCycle, a thermodynamic cycle analysis tool, to enable simultaneous optimization of hybrid propulsion system and vehicle parameters while taking full mission performance and constraints into account. This will improve current capabilities to assess impacts of electrified powertrain technologies on future aircraft designs.

Kenneth R. Lyons↗

Advancement of the General Aviation Synthesis Program Using Python to Enable Optimization-Based Hybrid-Propulsion Aircraft Design

In support of the Electrified Powertrain Flight Demonstrator and Advanced Air Transport Technologies projects at NASA, a new tool has been developed at NASA's Ames and Glenn Research Centers to enable coupled engine and airframe optimization and analysis. The new tool combines the engineering-level analysis methods and empirical models of the FORTRAN General Aviation Synthesis Program (GASP) with the Python-based OpenMDAO framework to provide a modular framework for efficient gradient-based optimization with the aim of incorporating new subsystem models for unconventional configurations. The tool has been verified against GASP analyses of several aircraft models and mission formulations. Preliminary efforts have been made to integrate pyCycle, a thermodynamic cycle analysis tool, to enable simultaneous optimization of hybrid propulsion system and vehicle parameters while taking full mission performance and constraints into account. This will improve current capabilities to assess impacts of electrified powertrain technologies on future aircraft designs.

Kenneth R Lyons↗

Validation of Shadowgraph Spectral Analysis using an SLS Block 2 Wind-Tunnel Model

A flow quantification method has previously been developed that uses high-speed shadowgraph images to extract frequency content of a transonic flow field around a wind-tunnel model. In this method, a classical spectral analysis is performed on a set of shadowgraph image pixels to identify their intensity fluctuation frequencies. The spatial change of the image intensity is related to the second derivative of the density gradient. Examining the change of the density gradient as a function of time can be used to understand the flow field pressure fluctuations. Pressure measurements collected during a Booster Obsolescence Life Extension Forward Attachment Structure Trade Study Test are used to calibrate the shadowgraph-based measurements. The test was performed at the NASA Ames Research Center Unitary Plan Wind Tunnel 11- by 11-foot test section and investigated aerodynamic environment downstream of various booster forward attachment hardware geometries of the Space Launch System Block 2 Cargo configuration. Time-correlated shadowgraph video and model surface pressure measurements were acquired during the test and served as an ideal opportunity to calibrate and validate the shadowgraph image intensity method.

Shadowgraph↗

A space-marching method for the computation of viscous internal flows

A space-marching method has been developed to compute 3-D viscous flows in internal geometries. The Navier-Stokes equations have been posed as an initial-value problem by neglecting the effects of streamwise diffusion and treating the streamwise pressure gradient as a known source term. The fully coupled system of equations has been solved by a noniterative algorithm at each streamwise step of the computation. A low Mach number formulation of the equations has been used to compute incompressible flow fields. A computer program has been written to implement all aspects of the space-marching algorithm. The program is modular and is easily adapted to the widely varying geometries of internal flows. The space-marching algorithm has been tested by computing simple flows with known analytical solutions. The method has been used to predict complex 3-D turbulent flows. The algorithm is stable and very economical. A single sweep of the flow field by the space-marching method is approximately equivalent to one time-step of the time-marching method.

Govindan, T. R.↗

High-Order Coronagraphic Wavefront Control With Algorithmic Differentiation: First Experimental Demonstration

Future space-based coronagraphs will rely critically on focal-plane wavefront sensing and control with deformable mirrors to reach deep contrast by mitigating optical aberrations in the primary beam path. Until now, most focal-plane wavefront control algorithms have been formulated in terms of Jacobian matrices, which encode the predicted effect of each deformable mirror actuator on the focal-plane electric field. A disadvantage of these methods is that Jacobian matrices can be cumbersome to compute and manipulate, particularly when the number of deformable mirror actuators is large. Recently, we proposed a new class of focal-plane wavefront control algorithms that utilize gradient-based optimization with algorithmic differentiation to compute wavefront control solutions while avoiding the explicit computation and manipulation of Jacobian matrices entirely. In simulations using a coronagraph design for the proposed Large UV/Optical/Infrared Surveyor (LUVOIR), we showed that our approach reduces overall CPU time and memory consumption compared to a Jacobian-based algorithm. Here, we expand on these results by implementing the proposed algorithm on the High Contrast Imager for Complex Aperture Telescopes (HiCAT) testbed at the Space Telescope Science Institute (STScI) and present initial experimental results, demonstrating contrast suppression capabilities equivalent to Jacobian-based methods.

wavefront control↗

Time dependent inflow-outflow boundary conditions for 2D acoustic systems

An analysis of the number and form of the required inflow-outflow boundary conditions for the full two-dimensional time-dependent nonlinear acoustic system in subsonic mean flow is performed. The explicit predictor-corrector method of MacCormack (1969) is used. The methodology is tested on both uniform and sheared mean flows with plane and nonplanar sources. Results show that the acoustic system requires three physical boundary conditions on the inflow and one on the outflow boundary. The most natural choice for the inflow boundary conditions is judged to be a specification of the vorticity, the normal acoustic impedance, and a pressure gradient-density gradient relationship normal to the boundary. Specification of the acoustic pressure at the outflow boundary along with these inflow boundary conditions is found to give consistent reliable results. A set of boundary conditions developed earlier, which were intended to be nonreflecting is tested using the current method and is shown to yield unstable results for nonplanar acoustic waves.

Watson, Willie R.↗

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

Global-Scale Consequences of Magnetic-Helicity Injection and Condensation on the Sun

In the recent paper of Antiochos, a new concept for the injection of magnetic helicity into the solar corona by small-scale convective motions and its condensation onto polarity inversion lines (PILs) has been developed. We investigate this concept through global simulations of the Sun's photospheric and coronal magnetic fields and compare the results with the hemispheric pattern of solar filaments. Assuming that the vorticity of the cells is predominately counter-clockwise/clockwise in the northern/southern hemisphere, the convective motions inject negative/positive helicity into each hemisphere. The simulations show that: (i) On a north-south orientated PIL, both differential rotation and convective motions inject the same sign of helicity which matches that required to reproduce the hemispheric pattern of filaments. (ii) On a high latitude east-west orientated polar crown or sub-polar crown PIL, the vorticity of the cells has to be approximately 2-3 times greater than the local differential rotation gradient in order to overcome the incorrect sign of helicity injection from differential rotation. (iii) In the declining phase of the cycle, as a bipole interacts with the polar field, in some cases helicity condensation can reverse the effect of differential rotation along the East-West lead arm, but not in all cases. The results show that this newly developed concept of magnetic helicity injection and condensation is a viable method to explain the hemispheric pattern of filaments in conjunction with the mechanisms used in Yeates et al. (2008). Future observational studies should focus on determining the vorticity component within convective motions to determine, both its magnitude and latitudinal variation relative to the differential rotation gradient on the Sun.

magnetic fields↗

System Sensitivity Analysis Applied to the Conceptual Design of a Dual-Fuel Rocket SSTO

This paper reports the results of initial efforts to apply the System Sensitivity Analysis (SSA) optimization method to the conceptual design of a single-stage-to-orbit (SSTO) launch vehicle. SSA is an efficient, calculus-based MDO technique for generating sensitivity derivatives in a highly multidisciplinary design environment. The method has been successfully applied to conceptual aircraft design and has been proven to have advantages over traditional direct optimization methods. The method is applied to the optimization of an advanced, piloted SSTO design similar to vehicles currently being analyzed by NASA as possible replacements for the Space Shuttle. Powered by a derivative of the Russian RD-701 rocket engine, the vehicle employs a combination of hydrocarbon, hydrogen, and oxygen propellants. Three primary disciplines are included in the design - propulsion, performance, and weights & sizing. A complete, converged vehicle analysis depends on the use of three standalone conceptual analysis computer codes. Efforts to minimize vehicle dry (empty) weight are reported in this paper. The problem consists of six system-level design variables and one system-level constraint. Using SSA in a 'manual' fashion to generate gradient information, six system-level iterations were performed from each of two different starting points. The results showed a good pattern of convergence for both starting points. A discussion of the advantages and disadvantages of the method, possible areas of improvement, and future work is included.

Olds, John R.↗

Radio jet refraction in galactic atmospheres with static pressure gradients

A theory of double radio sources which have a 'Z' or 'S' morphology is proposed, based on the refraction of radio jets in the extended atmosphere of an elliptical galaxy. The model describes a collimated jet of supersonic material bending self-consistently under the influence of external static pressure gradients. Gravity and magnetic fields are neglected in the simplest case except insofar as they determine the static pressure distribution. The calculation is a straightforward extension of a method used to calculate a ram-pressure model for twin radio trails ('C' morphology). It may also be described as a continuous-jet version of a buoyancy model proposed in 1973. The model has the added virtue of invoking a galactic atmosphere similar to those already indicated by X-ray measurements of some other radio galaxies and by models for the collimation of other radio jets.

Henriksen, R. N.↗

Propellant transfer: Tethered depot

Spacebasing of orbital transfer vehicles at a space station requires a depot that safely and efficiently stores and transfers the resupply propellants. In order to transfer propellants, a method effectively acquire only liquid and vent only gas must exist. A method that produces a low gravity to settle propellants would bypass these weaknesses, while allowing ground-like operations. This low gravity can be passively produced using gravity gradient techniques. A satellite with a large length to diameter ratio, such as a depot attached to a space station with a tether, stabilizes along on Earth radial because of an outward acceleration proportional to the distance from the satellite's center of gravity. Analysis indicates that liquid can be setteled with relatively short tether lengths. The feasibility, design requirements, and operational limitations of a tethered refueling depot were investigated with special emphasis on slosh control.

Kroll, K.↗

Investigations on entropy layer along hypersonic hyperboloids using a defect boundary layer

A defect approach coupled with matched asymptotic expansions is used to derive a new set of boundary layer equations. This method ensures a smooth matching of the boundary layer with the inviscid solution. These equations are solved to calculate boundary layers over hypersonic blunt bodies involving the entropy gradient effect. Systematic comparisons are made for both axisymmetric and plane flows in several cases with different Mach and Reynolds numbers. After a brief survey of the entropy layer characteristics, the defect boundary layer results are compared with standard boundary layer and full Navier-Stokes solutions. The entropy gradient effects are found to be more important in the axisymmetric case than in the plane one. The wall temperature has a great influence on the results through the displacement effect. Good predictions can be obtained with the defect approach over a cold wall in the nose region, with a first order solution. However, the defect approach gives less accurate results far from the nose on axisymmetric bodies because of the thinning of the entropy layer.

Brazier, J. P.↗

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗