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.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1971

Comparison of a discrete steepest ascent method with the continuous steepest ascent method for optimal programing

A discrete steepest ascent method which allows controls which are not piecewise constant (for example, it allows all continuous piecewise linear controls) was derived for the solution of optimal programming problems. This method is based on the continuous steepest ascent method of Bryson and Denham and new concepts introduced by Kelley and Denham in their development of compatible adjoints for taking into account the effects of numerical integration. The method is a generalization of the algorithm suggested by Canon, Cullum, and Polak with the details of the gradient computation given. The discrete method was compared with the continuous method for an aerodynamics problem for which an analytic solution is given by Pontryagin's maximum principle, and numerical results are presented. The discrete method converges more rapidly than the continuous method at first, but then for some undetermined reason, loses its exponential convergence rate. A comparsion was also made for the algorithm of Canon, Cullum, and Polak using piecewise constant controls. This algorithm is very competitive with the continuous algorithm.

Childs, A. G.↗

NASA NTRS1988

Digital active control law synthesis for aeroservoelastic systems

This paper presents a formulation for synthesis of digital active control laws for aeroservoelastic systems, which are typically modeled by large order equations in order to accurately represent the rigid and flexible body modes, unsteady aerodynamic forces, actuator dynamics, and gust spectra. The control law is expected to satisfy multiple design requirements on the dynamic loads, responses, actuator deflection and rate limitations, as well as maintain certain stability margins, yet should be simple enough to be implemented by an onboard digital microprocessor. The synthesis procedure minimizes a linear quadratic Gaussian type cost function, by updating selected free parameters of the control law, while satisfying a set of inequality constraints on the design loads, responses and stability margins. A stable classical control law or an estimator based full or reduced order control law can be used as an initial design starting point. The gradients of the cost function and the constraints, with respect to the digital control law design variables are derived analytically, to facilitate rapid convergence. Selected design responses can be treated as constraints instead of lumping them into the cost function, in order to satisfy individual root-mean-square load and response limitations. Constraints are also imposed on the minimum singular value requirements for stability robustness improvement.

Mukhopadhyay, Vivekananda↗

NASA NTRS1990

Viscosity and mass transport in nonuniform Keplerian disks

A quantitative formalism for Keplerian particulate disk dynamics having its basis in a heuristic description of viscous transport is presently applied to the case of a ring with optical thickness radial gradient. A steady-state velocity distribution solution directly yields both the radial mass transport and the viscosity. The analytical method employed involves the solution of a novel form of the Krook equation through the separation of the phase-space collisional-products distribution into a symmetrical component and a delta function-approachable remainder. While this model is simplified for the case treated, its general approach may be extended to less artificially restricted cases.

Ojakangas, G. W.↗

NASA NTRS1998

Study of Boundary Layer Development in a Two-Stage Low-Pressure Turbine

Experimental data from jet-engine tests have indicated that unsteady blade row (wake) interactions and separation can have a significant impact on the efficiency of turbine stages. The effects of these interactions can be intensified in low-pressure turbine stages because of the low Reynolds number operating environment. Measured turbine efficiencies at takeoff can be as much as two points higher than those at cruise conditions. Thus, during the last decade a significant amount of effort has been put into determining the effects of transition and turbulence on the performance of low pressure turbine stages. Experimental investigations have been performed, for example, by Hodson et al. and Halstead et al. These investigations have helped identify/clarify the roles that factors such as the Reynolds number, free stream turbulence intensity, pressure gradient and curvature have in the generation of losses. In parallel to the experimental investigations, there have been significant analytical efforts to improve the modeling of transition. Examples of such efforts include the works of Mayle and Gostelow et al. These newer models show promise of providing accurate transition predictions over a wide range of flow conditions, although they have yet to be implemented into the numerical flow analyses used by the turbine design community. Some recent computational investigations of interest include the works of Chernobrovkin and Lakshminarayana and Eulitz and Engel. The focus of the current effort has been to -use a viscous, unsteady quasi-three-dimensional Navier-Stokes analysis to study boundary layer development in a two-stage low-pressure turbine. A two-layer algebraic turbulence model, along with a natural transition model and a bubble transition model, have been used, The geometry used in the simulations has been the subject of extensive experiments. The predicted results have been compared with experimental data, including airfoil loadings and time-averaged/unsteady integral boundary layer quantities.

Dorney, Daniel J.↗

NASA NTRS2025

Photophoretic Propulsion Enabling Mesosphere Exploration NIAC Phase I Final Report

This Phase I report presents a comprehensive study on photophoretic flyers—innovative, ultralight, solar-powered vehicles that harness photophoretic forces generated via Knudsen pumping to achieve sustained flight in the mesosphere (50–80 km altitude). By integrating advanced materials such as nanocardboard— characterized by its extremely low areal density (~1 g/m²) and high bending stiffness—with ultrathin light-absorbing coatings, the project converts incident solar radiation directly into a directed thrust. Extensive experimental investigations, coupled with high-fidelity computational fluid dynamics (CFD) simulations using ANSYS Fluent, validate the concept across various three-dimensional geometries, including spherical, conical, and rocket-shaped configurations. These simulations bridge the gap between free-molecular and continuum flow regimes, demonstrating that optimized designs can generate lift forces sufficient to support kilogram-scale payloads even in low-pressure environments. At the heart of this innovation is the use of Knudsen pumping, where temperature gradients across porous surfaces induce directional gas flow, creating a modest overpressure that provides lift. The report introduces an analytical framework that interpolates between the well-known low-Reynolds number drag regime and the high-Reynolds number momentum theory. This model accurately predicts lift based on design parameters such as microchannel dimensions, porous wall geometry, areal density, and nozzle exit area. For instance, simulations indicate that 10-meter-scale structures with carefully engineered porous walls can achieve the necessary pressure differential to support scientifically significant payloads (~1 kg). The study also explores a hybrid propulsion approach that combines solar buoyancy with photophoretic lift. Initially, solar heating creates a buoyant force that elevates the flyer to mesospheric altitudes. Once in the optimal pressure range, the photophoretic mechanism—powered by Knudsen pumping—takes over as the primary source of lift, ensuring stable, long-duration flight. This dual-mode operation not only facilitates the deployment of photophoretic flyers but also broadens the potential applications for mesospheric exploration. In addition to propulsion, the report investigates the integration of photophoretic thrusters for trajectory control of existing research balloons in the upper stratosphere. By dynamically adjusting the nozzle orientation and controlling the flow-through velocity, these thrusters provide precise maneuverability, enabling the flyers to counteract atmospheric disturbances and adjust their flight paths in real time. For example, a photophoretic thruster approximately 7.5 by 7.5 meters in size could be unfolded below a payload gondola of a 60 million-cubic-foot zero-pressure balloon. Such a thruster can provide horizontal speed control of approximately 1 m/s using only sunlight and no moving parts (except those needed to track the Sun and control the jet direction). Importantly, photophoretic thrusters operate more efficiently at higher altitudes, which is complementary to known trajectory control techniques, such as propellers and tethered wings, which are more effective at lower altitudes. Finally, the report identifies three scientific research thrusts where mesospheric aircraft technology can have a profound impact: atmospheric tides, characterization of gravity waves, and investigation of mesospheric instabilities. Overall, the findings of this Phase I project represent a significant advancement in photophoretic propulsion technology. By demonstrating that large-scale, ultralight structures can be powered solely by solar radiation—via carefully engineered Knudsen pumping—this work lays a robust foundation for scalable, near-space flight architectures. Future refinements in material fabrication, structural optimization, and integrated trajectory control are expected to further enhance performance, paving the way for operational demonstrations that could revolutionize atmospheric science, remote sensing, and communication networks.

Knudsen Pump↗

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 NTRS2022

Low-Earth Orbit Trajectory Optimization in the Presence of Atmospheric Uncertainty

The previous 20 to 25 years have seen a tremendous increase in space exploration, and with that an increase in the level of logistics planning needed to ensure mission success. For spacecraft that are designed to be periodically re-supplied, a key logistics consumable is propellant, as it constitutes the greatest up-mass on re-supply vehicles. A trajectory design strategy is therefore desired that minimizes propellant usage in order to ease the demand for propellant re-supply missions. This thesis develops such a strategy in three stages, and uses the International Space Station (ISS) as its testbed, as no other LEO spacecraft is more challenging from a space logistics standpoint. First, the ISS trajectory planning problem is formulated as a constrained burn optimization problem assuming a deterministic atmosphere. The cost function is total ∆v, with constraints imposed on longitude of ascending viii node (LAN) and semi-major axis (SMA) altitude. Analytic derivatives are constructed for both the cost and constraints, which are necessary given the 6-week to 2-year time frames being considered. A gradient-based optimizer is then utilized to find locally-optimal solutions to real-world ISS trajectory planning problems. Second, atmospheric uncertainty is addressed by constructing a probabilistic model of space weather data using Gaussian Processes (GPs). Bayesian inference is performed using the GP model to generate mean and covariance estimates for space weather predictions, whose pedigree is assessed against test data. The predictions are then mapped into atmospheric density via the analytic Jacchia-Roberts density model, and the effect of space weather uncertainty on orbital lifetime is examined. Third, an ISS burn execution uncertainty model is developed. This model, along with the space weather uncertainty model, are deployed in a linear covariance analysis to ascertain their combined effect on LAN and SMA altitude dispersions. The deterministic constraints from the original problem are re-formulated as stochastic constraints, where now the constraint uncertainty interval is required to fall within specified bounds. An updated optimization framework is constructed using the original ∆v cost function along with the stochastic constraints to solve the trajectory optimization problem under atmospheric uncertainty. Finally, the complete architecture is summarized for deployment in an operational setting.

Trajectory Optimization↗

DOE OSTI2022

Focusing ions at atmospheric pressure using nonlinear DC voltage sequences applied to a stacked ring ion guide

Many modern ion mobility (IM) and mass spectrometers (MS) operate under low pressure (= 10 Torr) and employ high voltage radiofrequencies (RF) to provide ion confinement. Unfortunately, RF effectiveness drastically decreases as pressure increases, and few techniques for focusing ions at elevated pressures exist. Here we demonstrate a new approach for focusing ions at atmospheric pressure (AP) by applying nonlinear DC voltage sequences following quadratic and power (exponential) functions to a stacked ring ion guide. In this work, we used ion trajectory simulations to rigorously explore how ions react to nonlinear electric fields and validate the simulations with a set of ion current measurements performed at AP. Ion trajectory simulations show that ions initially defocus near the entrance of the device but then become intensely focused as they travel through the device. Contour plots for both nonlinear voltage sequences show electric field lines that increasingly curve inwards as a function of distance, resulting in spatial ion focusing. Experimental ion current and spot size measurements were performed at AP using a 10-cm stacked ring ion guide and a segmented Faraday cup detector. Quadratic sequences produced ~5% smaller spot sizes (~22.8 mm) and ~25% higher ion current compared to a linear voltage sequence (~24.0 mm). Alternatively, power sequences produced ~64% smaller spot sizes (~8.7 mm), albeit with ~10x lower ion current. However, both nonlinear voltage sequences produced similar ion currents at the center of the Faraday cup detector, indicating that higher ion densities are achieved when using nonlinear voltage gradients. These results demonstrate a new way to focus ions at AP, and the capabilities demonstrated here provide fundamental insights on how to keep ions inside analytical devices at elevated pressures without RF.

47 OTHER INSTRUMENTATION↗

NASA NTRS2004

Analytical Solution for the Critical Velocity of Pushing/Engulfment Transition

The distribution of ceramic particles in a metal matrix composite material depends primarily on the interaction of the particles with the solid/liquid interface during the solidification process. A numerical model that describes the evolution of the shape of the solid/liquid interface in the proximity of a foreign particle will presented in this paper. The model accounts for the influence of the temperature gradient and the Gibbs-Thomson and disjoining pressure effects. It shows that for the systems characterized by k(sub p) < k(sub L) the disjoining pressure causes the interface curvature to change its sign in the close-contact particle/interface region. It also shows that the increase of the temperature gradient diminishes the effect of the disjoining pressure. The analysis of the numerical results obtained for a large range of processing conditions and materials parameters has led to the development of an analytical solution for the critical velocity of pushing/engulfinent transition. The theoretical results will be discussed and compared with the experimental measurements performed under microgravity conditions.

Catalina, Adrian V.↗

DOE OSTI2024

Influence of Concentration Gradients on Electroconvection at a Cation-Exchange Membrane Surface

Membrane-based systems, such as electrodialysis, play an important role in desalination and industrial separation processes. Electrodialysis uses alternating anion- and cation-exchange membranes with a perpendicular electric field to generate concentrated and diluate streams from a feed solution. It is known that under overlimiting current conditions, reduced charge and mass transfer at the membrane interface leads to regions of high ion depletion generating instability and vortices termed electroconvection. While electroconvective mixing is known to directly impact the separation efficiency of electrodialysis, the influence of ion concentration gradients across the membrane experienced in a functional electrodialysis system is not known. Here, we report the influence of ion concentration gradients across a cation exchange membrane (Nafion) that is both aligned with and opposed to the applied electric field. Experiments were conducted by coflowing NaCl solutions of different concentrations (0.1–100 mM) on each side of the membrane, and electroconvection was visualized with a fluorescence dye (Rhodamine 6G). Here, we obtained concentration profiles from fluorescence image data and systematically measured the thickness of the depletion boundary layer d BL under different conditions. We found smaller d BL values at a higher flow rate both with and without concentration gradients. Our results show that electroconvection is enhanced when the electric field is opposite to the direction of the concentration gradient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

The effect of concentrated PV gradients on stationary waves

The author reexamines the Charney-Drazin problem with special attention to the concentration of potential vorticity gradient in the neighborhood of the tropopause. It is found that the degree of concentration has a profound effect on the response to stationary forcing, with greater concentration leading to greater response. Smoothing the concentration either analytically or numerically (by using coarser resolution) both lead to reduced responses, especially at higher wavenumbers. The results suggest a potentially important interaction between baroclinically unstable eddies and stationary waves. Insofar as the former act to mix potential vorticity in the troposphere while concentrating gradients at tropopause levels, they significantly condition the basic state for the latter.

Lindzen, Richard S.↗

DOE OSTI2021

Stochastic Gradient-Based Distributed Bayesian Estimation in Cooperative Sensor Networks

Distributed Bayesian inference provides a full quantification of uncertainty offering numerous advantages over point estimates that autonomous sensor networks are able to exploit. However, fully-decentralized Bayesian inference often requires large communication overheads and low network latency, resources that are not typically available in practical applications. In this paper, we propose a decentralized Bayesian inference approach based on stochastic gradient Langevin dynamics, which produces full posterior distributions at each of the nodes with significantly lower communication overhead. We provide analytical results on convergence of the proposed distributed algorithm to the centralized posterior, under typical network constraints. Finally, we also provide extensive simulation results to demonstrate the validity of the proposed approach.

42 ENGINEERING↗

NASA NTRS1972

Errors incurred in profile reconstruction and methods for increasing inversion accuracies for occultation type measurements

A method for augmenting the occultation experiment is described for slightly refractive media. This method which permits separation of the components of the gradient of refractivity, appears applicable to most of the planets for a major portion of their atmospheres and ionospheres. The analytic theory is given, and the results of numerical tests with a radially and angularly varying model of an ionosphere are discussed.

Gross, S. H.↗

NASA NTRS1973

Minimum impulse three-body trajectories.

A rapid and accurate method of calculating optimal impulsive transfers in the restricted problem of three bodies has been developed. The technique combines a multi-conic method of trajectory integration with primer vector theory and an accelerated gradient method of trajectory optimization. A unique feature is that the state transition matrix and the primer vector are found analytical without additional integrations or differentiations. The method has been applied to the determination of optimal two and three impulse transfers between the L2 libration point and circular orbits about both the earth and the moon.

D'Amario, L.↗

NASA NTRS1979

A prediction method for velocity and temperature profiles in a two-dimensional nominally steady turbulent boundary layer

This paper describes a recently developed boundary-layer prediction method for a variable property compressible flow, in which heat transfer takes place primarily by forced convection and for which the mainstream Mach number is small. The leading order terms, in asymptotic expansions for large Reynolds numbers, are obtained for the mean velocity and temperature distribution in both the inner and outer layer of the turbulent boundary layer. Closure in the inner layer is achieved using an analytical model for the mean profiles which is based on the observed coherent structure of the time-dependent inner layer flow. For the outer layer, simple eddy viscosity and conductivity models are developed without recourse to the Reynolds analogy. In the prediction method a numerical solution of the outer layer equations is matched to the analytical inner layer profiles as the computation procedes downstream. Calculations are presented for a range of adverse and favorable pressure gradient flows and the predicted results compare well with existing data.

Weigand, G. G.↗

NASA NTRS1985

Theoretical interpretation of the observed interplanetary magnetic field radial variation in the outer solar system

Observations of the azimuthal component of the IMF are evaluated through the use of an MHD model which shows the effect of magnetic flux tubes opening in the outer solar system. It is demonstrated that the inferred meridional transport of magnetic flux is consistent with predictions by the MHD model. The computed azimuthal and radial magnetic flux deficits are almost identical to the observations. It is suggested that the simplest interpretation of the observations is that meridional flows are created by a direct body force on the plasma. This is consistent with the analytic model of Nerney and Suess (1975), in which such flux deficits in the IMF arise naturally from the meridional gradient in the spiralling field.

Suess, S. T.↗

NASA NTRS1990

The baroeffect and an appropriate momentum boundary condition

As disparate molecular weight gases isothermally diffuse between two ends of a capillary tube, they can support a pressure gradient. The magnitude of this pressure gradient depends critically on viscous wall stress and becomes a measure of the boundary condition. Here, a one-dimensional analytical model is proposed that allows for specular gas reflection from the wall. It explains anomalous (4/3) correction factors required previously to match experiment to baroeffect models. It predicts a new physical phenomenon, a surface-driven baroeffect for equal molecular weight gases. Diffusive slip contributions exceed the order of convective diffusion for Peclet number Pe less than 1 and approximately equal convective diffusion for Pe values between 1 and 4. For binary gases, this model further extends the baroeffect experiment as a means to find momentum accommodation coefficients, slip lengths, and external friction coefficients. A corrected binary diffusion coefficient is defined.

Noever, David A.↗

NASA NTRS1993

Magnetohydrodynamic wave propagation in one-dimensional nonhomogeneous, self-gravitating clouds

We study the propagation of magnetohydrodynamic (MHD) waves through nonhomogeneous, selfgravitating, magnetic media representative of molecular cloud environments and focus on the issues of cloud support and line profiles. Since the general treatment of this topic is burdened by a complex mathematical formalism, we consider simplifying geometries which yield analytical solutions in the linear wave limit. In particular, we study both magnetoacoustic and Alfven wave propagation along the density gradient in a one-imensional slab. Of specific relevance to molecular clouds, we find that the back reaction of the Alfven waves can provide a pressure along the direction of the magnetic field lines; this pressure can help support a density enhancement against gravitational collapse. Furthermore, we find that the velocity amplitudes of these waves increase as the density decreases, in rough agreement with observational estimates of line width versus density relations.

Fatuzzo, Marco↗