Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic 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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

Development and application of a gradient method for solving differential games

A technique for solving n-dimensional games is developed and applied to two pursuit-evasion games. The first is a two-dimensional game similar to the homicidal chauffeur but modified to resemble an airplane-helicopter engagement. The second is a five-dimensional game of two airplanes at constant altitude and with thrust and turning controls. The performance function to be optimized by the pursuer and evader was the distance between the evader and a given target point in front of the pursuer. The analytic solution to the first game reveals that both unique and nonunique solutions exist. A comparison between the gradient results and the analytic solution shows a dependence on the nominal controls in regions where nonunique solutions exist. In the unique solution region, the results from the two methods agree closely. The results for the five-dimensional two-airplane game are also shown to be dependent on the nominal controls selected and indicate that initial conditions are in a region of nonunique solutions.

Roberts, D. A.↗
NASA NTRSDate not supplied

Fast Aircraft Separation Calculations for Gradient Based Optimization of Airspace Simulations

Simulations of airspace operational concepts can play a significant role in determining future paradigms that would allow for a safe increase in airspace density. In particular, airspace simulations which are capable of handling large numbers of aircraft act as an enabling capability for the testing of proposed airspace operational concepts. Simulations allowing for gradient based optimization methods are particularly attractive, since they would potentially allow for an efficient and empirical means to derive best operational practices. These could also allow for vehicle multidisciplinary design and optimization studies to include air traffic management considerations as a discipline. But any large scale simulation of airspace operations must include some methodology for addressing airspace separation requirements, which in the most direct sense would be tracked in a manner that computationally grows as a quadratic function of the number of simulated aircraft. Efficient indirect methods have been developed in certain contexts to address this limitation. However, any means of addressing separation requirements in a gradient based optimization context should be implemented by functions which provide analytic derivative information to maximize numerical precision and computational efficiency. In this paper, a fast and differentiable separation metric is described in application to gradient based optimization of airspace operations. Rather than computing the separation distance between every pair of aircraft in a simulation, this method effectively reduces the problem to a smaller relevant set using a geometric decomposition. This method guarantees that the smallest distance at all points in simulated time is determined exactly. When used in an optimization constraint context, this guarantees that a minimum separation is maintained between all pairs of aircraft. The presented metric has logarithmic computational growth with respect to the number of simulated aircraft, and is shown to perform well in a series of notional 2D airspace optimization problems when used to enforce specified airborne separation constraints. Results show that this is notably faster than a direct pairwise distance computing metric for optimizations involving both small and large numbers of aircraft, yet enforce separation requirements to the same tolerance. It is shown that this favorable scalability is an enabling capability for more sophisticated air traffic management conceptual studies.

Optimization↗
NASA NTRS1984

Consequences of a chromospheric temperature gradient on the width of H Alpha in late-type giants

An analytic expression for the integrated H alpha optical depth profile is derived for a one dimensional slab geometry model chromosphere, with electron temperature increasing as a power law with height. The formula predicts H alpha opacity and profile width to be sensitive functions of the thermal gradient. Application of the model to observation reveals that broad H alpha absorption widths in G and K giant stars are consistent with a mean H alpha chromospheric optical depth of 50, while narrower widths in M stars indicate slightly lower opacities. It is proposed that differences in H alpha width between late-type giants of similar spectral type may be due, in part, to differences in their chromospheric thermal gradient, and associated H alpha opacity.

Zarro, D. M.↗
NASA NTRS1981

Structural analysis of cylindrical thrust chambers, volume 2

Analytical results are presented which predict cumulative plastic deformation characteristic of damage observed in coolant channel walls of regeneratively cooled rocket thrust chambers. The damage consists of bulging and plastic flow which leads to thinout and rupture of the channel wall under the combined effects of high pressures, high temperatures and temperature gradients experienced during cyclic firings of actual test chambers. Analytical predictions correlate with test results.

Armstrong, W. H.↗
DOE OSTI2022

The evolution of composition and morphology during the initial growth of electrodeposited Ni-Fe films: Comparison between the potentiostatic mode and the pulse-reverse potential mode

The anomalous codeposition of Ni-Fe tends to develop a composition gradient zone along the growth direction within the near-substrate region. Although it has been shown that such a composition gradient can be mitigated by the pulsed deposition method, a fundamental understanding of this behavior based on the anomalous codeposition mechanism has not been achieved yet. In this work, the evolution of the composition and the morphology during the (pulsed) electrodeposition of Ni-Fe films was reported. Further, a theory based on the Fe deposition mechanism originally described by Bockris, Drazic, and Despic was proposed and described in detail, which was used to rationalize not only the occurrence and the decline of the composition gradient, but also the opposite signs of the composition gradient observed under the potentiostatic mode and the pulse-reverse potential mode.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 NTRS1993

Correlation of analytical and experimental hot structure vibration results

High surface temperatures and temperature gradients can affect the vibratory characteristics and stability of aircraft structures. Aircraft designers are relying more on finite-element model analysis methods to ensure sufficient vehicle structural dynamic stability throughout the desired flight envelope. Analysis codes that predict these thermal effects must be correlated and verified with experimental data. Experimental modal data for aluminum, titanium, and fiberglass plates heated at uniform, nonuniform, and transient heating conditions are presented. The data show the effect of heat on each plate's modal characteristics, a comparison of predicted and measured plate vibration frequencies, the measured modal damping, and the effect of modeling material property changes and thermal stresses on the accuracy of the analytical results at nonuniform and transient heating conditions.

Kehoe, Michael W.↗
NASA NTRS2005

Motion of a Drop on a Solid Surface Due to a Wettability Gradient

The hydrodynamic force experienced by a spherical-cap drop moving on a solid surface is obtained from two approximate analytical solutions and used to predict the quasi-steady speed of the drop in a wettability gradient. One solution is based on approximation of the shape of the drop as a collection of wedges, and the other is based on lubrication theory. Also, asymptotic results from both approximations for small contact angles, as well as an asymptotic result from lubrication theory that is good when the length scale of the drop is large compared with the slip length, are given. The results for the hydrodynamic force also can be used to predict the quasi-steady speed of a drop sliding down an incline.

Subramanian, R.↗
NASA NTRS2016

Optimization of Turbine Engine Cycle Analysis with Analytic Derivatives

A new engine cycle analysis tool, called Pycycle, was built using the OpenMDAO framework. Pycycle provides analytic derivatives allowing for an efficient use of gradient-based optimization methods on engine cycle models, without requiring the use of finite difference derivative approximation methods. To demonstrate this, a gradient-based design optimization was performed on a turbofan engine model. Results demonstrate very favorable performance compared to an optimization of an identical model using finite-difference approximated derivatives.

OpenMDAO↗
DOE OSTI2026

Onset of High-Harmonic Fast Wave Propagation in Annulus Geometry

Previous experiments on NSTX have reported high-harmonic fast wave (HHFW) power losses in the scrape-off layer (SOL), with dependencies on the wave parallel wavenumber, plasma density, and magnetic field. This study presents an analysis of HHFW propagation in the SOL plasma, focusing on parasitic coupling to the waveguide mode bounded by the steep gradient pedestal and the outer wall boundary. The analytic and finite element wave models are discussed using a 2D annulus geometry, which serves as a simplified representation of the NSTX-U plasma. The model shows that a minimum radial gap width is required to support wave propagation, analogous to the higher-order modes in coaxial geometries. The impact of several factors, including gap width, wave parallel phasing, and magnetic field strength, on parasitic wave coupling is examined, highlighting the role of the geometric dimension in relation to the wave perpendicular wavelength. The findings suggest that controlling the edge/SOL density profile (gap control) is crucial for mitigating parasitic coupling when the density in front of the antenna exceeds the density cutoff.

Baek, Seung-Gyou [Massachusetts Inst. of Technolog↗
DOE OSTI2022

Probing triaxial deformation of atomic nuclei in high-energy heavy ion collisions

Most atomic nuclei are deformed with a quadrupole shape described by its overall strength β 2 and triaxiality γ. The deformation can be accessed in high-energy heavy ion collisions by measuring the collective flow response of the produced quark-gluon plasma to the eccentricity ε 2 and the density gradient d ⊥ in the initial state. Using an analytical estimate and a Glauber model, I show that the variances $\langleε^{2}_{2}\rangle$ or $\langle$(δd ⊥ /d ⊥ ) 2 $\rangle$ and skewnesses $\langleε^{2}_{2}δd_{⊥}/d_{⊥}\rangle$ or $\langle(δd_{⊥}/d_{⊥})^{3}\rangle$ have a simple analytical form of a'+b'β$^{2}_{2}$ and a'+[b'+c' cos(3γ)]β$^{3}_{2}$, respectively. From these, I constructed several normalized skewnesses to isolate the γ dependence from that of β 2 and show that the correlations between a normalized skewness and a variance can constrain simultaneously β 2 and γ. Assuming a linear relation with elliptic flow v 2 and mean-transverse momentum [p T ] of final-state particles, v 2 ∝ ε 2 and δ [p T ] / [p T ] ∝ δd ⊥ /d ⊥ , similar conclusions are also expected for the variances and skewnesses of v 2 and [p T ], i.e., a + bβ$^{2}_{2}$ for $\langle v^{2}_{2}\rangle$ and $\langle(δ[p_{T}] / [p_{T}])^{2}\rangle$ and a + [b + c cos (3γ)] β$^{3}_{2}$ for $\langle v^{2}_{2}δ [p_{T}] / [p_{T}]\rangle$ or $\langle(δ[p_{T}] / [p_{T}])^{3}\rangle$. My findings motivate a dedicated system scan of high-energy heavy ion collisions at RHIC and LHC to measure triaxiality of atomic nuclei: one first determines the coefficients b and c by collisions of isobaric near prolate nuclei, cos (3γ) ≈ 1, and near oblate nuclei, cos (3γ) ≈ –1 , with known β 2 values, followed by collisions of other species of interest with similar mass number. The (β 2 , γ) values for this species can be inferred directly from the measured variance and skewness observables from these collisions. Finally, the results demonstrate the unique opportunities offered by high-energy collisions as a tool to perform interdisciplinary nuclear physics studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Robust Optimal Experimental Design of Infinite-Dimensional Bayesian Nonlinear Inverse Problems

Abstract. We consider robust optimal experimental design (ROED) for nonlinear Bayesian inverse problems governed by partial differential equations (PDEs). An optimal design is one that maximizes some utility quantifying the quality of the solution of an inverse problem. However, the optimal design is dependent on elements of the inverse problem such as the simulation model, the prior, or the measurement error model. ROED aims to produce an optimal design that is aware of the additional uncertainties encoded in the inverse problem and remains optimal even after variations in them. We follow a worst-case scenario approach to develop a new framework for robust optimal design of nonlinear Bayesian inverse problems. The proposed framework (a) is scalable and designed for infinite-dimensional Bayesian nonlinear inverse problems constrained by PDEs; (b) develops efficient approximations of the utility, namely the expected information gain; (c) employs eigenvalue sensitivity techniques to develop analytical forms and efficient evaluation methods of the gradient of the utility with respect to the uncertainties against which we wish to be robust; and (d) employs a probabilistic optimization paradigm that properly defines and efficiently solves the resulting combinatorial max-min optimization problem. The effectiveness of the proposed approach is illustrated for optimal sensor placement problem in an inverse problem governed by an elliptic PDE.

Chowdhary, Abhijit↗
NASA NTRS1984

Consequences of a chromospheric temperature gradient on the width of H-alpha in late-type giants

An analytic expression for the integrated H alpha optical depth profile is derived for a one dimensional slab geometry model chromosphere, with electron temperature increasing as a power law with height. The formula predicts H alpha opacity and profile width to be sensitive functions of the thermal gradient. Application of the model to observation reveals that broad H alpha absorption widths in G and K giant stars are consistent with a mean H alpha chromospheric optical depth of 50, while narrower widths in M stars indicate slightly lower opacities. It is proposed that differences in H alpha width between late-type giants of similar spectral type may be due, in part, to differences in their chromospheric thermal gradient, and associated H alpha opacity.

Zarro, D. M.↗
NASA NTRS1985

Least squares inversion of Stokes profiles in the presence of velocity gradients

The Auer, Heasley and House Stokes inversion procedure in use at High Altitude Observatory is based on the analytic solution of the equation of transfer for polarized light where the representation of the thermodynamic and magnetic structure of the atmosphere is assumed to have a high degree of invariance, namely, a Milne-Eddington (ME) structure with a constant magnetic field. In the presence of invariance breaking gradients the resultant Stokes profiles are represented only approximately, if at all, by analytic forms. The accuracy of the inversion parameters and their significance as measures of actual structure are explored for the ME and the Landman-Finn sunspot models under the effects of velocity gradients. The resulting field parameters are good to a few percent and prove to be insensitive to the errors committed by the use of a ME-representation, but the resulting ME parameters yield a less precise measure of thermal structure.

Skumanich, A.↗
DOE OSTI2021

Gradient sensing via cell communication

Experimental evidence lends support to the conjecture that cell-to-cell communication plays a role in the gradient sensing of chemical species by certain chains of cells. Models have been formulated to explore this idea. For cells with no identifiable sensing structure, Mugler et al. [Proc. Natl. Acad. Sci. (U.S.A.) 113, E689 (2016)] have defined a particular local excitation, global inhibition (LEGI) model that pits nearest-neighbor communication against local reactions in a noisy environment to suggest how this sensing capability might arise in a physical system. In this study, we generalize the nearest-neighbor communication mechanism in the aforementioned LEGI model in order to explore the extent to which the gradient sensing characteristics depend on the parametrization of the communication itself, as well as on the cell size, the radius of influence of neighboring cells, and the influence of the background noise. Using our generalization and a collection of particular candidate communication models, we find that the precision of gradient sensing is indeed sensitive to the particular communication model, and we derive physical and analytic explanations for these results. The framework established and the associated results should prove useful in understanding the appropriateness of particular cell-to-cell communication models in gradient sensing studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Nonlinear sparse Bayesian learning for physics-based models

This paper addresses the issue of overfitting while calibrating unknown parameters of over-parameterized physics-based models with noisy and incomplete observations. Here, a semi-analytical Bayesian framework of nonlinear sparse Bayesian learning (NSBL) is proposed to identify sparsity among model parameters during Bayesian inversion. NSBL offers significant advantages over machine learning algorithm of sparse Bayesian learning (SBL) for physics-based models, such as 1) the likelihood function or the posterior parameter distribution is not required to be Gaussian, and 2) prior parameter knowledge is incorporated into sparse learning (i.e. not all parameters are treated as questionable). NSBL employs the concept of automatic relevance determination (ARD) to facilitate sparsity among questionable parameters through parameterized prior distributions. The analytical tractability of NSBL is enabled by employing Gaussian ARD priors and by building a Gaussian mixture-model approximation of the posterior parameter distribution that excludes the contribution of ARD priors. Subsequently, type-II maximum likelihood is executed using Newton's method whereby the evidence and its gradient and Hessian information are computed in a semi-analytical fashion. We show numerically and analytically that SBL is a special case of NSBL for linear regression models. Subsequently, a linear regression example involving multimodality in both parameter posterior pdf and model evidence is considered to demonstrate the performance of NSBL in cases where SBL is inapplicable. Next, NSBL is applied to identify sparsity among the damping coefficients of a mass-spring-damper model of a shear building frame. These numerical studies demonstrate the robustness and efficiency of NSBL in alleviating overfitting during Bayesian inversion of nonlinear physics-based models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Analytical Partial Derivatives of the Q-Law Guidance Algorithm

One method of producing low-thrust trajectories is the closed-loop Q-Law guidance algorithm. This work wraps Q-Law inside a nonlinear programming problem to search for locally optimal gain combinations and enforce nonlinear constraints on the initial state. Gradient-based trajectory optimization has been shown to benefit greatly when analytical partial derivatives are supplied to the optimizer. This paper presents the Q-Law optimization problem setup and derives the Q-Law thrust vector partial derivatives for the State Transition Matrix required by standard nonlinear optimization software.

Jackson L Shannon↗
NASA NTRS2018

Predicting Global Ground Geoelectric Field With Coupled Geospace and Three-Dimensional Geomagnetic Induction Models

We forecast the global effects of space weather on the geoelectric and geomagnetic fields using a novel combination of methods. We use a realistic three-dimensional (3-D) model of Earth's electrical conductivity and a realistic representation of magnetospheric and ionospheric current systems. Our scheme involves the following steps: (1) We run a global magnetohydrodynamic model of the magnetosphere coupled to an electrostatic model of the ionosphere. (2) We calculate a global time series of the ground magnetic field resulting from the ionospheric, field-aligned, and magnetospheric currents of the global magnetohydrodynamic model. (3) We approximate this external field by an equivalent source current flowing in a thin shell above Earth. (4) We calculate a global time series of geoelectric and geomagnetic fields from the equivalent current and a 3-D conductivity model of Earth that also takes into account the coast effect due to large horizontal conductivity gradient. We verify our implementation by comparing the results against known analytic and numeric solutions, and then apply our scheme to the geomagnetic storm of 14 and 15 December 2006. In particular, we show that accounting for 3-D structure of Earth's conductivity results in significantly enhanced geoelectric field at large lateral gradients of conductivity, especially in coastal regions, both at middle and high latitudes. In the studied geomagnetic storm the largest values of 3-D geoelectric field are detected at high latitudes reaching 2.5 volts per kilometer and the 3-D effect extends inland by a few hundred kilometers.

Honkonen, I.↗