Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adjoint 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ecosystem biogeochemistry model parameterization: Do more flux data result in a better model in predicting carbon flux?

Reliability of terrestrial ecosystem models highly depends on the quantity and quality of thedata that have been used to calibrate the models. Nowadays, in situ observations of carbon fluxes areabundant. However, the knowledge of how much data (data length) and which subset of the time seriesdata (data period) should be used to effectively calibrate the model is still lacking. This study uses theAmeriFlux carbon flux data to parameterize the Terrestrial Ecosystem Model (TEM) with an adjoint-baseddata assimilation technique for various ecosystem types. Parameterization experiments are thus conductedto explore the impact of both data length and data period on the uncertainty reduction of the posteriormodel parameters and the quantification of site and regional carbon dynamics. We find that: the modelis better constrained when it uses two-year data comparing to using one-year data. Further, two-year datais sufficient in calibrating TEM’s carbon dynamics, since using three-year data could only marginallyimprove the model performance at our study sites; the model is better constrained with the data thathave a higher‘‘climate variability’’than that having a lower one. The climate variability is used to measurethe overall possibility of the ecosystem to experience all climatic conditions including drought and extremeair temperatures and radiation; the U.S. regional simulations indicate that the effect of calibration datalength on carbon dynamics is amplified at regional and temporal scales, leading to large discrepanciesamong different parameterization experiments, especially in July and August. Our findings areconditioned on the specific model we used and the calibration sites we selected. The optimal calibrationdata length may not be suitable for other models. However, this study demonstrates that there may exist athreshold for calibration data length and simply using more data would not guarantee a better modelparameterization and prediction. More importantly, climate variability might be an effective indicator ofinformation within the data, which could help data selection for model parameterization. As a result, we believe ourfindings will benefit the ecosystem modeling community in using multiple-year data to improve modelpredictability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Nonlinear Optimal Control of Electron Dynamics Within Hartree-Fock Theory

Consider the problem of determining the optimal applied electric field to drive a molecule from an initial state to a desired target state. For even moderately sized molecules, solving this problem directly using the exact equations of motion—the time-dependent Schrödinger equation (TDSE)—is numerically intractable. Here, we present a solution of this problem within time-dependent Hartree-Fock (TDHF) theory, a mean field approximation of the TDSE. Optimality is defined in terms of minimizing the total control effort while maximizing the overlap between desired and achieved target states. We frame this problem as an optimization problem constrained by the nonlinear TDHF equations; we solve it using trust region optimization with gradients computed via a custom-built adjoint state method. For three molecular systems, we show that with very small neural network parametrizations of the control, our method yields solutions that achieve desired targets within acceptable constraints and tolerances.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Sensitivity derivatives for static test loading boundary conditions

Attention is given to the application of the adjoint variable method under the special combination of displacement and stress boundary conditions which occur when static test conditions are simulated. In the two-step problem-solving method employed, calculation of internal stresses, displacements and total reaction force at the boundary is followed by the calculation of the applied load and the total reaction force due to unit displacements.

Gurdal, Z.↗
NASA NTRS1987

Shape design sensitivity analysis and optimal design of structural systems

The material derivative concept of continuum mechanics and an adjoint variable method of design sensitivity analysis are used to relate variations in structural shape to measures of structural performance. A domain method of shape design sensitivity analysis is used to best utilize the basic character of the finite element method that gives accurate information not on the boundary but in the domain. Implementation of shape design sensitivty analysis using finite element computer codes is discussed. Recent numerical results are used to demonstrate the accuracy obtainable using the method. Result of design sensitivity analysis is used to carry out design optimization of a built-up structure.

Choi, Kyung K.↗
NASA NTRS1989

On the accuracy of shape sensitivity derivatives

Six methods for calculating shape sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the finite-element mesh is refined, displacement sensitivity derivatives converge more slowly than the displacement themselves. Five of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer, rather than reduced accuracy.

Haftka, Raphael T.↗
NASA NTRS1991

On the accuracy of shape sensitivity

The calculation of sensitivity of the response of a structure modeled by finite elements to shape variation is known to be subject to numerical difficulties. The accuracy of a given method is typically measured against the yard stick of finite-difference sensitivity calculation. The present paper demonstrates with a simple example that this approach may be flawed because of discretization errors associated with the finite element mesh. Seven methods for calculating sensitivity derivatives are compared for a two-material beam problem with a moving interface. It is found that as the mesh is refined, displacement sensitivity derivatives converge more slowly than the displacements. Six of the methods agree fairly well, but the adjoint variational surface method provides substantially different results. However, the difference is found to reflect convergence from another direction to the same answer rather than reduced accuracy. Additionally, it is observed that small derivatives are particularly prone to accuracy problems.

Haftka, R. T.↗
NASA NTRS1991

Discrete shape sensitivity equations for aerodynamic problems

Intensive research and development in Computational Fluid Dynamics (CFD) has recently produced many powerful CFD codes to simulate complex aerodynamic phenomena. However, in order to enhance the usefulness of these CFD codes for design practice, development of design sensitivity equations compatible to these codes becomes very important. This paper represents a part of such an effort to develop a sensitivity analysis methodology that enables the sensitivity equations to be implemented into existing CFD codes with minimal code modification. The methodology is based upon a preelimination procedure which accounts for consistently linearized boundary conditions. Formulations of both the direct differentiation and the adjoint variable methods will be presented in the paper.

Hou, Gene J.-W.↗
NASA NTRS1993

Sensitivity analysis of Lyapunov and Riccati equations with application to controls-structures integrated design

Controls-structure integrated design is a complicated multidisciplinary design optimization problem which involves the state equations pertaining to open-loop eigenvalues and control laws. In order to alleviate the intensity of the computation, this study uses the adjoint variable method to derive sensitivity equations for the eigenvalue, Liapunov, and Riccati equations. These individual sensitivity equations are then combined together to form the multidisciplinary sensitivity equations for the control structure integrated design problems. A set of linear sensitivity equations, proportional in number to the number of performance functions involved in the optimization process, are solved. This proposed approach may provide a great saving in computer resources. The validity of the newly developed sensitivity equations is verified by numerical examples.

Hou, Gene↗
NASA NTRS1993

Formulation for Simultaneous Aerodynamic Analysis and Design Optimization

An efficient approach for simultaneous aerodynamic analysis and design optimization is presented. This approach does not require the performance of many flow analyses at each design optimization step, which can be an expensive procedure. Thus, this approach brings us one step closer to meeting the challenge of incorporating computational fluid dynamic codes into gradient-based optimization techniques for aerodynamic design. An adjoint-variable method is introduced to nullify the effect of the increased number of design variables in the problem formulation. The method has been successfully tested on one-dimensional nozzle flow problems, including a sample problem with a normal shock. Implementations of the above algorithm are also presented that incorporate Newton iterations to secure a high-quality flow solution at the end of the design process. Implementations with iterative flow solvers are possible and will be required for large, multidimensional flow problems.

Hou, G. W.↗
NASA NTRS2001

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗
NASA NTRS2016

Challenges and Progress in Aerodynamic Design of Hybrid Wingbody Aircraft with Embedded Engines

We summarize the contributions to high-fidelity capabilities for analysis and design of hybrid wingbody (HWB) configurations considered by NASA. Specifically, we focus on the embedded propulsion concepts of the N2-B and N3-X configurations, some of the future concepts seriously investigated by the NASA Fixed Wing Project. The objective is to develop the capability to compute the integrated propulsion and airframe system realistically in geometry and accurately in flow physics. In particular, the propulsion system (including the entire engine core-compressor, combustor, and turbine stages) is vastly more difficult and costly to simulate with the same level of fidelity as the external aerodynamics. Hence, we develop an accurate modeling approach that retains important physical parameters relevant to aerodynamic and propulsion analyses for evaluating the HWB concepts. Having the analytical capabilities at our disposal, concerns and issues that were considered to be critical for the HWB concepts can now be assessed reliably and systematically; assumptions invoked by previous studies were found to have serious consequences in our study. During this task, we establish firmly that aerodynamic analysis of a HWB concept without including installation of the propulsion system is far from realistic and can be misleading. Challenges in delivering the often-cited advantages that belong to the HWB are the focus of our study and are emphasized in this report. We have attempted to address these challenges and have had successes, which are summarized here. Some can have broad implications, such as the concept of flow conditioning for reducing flow distortion and the modeling of fan stages. The design optimization capability developed for improving the aerodynamic characteristics of the baseline HWB configurations is general and can be employed for other applications. Further improvement of the N3-X configuration can be expected by expanding the design space. Finally, the support of the System Analysis and Integration Element under the NASA Fixed Wing Project has enabled the development and helped deployment of the capabilities shown in this report.

design optimization↗
NASA NTRS2017

Aerodynamic Design of the Hybrid Wing Body Propulsion-Airframe Integration

A hybrid wingbody (HWB) concept is being considered by NASA as a potential subsonic transport aircraft that meets aerodynamic, fuel, emission, and noise goals in the time frame of the 2030s. While the concept promises advantages over conventional wing-and-tube aircraft, it poses unknowns and risks, thus requiring in-depth and broad assessments. Specifically, the configuration entails a tight integration of the airframe and propulsion geometries; the aerodynamic impact has to be carefully evaluated. With the propulsion nacelle installed on the (upper) body, the lift and drag are affected by the mutual interference effects between the airframe and nacelle. The static margin for longitudinal stability is also adversely changed. We develop a design approach in which the integrated geometry of airframe (HWB) and propulsion is accounted for simultaneously in a simple algebraic manner, via parameterization of the planform and airfoils at control sections of the wingbody. In this paper, we present the design of a 300-passenger transport that employs distributed electric fans for propulsion. The trim for stability is achieved through the use of the wingtip twist angle. The geometric shape variables are determined through the adjoint optimization method by minimizing the drag while subject to lift, pitch moment, and geometry constraints. The design results clearly show the influence on the aerodynamic characteristics of the installed nacelle and trimming for stability. A drag minimization with the trim constraint yields a reduction of 10 counts in the drag coefficient.

hybrid wing↗
NASA NTRS2018

Aerodynamic Design and Optimization of Fan Stage for Boundary Layer Ingestion Propulsion System

The present paper addresses the process of preliminary design of a low-pressure fan and outlet guide vane (OGV) of a boundary layer ingestion (BLI) propulsion system. The tail-cone thruster systems of NASA's STARC_ABL (Single-aisle Turboelectric Aircraft with an Aft Boundary-Layer propulsor) adopts an axi-symmetric BLI type inlet as opposed to other embedded engine systems. Thus, the focus of the present work is placed on maximizing the efficiency of the fan and OGV stages under a significant radial distortion. A parameterization with B-spline function for camber line angles, metal chord, thickness distribution and stacking axis of blades is presented. The flowpath lines are also parameterized by B-spline function and aggregated in the design system of blades. The design optimization with evolutionary algorithm is performed with constraints of fan pressure ratio, OGV exit swirl angle and nozzle exit properties. The inlet conditions for the turbo-machinery CFD (Computational Fluid Dynamics) domain and the design goal of the fan stage are driven by a propulsion airframe integration (PAI) model that uses a 3-D unstructured RANS (Reynolds Average Navier Stokes) solver and actuator disk model. The expected power saving of the BLI propulsor is quantified via PAI analysis and the resulting preliminary design of the fan stages is compared with a clean-inlet flow propulsor.

Turbo-machinery↗
NASA NTRSDate not supplied

Goal-Oriented Discretization Error Control in Coupled Nearfield-Farfield Low-Boom Simulations

The method of adjoint weighted residuals is used to determine the level of discretization error in loudness predictions of sonic booms on the ground. We analyze the standard nearfied-farfield domain decomposition approach. In the nearfield domain, the three-dimensional Euler equations are solved to obtain a pressure signature generated by the aircraft. In the farfield, this waveform is propagated through the atmosphere to the ground by solving the augmented Burgers’ equation. Loudness is characterized using weighted sound-exposure-level metrics. We formulate discretization error estimates for the ground signature and loudness metrics for this one-way coupled system. Although the nearfield solution is independent of the farfield, the adjoint formulation for the coupled system provides feedback from the farfield to identify high-error regions in the nearfield. The results demonstrate that the discrete adjoint implementation is asymptotically consistent and provides reliable error estimates. Furthermore, we show how the error can be controlled through adaptive refinement of the nearfield mesh. The approach is evaluated on two- and three-dimensional problems, including the X-59 flight demonstration aircraft.

CST↗
DOE OSTI2024

Exploring gauge-fixing conditions with gradient-based optimization

Lattice gauge fixing is required to compute gauge-variant quantities, for example those used in RI-MOM renormalization schemes or as objects of comparison for model calculations. Recently, gauge-variant quantities have also been found to be more amenable to signal-to-noise optimization using contour deformations. These applications motivate systematic parameterization and exploration of gauge-fixing schemes. This work introduces a differentiable parameterization of gauge fixing which is broad enough to cover Landau gauge, Coulomb gauge, and maximal tree gauges. The adjoint state method allows gradient-based optimization to select gauge-fixing schemes that minimize an arbitrary target loss function.

Detmold, William↗