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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An internal variable constitutive model for the large deformation of metals at high temperatures

The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.

Brown, Stuart↗
NASA NTRS1988

Analytical study of the liquid phase transient behavior of a high temperature heat pipe

The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.

Roche, Gregory Lawrence↗
NASA NTRS1979

Propagation of waves in a medium with high radiation pressure

The propagation and mutual transformation of acoustic and thermal waves are investigated in media with a high radiative pressure. The equations of hydrodynamics for matter and the radiative transfer equations in a moving medium in the Eddington approximation are used in the investigation. Model problems of waves in a homogeneous medium with an abrupt jump in opacity and in a medium of variable opacity are presented. The characteristic and the times of variability are discussed. Amplitude for the brightness fluctuations for very massive stars are discussed.

Bisnovatyy-Kogan, G. S.↗
NASA NTRS2002

Prediction of Shock Wave Structure in Weakly Ionized Gas Flow by Solving MGD Equation

This paper reports the recent research results of shockwave structure predictions using a new developed code. The modified Rankine-Hugoniot relations across a standing normal shock wave are discussed and adopted to obtain jump conditions. Coupling a electrostatic body force to the Burnett equations, the weakly ionized flow field across the shock wave was solved. Results indicated that the Modified Rankine-Hugoniot equations for shock wave are valid for a wide range of ionization fraction. However, this model breaks down with small free stream Mach number and with large ionization fraction. The jump conditions also depend on the value of free stream pressure, temperature and density. The computed shock wave structure with ionization provides results, which indicated that shock wave strength may be reduced by existence of weakly ionized gas.

Deng, Z. T.↗
NASA NTRS1992

A model for the radio brightness of the supernova remnant 1987A

The recently observed rise in radio flux from SN1987A is earlier than predicted in models of the interaction with the observed line-emitting gas. The observed size of the radio source is roughly consistent with the expected position of the SN shock wave and the position of the termination shock of the blue supergiant progenitor wind. The radio increase is attributed to the interaction of the SN shock with this density jump. The required efficiency of electron acceleration and magnetic field amplification is larger than usually occurs in SNe and SNR. Possible reasons for the enhanced radio emission are particle acceleration by the weak reflected shock that travels through the interaction region, and the action of the Richtmyer-Meshkov instability near the initial termination-shock density jump. The action of these mechanisms is time-limited; the prediction of this model is that the radio flux should begin to decline within the next year.

Chevalier, Roger A.↗
NASA NTRS1987

The effect of nonlinearities on the dynamic response of a large shuttle payload

The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed under development when the program was cancelled. The first version, designated G-prime, was for planetary missions. The second version, designated G, was to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic jump that occurs in systems with nonlinear stiffness.

Sullivan, Timothy L.↗
NASA NTRS1987

The effect of nonlinearities on the dynamic response of a large Shuttle payload

The STS Centaur was designed to be a high energy upper stage for use with the Space Shuttle. Two versions were designed and under development when the program was cancelled. The first version, designated G-prime, was designed for planetary missions. The second version, designated G, was designed to place spacecraft in geosynchronous orbit. As a part of the STS Centaur finite-element model verification effort, test articles of both versions were subjected to a series of static tests. In addition the Centaur G-prime test article was subjected to a series of dynamic tests including a modal survey. Both the static and dynamic tests showed that nonlinearities existed in the Centaur and its support system. The support system included flight-like latches. The nonlinearities were particularly apparent in tests that loaded the forward support structure of the Centaur. These test results were used to aid in the development of two improved finite-element models. The first was a linear model, while the second contained nonlinear elements at the boundaries. Results from both models were compared with the transient response obtained from a step-relaxation or twang test. The linear model was able to accurately match the low frequency response found in the test data. However, only the nonlinear model was able to match higher frequency response that was present in some of the test data. In addition the nonlinear model was able to predict other nonlinear behavior such as the dynamic 'jump' that occurs in systems with nonlinear stiffness.

Sullivan, Timothy L.↗
NASA NTRS2023

JWST NIRCam Defocused Imaging Photometric Stability Performance and How It Can Sense Mirror Tilts

We use JWST NIRCam short-wavelength photometry to capture a transit lightcurve of the exoplanet HAT-P-14 b to assess performance as part of instrument commissioning. The short-wavelength precision is 152 ppm per 27 s integration as measured over the full time series compared to a theoretical limit of 107 ppm, after corrections to spatially correlated 1/f noise. Persistence effects from charge trapping are well fit by an exponential function with short characteristic timescales, settling on the order of 5–15 minutes. The short-wavelength defocused photometry is also uniquely well suited to measure the real-time wave-front error of JWST. Analysis of the images and reconstructed wave-front maps indicates that two different hexagonal primary mirror segments exhibited “tilt events,” where they changed orientation rapidly in less than ∼1.4 s. In some cases, the magnitude and timing of the flux jumps caused by tilt events can be accurately predicted with a telescope model. These tilt events can be sensed by simultaneous longer-wavelength NIRCam grism spectral images alone in the form of changes to the point-spread function, diagnosed from the full width at half maximum. They can also be sensed with the fine guidance sensor instrument from difference images. Tilt events possibly from sudden releases of stress in the backplane structure behind the mirrors were expected during the commissioning period because they were found in ground-based testing. Tilt events have shown signs of decreasing in frequency but have not disappeared completely. The detectors exhibit some minor (less than 1%) deviations from linear behavior in the first few groups of each integration, potentially impacting absolute fluxes and transit depths on bright targets, where only a handful of groups are possible. Overall, the noise is within 50% of the theoretical photon noise and read noise. This bodes well for high-precision measurements of transiting exoplanets and other time variable targets.

Everett Schlawin↗
NASA NTRS1946

Free-Spinning-Tunnel Tests of a 1/18-Scale Model of the Fairchild XNQ-1 Airplane, TED No. NACA 2398

Spin tests have been performed in the Langley 20-foot free-spinning tunnel on a 1/18-scale model of the Fairchild XNQ-1 airplane. The spin and recovery characteristics of the model were determined for the normal gross-weight loading and for two variations from this loading - center of gravity moved rearward and relative mass distribution increased along the fuselage. These tests were performed for two vertical-tail plan forms. The investigation also included simulated pilot-escape tests and rudder-force tests. The recovery characteristics of the model were satisfactory for all conditions tested by full reversal of the rudder and by simultaneous neutralization of the rudder and elevator. It was indicated that if necessary to escape from the spinning airplane, the pilot should jump from the outboard side of the fuselage and as far rearward as possible. Aa determined from spin model tests, the rudder pedal force required to reverse the rudder for recovery from the spin will be light.

Daughtridge, Lee T., Jr.↗
NASA NTRS1984

The adiabatic energy change of plasma electrons and the frame dependence of the cross-shock potential at collisionless magnetosonic shock waves

The adiabatic energy gain of electrons in the stationary electric and magnetic field structure of collisionless shock waves was examined analytically in reference to conditions of the earth's bow shock. The study was performed to characterize the behavior of electrons interacting with the cross-shock potential. A normal incidence frame (NIF) was adopted in order to calculate the reversible energy change across a time stationary shock, and comparisons were made with predictions made by the de Hoffman-Teller (HT) model (1950). The electron energy gain, about 20-50 eV, is demonstrated to be consistent with a 200-500 eV potential jump in the bow shock quasi-perpendicular geometry. The electrons lose energy working against the solar wind motional electric field. The reversible energy process is close to that modeled by HT, which predicts that the motional electric field vanishes and the electron energy gain from the electric potential is equated to the ion energy loss to the potential.

Goodrich, C. C.↗
NASA NTRS2014

The Hot Horizontal-Branch Stars in NGC288 - Effects of Diffusion and Stratification on Their Atmospheric Parameters*

Context. NGC288 is a globular cluster with a well developed blue horizontal branch covering the so-called u-jump which indicates the onset of diffusion. It is therefore well suited to study the effects of diffusion in blue horizontal branch (HB) stars. Aims. We compare observed abundances to predictions from stellar evolution models calculated with diffusion and from stratified atmospheric models. We verify the effect of using stratified model spectra to derive atmospheric parameters. In addition we investigate the nature of the overluminous blue HB stars around the u-jump. Methods. We define a new photometric index sz from uvby measurements that is gravity sensitive between 8 000K and 12 000 K. Using medium-resolution spectra and Stroemgren photometry we determine atmospheric parameters (Teff, logg) and abundances for the blue HB stars. We use both homogeneous and stratified model spectra for our spectroscopic analyses. Results. The atmospheric parameters and masses of the hot HB stars in NGC288 show a behaviour seen also in other clusters for temperatures between 9 000K and 14 000 K. Outside this temperature range, however, they follow rather the results found for such stars in (omega)Cen. The abundances derived from our observations are for most elements (except He and P) within the abundance range expected from evolutionary models that include the effects of atomic diffusion and assume a surface mixed mass of 10(exp −7) M. The abundances predicted by stratified model atmospheres are generally significantly more extreme than observed, except for Mg. The use of stratified model spectra to determine effective temperatures, surface gravities and masses moves the hotter stars to a closer agreement with canonical evolutionary predictions. Conclusions. Our results show definite promise towards solving the long-standing issue of surface gravity and mass discrepancies for hot HB stars, but there is still much work needed to arrive at a self-consistent solution.

horizontal branch (HB)↗
NASA NTRS2003

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗
NASA NTRS2008

Nonstandard Analysis and Jump Conditions for Converging Shock Waves

Nonstandard analysis is an area of modern mathematics which studies abstract number systems containing both infinitesimal and infinite numbers. This article applies nonstandard analysis to derive jump conditions for one-dimensional, converging shock waves in a compressible, inviscid, perfect gas. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions in the thermodynamic and fluid dynamic parameters occur smoothly across this interval. Predistributions of the Heaviside function and the Dirac delta measure are introduced to model the flow parameters across a shock wave. The equations of motion expressed in nonconservative form are then applied to derive unambiguous relationships between the jump functions for the flow parameters.

Baty, Roy S.↗
NASA NTRS2011

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 0.4 - 1.8 x 10(exp 3) s. The first five burst spectra are well modeled by blackbodies, with temperatures kT approx 2 - 9 keV. However, the sixth burst had a complicated spectrum that is well characterized by a blackbody plus two 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) x 10(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. Within the framework of the magnetar model, the net spin-down of the star could be explained by regions of the superfluid that rotate. slower than the rest. The bursts, flux enhancements, and pulse morphology changes can be explained as arising from crustal deformations due to stresses imposed by the highly twisted internal magnetic field. However, unlike other AXP outbursts, we cannot account for a major twist being implanted in the magnetosphere.

Gavriil, Fotis P.↗
NASA NTRS1990

Properties of large scale plasma flow during the early stage of the plasmaspheric refilling

The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.

Singh, Nagendra↗
NASA NTRS2019

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 duringgeomagnetically 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 thephotoelectron 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, potentialdrop of 3–4 V and a potential jump. The value of the potential jump depends on the hydrogen ion tophotoelectron flux ratio and is also controlled by polar rain electrons. With quasineutrality required at itsupper boundary, the jump only occurs in the presence of polar rain and its location depends on the polarrain flux. Model predictions compare well with FAST observations presented by Kitamura et al.(2012, https://doi.org/10.1029/2011JA017459).

Khazanov, George V.↗
NASA NTRS1995

Analytic MHD Theory for Earth's Bow Shock at Low Mach Numbers

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven M(A) and sonic M(s) Mach numbers are both much greater than 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first-order perturbation solution is identical to the solution found in the earlier theory. The second-order perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle theta between the normal to the bow shock and magnetic field is not close to 0 deg or 180 deg (the solutions are symmetric about 90 deg). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 deg. In this exceptional case the new solution does not improve on the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta approx. = 90 deg another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) is less than or = 2. MHD and gasdynamic simulations have produced empirical models in which the shock's standoff distance a(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between a(s) and X obtained from MHD simulations, a(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots of X and a(s) predicted for the Earth's bow shock. The plots show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the corresponding phenomenological relations in comparisons with data. Furthermore, significant differences exist between the standoff distances predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.

Grabbe, Crockett L.↗
NASA NTRS1993

Dynamics and control of coherent structures in the turbulent wall layer: An overview

We expand the velocity field in the vicinity of the wall in empirical eigenfunctions obtained from experiment. Truncating our system and using Galerkin projection, we obtain a closed set of non-linear ordinary differential equations with ten degrees of freedom. We find a rich dynamical behavior, including in particular a heteroclinic attracting orbit giving rise to intermittency. The intermittent jump from one attracting point to the other resembles in many respects the bursts observed in experiments. Specifically, the time between jumps and the duration of the jumps, is approximately that observed in a burst; the jump begins with the formation of a narrowed and intensified updraft, like the ejection phase of a burst, and is followed by a gentle, diffuse downdraft like the sweep phase of a burst. The magnitude of the Reynolds stress spike produced during a burst is limited by our truncation. The behavior is quite robust, much of it being due to the symmetries present (Aubry's group has examined dimensions up to 128 with persistence of the global behavior). We have examined eigenvalues and coefficients obtained from experiment, and from exact simulation, which differ in magnitude. Similar behavior is obtained in both cases; in the latter case, the heteroclinic orbits connect limit cycles instead of fixed points, corresponding to cross-stream waving of the streamwise rolls. The bifurcation diagram remains structurally similar, but somewhat distorted. The role of the pressure term is made clear - it triggers the intermittent jumps, which otherwise would occur at longer and longer intervals, as the system trajectory is attracted closer and closer to the heteroclinic cycle. The pressure term results in the jumps occurring at essentially random times, and the magnitude of the signal determines the average timing. Stretching of the wall region shows that the model is consistent with observations of polymer drag reduction. Change of the third order coefficients, corresponding to acceleration or deceleration of the mean flow, changes the heteroclinic cycles from attracting to repelling, increasing or decreasing the stability, in agreement with observations. The existence of fixed points is an artifact introduced by the projection; however, a decoupled model still displays the rich dynamics. Numerous assumptions made in Aubry et al. (1988) can now be proved exactly. Feeding back eigenfuncitons with the proper phase can delay the bursting, (the heteroclinic jump to the other fixed point), decreasing the drag. It is also possible to speed up the bursting, increasing mixing to control separation. Our approach is optimal for short time tracking in control.

Berkooz, Gal↗