Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JuMP model”

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 235 records · Page 13

Statistical detection of the mid-Pleistocene transition

Statistical methods have been used to show quantitatively that the transition in mean and variance observed in delta O-18 records during the middle of the Pleistocene was abrupt. By applying these methods to all of the available records spanning the entire Pleistocene, it appears that this jump was global and primarily represents an increase in ice mass. At roughly the same time an abrupt decrease in sea surface temperature also occurred, indicative of sudden global cooling. This kind of evidence suggests a possible bifurcation of the climate system that must be accounted for in a complete explanation of the ice ages. Theoretical models including internal dynamics are capable of exhibiting this kind of rapid transition.

Maasch, K. A.↗

The relation between ferroelasticity and superconductivity

The high-temperature superconductivity is explained widely by the layered crystal structure. The one- and two-dimensional subsystems and their interaction are investigated here. It is assumed that the high-T(sub c) superconductivity takes place in the two-dimensional subsystem and the increase of the phase transition temperature from 60 K up to 90 K is the consequence of turning on the influence of one-dimensional chains. The interaction between the two subsystems is transferred along the c axis by the phonons of breathing mode, which causes the hybridization of the electronic bonds between these subsystems. The experimental works indicate that the existence of both the chains Cu(1)-O and their interaction with the superconducting plane of Cu(2)-O modify the temperature of the transition to the superconducting state. It is seen from the neutron scattering data that the rates of the interatomic distance dependencies on temperature are changed around 240 K and 90 K. The 'zig-zag' order in Cu(1)-O chains has been postulated but, on the other hand, the vibrations with a large amplitude only were reported. The bi-stabilized situation of the oxygen ions can be caused by the change of distance between these ions and the Ba ions. It leads to the appearance of a two-well potential. Its parameters depend on temperature and the dynamics of the oxygen ions' movement. They can induce the antipolar order, which can be, however, more or less chaotic. The investigation of the ferroelastic properties of Y-Ba-Cu-O samples lead to the conclusion that they are related to jumps of ions inside the given chain and not to a diffusion between different sites in the ab plane. Researchers deduce, thus, that the fluctuating oxygen ions from these chains create dipoles in the ab plane. They can be described with the pseudo-spin formalism (- Pauli matrices). The system can be described with the Ising model. The pseudo-spins interact with phonons and influence the superconductivity in the second subsystem.

Molak, A.↗

Spherically symmetric, expanding, non-LTE model atmospheres for novae during their early stages

In the continuum and line-blanketed models presented here, nova atmospheres are characterized by a very slow decrease of density with increasing radius. This feature leads to very large geometrical extensions so that there are large temperature differences between the inner and outer parts of the line-forming regions. The theoretical spectra show a large IR excess and a small Balmer jump which may be either in absorption or in emission. For the parameters considered (effective temperature of about 10 exp 4 K, L = 2 x 10 exp 4 solar luminosities, outer boundary density of about 3 x 10 exp -15 g cm exp -3, mass-loss rate of 10 exp -5 solar masses/yr), most lines are in absorption. The effects of changes in the abundances of the heavy elements on the emergent spectra are discussed. The strong unidentified features observed in ultraviolet spectra of novae are found in actuality to be regions of transparency within the Fe 'forest'. Ultraviolet spectra obtained from the IUE archives are displayed, and spectral synthesis of these spectra is done using the theoretical atmospheres.

Hauschildt, P. H.↗

Supercomputer optimizations for stochastic optimal control applications

Supercomputer optimizations for a computational method of solving stochastic, multibody, dynamic programming problems are presented. The computational method is valid for a general class of optimal control problems that are nonlinear, multibody dynamical systems, perturbed by general Markov noise in continuous time, i.e., nonsmooth Gaussian as well as jump Poisson random white noise. Optimization techniques for vector multiprocessors or vectorizing supercomputers include advanced data structures, loop restructuring, loop collapsing, blocking, and compiler directives. These advanced computing techniques and superconducting hardware help alleviate Bellman's curse of dimensionality in dynamic programming computations, by permitting the solution of large multibody problems. Possible applications include lumped flight dynamics models for uncertain environments, such as large scale and background random aerospace fluctuations.

Chung, Siu-Leung↗

A numerical method for acoustic oscillations in tubes

A numerical method to obtain the neutral curve for the onset of acoustic oscillations in a helium-filled tube is described. Such oscillations can cause a serious heat loss in the plumbing associated with liquid helium dewars. The problem is modelled by a second-order, ordinary differential eigenvalue problem for the pressure perturbation. The numerical method to find the eigenvalues and track the resulting points along the neutral curve is tailored to this problem. The results show that a tube with a uniform temperature gradient along it is much more stable than one where the temperature suddenly jumps from the cold to the hot value in the middle of the tube.

Gary, John M.↗

Theoretical analysis of a rotating two-phase detonation in liquid rocket motors.

A nonlinear analysis to study tangential mode shock instabilities in a thin annular chamber is carried out by employing a one dimensional two phase detonation wave as a reaction model. It is assumed that phase change and reaction take place only within the wave, which is treated as a discontinuity. The annulus is unrolled and the flow is considered as two dimensional with the coordinate system fixed on the wave front. Between waves, the flow is assumed to be isentropic with no interaction between droplets and burned gases. Jump conditions across the wave are solved for two phase flow. The average pressure along the injection plate is related to the design chamber pressure by the use of overall conservation equations. The wave strength is written in terms of the design parameters of the chamber. The results compare favorably with existing experiments. Finally, the effects of drop size on the wave strength are discussed and a simple criterion which sets the lower limit of validity for this strong wave analysis, is presented.

Shen, P. I.-W.↗

Spline-Based Smoothing of Airfoil Curvatures

Constrained fitting for airfoil curvature smoothing (CFACS) is a splinebased method of interpolating airfoil surface coordinates (and, concomitantly, airfoil thicknesses) between specified discrete design points so as to obtain smoothing of surface-curvature profiles in addition to basic smoothing of surfaces. CFACS was developed in recognition of the fact that the performance of a transonic airfoil is directly related to both the curvature profile and the smoothness of the airfoil surface. Older methods of interpolation of airfoil surfaces involve various compromises between smoothing of surfaces and exact fitting of surfaces to specified discrete design points. While some of the older methods take curvature profiles into account, they nevertheless sometimes yield unfavorable results, including curvature oscillations near end points and substantial deviations from desired leading-edge shapes. In CFACS as in most of the older methods, one seeks a compromise between smoothing and exact fitting. Unlike in the older methods, the airfoil surface is modified as little as possible from its original specified form and, instead, is smoothed in such a way that the curvature profile becomes a smooth fit of the curvature profile of the original airfoil specification. CFACS involves a combination of rigorous mathematical modeling and knowledge-based heuristics. Rigorous mathematical formulation provides assurance of removal of undesirable curvature oscillations with minimum modification of the airfoil geometry. Knowledge-based heuristics bridge the gap between theory and designers best practices. In CFACS, one of the measures of the deviation of an airfoil surface from smoothness is the sum of squares of the jumps in the third derivatives of a cubicspline interpolation of the airfoil data. This measure is incorporated into a formulation for minimizing an overall deviation- from-smoothness measure of the airfoil data within a specified fitting error tolerance. CFACS has been extensively tested on a number of supercritical airfoil data sets generated by inverse design and optimization computer programs. All of the smoothing results show that CFACS is able to generate unbiased smooth fits of curvature profiles, trading small modifications of geometry for increasing curvature smoothness by eliminating curvature oscillations and bumps (see figure).

Li, W.↗

Turbulence in electrostatic ion acoustic shocks

Three types of collisionless electrostatic ion acoustic shocks are investigated using a double plasma (DP) device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distrubution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.

Means, R. W.↗

Turbulence in electrostatic ion-acoustic shocks

Three types of collisionless electrostatic ion-acoustic shocks are investigated using the University of California, Los Angeles, double plasma device: (1) laminar shocks; (2) small amplitude turbulent shocks in which the turbulence is confined to be upstream of the shock potential jump; and (3) large amplitude turbulent shocks in which the wave turbulence occurs throughout the shock transition. The wave turbulence is generated by ions which are reflected from the shock potential; linear theory spatial growth increments agree with experimental values. The experimental relationship between the shock Mach number and the shock potential is shown to be inconsistent with theoretical shock models which assume that the electrons are isothermal. Theoretical calculations which assume a trapped electron equation of a state and a turbulently flattened velocity distribution function for the reflected ions yields a Mach number vs potential relationship in agreement with experiment.

Means, R. W.↗

Infrared emission and tidal interactions of spiral galaxies

Computer simulations of tidal interactions of spiral galaxies are used to attempt to understand recent discoveries about infrared (IR) emitting galaxies. It is found that the stronger tidal perturbation by a companion the more disk gas clouds are thrown into nucleus crossing orbits and the greater the velocity jumps crossing spiral arms. Both these tidally created characteristics would create more IR emission by high speed cloud collisions and more IR via effects of recently formed stars. This expectation at greater tidal perturbation matches the observation of greater IR emission for spiral galaxies with closer and/or more massive companions. The greater collision velocities found at stronger perturbations on the models will also result in higher dust temperature in the colliding clouds. In the IR pairs examined, most have only one member, the larger, detected and when both are detected, the larger is always the more luminous. In simulations and in a simple analytic description of the strong distance dependence of the tidal force, it is found that the big galaxy of a pair is more strongly affected than the small.

Byrd, Gene G.↗

Evaporation, ignition and combustion of nondilute clusters of drops

A theory of evaporation, ignition, and burning of moderately dense spherical drop clusters has been developed. The theory takes into account burning of premixed air and fuel internal to the cluster at ignition and subsequent burning of fuel emitted from the cluster by a flame sheet surrounding it. The model considers interdrop interaction, momentum exchange between drops and gas, and turbulent exchange processes between the cluster and its surroundings. Calcualtions are performed for varying initial air to fuel mass ratios, initial cluster radii, ambient gas temperatures, and initial drop temperatures. Results are presented for fuel burn fractions at ignition and at the moment of drop disappearance, as well as jump conditions at ignition.

Bellan, J.↗

Unusual Bow Shock and Magnetopause Locations During May 11-12, 1999

On May 11 and 12, 1999, when the solar wind proton density dropped to only about 3% of its nominal value, WIND and IMP 8 (and possibly Lunar Prospector) encountered the Earth's bow shock at unusually large distances. These observations, at 57-63 Re, correspond to the most distant crossings of the bow shock ever recorded. This paper will present observations from Interball-1 and Geotail, showing the bow shock near its nominal position, and from IMP 8 and WIND (and perhaps Lunar Prospector) crossing the bow shock at increasing radial distances. The bow shock is extremely weak by the time of its observation by WIND. However, a clear jump in both the magnetic field strength and solar wind bulk velocity is quite apparent, although with no significant changes in solar wind proton density. The nature of this dissipating shock will be discussed. Also, interball-1 has observed the magnetopause at an unusually large standoff distance. The observed motion of the magnetopause and bow shock will be compared with model predictions. This set of observations provide an excellent test case for the behavior of gas dynamic and MHD models under extreme conditions.

Szabo, Adam↗

The Interaction of the Solar Wind with Solar Probe Plus - 3D Hybrid Simulation: The Study for the Distance 4.5Rs - Report 1

Our report devotes a 3D numerical hybrid model of the interaction of the solar wind with the Solar Probe spacecraft. The SPP model includes 3 main parts, namely, a non-conducting heat shield, a support system, and cylindrical section or spacecraft bus that contains the particle analysis devices and antenna. One observes an excitation of the low frequency Alfven and whistler type wave directed by the magnetic field with an amplitude of about (0.06-0.6) V/m. The compression waves and the jumps in an electric field with an amplitude of about (0.15-0.7) V/m were also observed. The wave amplitudes are comparable to or greater than previously estimated max wave amplitudes that SPP is expected to measure. The results of our hybrid simulation will be useful for understanding the plasma environment near the SPP spacecraft at the distance 4.5 Rs. Future simulation will take into account the charging of the spacecraft, the charge separation effects, an outgassing from heat shield, a photoionization and an electron impact ionization effects near the spacecraft.

Lipatov, Alexander S.↗

The Interaction of the Solar Wind with Solar Probe Plus - 3D Hybrid Simulation: The Study for the Distance 4.5Rs - Report 1

Our report devotes a 3D numerical hybrid model of the interaction of the solar wind with the Solar Probe spacecraft. The Solar Probe Plus (SPP) model includes 3 main parts, namely, a non-conducting heat shield, a support system, and cylindrical section or spacecraft bus that contains the particle analysis devices and antenna. One observes an excitation of the low frequency Alfven and whistler type wave directed by the magnetic field with an amplitude of about (0.06-0.6) V/m. The compression waves and the jumps in an electric field with an amplitude of about (0.15-0.7) V/m were also observed. The wave amplitudes are comparable to or greater than previously estimated max wave amplitudes that SPP is expected to measure. The results of our hybrid simulation will be useful for understanding the plasma environment near the SPP spacecraft at the distance 4.5 Rs. Future simulation will take into account the charging of the spacecraft, the charge separation effects, an outgassing from heat shield, a photoionization and an electron impact ionization effects near the spacecraft.

Lipatov, Alexander S.↗

Advances in Integrally Formed Metallic Structures

NASA Langley Research Center (LaRC) is applying flow forming technology to aerospace aluminum alloys for the manufacture of integrally stiffened cylinders that meet aircraft fuselage or launch vehicle specifications. High-strength 2xxx series alloys, in various temper/product form conditions, are being compared with medium-strength, very formable AA 6061 as the benchmark. Screening of candidate materials traditionally relies on empirical forming trials, but more recently on metrics derived from constant nominal strain rate tensile testing. In this work, the established hierarchy of formability metrics is adopted as the baseline; tensile area reduction, work hardening exponent, and modulus of resilience. The suite of metrics can be expanded based on phenomenological and mechanistic considerations of plastic flow instabilities to incorporate anisotropy ratio, rate of work hardening, and dynamic strain aging effects. Separate research thrusts at LaRC are aimed at formulating a protocol for scalable production of fully formed, defect-free components: 1) quantifying the extended set of formability metrics from constant and jump strain rate tensile tests with digital image correlation (DIC); 2) using statistically designed experiments on a vertical spin/flow forming facility (WF VUD-600®) to isolate the effect of convenient process variables, such as axial feed rate and thickness reduction; and 3) predicting the effect of feed ratio (roller translation relative to mandrel rotation speeds) on material response via process modeling using the DEFORM® simulation software. Current results will be discussed in terms of correlating intrinsic material properties with the bulk workability characteristics to successfully fabricate cylinders with complex cross-sections.

Aluminum↗

Multivalued velocity field model of turbulence

The considered study is based on new theoretical concepts regarding a post-instability model of a fluid discussed by Zak (1980). The model permits the completion of the governing equations of turbulence by introducing multivalued fields of velocities. Attention is given to the mechanism of energy dissipation, the characteristic wave propagation, a simplified model, the formation of turbulence around stagnation points, the formulation of boundary conditions, and the mechanism of turbulence formation. The mechanism of turbulence formation can be understood as propagation of initial discontinuities from the boundaries into a flow with the characteristic velocity which is defined by the normal (to the boundary) velocity components. These components emerge at the boundary as a result of jumps in the tangential components due to the continuity equation.

Zak, M.↗

Application of dynamical systems theory to the high angle of attack dynamics of the F-14

Dynamical systems theory has been used to study the nonlinear dynamics of the F-14. An eight degree of freedom model that does not include the control system present in operational F-14s has been analyzed. The aerodynamic model, supplied by NASA, includes nonlinearities as functions of the angles of attack and sideslip, the rotation rate, and the elevator deflection. A continuation method has been used to calculate the steady states of the F-14 as continuous functions of the control surface deflections. Bifurcations of these steady states have been used to predict the onset of wing rock, spiral divergence, and jump phenomena which cause the aircraft to enter a spin. A simple feedback control system was designed to eliminate the wing rock and spiral divergence instabilities. The predictions were verified with numerical simulations.

Jahnke, Craig C.↗

Boundary-value problem of configurations with compressible free vortex flow

A self-consistent version of the compressible boundary-value problem for configurations with leading-edge vortex separation is formulated, based on the assumption that the compressible flow field is controlled by the linearized potential equation. The stream surface boundary condition and the zero pressure jump condition of the compressible free vortex flows are analyzed; application of the Goethert rule permits the compressible nonlinear boundary-value problem for the subsonic flow domain to be transformed into an equivalent nonlinear incompressible problem. The compressibility corrections developed are used in numerical calculations of subsonic leading-edge vortex flows about planar wing geometries. The sample calculations, employing an inviscid flow model in which the wing and vortex sheets are represented by piecewise continuous quadratic doublet sheet distributions, are applicable to high subsonic Mach numbers.

Brune, G. W.↗