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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A Multigrid Algorithm for Immersed Interface Problems

Many physical problems involve interior interfaces across which the coefficients in the problem, the solution, its derivatives, the flux, or the source term may have jumps. These interior interfaces may or may not align with a underlying Cartesian grid. Zhilin Li, in his dissertation, showed how to discretize such elliptic problems using only a Cartesian grid and the known jump conditions to second order accuracy. In this paper, we describe how to apply the full multigrid algorithm in this context. In particular, the restriction, interpolation, and coarse grid problem will be described. Numerical results for several model problems are given to demonstrate that good rates can be obtained even when jumps in the coefficients are large and do not align with the grid.

Adams, Loyce↗
NASA NTRS1999

Numerical and Analytical Solutions of Hypersonic Interactions Involving Surface Property Discontinuities

The local viscous-inviscid interaction field generated by a wall temperature jump on a flat plate in supersonic flow and on the windside of a Reusable Launch Vehicle in hypersonic flow is studied in detail by both a Navier-Stokes numerical code and an analytical triple-deck model. Treatment of the rapid heat transfer changes both upstream and downstream of the jump is included. Closed form relationships derived from the triple-deck theory are presented. The analytically predicted pressure and heating variations including upstream influence are found to be in generally good agreement with the Computational Fluid Dynamic (CFD) predictions. These analyses not only clarify the interactive physics involved but also are useful in preliminary design of thermal protection systems and as an insertable module to improve CFD code efficiency when applied to such small-scale interaction problems. The analyses only require conditions at the wall and boundary-layer edge which are easily extracted from a baseline, constant wall temperature, CFD solution.

Gnoffo, Peter A.↗
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↗
DOE OSTI2021

Unfolding bovine α-lactalbumin with T-jump: Characterizing disordered intermediates via time-resolved x-ray solution scattering and molecular dynamics simulations

The protein folding process often proceeds through partially folded transient states. Therefore, a structural understanding of these disordered states is crucial for developing mechanistic models of the folding process. Characterization of unfolded states remains challenging due to their disordered nature, and incorporating multiple methods is necessary. Combining the time-resolved x-ray solution scattering (TRXSS) signal with molecular dynamics (MD), we are able to characterize transient partially folded states of bovine α-lactalbumin, a model system widely used for investigation of molten globule states, during its unfolding triggered by a temperature jump. We track the unfolding process between 20 μs and 70 ms and demonstrate that it passes through three distinct kinetic states. The scattering signals associated with these transient species are then analyzed with TRXSS constrained MD simulations to produce protein structures that are compatible with the input signals. Without utilizing any experimentally extracted kinetic information, the constrained MD simulation successfully drove the protein to an intermediate molten globule state; signals for two later disordered states are refined to terminal unfolded states. From our examination of the structural characteristics of these disordered states, we discuss the implications disordered states have on the folding process, especially on the folding pathway. Finally, we discuss the potential applications and limitations of this method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
DOE OSTI2020

Thermal characterization of the build chamber in electron beam melting

Electron beam powder bed fusion, commonly termed electron beam melting (EBM), offers great versatility in multiple-part processes and can produce high quality as-built components due to low residual stresses. The EBM process is complex and requires careful thermal management, including uniform and consistent preheating of the powder bed, in order to ensure quality and consistency of the product. However, most of the simulations in the literature focus on the selective melting stage of the process. As of today, optimal conditions for initial pre-heating temperatures are only available for specific materials. Thus, in order to extend the EBM technology to other desired build materials, a much better understanding of the pre-heating stages is required. In this work, numerical and experimental approaches are combined in order to investigate the effects and sensitivities of heat removal from the build plate during EBM pre-heating stages. For this purpose, a carefully reconstructed numerical model of the build chamber of an ARCAM Q20 + machine is developed. It includes all main parts of the chamber and all relevant heat transfer mechanisms, whereas special attention is paid to radiation heat exchange between various bodies. In order to validate the model, dedicated experiments are performed, in which a system of thermocouples is installed in the build chamber, allowing direct measurement of the local temperatures of the start-plate and heat shields. A good agreement between the simulation and experimental findings is achieved, leading to a better understanding of the thermal processes characteristic to the pre-heating stages. This basic analysis is followed by a representative pre-heating case, where a powder bed is also considered. The energy required to obtain the desired pre-heating temperatures is evaluated, and the role of the powder bed in heat transfer within the chamber is assessed. The pre-heating stage, simulated in the present work, is supposed to create proper conditions for sintering and consequent melting of the powder, which are highly dependent on the local temperatures and heat transfer features. Thus, the findings of the reported work present a step towards a better understanding of the thermal processes that characterize EBM. The reported model can be further used to provide realistic boundary condition inputs for other meso- or macro-scale models as a function of time and geometry. The model can serve also for verification of machine settings, i.e., jump-safe and melt-safe ones which actually provide desired preheat, and for development of settings for new powders.

36 MATERIALS SCIENCE↗
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.↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
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.↗