Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shape from gradient”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Effect of modulator sorption on gradient shape in ion-exchange chromatography

Mobile phase additives, or modulators, are used in gradient elution chromatography to facilitate separation and reduce separation time. The modulators are usually assumed to be linearly adsorbed or unadsorbed. Here, the consequences of nonlinear modulator adsorption are examined for ion-exchange gradient elution through a series of simulations. Even when the buffer salt is identical to the modulator salt, gradient deformation is observed; the extent of deformation increases as the volume of the feed is increased. When the modulator salt is different from the buffer salt, unusual effects are observed, and the chromatograms are quite different from those predicted by classical gradient elution theory. In particular, local increases in the buffer concentration are found between feed bands, and serve to improve the separation. These effects become more pronounced as the feed volume increases, and could therefore prove valuable in preparative applications.

NASA Discipline Life Support Systems↗
NASA NTRS1997

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗
NASA NTRS1997

Methods for Sensitivity Analysis in Gradient-Based Shape Optimization: A Review

CFD is now being widely used in the design of mechanical and aerospace systems and their components. Improving the accuracy and efficiency of this process can reduce its cost and increase its effectiveness. Other motivating factors that push and pull CFD into design optimization are; (1) the need for design methodologies (single- and multi-disciplinary) to reduce cost and time in developing new fluid and aero systems for fast response to market changes; (2) new design requirements outside existing database; (3) new CFD algorithms and growth in computing power; (4) need for increased product quality and reliability; (5) lack of test facilities; (6) limitations and cost of testing techniques; (7) limitations of existing design practices; and (8) prohibitive cost of and conflicting designs derived from isolated component analyses and design, which are repetitive in nature and omit mutual interference.

Baysal, Oktay↗
NASA NTRS1997

Gradient-Based Aerodynamic Shape Optimization Using ADI Method for Large-Scale Problems

A gradient-based shape optimization methodology, that is intended for practical three-dimensional aerodynamic applications, has been developed. It is based on the quasi-analytical sensitivities. The flow analysis is rendered by a fully implicit, finite volume formulation of the Euler equations.The aerodynamic sensitivity equation is solved using the alternating-direction-implicit (ADI) algorithm for memory efficiency. A flexible wing geometry model, that is based on surface parameterization and platform schedules, is utilized. The present methodology and its components have been tested via several comparisons. Initially, the flow analysis for for a wing is compared with those obtained using an unfactored, preconditioned conjugate gradient approach (PCG), and an extensively validated CFD code. Then, the sensitivities computed with the present method have been compared with those obtained using the finite-difference and the PCG approaches. Effects of grid refinement and convergence tolerance on the analysis and shape optimization have been explored. Finally the new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4. Despite the expected increase in the computational time, the results indicate that shape optimization, which require large numbers of grid points can be resolved with a gradient-based approach.

Pandya, Mohagna J.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Bune, Andris V.↗
NASA NTRS1997

Three-Dimensional Viscous Alternating Direction Implicit Algorithm and Strategies for Shape Optimization

A gradient-based shape optimization based on quasi-analytical sensitivities has been extended for practical three-dimensional aerodynamic applications. The flow analysis has been rendered by a fully implicit, finite-volume formulation of the Euler and Thin-Layer Navier-Stokes (TLNS) equations. Initially, the viscous laminar flow analysis for a wing has been compared with an independent computational fluid dynamics (CFD) code which has been extensively validated. The new procedure has been demonstrated in the design of a cranked arrow wing at Mach 2.4 with coarse- and fine-grid based computations performed with Euler and TLNS equations. The influence of the initial constraints on the geometry and aerodynamics of the optimized shape has been explored. Various final shapes generated for an identical initial problem formulation but with different optimization path options (coarse or fine grid, Euler or TLNS), have been aerodynamically evaluated via a common fine-grid TLNS-based analysis. The initial constraint conditions show significant bearing on the optimization results. Also, the results demonstrate that to produce an aerodynamically efficient design, it is imperative to include the viscous physics in the optimization procedure with the proper resolution. Based upon the present results, to better utilize the scarce computational resources, it is recommended that, a number of viscous coarse grid cases using either a preconditioned bi-conjugate gradient (PbCG) or an alternating-direction-implicit (ADI) method, should initially be employed to improve the optimization problem definition, the design space and initial shape. Optimized shapes should subsequently be analyzed using a high fidelity (viscous with fine-grid resolution) flow analysis to evaluate their true performance potential. Finally, a viscous fine-grid-based shape optimization should be conducted, using an ADI method, to accurately obtain the final optimized shape.

Pandya, Mohagna J.↗
NASA NTRS1977

Shape optimization of pressure gradient microphones

Recently developed finite element computer programs were utilized to investigate the influence of the shape of a body on its scattering field with the aim of determining the optimal shape for a Pressure Gradient Microphone (PGM). Circular cylinders of various aspect ratios were evaluated to choose the length to diameter ratio best suited for a dual element PGM application. Alterations of the basic cylindrical shape by rounding the edges and recessing at the centerline were also studied. It was found that for a + or - 1 db deviation from a linear pressure gradient response, a circular cylinder of aspect ratio near 0.5 was most suitable, yielding a useful upper frequency corresponding to ka = 1.8. The maximum increase in this upper frequency limit obtained through a number of shape alterations was only about 20 percent. An initial experimental evaluation of a single element cylindrical PGM of aspect ratio 0.18 utilizing a piezoresistive type sensor was also performed and is compared to the analytical results.

Norum, T. D.↗
NASA NTRS1999

Demonstration of Automatically-Generated Adjoint Code for Use in Aerodynamic Shape Optimization

Gradient-based optimization requires accurate derivatives of the objective function and constraints. These gradients may have previously been obtained by manual differentiation of analysis codes, symbolic manipulators, finite-difference approximations, or existing automatic differentiation (AD) tools such as ADIFOR (Automatic Differentiation in FORTRAN). Each of these methods has certain deficiencies, particularly when applied to complex, coupled analyses with many design variables. Recently, a new AD tool called ADJIFOR (Automatic Adjoint Generation in FORTRAN), based upon ADIFOR, was developed and demonstrated. Whereas ADIFOR implements forward-mode (direct) differentiation throughout an analysis program to obtain exact derivatives via the chain rule of calculus, ADJIFOR implements the reverse-mode counterpart of the chain rule to obtain exact adjoint form derivatives from FORTRAN code. Automatically-generated adjoint versions of the widely-used CFL3D computational fluid dynamics (CFD) code and an algebraic wing grid generation code were obtained with just a few hours processing time using the ADJIFOR tool. The codes were verified for accuracy and were shown to compute the exact gradient of the wing lift-to-drag ratio, with respect to any number of shape parameters, in about the time required for 7 to 20 function evaluations. The codes have now been executed on various computers with typical memory and disk space for problems with up to 129 x 65 x 33 grid points, and for hundreds to thousands of independent variables. These adjoint codes are now used in a gradient-based aerodynamic shape optimization problem for a swept, tapered wing. For each design iteration, the optimization package constructs an approximate, linear optimization problem, based upon the current objective function, constraints, and gradient values. The optimizer subroutines are called within a design loop employing the approximate linear problem until an optimum shape is found, the design loop limit is reached, or no further design improvement is possible due to active design variable bounds and/or constraints. The resulting shape parameters are then used by the grid generation code to define a new wing surface and computational grid. The lift-to-drag ratio and its gradient are computed for the new design by the automatically-generated adjoint codes. Several optimization iterations may be required to find an optimum wing shape. Results from two sample cases will be discussed. The reader should note that this work primarily represents a demonstration of use of automatically- generated adjoint code within an aerodynamic shape optimization. As such, little significance is placed upon the actual optimization results, relative to the method for obtaining the results.

Green, Lawrence↗
NASA NTRS1993

Aerodynamic shape optimization using preconditioned conjugate gradient methods

In an effort to further improve upon the latest advancements made in aerodynamic shape optimization procedures, a systematic study is performed to examine several current solution methodologies as applied to various aspects of the optimization procedure. It is demonstrated that preconditioned conjugate gradient-like methodologies dramatically decrease the computational efforts required for such procedures. The design problem investigated is the shape optimization of the upper and lower surfaces of an initially symmetric (NACA-012) airfoil in inviscid transonic flow and at zero degree angle-of-attack. The complete surface shape is represented using a Bezier-Bernstein polynomial. The present optimization method then automatically obtains supercritical airfoil shapes over a variety of freestream Mach numbers. Furthermore, the best optimization strategy examined resulted in a factor of 8 decrease in computational time as well as a factor of 4 decrease in memory over the most efficient strategies in current use.

Burgreen, Greg W.↗
NASA NTRS2004

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗
NASA NTRS2006

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗
NASA NTRS2015

Seasonal Variability of Saturn's Tropospheric Temperatures, Winds and Para-H2 from Cassini Far-IR Spectroscopy

Far-IR 16-1000 micrometer spectra of Saturn's hydrogen-helium continuum measured by Cassini's Composite Infrared Spectrometer (CIRS) are inverted to construct a near-continuous record of upper tropospheric (70-700 mbar) temperatures and para-H2 fraction as a function of latitude, pressure and time for a third of a saturnian year (2004-2014, from northern winter to northern spring). The thermal field reveals evidence of reversing summertime asymmetries superimposed onto the belt/zone structure. The temperature structure is almost symmetric about the equator by 2014, with seasonal lag times that increase with depth and are qualitatively consistent with radiative climate models. Localised heating of the tropospheric hazes (100-250 mbar) create a distinct perturbation to the temperature profile that shifts in magnitude and location, declining in the autumn hemisphere and growing in the spring. Changes in the para-H2 (f(sub p)) distribution are subtle, with a 0.02-0.03 rise over the spring hemisphere (200-500 mbar) perturbed by (i) low-f(sub p) air advected by both the springtime storm of 2010 and equatorial upwelling; and (ii) subsidence of high-f(sub p) air at northern high latitudes, responsible for a developing north-south asymmetry in f(sub p). Conversely, the shifting asymmetry in the para-H2 disequilibrium primarily reflects the changing temperature structure (and hence the equilibrium distribution of f(sub p)), rather than actual changes in f(sub p) induced by chemical conversion or transport. CIRS results interpolated to the same point in the seasonal cycle as re-analysed Voyager-1 observations (early northern spring) show qualitative consistency from year to year (i.e., the same tropospheric asymmetries in temperature and f(sub p)), with the exception of the tropical tropopause near the equatorial zones and belts, where downward propagation of a cool temperature anomaly associated with Saturn's stratospheric oscillation could potentially perturb tropopause temperatures, para-H2 and winds. Quantitative differences between the Cassini and Voyager epochs suggest that the oscillation is not in phase with the seasonal cycle at these tropospheric depths (i.e., it should be described as quasi-periodic rather than 'semi annual'). Variability in the zonal wind field derived from latitudinal thermal gradients is small (less than 10 m/s per scale height near the tropopause) and mostly affects the broad retrograde jets, with the notable exception of large variability on the northern flank of the equatorial jet. The meridional potential vorticity (PV) gradient, and hence the 'staircase of PV' associated with spatial variations in the vigour of vertical mixing, has varied over the course of the mission but maintained its overall shape. PV gradients in latitude and altitude are used to estimate the atmospheric refractive index for the propagation of stationary planetary (Rossby) waves, predicting that such wave activity would be confined to regions of real refractivity (tropical regions plus bands at 35-45 in both hemispheres). The penetration depth of these regions into the upper troposphere is temporally variable (potentially associated with stratification changes), whereas the latitudinal structure is largely unchanged over time (associated with the zonal jet system).

Fletcher, Leigh N.↗
NASA NTRS1980

The albedo of snow for partially cloudy skies

The input parameters of the model are atmospheric precipitable water, ozone content, turbidity, cloud optical thickness, size and shape of ice crystal of snow and surface pressure. The model outputs spectral and integrated solar flux snow reflectance as a function of solar elevation and fractional cloudcover. The model is illustrated using representative parameters for the Antarctic coastal regions. The albedo for a clear sky depends inversely on the solar elevation. At high elevation the albedo depends primarily upon the grain size; at low elevation this dependence is on grain size and shape. The gradient of the albedo-elevation curve increases as the grains get larger and faceted. The albedo for a dense overcast is a few percent higher than the clear sky albedo at high elevations. A simple relation between the grain size and the overcast albedo is obtained. For a set of grain size and shape, the albedo matrices (the albedo as a function of solar elevation and fractional cloudcover) are tabulated.

Choudhury, B. J.↗
NASA NTRS1973

Automated Shell Theory for Rotating Structures (ASTROS)

Computer program can be used to analyze any disk or shell of revolution of arbitrary cross section under inertial loads caused by rotation about shell axis and under various static loads, including thermal gradients. Geometric shapes incorporated in program are ellipsoidal, spherical, ogival, toroidal, conical, circular plate, cylindrical, and parabolic.

Foster, B. J.↗
NASA NTRS1982

Analytical electron microscope study of eight ataxites

Optical and electron optical (SEM, TEM, AEM) techniques were employed to investigate the fine structure of eight ataxite-iron meteorites. Structural studies indicated that the ataxites can be divided into two groups: a Widmanstaetten decomposition group and a martensite decomposition group. The Widmanstaetten decomposition group has a Type I plessite microstructure and the central taenite regions contain highly dislocated lath martensite. The steep M shaped Ni gradients in the taenite are consistent with the fast cooling rates, of not less than 500 C/my, observed for this group. The martensite decomposition group has a Type III plessite microstructure and contains all the chemical group IVB ataxites. The maximum taenite Ni contents vary from 47.5 to 52.7 wt % and are consistent with slow cooling to low temperatures of not greater than 350 C at cooling rates of not greater than 25 C/my.

Novotny, P. M.↗
NASA NTRS1941

Tests of a Gust-Alleviating Wing in the Gust Tunnel

Tests were made in the NACA gust tunnel to determine the effectiveness of a torsionally flexible wing with the torsion axis ahead of the locus of the section aerodynamic centers in reducing airplane accelerations due to atmospheric gusts. For three gust shapes, a series of flights was made with the airplane model equipped with either a torsionally flexible or a rigid wing. The results indicated that the torsionally flexible wing reduced the maximum acceleration increment 5 percent for the sharp-edge gust and about 17 percent for gust shapes with gradient distances of 6.8 and 15 chord lengths. The analysis indicated that the effectiveness of this method of gust alleviation was independent of the gust velocity and that, for the same total load increment, the torsionally flexible wing would have 10 percent less bending-moment increment at the root section of the wing than a rigid wing in all but the sharpest gusts. The results also indicated that the torsionally flexible wing slightly increased the longitudinal stability of the airplane model in a gust.

Shufflebarger, C C↗