Thermally induced electroconvection
Conducting fluid steady convection due to combined effects of thermal gradient and DC electric field predicted by analytical model
Melcher, J. R.↗
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.
An investigation of gap heating due to stepped tiles in zero pressure gradient regions of the Shuttle Orbiter Thermal Protection System
An analytical study is presented which investigates the cause of the excessive heating in the tile-to-tile gaps of the Shuttle Orbiter Thermal Protection System, as evidenced by the visible discoloration and charring of the filler bar and strain isolation pad used in the attachment of tiles to the aluminum substrate. Techniques are developed to estimate the pressure disturbances due to a stepped tile and to calculate the disturbance-induced mass flow rates in the tile-to-tile gaps, filler bar, strain isolation pad, and the tile itself. A thermal analysis in the tile-to-tile gap is then performed in order to determine the temperature response to the hot gas flow. Calculations are performed at locations on the fuselage and the left wing where damaged filler bars were observed on the first flight of the Shuttle. It is concluded that the steps and gaps must be controlled within tight tolerances during tile installation and the tolerances must be maintained in flight. If the tolerances cannot be maintained, tile-to-tile gap filler may be an alternative.
Smith, D. M.↗
Double diffusive convection during directional solidification - Nonlinear analysis
The dynamics of doubly diffusive convective flow during the constant-velocity vertical solidification of binary melts to form single-phase solids is examined, summarizing the results of recent experimental and analytical investigations. The effects of horizontal temperature gradients are considered; the nonlinearities of the system and their numerical treatment are discussed; and numerical results are presented graphically for the solidification of a Pb melt containing 0.0015 wt pct Sn, with growth velocity 0.0002 cm/s and a temperature gradient of 200 K/cm in the liquid.
Coriell, S. R.↗
Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model
In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.
Wu, Jie↗
A Robust Initialization Scheme for a Lateral Trajectory Optimization Problem with Time of Arrival Windows
We present a robust initialization scheme that estimates parameter values for the numerical solution of a two-point boundary value problem. The two-point boundary value problem formulation stems from the optimization of a cost functional subject to the dynamics of a simplified lateral aircraft model and other constraints. Leveraging regular perturbation methods, initial parameter estimates are analytically determined and used to initialize a gradient descent optimization routine which is shown to rapidly converge over a range of initial aircraft positions and heading angles. Additionally, the velocity of the aircraft is optimized to ensure the trajectory of the aircraft terminates within a desired region in both time and space.
lateral trajectory optimization↗
A self-consistent two-dimensional resistive fluid theory of field-aligned potential structures including charge separation and magnetic and velocity shear
A self-consistent two-fluid theory that includes the magnetic field and shear patterns is developed to model stationary electrostatic structures with field-aligned potential drops. Shear flow is also included in the theory since this seems to be a prominent feature of the structures of interest. In addition, Ohmic dissipation, a Hall term, and pressure gradients in a generalized Ohm's law, modified for cases without quasi-neutrality, are included. In the analytic theory, the electrostatic force is balanced by field-aligned pressure gradients (i.e., thermal effects in the direction of the magnetic field) and by pressure gradients and magnetic stresses in the perpendicular direction. Within this theory, simple examples of applications are presented to demonstrate the kind of solutions resulting from the model. The results show how the effects of charge separation and shear in the magnetic field and the velocity can be combined to form self-consistent structures such as are found to exist above the aurora, suggested also in association with solar flares.
Hesse, Michael↗
The Hadley and Rossby regimes in a spherical atmosphere
The properties of the steady Hadley and Rossby regimes for a thermally forced rotating fluid on a sphere are studied. The two layer modified geostrophic model is employed which allows for thermal advection by the divergent wind and time dependent static stability. Heating processes are parameterized using the Newtonian approximation and Rayleigh friction is accounted for. The equations are transformed to spectral form using spherical harmonics and then truncated retaining a simple axisymmetric state and initial, one wave. A time independent Hadley circulation is obtained which is neutral to axisymmetric disturbances but unstable to wave like perturbations for intermediate values of the meridional temperature gradient, indicating the existence of both an upper and lower symmetric Hadley regime. An analytical solution for the steady Rossby circulation is determined for values of the meridional temperature gradient where the Hadley regime is unstable. Linear perturbation theory is used to show that within the steady Rossby regime two or more waves cannot exist simultaneously.
Feldstein, S. B.↗
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.↗
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↗
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.↗
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.↗
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.↗
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.↗
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.↗
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↗
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.↗
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.↗
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↗