Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JuMP model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

Energy management during the space shuttle transition

An approach to calculating optimal, gliding flight paths of the type associated with the space shuttle's transition from entry to cruising flight is presented. Kinetic energy and total energy (per unit weight) replace velocity and time in the dynamic equations, reducing the dimension and complexity of the problem. The capability for treating integral and terminal penalties (as well as Mach number effects) is retained in the numerical optimization; hence, stability and control boundaries can be observed as trajectories to the desired final energy, flight path angle, and range are determined. Numerical results show that the jump to the front-side of the L/D curve need not be made until the end of the transition and that the dynamic model provides a conservative range estimate. Alternatives for real time trajectory control are discussed.

Stengel, R. F.↗

Interstellar shock waves

The structure of interstellar shocks driven by supernova remnants and by expanding H II regions around early-type stars is discussed. Jump conditions are examined, along with shock fronts, post-shock relaxation layers, collisional shocks, collisionless shocks, nonradiative shocks, radiative atomic shocks, and shock models of observed nebulae. Effects of shock waves on interstellar molecules are examined, with reference to the chemistry behind shock fronts, infrared and vibrational-rotational cooling by molecules, and observations of shocked molecules. Some current problems and applications of the study of interstellar shocks are summarized, including the initiation of star formation by radiative shock waves, interstellar masers, the stability of shocks, particle acceleration in shocks, and shocks in galactic nuclei.

Mckee, C. F.↗

Wave particle interactions in the foot of the Saturnian bow shock

Voyager 1 plasma wave data show that the quasi-perpendicular, supercritical bow shock of Saturn exhibits the same plasma wave phenomenology observed near Jupiter. Using a quasi-linear model developed for the Jovian bow shock, it is shown that plasma waves at Saturn can generate a significant portion of the totalelectron temperature jump measured across the shock. In this respect, Saturn's bow shock more closely resembles Jupiter's than the earth's, where plasma waves contribute negligibly to the total electron temperature jump.

Moses, S. L.↗

Yang-Lee edge singularity triggered entanglement transition

Here we show that a class of $\mathscr{PT}$ symmetric non-Hermitian Hamiltonians realizing the Yang-Lee edge singularity exhibits an entanglement transition in the long-time steady state evolved under the Hamiltonian. Such a transition is induced by a level crossing triggered by the critical point associated with the Yang-Lee singularity and hence is first order in nature. At the transition, the entanglement entropy of the steady state jumps discontinuously from a volume-law to an area-law scaling. We exemplify this mechanism using a one-dimensional transverse field Ising model with additional imaginary fields, as well as the spin-1 Blume-Capel model and the three-state Potts model. We further make a connection to the forced-measurement induced entanglement transition in a Floquet nonunitary circuit subject to continuous measurements followed by post-selections. Our results demonstrate a new mechanism for entanglement transitions in non-Hermitian systems harboring a critical point.

36 MATERIALS SCIENCE↗

Stochastic Stability of Nonlinear Sampled Data Systems with a Jump Linear Controller

This paper analyzes the stability of a sampled- data system consisting of a deterministic, nonlinear, time- invariant, continuous-time plant and a stochastic, discrete- time, jump linear controller. The jump linear controller mod- els, for example, computer systems and communication net- works that are subject to stochastic upsets or disruptions. This sampled-data model has been used in the analysis and design of fault-tolerant systems and computer-control systems with random communication delays without taking into account the inter-sample response. To analyze stability, appropriate topologies are introduced for the signal spaces of the sampled- data system. With these topologies, the ideal sampling and zero-order-hold operators are shown to be measurable maps. This paper shows that the known equivalence between the stability of a deterministic, linear sampled-data system and its associated discrete-time representation as well as between a nonlinear sampled-data system and a linearized representation holds even in a stochastic framework.

Gonzalez, Oscar R.↗

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗

Defensive Cybersecurity Architecture Design Using Force-on-Force Cyber-Physical Modeling

Currently, nuclear power plant physical security systems are highly dependent on air-gaps as a protective measure against cyber-threats. Cyber-physical threats become more likely as advanced cyber-threat capabilities to jump air-gaps transition into common use. Defending against the emerging threat of cyber-enabled physical intrusions is poorly understood. The consequence of these cyber-physical attacks has no quantitative analysis method to inform risk-informed, performance-based cybersecurity approaches. By modifying the physical security simulation tool Dante, cyber-physical threat consequence was able to be analyzed on a notional facility. The results of this analysis are used to design a Defensive Cybersecurity Architecture (DCSA) for the physical security system to produce example resilience measures for this notional facility. A DCSA defines security levels to provide a graded approach for defending plant functions, and security zones for trusted communication between systems. This approach can be applied to real world systems to produce physical protection systems and response measures that are resilient to cyber-physical threats.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The computational treatment of supersonic wake flows in INCOREL

Supersonic wake flows are computed in both conical and three-dimensional full potential flows by modeling the wake as a planar cut with a potential discontinuity. Continuity of pressure is satisfied exactly on the wake cut by computing the necessary potential jump at each wake cut grid point. In the conical wake problem, the jump in potential vanishes in the crossflow symmetry plane. Generally, good correlation is achieved for fully three-dimensional wake flows with the planar wake approximation. Matching pressures in the isolated wake flow leads to small variations in the spanwise jump in potential. On the other hand, a small variation in the jump in potential near the wing tip flow can lead to a significant loading in the wake flow if the pressures are not matched.

Siclari, M. J.↗

Analysis of suprathermal electron properties at the magnetic pile-up boundary of Comet P/Halley

Among the plasma discontinuities detected by the Giotto spacecraft around Comet P/Halley, the magnetic pile-up boundary, located at about 135,000 km from the nucleus, has a sharpness which was not foreseen by theoretical models. At this boundary, which marks the beginning of the region where the field lines draped around the nucleus have been piled up, the magnetic field jumps sharply. Electron measurements provided by the RPA experiment show that a clear plasma discontinuity coincides with this magnetic feature. Significant changes occur here in the suprathermal electron distribution function. A magneto-plasma sheet is clearly defined after the boundary. Inside this sheet, close correlations exist between the parameters describing the magnetic field and the electron population. The polytropic equation of state governing the suprathermal electrons in the sheet has been deduced from RPA measurements. Some implications of this law are discussed.

Mazelle, C.↗

Termination shock - Solar cycle variations of location and speed

Plasma and magnetic field data from Voyagers 1 and 2 over an 11-year period (1978-1988) together with a one-fluid MHD model are used to study the motion of the termination shock. The observed large fluctuations in solar wind parameters cause variations in the speed, in the jump conditions, and in the location of the termination shock. If the interstellar stagnation pressure is equivalent to a magnetic field of 0.5 nT plus draping enhancements, the termination shock moves between 88 AU and 102 AU. The location of the termination shock is anticorrelated with the sunspot number. The termination shock is very strong with pressure ratios of a few thousands, the fast Mach number is about 12, and the plasma beta ratio jumps from about 0.1 to about 14. During the declining phase of the solar cycle, the shock speed is often greater than or comparable to the spacecraft speed. If the crossing of the termination shock by a spacecraft occurs during the declining phase, the spacecraft will cross the termination shock 3 or more times over a period of about 4 years.

Whang, Y. C.↗

Altered astronaut lower limb and mass center kinematics in downward jumping following space flight

Astronauts exposed to the microgravity conditions encountered during space flight exhibit postural and gait instabilities upon return to earth that could impair critical postflight performance. The aim of the present study was to determine the effects of microgravity exposure on astronauts' performance of two-footed jump landings. Nine astronauts from several Space Shuttle missions were tested both preflight and postflight with a series of voluntary, two-footed downward hops from a 30-cm-high step. A video-based, three-dimensional motion-analysis system permitted calculation of body segment positions and joint angular displacements. Phase-plane plots of knee, hip, and ankle angular velocities compared with the corresponding joint angles were used to describe the lower limb kinematics during jump landings. The position of the whole-body center of mass (COM) was also estimated in the sagittal plane using an eight-segment body model. Four of nine subjects exhibited expanded phase-plane portraits postflight, with significant increases in peak joint flexion angles and flexion rates following space flight. In contrast, two subjects showed significant contractions of their phase-plane portraits postflight and three subjects showed insignificant overall changes after space flight. Analysis of the vertical COM motion generally supported the joint angle results. Subjects with expanded joint angle phase-plane portraits postflight exhibited larger downward deviations of the COM and longer times from impact to peak deflection, as well as lower upward recovery velocities. Subjects with postflight joint angle phase-plane contraction demonstrated opposite effects in the COM motion. The joint kinematics results indicated the existence of two contrasting response modes due to microgravity exposure. Most subjects exhibited "compliant" impact absorption postflight, consistent with decreased limb stiffness and damping, and a reduction in the bandwidth of the postural control system. Fewer subjects showed "stiff" behavior after space flight, where contractions in the phase-plane portraits pointed to an increase in control bandwidth. The changes appeared to result from adaptive modifications in the control of lower limb impedance. A simple 2nd-order model of the vertical COM motion indicated that changes in the effective vertical stiffness of the legs can predict key features of the postflight performance. Compliant responses may reflect inflight adaptation due to altered demands on the postural control system in microgravity, while stiff behavior may result from overcompensation postflight for the presumed reduction in limb stiffness inflight.

Clinical Trial↗

Local spin polarizations in relativistic heavy-ion collisions

Here, based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total-angular-momentum-conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sheet Flows, Avalanches, and Dune Evolution on Earth and Mars

This investigation is a collaboration between researchers at Cornell University, the University of Florida, and the University of Rennes 1, France. Flow modeling at Cornell University focused on mechanisms for the suspension and transport of wind-blown sand that are important in both terrestrial and Martian environments. These mechanisms include the saltation (or jumping) of grains, collisions between grains, and the interaction of grains with the velocity fluctuations of the turbulent wind. Of particular interest are sheet flows; these are relatively thin, highly concentrated regions of grains flowing near the ground under the influence of a strong turbulent wind. In them, the grains are suspended by interparticle collisions. Sheet flows may be relatively rare events, but they have the capacity to move great amounts of sand. In order to describe sheet flows, a turbulent mixture theory was formulated for particles in a fluid in which fluctuations in the volume fiaction of the particles take place on the scale of the turbulent eddies. Ensemble averaged equations for particle and fluid mass, momentum, and energy and fluid rate of dissipation were expressed in terms of Farve (concentration) averaged velocities and the associated velocity fluctuations. Correlations that describe the turbulent suspension of particles and dissipation of turbulent energy of both phases due to fluid particle interactions were modeled and boundary conditions at the bed and at the upper surface of the collisional flow were formulated. The boundary conditions at the upper surface were tested in a numerical simulation developed at the University of Florida. Steady and unsteady solutions for steady and unsteady fully-developed flows were obtained over a range of wind speeds fiom the lowest for which collisional between particles occurred to at which turbulent suspension is found to dominate collisional suspension. Below the value of the wind speed at which collisions between particles were unimportant, numerical solutions were obtained for the velocity distribution function and the resulting fields of concentration, particle and gas mean velocity, and particle shear stress for the steady two-dimensional saltation of spherical sand particles driven by a turbulent wind over a bed characterized by a simple relationship (the splash function) between the properties of incoming particles and those of the rebounding particles and other particles ejected fiom the bed. At the University of Rennes 1, experiments devoted to the characterization of the splash function for beds consisting of either random or ordered arrays of spheres in two- dimensions were completed. These indicated the role played by the packing geometry in the rebound and ejection of grains. Preliminary experiments on response of a three- dimensional collision bed to a collision with a single particle were performed. Data was taken with a single camera focused on the plane of collision. Here, for example, the decrease of the effective coefficient of restitution of the bed with an increase of the angle of incidence of the incoming particle has been measured. Other experiments on avalanches at Rennes studied the properties of the flows of particles that are responsible for the motion of the leeward side of a dune. In these, the dependence of the initiation of avalanches on the packing and depth of the particles was measured. Particle migration was studied in inclined flows of a binary mixture of disks and the mechanisms of diffision and segregation were isolated and characterized. The influence of side wall on dense, rapid inclined flows was measured and shown to be the reason why the angle of the free surface in such flows can exceed the static angle of repose. Future research will be devoted to a better understanding the transition between saltating (collisionless) and collisional flows as the wind speed the increases. This will involve the understanding of the evolution of the splash function as clisions with the bed become more numerous, more frequent, and more violent.

Source record↗

A-LEAF

A-LEAF is an electricity market simulation framework. It includes models for short-term market operations and long-term investment planning of generation, storage, and transmission. Model formulations include least cost and two-stage market equilibrium. Simulations can be flexibly structured over one or more years, including options for periodic samples or jumping to a single future year. Years can be simulated as 365 consecutive days at hourly or 5-minute sample rates, or they can be subsampled using a representative day group selection algorithm.

MANN, WILLIAM↗

Quantum crystals in neutron stars

Using the many-body techniques appropriate for quantum crystals it is shown that the deep interior of a neutron star is most likely an orderly arrangement of neutrons, protons and hyperons forming a solid. It is shown that a liquid or gas arrangement would produce higher energy. If so, a neutron star can be viewed as two solids (crust and core) permeated by a layer of ordinary or (perhaps) superfluid liquid. Astronomical evidence is in favor of such a structure: the sudden jumps in the periods of the Crab and Vela pulsars that differ by a factor of about 100 can be easily explained by the star-quake model. If the Crab is less massive than Vela (i.e., if it is not dense enough to have a solid core), the star-quakes take place in the crust whereas for Vela they occur in the core.

Canuto, V.↗

Hypersonic merged layer blunt body flows with wakes

An attempt is made here to understand the basic physics of the flowfield with wake on a blunt body of revolution under hypersonic rarefied conditions. A merged layer model of flow is envisioned. Full steady-state Navier-Stokes equations in spherical polar coordinate system are computed from the surface with slip and temperature jump conditions to the free stream by the Accelerated Successive Replacement method of numerical integration. Analysis is developed for bodies of arbitrary shape, but actual computations have been carried out for a sphere and sphere-cone body. Particular attention is paid to set the limit of the onset of separation, wake closure, shear-layer impingement, formation and dissipation of the shocks in the flowfield. Validity of the results is established by comparing the present results for sphere with the corresponding results of the SOFIA code in the common region of their validity and with the experimental data.

Jain, Amolak C.↗

Automatic control of cryogenic wind tunnels

Inadequate Reynolds number similarity in testing of scaled models affects the quality of aerodynamic data from wind tunnels. This is due to scale effects of boundary-layer shock wave interaction which is likely to be severe at transonic speeds. The idea of operation of wind tunnels using test gas cooled to cryogenic temperatures has yielded a quantrum jump in the ability to realize full scale Reynolds number flow similarity in small transonic tunnels. In such tunnels, the basic flow control problem consists of obtaining and maintaining the desired test section flow parameters. Mach number, Reynolds number, and dynamic pressure are the three flow parameters that are usually required to be kept constant during the period of model aerodynamic data acquisition. The series of activity involved in modeling, control law development, mechanization of the control laws on a microcomputer, and the performance of a globally stable automatic control system for the 0.3-m Transonic Cryogenic Tunnel (TCT) are discussed. A lumped multi-variable nonlinear dynamic model of the cryogenic tunnel, generation of a set of linear control laws for small perturbation, and nonlinear control strategy for large set point changes including tunnel trajectory control are described. The details of mechanization of the control laws on a 16 bit microcomputer system, the software features, operator interface, the display and safety are discussed. The controller is shown to provide globally stable and reliable temperature control to + or - 0.2 K, pressure to + or - 0.07 psi and Mach number to + or - 0.002 of the set point value. This performance is obtained both during large set point commands as for a tunnel cooldown, and during aerodynamic data acquisition with intrusive activity like geometrical changes in the test section such as angle of attack changes, drag rake movements, wall adaptation and sidewall boundary-layer removal. Feasibility of the use of an automatic Reynolds number control mode with fixed Mach number control is demonstrated.

Balakrishna, S.↗

Representing Matrix Cracks Through Decomposition of the Deformation Gradient Tensor in Continuum Damage Mechanics Methods

A method is presented to represent the large-deformation kinematics of intraply matrix cracks and delaminations in continuum damage mechanics (CDM) constitutive material models. The method involves the additive decomposition of the deformation gradient tensor into 'crack' and 'bulk material' components. The response of the intact bulk material is represented by a reduced deformation gradient tensor, and the opening of an embedded cohesive interface is represented by a normalized cohesive displacement-jump vector. The rotation of the embedded interface is tracked as the material deforms and as the crack opens. The distribution of the total local deformation between the bulk material and the cohesive interface components is determined by minimizing the difference between the cohesive stress and the bulk material stress projected onto the cohesive interface. The improvements to the accuracy of CDM models that incorporate the presented method over existing approaches are demonstrated for a single element subjected to simple shear deformation and for a finite element model of a unidirectional open-hole tension specimen. The material model is implemented as a VUMAT user subroutine for the Abaqus/Explicit finite element software. The presented deformation gradient decomposition method reduces the artificial load transfer across matrix cracks subjected to large shearing deformations, and avoids the spurious secondary failure modes that often occur in analyses based on conventional progressive damage models.

Leone, Frank A., Jr.↗