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 127 records · Page 7

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings.

Vijayaraghavan, D.↗

Formation of slow shock pairs associated with coronal mass ejections

The formation of a forward-reverse slow shock pair in the solar corona is presently simulated by an MHD model that uses the Rankine-Hugoniot solution to calculate the flow-property jumps at all shock crossings. The shocks divide the solution-domain into several continuous flow regions whose respective governing characteristics are solved by the method of characteristics. The plasma impact compresses the plasma near the front of the coronal mass ejection (CME); as the CME-associated slow shock pair moves outwards in interplanetary space, it evolves into a pair of fast shocks. All three phenomena are eventually manifested in interplanetary space as a magnetic cloud accompanied by a fast shock pair, with a forward shock preceding the cloud and a reverse shock appearing either within or behind the cloud.

Whang, Y. V.↗

Episodes of terrestrial geologic activity during the past 260 million years - A quantitative approach

Records of major geologic events of the past about 260 Myr including: biologic extinction events, ocean-anoxic and black-shale events, major changes in sea level, major evaporite (salt) deposits, continental flood-basalt eruptions, first-order discontinuities in sea-floor spreading, and major mountain building events, have been aggregated and analyzed with moving-window and spectral techniques that facilitate recognition of clustering and possible cyclicity. Significant clustering of events suggests a model in which changes in rates and directions of sea-floor spreading ('ridge jumps') are associated with episodic rifting, volcanism, mountain building, global sea level and changes in the composition of the earth's atmosphere via the carbon cycle. The geologic data formally show a statistically significant underlying periodicity of 26.6 Myr for the Mesozoic and Cenozoic. Phase information suggests that the most recent maximum of the cycle occurred with the last 9 Myr, and may be close to the present time.

Rampino, Michael R.↗

A Physical Relationship Between Electron-Proton Temperature Equilibration and Mach Number in Fast Collisionless Shocks

The analysis of Balmer-dominated optical spectra from nonradiative (adiabatic) SNRs has shown that the ratio of the electron to proton temperature at the blast wave is close to unity at v(s) approximate or smaller than 400 km s(-1) but declines sharply down to the minimum value of m(e) /m(p) dictated by the jump conditions at shock speeds exceeding 2000 km s(-1). We propose a physical model for the heating of electrons and ions in non-cosmic-ray-dominated, strong shocks (v(s) 400 km s(-1)) wherein the electrons are heated by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of electron heating over a wide range of shock speeds, producing a relationship (T(e)/T(p))sub 0 proportional to v(-2/s) (proportional to M(-2)) that is fully consistent with the observations.

SNR(ONRADIATIVE SUPERNOVA REMNANTS)↗

Formation of the Potential Jump over the Geomagnetically Quiet Sunlit Polar Cap Region

We consider the formation of a potential drop over the Earth's polar cap during geomagnetically quiet daytime. The observed potential drop is primarily defined by the hydrogen, photoelectron, and polar rain fluxes ratios and depends strongly on the energy distribution of the photoelectron flux. Polar rain is an essential component of the model required for plasma quasineutrality. The potential distribution along the magnetic field line has two regions, with a small, gradual, potential drop of 3 to 4 volts and a potential jump. The value of the potential jump depends on the hydrogen ion to photoelectron flux ratio and is also controlled by polar rain electrons. With quasineutrality required at its upper boundary, the jump only occurs in the presence of polar rain and its location depends on the polar rain flux. Model predictions compare well with FAST (Fast Auroral SnapshoT) observations presented by Kitamura et al.(2012, https://doi.org/10.1029/2011JA017459).

Khazanov, G. V.↗

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗

Hidden Statistics Approach to Quantum Simulations

Recent advances in quantum information theory have inspired an explosion of interest in new quantum algorithms for solving hard computational (quantum and non-quantum) problems. The basic principle of quantum computation is that the quantum properties can be used to represent structure data, and that quantum mechanisms can be devised and built to perform operations with this data. Three basic non-classical properties of quantum mechanics superposition, entanglement, and direct-product decomposability were main reasons for optimism about capabilities of quantum computers that promised simultaneous processing of large massifs of highly correlated data. Unfortunately, these advantages of quantum mechanics came with a high price. One major problem is keeping the components of the computer in a coherent state, as the slightest interaction with the external world would cause the system to decohere. That is why the hardware implementation of a quantum computer is still unsolved. The basic idea of this work is to create a new kind of dynamical system that would preserve the main three properties of quantum physics superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. In other words, such a system would reinforce the advantages and minimize limitations of both quantum and classical aspects. Based upon a concept of hidden statistics, a new kind of dynamical system for simulation of Schroedinger equation is proposed. The system represents a modified Madelung version of Schroedinger equation. It preserves superposition, entanglement, and direct-product decomposability while allowing one to measure its state variables using classical methods. Such an optimal combination of characteristics is a perfect match for simulating quantum systems. The model includes a transitional component of quantum potential (that has been overlooked in previous treatment of the Madelung equation). The role of the transitional potential is to provide a jump from a deterministic state to a random state with prescribed probability density. This jump is triggered by blowup instability due to violation of Lipschitz condition generated by the quantum potential. As a result, the dynamics attains quantum properties on a classical scale. The model can be implemented physically as an analog VLSI-based (very-large-scale integration-based) computer, or numerically on a digital computer. This work opens a way of developing fundamentally new algorithms for quantum simulations of exponentially complex problems that expand NASA capabilities in conducting space activities. It has been illustrated that the complexity of simulations of particle interaction can be reduced from an exponential one to a polynomial one.

Zak, Michail↗

Thermal catastrophe in the plasma sheet boundary layer

This letter presents a first step towards a substorm model including particle heating and transport in the plasma sheet boundary layer (PSBL). The heating mechanism discussed is resonant absorption of Alfven waves. For some assumed MHD perturbation incident from the tail lobes onto the plasma sheet, the local heating rate in the PSBL has the form of a resonance function of the one-fluid plasma temperature. Balancing the local heating by convective transport of the heated plasma toward the central plasma sheet, an 'equation of state" is found for the steady-state PSBL whose solution has the form of a mathematical catastrophe: at a critical value of a parameter containing the incident power flux, the local density, and the convection velocity, the equilibrium temperature jumps discontinuously. Associating this temperature increase with the abrupt onset of the substorm expansion phase, the catastrophe model indicates at least three ways in which the onset may be triggered. Several other consequences related to substorm dynamics are suggested by the simple catastrophe model.

Smith, Robert A.↗

Unsupervised Detection of Planetary Craters by a Marked Point Process

With the launch of several planetary missions in the last decade, a large amount of planetary images is being acquired. Preferably, automatic and robust processing techniques need to be used for data analysis because of the huge amount of the acquired data. Here, the aim is to achieve a robust and general methodology for crater detection. A novel technique based on a marked point process is proposed. First, the contours in the image are extracted. The object boundaries are modeled as a configuration of an unknown number of random ellipses, i.e., the contour image is considered as a realization of a marked point process. Then, an energy function is defined, containing both an a priori energy and a likelihood term. The global minimum of this function is estimated by using reversible jump Monte-Carlo Markov chain dynamics and a simulated annealing scheme. The main idea behind marked point processes is to model objects within a stochastic framework: Marked point processes represent a very promising current approach in the stochastic image modeling and provide a powerful and methodologically rigorous framework to efficiently map and detect objects and structures in an image with an excellent robustness to noise. The proposed method for crater detection has several feasible applications. One such application area is image registration by matching the extracted features.

Troglio, G.↗

Development and Evaluation of Fault-Tolerant Flight Control Systems

The research is concerned with developing a new approach to enhancing fault tolerance of flight control systems. The original motivation for fault-tolerant control comes from the need for safe operation of control elements (e.g. actuators) in the event of hardware failures in high reliability systems. One such example is modem space vehicle subjected to actuator/sensor impairments. A major task in flight control is to revise the control policy to balance impairment detectability and to achieve sufficient robustness. This involves careful selection of types and parameters of the controllers and the impairment detecting filters used. It also involves a decision, upon the identification of some failures, on whether and how a control reconfiguration should take place in order to maintain a certain system performance level. In this project new flight dynamic model under uncertain flight conditions is considered, in which the effects of both ramp and jump faults are reflected. Stabilization algorithms based on neural network and adaptive method are derived. The control algorithms are shown to be effective in dealing with uncertain dynamics due to external disturbances and unpredictable faults. The overall strategy is easy to set up and the computation involved is much less as compared with other strategies. Computer simulation software is developed. A serious of simulation studies have been conducted with varying flight conditions.

Song, Yong D.↗

The 2006-2007 Active Phase Of Anomalous X-Ray Pulsar 4U 0142+61: Radiative and Timing Changes, Bursts, and Burst Spectral Features

After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March that lasted several months and included six X-ray bursts as well as many changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 2007 February 7. The burst durations ranged from 8-3x10(exp 3)s. The first five burst spectra are well modeled by blackbodies, with temperatures kT approx. 2 - 6 keV. However, the sixth burst had a complicated spectrum that is well characterized by a blackbody plus three emission features whose amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry into the active phase the pulsar showed a significant change in pulse morphology and a likely timing glitch. The glitch had a total frequency jump of (1.9+/-0.4)x10(exp -7) Hz, which recovered with a decay time of 17+/-2 days by more than the initial jump, implying a net spin-down of the pulsar. We discuss these events in the context of the magnetar model.

Gavril, Fotis P.↗

Model simulation studies to clarify the effect on saccadic eye movements of initial condition velocities set by the Vestibular Ocular Reflex (VOR)

Voluntary active head rotations produced vestibulo-ocular reflex eye movements (VOR) with the subject viewing a fixation target. When this target jumped, the size of the refixation saccades were a function of the ongoing initial velocity of the eye. Saccades made against the VOR were larger in magnitude. Simulation of a reciprocally innervated model eye movement provided results comparable to the experimental data. Most of the experimental effect appeared to be due to linear summation for saccades of 5 and 10 degree magnitude. For small saccades of 2.5 degrees, peripheral nonlinear interaction of state variables in the neuromuscular plant also played a role as proven by comparable behavior in the simulated model with known controller signals.

Nam, M. H.↗

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗

Hard X-ray time profiles and acceleration processes in large solar flares

The hard X-ray time profiles of the (1972) August 4 and 7 flares are investigated, taking into account a comparison of the time profiles of different energy channels. It is shown that for these flares the temporal features of the intensity profiles of higher energy channels are delayed with respect to those of channel 1. The delay time gradually increases to approximately 5 sec as the channel number increases from 1 to 5, and it jumps to approximately 15 sec for channels 6 and 7. A description is presented of a model in which the delay and other characteristics of the observed time profiles in channels 1-5 are self-consistently explained by the increase of the electron energy loss time with electron energy.

Bai, T.↗

Modulational development of nonlinear gravity-wave groups

Observations of the development of nonlinear surface gravity-wave groups are presented, and the amplitude and phase modulations are calculated using Hilbert-transform techniques. With increasing propagation distance and wave steepness, the phase modulation develops local phase reversals whose locations correspond to amplitude minima or nodes. The concomitant frequency modulation develops jumps or discontinuities. The observations are compared with recent similar results for wavetrains. The observations are modelled numerically using the cubic nonlinear Schroedinger equation. The motivation is twofold: to examine quantitatively the evolution of phase as well as amplitude modulation, and to test the inviscid predictions for the asymptotic behavior of groups versus long-time observations. Although dissipation rules out the recurrence, there is a long-time coherence of the groups. The phase modulation is found to distinguish between dispersive and soliton behavior.

Chereskin, T. K.↗

Direct simulation of transitional flow for hypersonic reentry conditions

This paper presents results of flowfield calculations for typical hypersonic reentry conditions encountered by the nose region of the Space Shuttle Orbiter. Most of the transitional flow regime is covered by the altitude range of 150 to 92 km. Calculations were made with the Direct Simulation Monte Carlo (DSMC) method that accounts for translational, rotational, vibrational, and chemical nonequilibrium effects. Comparison of the DSMC heating results with both Shuttle flight data and continuum predictions showed good agreement at the lowest altitude considered. However, as the altitude increased, the continuum predictions, which did not include slip effects, departed rapidly from the DSMC results by overpredicting both heating and drag. The results demonstrate the effects of rarefaction on the shock and the shock layer, along with the extent of the slip and temperature jump at the surface. Also, the sensitivity of the flow structure to the gas-surface interaction model, thermal accommodation, and surface catalysis are studied.

Moss, J. N.↗

Adaptive control and noise suppression by a variable-gain gradient algorithm

An adaptive control system based on normalized LMS filters is investigated. The finite impulse response of the nonparametric controller is adaptively estimated using a given reference model. Specifically, the following issues are addressed: The stability of the closed loop system is analyzed and heuristically established. Next, the adaptation process is studied for piecewise constant plant parameters. It is shown that by introducing a variable-gain in the gradient algorithm, a substantial reduction in the LMS adaptation rate can be achieved. Finally, process noise at the plant output generally causes a biased estimate of the controller. By introducing a noise suppression scheme, this bias can be substantially reduced and the response of the adapted system becomes very close to that of the reference model. Extensive computer simulations validate these and demonstrate assertions that the system can rapidly adapt to random jumps in plant parameters.

Merhav, S. J.↗