SEARCH · Engineering Papers

Results for “analytical gradients”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Active control law synthesis for flexible aircraft

An application of an active control law synthesis procedure is presented, to meet multiple design requirements for a flexible aircraft modeled by a large order state space system. In this synthesis procedure, a linear quadratic Gaussian type cost function is minimized by updating the free parameters of the control law, while satisfying a set of constraints on the design loads, responses and stability margins. Analytical expressions for gradients of the cost function and the constraints, with respect to the control law design variables are used to facilitate rapid numerical convergence. These gradients can also be used for sensitivity study. A stable classical control law as well as an estimator-based full or reduced order control law can be modified, in order to meet individual root-mean-square response limitations as well as minimum singular value restrictions. Both analog and digital control laws can be optimized. Low order, robust control laws were synthesized for flutter suppression of a flexible aircraft.

Mukhopadhyay, V.↗

NASA NTRS1992

Metastable State Relaxation in a Gravitational Field

A metastable state relaxation equation for a physical system placed into a gravitational field is constructed for non-critical supersaturated solutions which arc in the immediate neighborhood of the coexistence line. Solutions of this equation are obtained in two different regimes: stationary and dynamic. The sedimentation time which can be defined as the time of the subcritical solute cluster redistribution corresponding to the final steady state in the gravitational field is found. The formation of the concentration gradient is proved analytically and its expression through the model parameters is obtained. The following analysis gives the expression for the sedimentation time which does not depend on the column height. The law of the concentration change with respect to the column height is also found and analyzed.

Izmailov, Alexander F.↗

NASA NTRS2001

Quantitation of promethazine and metabolites in urine samples using on-line solid-phase extraction and column-switching

A chromatographic method for the quantitation of promethazine (PMZ) and its three metabolites in urine employing on-line solid-phase extraction and column-switching has been developed. The column-switching system described here uses an extraction column for the purification of PMZ and its metabolites from a urine matrix. The extraneous matrix interference was removed by flushing the extraction column with a gradient elution. The analytes of interest were then eluted onto an analytical column for further chromatographic separation using a mobile phase of greater solvent strength. This method is specific and sensitive with a range of 3.75-1400 ng/ml for PMZ and 2.5-1400 ng/ml for the metabolites promethazine sulfoxide, monodesmethyl promethazine sulfoxide and monodesmethyl promethazine. The lower limits of quantitation (LLOQ) were 3.75 ng/ml with less than 6.2% C.V. for PMZ and 2.50 ng/ml with less than 11.5% C.V. for metabolites based on a signal-to-noise ratio of 10:1 or greater. The accuracy and precision were within +/- 11.8% in bias and not greater than 5.5% C.V. in intra- and inter-assay precision for PMZ and metabolites. Method robustness was investigated using a Plackett-Burman experimental design. The applicability of the analytical method for pharmacokinetic studies in humans is illustrated.

NASA Center JSC↗

DOE OSTI2021

Correspondence between neuroevolution and gradient descent

Abstract We show analytically that training a neural network by conditioned stochastic mutation or neuroevolution of its weights is equivalent, in the limit of small mutations, to gradient descent on the loss function in the presence of Gaussian white noise. Averaged over independent realizations of the learning process, neuroevolution is equivalent to gradient descent on the loss function. We use numerical simulation to show that this correspondence can be observed for finite mutations, for shallow and deep neural networks. Our results provide a connection between two families of neural-network training methods that are usually considered to be fundamentally different.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Absence of Barren Plateaus and Scaling of Gradients in the Energy Optimization of Isometric Tensor Network States

Abstract Vanishing gradients can pose substantial obstacles for high-dimensional optimization problems. Here we consider energy minimization problems for quantum many-body systems with extensive Hamiltonians and finite-range interactions, which can be studied on classical computers or in the form of variational quantum eigensolvers on quantum computers. Barren plateaus correspond to scenarios where the average amplitude of the energy gradient decreases exponentially with increasing system size. This occurs, for example, for quantum neural networks and for brickwall quantum circuits when the depth increases polynomially in the system size. Here we prove that the variational optimization problems for matrix product states, tree tensor networks, and the multiscale entanglement renormalization ansatz are free of barren plateaus. The derived scaling properties for the gradient variance provide an analytical guarantee for the trainability of randomly initialized tensor network states (TNS) and motivate certain initialization schemes. In a suitable representation, unitary tensors that parametrize the TNS are sampled according to the uniform Haar measure. We employ a Riemannian formulation of the gradient based optimizations which simplifies the analytical evaluation.

Barthel, Thomas↗

NASA NTRS1993

Phase-retrieval algorithms for a complicated optical system

Phase-retrieval algorithms have been developed that handle a complicated optical system that requires multiple Fresnellike transforms to propagate from one end of the system to the other including the absorption by apertures in more than one plane and allowance for bad detector pixels. Gradient-search algorithms and generalizations of the iterative-transform phase-retrieval algorithms are derived. Analytic expressions for the gradient of an error metric, with respect to polynomial coefficients and with respect to point-by-point phase descriptions, are given. The entire gradient can be computed with the number of transforms required to propagate a wave front from one end of the optical system to the other and back again, independent of the number of coefficients or phase points. This greatly speeds the computation. The reconstruction of pupil amplitude is also given. A convergence proof of the generalized iterative transform algorithm is given. These improved algorithms permit a more accurate characterization of complicated optical systems from their point spread functions.

Fienup, J. R.↗

DOE OSTI2021

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Corrections to fringe distortion due to flow density gradients in optical interferometry

An analytical method is formulated to account for distortions of optical interferograms used for studies of flow over airfoils experiencing dynamic stall. It is shown that such distortions are generated primarily due to optical path deflections in the test flow, caused by large density gradients. Such in-flow optical path deflections are neglected in conventional optical techniques for flow studies. The present method employs a ray analysis to determine these in-flow optical path deflections, and accurately predicts the interferogram distortions.

Cho, Y. C.↗

DOE OSTI2020

Formulation and solution approach for calibrating activity-based travel demand model-system via microsimulation

This study addresses the problem of calibrating utility-maximizing nested logit activity-based travel demand model-systems. After estimation, it is common practice to use aggregate measurements to calibrate the estimated model-system’s parameters prior to their application in transportation planning, policy making, and operations. However, calibration of activity-based model-systems has received much less attention. Existing calibration approaches are myopic heuristics in the sense that they do not consider the fundamental inter-dependencies among choice-models and do not have a systematic way to adjust model parameters. Also, other purely simulation-based approaches do not perform well in large-scale applications. In this study, we focus on utility-maximizing nested logit activity-based model-systems and calibrating aggregate statistics such as activity shares, mode shares, time-dependent & mode-specific OD flows, and time-dependent & mode-specific sensor counts. We formulate the calibration problem as a simulation-based optimization problem and propose a stochastic gradient-based solution procedure to solve it. The solution procedure relies on microsimulation to calculate expectations of the aggregate statistics of interest to the calibration problem. Additionally, we derive approximate analytical expressions for the gradient of the objective function —that are evaluated through microsimulation on mini-batches of the population. The proposed solution procedure is sensitive to the fundamental structure of the activity-based model-system and is non-myopic in considering the dependencies across its model components. The formulated optimization problem is non-convex, highly nonlinear, and potentially has multiple-minima. Lastly, we show —through a real-world application— that the proposed solution procedure outperforms other state-of-the-art purely simulation-based optimization approaches in terms of computational efficiency, stability, and convergence. We also compare various gradient-based solution algorithms to determine the best algorithm to update the parameters. This work has the potential to facilitate wider and easier application of activity-based model-systems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1973

Transonic laminar boundary layers with surface curvature.

The effect of surface curvature (both longitudinal and transverse) and the associated pressure gradient across the flow is investigated analytically for a laminar boundary layer subjected to pressure gradients along the flow. Property variation which results from heat transfer and compressibility is taken into account. Numerical solutions of the boundary layer equations are obtained for locally similar sonic flow through the throat of a nozzle for a range of flow conditions and for various shaped nozzle surfaces with different amounts of wall cooling. A few solutions were also obtained for the analogous flow around the shoulder of a flat-faced body in a supersonic flow. The effect of various parameters that arise in the equations upon application of the Levy-Mangler transformation are investigated and discussed with respect to their influence on the velocity and total enthalpy profiles and the corresponding profile slopes at the surface to which the shear stress and heat transfer are related. An important finding is that at throat Reynolds numbers less than 100,000 the heat transfer parameter at a nozzle throat decreases as the throat radius of curvature decreases.

Back, L. H.↗

NASA NTRSDate not supplied

Gradient-Based Optimization of the Common Research Model Wing Subject to CFD-Based Gust and Flutter Constraints

The linearized frequency-domain method was recently implemented in the stabilized finite element solver in NASA’s FUN3D code. Previous work by the authors used this method for enforcing flutter constraints during gradient-based optimizations. More recently, the solver was expanded to account for continuous (also known as stochastic) gust responses. This paper expands on recent Common Research Model wing optimization work, which demonstrated gradient-based optimization with flutter and stochastic gust constraints, among others. While that work utilized FUN3D for static aeroelastic solutions but relied on doublet lattice aerodynamics for gust and flutter responses, the present work replaces these unsteady aerodynamic analyses with those of FUN3D’s linearized frequency-domain solver. With analytic derivatives available, gradient-based optimization is performed through the use of the OpenMDAO/MPhys libraries with over 700 shape, structural, and aerodynamic design variables and over 10 nonlinear constraints. Comparisons of analysis results and optimized designs are made between doublet lattice and linearized frequency-domain solutions.

aeroelasticity↗

NASA NTRS2016

Parallel Aircraft Trajectory Optimization with Analytic Derivatives

Trajectory optimization is an integral component for the design of aerospace vehicles, but emerging aircraft technologies have introduced new demands on trajectory analysis that current tools are not well suited to address. Designing aircraft with technologies such as hybrid electric propulsion and morphing wings requires consideration of the operational behavior as well as the physical design characteristics of the aircraft. The addition of operational variables can dramatically increase the number of design variables which motivates the use of gradient based optimization with analytic derivatives to solve the larger optimization problems. In this work we develop an aircraft trajectory analysis tool using a Legendre-Gauss-Lobatto based collocation scheme, providing analytic derivatives via the OpenMDAO multidisciplinary optimization framework. This collocation method uses an implicit time integration scheme that provides a high degree of sparsity and thus several potential options for parallelization. The performance of the new implementation was investigated via a series of single and multi-trajectory optimizations using a combination of parallel computing and constraint aggregation. The computational performance results show that in order to take full advantage of the sparsity in the problem it is vital to parallelize both the non-linear analysis evaluations and the derivative computations themselves. The constraint aggregation results showed a significant numerical challenge due to difficulty in achieving tight convergence tolerances. Overall, the results demonstrate the value of applying analytic derivatives to trajectory optimization problems and lay the foundation for future application of this collocation based method to the design of aircraft with where operational scheduling of technologies is key to achieving good performance.

aircraft↗

NASA NTRS1995

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗

NASA NTRS2017

Constrained Burn Optimization for the International Space Station

In long-term trajectory planning for the International Space Station (ISS), translational burns are currently targeted sequentially to meet the immediate trajectory constraints, rather than simultaneously to meet all constraints, do not employ gradient-based search techniques, and are not optimized for a minimum total deltav (v) solution. An analytic formulation of the constraint gradients is developed and used in an optimization solver to overcome these obstacles. Two trajectory examples are explored, highlighting the advantage of the proposed method over the current approach, as well as the potential v and propellant savings in the event of propellant shortages.

Brown, Aaron J.↗

DOE OSTI2022

Accurate parameterization of the kinetic energy functional

The absence of a reliable formulation of kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions which have been proposed in the literature and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals - a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since, the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H 2 , LiH, LiF and a one-dimensional chain of 8 hydrogens using data from Kohn-Sham density functional theory calculations for only a few structures. For the ease of reproduction, codes used to generate the models are provided in the supporting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1976

Solving the optimal attention allocation problem in manual control

Within the context of the optimal control model of human response, analytic expressions for the gradients of closed-loop performance metrics with respect to human operator attention allocation are derived. These derivatives serve as the basis for a gradient algorithm that determines the optimal attention that a human should allocate among several display indicators in a steady-state manual control task. Application of the human modeling techniques are made to study the hover control task for a CH-46 VTOL flight tested by NASA.

Kleinman, D. L.↗

NASA NTRS1989

Digital robust control law synthesis using constrained optimization

Development of digital robust control laws for active control of high performance flexible aircraft and large space structures is a research area of significant practical importance. The flexible system is typically modeled by a large order state space system of equations in order to accurately represent the dynamics. The active control law must satisy multiple conflicting design requirements and maintain certain stability margins, yet should be simple enough to be implementable on an onboard digital computer. Described here is an application of a generic digital control law synthesis procedure for such a system, using optimal control theory and constrained optimization technique. A linear quadratic Gaussian type cost function is minimized by updating the free parameters of the digital control law, while trying to satisfy a set of constraints on the design loads, responses and stability margins. Analytical expressions for the gradients of the cost function and the constraints with respect to the control law design variables are used to facilitate rapid numerical convergence. These gradients can be used for sensitivity study and may be integrated into a simultaneous structure and control optimization scheme.

Mukhopadhyay, Vivekananda↗

NASA NTRS1988

Digital robust active control law synthesis for large order flexible structure using parameter optimization

A generic procedure for the parameter optimization of a digital control law for a large-order flexible flight vehicle or large space structure modeled as a sampled data system is presented. A linear quadratic Guassian type cost function was minimized, while satisfying a set of constraints on the steady-state rms values of selected design responses, using a constrained optimization technique to meet multiple design requirements. Analytical expressions for the gradients of the cost function and the design constraints on mean square responses with respect to the control law design variables are presented.

Mukhopadhyay, V.↗