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 505 records · Page 28

Shape optimization of single- and two-element airfoils on multiblock grids

A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Lacasse, James M.↗

Multigrid properties of upwind-biased data reconstructions

The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.

Warren, Gary P.↗

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t 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 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.

Nemec, Marian↗

Mode Propagation in Nonuniform Circular Ducts with Potential Flow

A previously reported closed form solution is expanded to determine effects of isentropic mean flow on mode propagation in a slowly converging-diverging duct, a circular cosh duct. On the assumption of uniform steady fluid density, the mean flow increases the power transmission coefficient. The increase is directly related to the increase of the cutoff ratio at the duct throat. With the negligible transverse gradients of the steady fluid variables, the conversion from one mode to another is negligible, and the power transmission coefficient remains unchanged with the mean flow direction reversed. With a proper choice of frequency parameter, many different modes can be made subject to a single value of the power transmission loss. A systematic method to include the effects of the gradients of the steady fluid variables is also described.

Cho, Y. C.↗

An implementation of the look-ahead Lanczos algorithm for non-Hermitian matrices, part 2

It is shown how the look-ahead Lanczos process (combined with a quasi-minimal residual QMR) approach) can be used to develop a robust black box solver for large sparse non-Hermitian linear systems. Details of an implementation of the resulting QMR algorithm are presented. It is demonstrated that the QMR method is closely related to the biconjugate gradient (BCG) algorithm; however, unlike BCG, the QMR algorithm has smooth convergence curves and good numerical properties. We report numerical experiments with our implementation of the look-ahead Lanczos algorithm, both for eigenvalue problem and linear systems. Also, program listings of FORTRAN implementations of the look-ahead algorithm and the QMR method are included.

Freund, Roland W.↗

Characterization of SiC/SiC Composites Subjected to Sequential Tensile Creep and Dwell Fatigue in Air under Thermal Gradient Conditions

Uncoated SiC/SiC composites, with 2D and 3D fiber architectures, were investigated under steady state thermal gradient conditions with sequential tensile creep and dwell fatigue. Melt infiltration (MI), chemical vapor infiltration (CVI), polymer impregnation and pyrolysis (PIP), and hybrid (CVI and PIP) techniques were used to fabricate ceramic matrix composites (CMCs) with a 2D balanced fiber architecture. CMCs with 3D fiber weaves were also fabricated using MI and hybrid methods. Steady state through the thickness thermal gradients were generated for the tensile specimens using laser heating on the front (hot side) and forced air cooling on the back (cold side). Hot side and cold side temperatures of 2400 and 2200 °F, respectively, were used for testing CMCs fabricated with the MI process, whereas hot side and cold side temperatures of 2700 and 2400 °F, respectively, were used for CMCs fabricated with CVI, PIP, and hybrid processes. Each CMC specimen was initially tested under tensile creep at 69 MPa (10 ksi) for 80 hours, which was then followed by tensile dwell fatigue testing at a maximum stress of 69 MPa (10 ksi) for another 80 hours for a total of 160 hours. This type of sequential tensile creep followed by tensile dwell fatigue testing was performed to ascertain the integrity of uncoated CMCs under steady state thermal gradient conditions. Tensile creep and dwell fatigue test results, generated under steady state thermal gradient conditions on different types of SiC/SiC composites, are compared with corresponding creep behaviors observed under isothermal conditions to determine the influence of both monotonic creep (static fatigue) and dwell fatigue on accumulated strains and deformation behaviors of the CMCs.

Kalluri, Sreeramesh↗

Mode propagation in nonuniform circular ducts with potential flow

A previously reported closed form solution is expanded to determine effects of isentropic mean flow on mode propagation in a slowly converging-diverging duct - a circular cosh duct. On the assumption of uniform steady fluid density, the mean flow increases the power transmission coefficient. The increase is directly related to the increase of the cutoff ratio at the duct throat. With the negligible transverse gradients of the steady fluid variables, the conversion from one mode to another is negligible, and the power transmission coefficient remains unchanged with the mean flow direction reversed. With a proper choice of frequency parameter, many different modes can be made subject to a single value of the power transmission loss. The paper also describes a systematic method to include the effects of the gradients of the steady fluid variables.

Cho, Y. C.↗

Noise Reduction Trajectory Analysis of a Supersonic Business Jet using Novel Optimization Tools

Proposals to reduce airport noise during takeoff and landing for supersonic aircraft using methods such as variable noise reduction systems add complexity to conceptual flight models. To better optimize these takeoff and landing profiles for noise certification, new modeling methods are explored in this paper with the end goal of more effectively determining the sensitivities of airframe, propulsion, and mission design variables on overall airport noise. Takeoff and landing trajectories are modeled for a notional supersonic business jet concept developed by NASA for use in environmental impact studies conducted by the International Civil Aviation Organization. Optimization tools capable of gradient-based optimal control and collocation solving methods are examined with the aim of achieving faster solutions in a more comprehensive design space. The benefits and limitations of these new methods are compared with existing methods used in previous studies of the airplane concept. It is found that the new modeling methods match closely when compared with existing tools for a standard takeoff and landing case, with a cumulative effective perceived noise difference of 0.1 EPNdB. Comparisons between takeoff trajectories with a variable noise reduction system match within 0.8 EPNdB and 1.3 EPNdB, due to highlighted differences in the modeling approaches.

supersonic↗

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗

Application of static and dynamic enclosures for determining dimethyl sulfide and carbonyl sulfide exchange in Sphagnum peatlands: Implications for the magnitude and direction of flux

A static enclosure method was applied to determine the exchange of dimethyl sulfide (DMS) and carbonyl sulfide (OCS) between the surface of Sphagnum peatlands and the atmosphere. Measurements were performed concurrently with dynamic (flow through) enclosure measurements with sulfur-free air used as sweep gas. This latter technique has been used to acquire the majority of available data on the exchange of S gases between the atmosphere and the continental surfaces and has been criticized because it is thought to overestimate the true flux of gases by disrupting natural S gas gradients. DMS emission rates determined by both methods were not statistically different between 4 and greater than 400 nmol/sq m/h, indicating that previous data on emissions of at least DMS are probably valid. However, the increase in DMS in static enclosures was not linear, indicating the potential for a negative feedback of enlosure DMS concentrations on efflux. The dynamic enclosure method measured positive OCS flux rates (emission) at all sites, while data using static enclosures indicated that OCS was consumed from the atmosphere at these same sites at rates of 3.7 to 55 nmol/sq m/h. Measurements using both enclosure techniques at a site devoid of vegetation showed that peat was a source of both DMS and OCS. However, the rate of OCS efflux from decomposing peat was more than counterbalanced by OCS consumption by vegetation, including Sphagnum mosses, and net OCS uptake occurred at all sites. We propose that all wetlands are net sinks for OCS.

De Mello, William Z.↗

Fronts propagating with curvature dependent speed: Algorithms based on Hamilton-Jacobi formulations

New numerical algorithms are devised (PSC algorithms) for following fronts propagating with curvature-dependent speed. The speed may be an arbitrary function of curvature, and the front can also be passively advected by an underlying flow. These algorithms approximate the equations of motion, which resemble Hamilton-Jacobi equations with parabolic right-hand-sides, by using techniques from the hyperbolic conservation laws. Non-oscillatory schemes of various orders of accuracy are used to solve the equations, providing methods that accurately capture the formation of sharp gradients and cusps in the moving fronts. The algorithms handle topological merging and breaking naturally, work in any number of space dimensions, and do not require that the moving surface be written as a function. The methods can be used also for more general Hamilton-Jacobi-type problems. The algorithms are demonstrated by computing the solution to a variety of surface motion problems.

Osher, Stanley↗

The use of singular value gradients and optimization techniques to design robust controllers for multiloop systems

A method for designing robust feedback controllers for multiloop systems is presented. Robustness is characterized in terms of the minimum singular value of the system return difference matrix at the plant input. Analytical gradients of the singular values with respect to design variables in the controller are derived. A cumulative measure of the singular values and their gradients with respect to the design variables is used with a numerical optimization technique to increase the system's robustness. Both unconstrained and constrained optimization techniques are evaluated. Numerical results are presented for a two-input/two-output drone flight control system.

Newsom, J. R.↗

The use of singular value gradients and optimization techniques to design robust controllers for multiloop systems

A method for designing robust feedback controllers for multiloop systems is presented. Robustness is characterized in terms of the minimum singular value of the system return difference matrix at the plant input. Analytical gradients of the singular values with respect to design variables in the controller are derived. A cumulative measure of the singular values and their gradients with respect to the design variables is used with a numerical optimization technique to increase the system's robustness. Both unconstrained and constrained optimization techniques are evaluated. Numerical results are presented for a two output drone flight control system.

Newsom, J. R.↗

Airfoil design by numerical optimization using a minicomputer

A computer program developed for the automated design of low speed airfoils utilizes a generalized Joukowski method for aerodynamic analysis coupled with a conjugate gradient, penalty function, numerical optimization algorithm to give an efficient calculation technique for use with minicomputers. The program designs airfoils with a prescribed pressure distribution as well as those which minimize or maximize some aerodynamic force coefficient. At present the method is restricted to inviscid, incompressible flow. A typical design problem will execute in 4.5 hr on an HP 9830 minicomputer.

Hicks, R. M.↗

Comparison of four stable numerical methods for Abel's integral equation

The 3-D image reconstruction from cone-beam projections in computerized tomography leads naturally, in the case of radial symmetry, to the study of Abel-type integral equations. If the experimental information is obtained from measured data, on a discrete set of points, special methods are needed in order to restore continuity with respect to the data. A new combined Regularized-Adjoint-Conjugate Gradient algorithm, together with two different implementations of the Mollification Method (one based on a data filtering technique and the other on the mollification of the kernal function) and a regularization by truncation method (initially proposed for 2-D ray sample schemes and more recently extended to 3-D cone-beam image reconstruction) are extensively tested and compared for accuracy and numerical stability as functions of the level of noise in the data.

Murio, Diego A.↗

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

An analytical and experimental study of sound propagation and attenuation in variable-area ducts

The performance of sound suppression techniques in ducts that produce refraction effects due to axial velocity gradients was evaluated. A computer code based on the method of multiple scales was used to calculate the influence of axial variations due to slow changes in the cross-sectional area as well as transverse gradients due to the wall boundary layers. An attempt was made to verify the analytical model through direct comparison of experimental and computational results and the analytical determination of the influence of axial gradients on optimum liner properties. However, the analytical studies were unable to examine the influence of non-parallel ducts on the optimum linear conditions. For liner properties not close to optimum, the analytical predictions and the experimental measurements were compared. The circumferential variations of pressure amplitudes and phases at several axial positions were examined in straight and variable-area ducts, hard-wall and lined sections with and without a mean flow. Reasonable agreement between the theoretical and experimental results was obtained.

Nayfeh, A. H.↗

Test of Shi et al. Method to Infer the Magnetic Reconnection Geometry from Spacecraft Data: MHD Simulation with Guide Field and Antiparallel Kinetic Simulation

When analyzing data from an array of spacecraft (such as Cluster or MMS) crossing a site of magnetic reconnection, it is desirable to be able to accurately determine the orientation of the reconnection site. If the reconnection is quasi-two dimensional, there are three key directions, the direction of maximum inhomogeneity (the direction across the reconnection site), the direction of the reconnecting component of the magnetic field, and the direction of rough invariance (the "out of plane" direction). Using simulated spacecraft observations of magnetic reconnection in the geomagnetic tail, we extend our previous tests of the direction-finding method developed by Shi et al. (2005) and the method to determine the structure velocity relative to the spacecraft Vstr. These methods require data from four proximate spacecraft. We add artificial noise and calibration errors to the simulation fields, and then use the perturbed gradient of the magnetic field B and perturbed time derivative dB/dt, as described by Denton et al. (2010). Three new simulations are examined: a weakly three-dimensional, i.e., quasi-two-dimensional, MHD simulation without a guide field, a quasi-two-dimensional MHD simulation with a guide field, and a two-dimensional full dynamics kinetic simulation with inherent noise so that the apparent minimum gradient was not exactly zero, even without added artificial errors. We also examined variations of the spacecraft trajectory for the kinetic simulation. The accuracy of the directions found varied depending on the simulation and spacecraft trajectory, but all the directions could be found within about 10 for all cases. Various aspects of the method were examined, including how to choose averaging intervals and the best intervals for determining the directions and velocity. For the kinetic simulation, we also investigated in detail how the errors in the inferred gradient directions from the unmodified Shi et al. method (using the unperturbed gradient) depended on the amplitude of the calibration errors. For an accuracy of 3 for the maximum gradient direction, the calibration errors could be as large as 3% of reconnection magnetic field, while for the same accuracy for the minimum gradient direction, the calibration errors could only be as large as 0.03% of the reconnection magnetic field. These results suggest that the maximum gradient direction can normally be determined by the unmodified Shi et al. method, while the modified method or some other method must be used to accurately determine the minimum gradient direction. The structure velocity was found with magnitude accurate to 2% and direction accurate to within 5%.

Denton, R.↗