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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

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 140 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.↗
NASA NTRS1990

Effects of material non-linearity on the residual stresses in a dendritic silicon crystal ribbon

Thermal stresses developed in a dendritic silicon crystal ribbon have been shown to cause plastic deformation and residual stresses in the ribbon. This paper presents an implementation of a numerical model proposed for thermoelastoplastic behavior of a material. The model has been used to study the effects of plasticity of silicon on the residual stresses. The material properties required to implement this model are all assumed, and the response of the material to the variations in these assumed parameters of the constitutive law and in the finite element mesh is investigated. The steady state growth process is observed to be periodic with nonzero residual stresses. Numerical difficulties are also encountered in the computer solution process, resulting in sharp jumps and large oscillations in the stress responses.

Ray, Sujit K.↗
NASA NTRS1988

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.↗
NASA NTRS1991

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.↗
NASA NTRS1992

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.↗
DOE OSTI2024

Vacancy diffusion barrier spectrum and diffusion correlation in multicomponent alloys

Vacancy diffusion serves a crucial role in many important kinetic behaviors and properties of multicomponent alloys. Essential questions, however, persist regarding how chemical complexity affects diffusion and what unique characteristics, if any, set these alloys apart from traditional metals. Using neural network kinetics model, we study vacancy diffusion in NbMoTa alloy across a broad temperature range (2600 to 800 K). Unlike pure metals, the two key diffusion parameters—diffusion correlation factor f and activation energy ΔG m —are not constant in alloys, but instead substantially decrease with decreasing temperature. This temperature dependence arises from a reduced number of active vacancy jump pathways at lower temperatures, leading to more correlated diffusion. Upon examining vacancy diffusion throughout the entire compositional space of the Nb-Mo-Ta system, we discover that the slowest vacancy diffusion surprisingly occurs in the non-equimolar region, rather than the equimolar concentration where the configurational entropy is highest. The diffusion barrier spectrum, characterizing the diffusion energy landscape, is an intrinsic material characteristic, which controls both f and ΔG m and, thereby, the diffusivity. Lastly, we find that the vacancy diffusion rate drops noticeably in the presence of local chemical order in the NbMoTa system, particularly for MoTa alloys with long-range B2 order.

Defects↗
DOE OSTI2023

Dynamic Heterogeneity and Kovacs’ Memory Effects in Supercooled Water

Understanding the properties of supercooled water is important for developing a comprehensive theory for liquid water and amorphous ices. Because of rapid crystallization for deeply supercooled water, experiments on it are typically carried out under conditions in which the temperature and/or pressure are rapidly changing. As a result, information on the structural relaxation kinetics of supercooled water as it approaches (metastable) equilibrium is useful for interpreting results obtained in this experimentally challenging region of phase space. Here, we used infrared spectroscopy and the fast time resolution obtained by transiently heating nanoscale water films to investigate relaxation kinetics (aging) in supercooled water. When the structural relaxation of the water films was followed using a temperature jump protocol analogous to the classic experiments of Kovacs, similar memory effects were observed. In particular, after suitable aging at one temperature, water’s structure displayed an extremum versus the number of heat pulses upon changing to a second temperature before eventually relaxing to a steady-state structure characteristic of that temperature. A random double well model based on the idea of dynamic heterogeneity in supercooled water accounts for the observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

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↗
DOE OSTI2024

Disorder Enhanced Thermalization in Interacting Many-Particle System

We introduce an extension of the non-equilibrium dynamical mean field theory to incorporate the effects of static random disorder in the dynamics of a many-particle system by integrating out different disorder configurations resulting in an effective time-dependent density-density interaction. We use this method to study the non-equilibrium transient dynamics of a system described by the Fermi Anderson-Hubbard model following an interaction and disorder quench. The method recovers the solution of the disorder-free case for which the system exhibits qualitatively distinct dynamical behaviors in the weak-coupling (prethermalization) and strong-coupling regimes (collapse-and-revival oscillations). However, we find that weak random disorder promotes thermalization. In the weak coupling regime, the jump in the quasiparticle weight in the prethermal regime is suppressed by random disorder while in the strong-coupling regime, random disorder reduces the amplitude of the quasiparticle weight oscillations. These results highlight the importance of disorder in the dynamics of realistic many-particle systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1988

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.↗
DOE OSTI2020

Adhesion-force micro-scale study of desiccating granular material

Experiments on five-, four-, three- and two-wet-hydrophilic-grain clusters were performed to investigate evolution of adhesion of granular media during drying on the micro-scale. The experiments show that the adhesion-force of a cluster initially grows at most to three times the original value before decreasing to zero by the end of evaporation. The adhesion-force is composed of capillary pressure force acting over the liquid/solid contact surface area, and surface tension forces acting over the three-phase contact perimeter length. This is in contrast with most macro-scale phenomenological models, in which the only desaturation process variables affecting strength are suction and saturation. Both the contact surface area and contact perimeter length are reduced to zero upon complete liquid evaporation. The morphology of an evaporating water body evolves through slow flow controlled by evaporation rate, interrupted by various modes of fast air entry, which are non-equilibrium jumps of liquid/gas interfaces (Haines jumps). The instabilities involve large adhesion force discontinuities and substantial water mass reconfiguration with water flow in an extremely short time, which makes the process transient. The reconfigurations can reduce the original multi-grain water clusters to four-, three- and two-grain clusters by way of three different instability modes: of thin-sheet instability, or meniscus snap-through instability, depending on the sign of the Gauss curvature of the liquid surface, or finally, for two-grain bridges only, a liquid wire pinch-off. For larger meso-scale assemblies, however, the global adhesion-force evolution is little affected by the jumps. The air entries are potential sites for drying cracks. The (approximately) calculated capillary pressure for two- and three-grain clusters, in no cases is seen to reach high values, predicted from water retention curves.

Engineering↗
DOE OSTI2020

Tuning the critical Li intercalation concentrations for MoX 2 bilayer phase transitions using classical and machine learning approaches

Transition metal dichalcogenides (TMDs) such as MoX 2 are known to undergo a structural phase transformation as well as a change in the electronic conductivity upon Li intercalation. These properties make them candidates for charge tunable ion-insertion materials that could be used in electro-chemical devices for neuromorphic computing applications. In this work we study the phase stability and electronic structure of Li-intercalated bilayer MoX 2 with X=S, Se or Te. Using first-principles calculations in combination with classical and machine learning modeling approaches we find that the energy needed to stabilize the conductive phase decreases with increasing atomic mass of the chalcogen atom X. A similar decreasing trend is found in the threshold Li concentration where the structural phase transition takes place. While the electronic conductivity increases with increasing ion concentration at low concentrations, we do not observe a conductivity jump at the phase transition point.

36 MATERIALS SCIENCE↗
NASA NTRS1972

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.↗
NASA NTRS2008

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.↗
NASA NTRS1973

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.↗
NASA NTRS1973

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.↗
NASA NTRS1987

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.↗
NASA NTRS1988

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.↗