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 379 records · Page 21

JPL pulsar timing observations. I - The Vela pulsar

It is found through an analysis of 12 years of arrival time data for PSR 0833-45, by fitting a model of the pulse period to the data in a least-squares sense, that there exist regions of exceptionally stable behavior in which rms residuals are typically about 40 microsec. In addition, each residual exhibits a significant period parameter second derivative which is found to change little within a given interjump era. Parameter changes in the period P, and the derivatives introduced by each of the four large observed jumps in period, are recomputed to demonstrate that the post-jump decay is nearly exponential in form. It is found that random fluctuations in P are an order of magnitude less than the systematic variations caused by the seond derivative period parameter, whose component associated with the general slowdown of the pulsar can be reliably measured when the observations span 55 years or more.

Downs, G. S.↗

The Dryline on 22 May 2002 during IHOP_2002: Convective-Scale Measurements at the Profiling Site

A detailed analysis of the structure of a double dryline observed over the Oklahoma panhandle during the first International H2O Project (IHOP_2002) convective initiation (CI) mission on 22 May 2002 is presented. A unique and unprecedented set of high temporal and spatial resolution measurements of water vapor mixing ratio, wind, and boundary layer structure parameters were acquired using the National Aeronautics and Space Administration (NASA) scanning Raman lidar (SRL), the Goddard Lidar Observatory for Winds (GLOW), and the Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE), respectively. These measurements are combined with the vertical velocity measurements derived from the National Center for Atmospheric Research (NCAR) Multiple Antenna Profiler Radar (MAPR) and radar structure function from the high-resolution University of Massachusetts frequency-modulated continuous-wave (FMCW) radar to reveal the evolution and structure of the late afternoon double-dryline boundary layer. The eastern dryline advanced and then retreated over the Homestead profiling site in the Oklahoma panhandle, providing conditions ripe for a detailed observation of the small-scale variability within the boundary layer and the dryline. In situ aircraft data, dropsonde and radiosonde data, along with NCAR S-band dual-polarization Doppler radar (S-Pol) measurements, are also used to provide the larger-scale picture of the double-dryline environment. Moisture and temperature jumps of about 3 g kg(sup -1) and 1 -2 K, respectively, were observed across the eastern radar fine line (dryline), more than the moisture jumps (1-2 g kg(sup -1)) observed across the western radar fine line (secondary dryline). Most updraft plumes observed were located on the moist side of the eastern dryline with vertical velocities exceeding 3 m s(sup -1) and variable horizontal widths of 2-5 km, although some were as wide as 7-8 km. These updrafts were up to 1.5 g kg(sup -1) moister than the surrounding environment. Although models suggested deep convection over the Oklahoma panhandle and several cloud lines were observed near the dryline, the dryline itself did not initiate any storms over the intensive observation region (IOR). Possible reasons for this lack of convection are discussed. Strong capping inversion and moisture detrainment between the lifting condensation level and the level of free convection related to an overriding drier air, together with the relatively small near-surface moisture values (less than 10 g kg(sup -1)), were detrimental to CI in this case.

Demoz, Belay↗

Bounded state space

This investigation is divided functionally into three different areas: (1) study of bounded state space, (2) nonlinear smoothing theory, and (3) system identification. (1) Study of bounded state space: necessary and sufficient conditions for an optimal control are obtained for a bounded state space optimal control problem. The difficulty of determining the so-called jump conditions is eliminated; however, the problem of determining the points where the response either enters or leaves the boundary still remains unsolved. (2) Nonlinear smoothing theory: nonlinear fixed-interval, fixed-point and fixed-lag smoothing of a random signal generated by a stochastic differential equation are investigated. Results on the asymptotic stability of a linear constant-parameter fixed-interval smoothing filter are obtained. (3) System identification: a particular stochastic modelling problem is solved. An Ito stochastic integral equation is used to mathematically model a black box having multiple inputs and multiple outputs. A new method for identifying system parameters is presented.

Eyman, E. D.↗

A Strong Merger Shock in Abell 665

Deep (103 ks) Chandra observations of Abell 665 have revealed rich structures in this merging galaxy cluster, including a strong shock and two cold fronts. The newly discovered shock has a Mach number of M =3.0 +/- 0.6, propagating in front of a cold disrupted cloud. This makes Abell 665 the second cluster, after the Bullet cluster, where a strong merger shock of M is approximately 3 has been detected. The shock velocity from jump conditions is consistent with (2.7 +/- 0.7) × 10(exp 3) km s(exp −1). The new data also reveal a prominent southern cold front with potentially heated gas ahead of it. Abell 665 also hosts a giant radio halo. There is a hint of diffuse radio emission extending to the shock at the north, which needs to be examined with better radio data. This new strong shock provides a great opportunity to study the reacceleration model with the X-ray and radio data combined.

galaxies: clusters: general- galaxies: clisters; i↗

The Origin of Hot Subluminous Horizontal-Branch Stars in (omega) Centauri and NGC 2808

Hot subluminous stars lying up to 0.7 mag below the extreme horizontal branch (EHB) are found in the ultraviolet (UV) color magnitude diagrams of both (omega) Cen and NGC 2808. In order to explore the evolutionary status of these subluminous stars, we have evolved a set of low-mass stars continuously from the main sequence through the helium-core flash to the HB (horizontal branch) for a wide range in the mass loss along the red-giant branch (RGB). Stars with the largest mass loss evolve off the RGB to high effective temperatures before igniting helium in their cores. Our results indicate that the subluminous EHB stars, as well as the gap within the EHB of NGC 2808, can be explained if these stars undergo a late helium-core flash while descending the white-dwarf cooling curve. Under these conditions the convection zone produced by the helium flash will penetrate into the stellar envelope, thereby mixing most, if not all, of the envelope hydrogen into the hot helium-burning interior, where it is rapidly consumed. This phenomenon is analogous to the 'born-again' scenario for producing hydrogen-deficient stars following a very late helium-shell flash. This 'flash mixing' of the stellar envelope greatly enhances the envelope helium and carbon abundances and, as a result, leads to a discontinuous jump in the HB effective temperature. We argue that the EHB gap in NGC 2808 is associated with this theoretically predicted dichotomy in the HB morphology. Using new helium- and carbon-rich stellar atmospheres, we show that these changes in the envelope abundances of the flash-mixed stars will suppress the UV flux by the amount needed to explain the hot subluminous EHB stars in (omega) Cen and NGC 2808. Moreover, we demonstrate that models without flash mixing lie, at most, only approximately 0.1 mag below the EHB, and hence fail to explain the observations. Flash mixing may also provide a new evolutionary channel for producing the high gravity, helium-rich sdO and sdB stars.

Sweigart, Allen V.↗

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. Calculations are performed for varying initial air-to-fuel-mass ratios, initial cluster radii, ambient gas temperatures and initial drop temperatures. Results are presented for ratios of fuel mass burned to fuel mass lost from the cluster between drop ignition and drop disappearance, fuel burned fractions at ignition and at the moment of drop disappearance, and jump conditions at ignition.

Bellan, J.↗

Extended behavioral decomposition for estimating ultrahigh reliability

Limitations of some analytic techniques in approximating the reliability of life-critical electronic systems are discussed, and a framework for the specification of recovery and fault-handling submodels is suggested. The framework makes full use of the instantaneous jump theorem by viewing the collection of interfering, premature exits from any fault handling and recovery submodel as defining a new, competing process submodel. This approach allows a greater flexibility in submodel representation, since submodels may contain arbitrary entrance arcs, exit arcs, and competing, interfering transitions with arbitrary destinations. Since the effects of near-coincident faults need not be represented as system failure events, the reliability estimates produced need not be unduly conservative. Comparisons on small models, where exact results can be computed, show substantial improvement in accuracy over earlier techniques. Implementation of the technique in an X-windows-based system, XHARP, is described.

Geist, Robert↗

The flow of a thin liquid film on a stationary and rotating disk. II - Theoretical prediction

The existing theoretical models are improved and a systematic procedure to compute the free surface flow of a thin liquid film is suggested. The solutions for axisymmetric radial flow on a stationary horizontal disk and for the disk rotating around its axis are presented. The theoretical predictions are compared with the experimental data presented in Part I of this report. The analysis shows results for both supercritical and subcritical flows and the flow structure in the vicinity of a hydraulic jump which isolates these two flow types. The detailed flow structure in a hydraulic jump was computed and shown to contain regions of separation including a 'surface roller'. The effects of surface tension are found to be important near the outer edge of the disk where the fluid experiences a free fall. At other locations, the surface tension is negligible. For a rotating disk, the frictional resistance in the angular direction is found to be as important as that in the radial direction.

Rahman, M. M.↗

A New Look into the Effect of Large Drops on Radiative Transfer Process

Recent studies indicate that a cloudy atmosphere absorbs more solar radiation than any current 1D or 3D radiation model can predict. The excess absorption is not large, perhaps 10-15 W/sq m or less, but any such systematic bias is of concern since radiative transfer models are assumed to be sufficiently accurate for remote sensing applications and climate modeling. The most natural explanation would be that models do not capture real 3D cloud structure and, as a consequence, their photon path lengths are too short. However, extensive calculations, using increasingly realistic 3D cloud structures, failed to produce photon paths long enough to explain the excess absorption. Other possible explanations have also been unsuccessful so, at this point, conventional models seem to offer no solution to this puzzle. The weakest link in conventional models is the way a size distribution of cloud particles is mathematically handled. Basically, real particles are replaced with a single average particle. This "ensemble assumption" assumes that all particle sizes are well represented in any given elementary volume. But the concentration of larger particles can be so low that this assumption is significantly violated. We show how a different mathematical route, using the concept of a cumulative distribution, avoids the ensemble assumption. The cumulative distribution has jumps, or steps, corresponding to the rarer sizes. These jumps result in an additional term, a kind of Green's function, in the solution of the radiative transfer equation. Solving the cloud radiative transfer equation with the measured particle distributions, described in a cumulative rather than an ensemble fashion, may lead to increased cloud absorption of the magnitude observed.

Marshak, Alexander↗

Isolator Shock Train Response to Dual-Mode Scramjet Throttling-Like Forcing

Typical dual-mode scramjet designs employ an isolator component situated between the inlet and combustor. The isolator fluidically buffers the inlet from the combustor via the formation of a shock train, an interconnected series of shock-boundary layer interactions. The isolator is required during ramjet-engine operation and at the lower edges of the scramjet operational envelope. In an operational engine, transient behavior imposes fluidic dynamics due to attitude adjustments, combustor throttling, aberrant resonances within the combustor, etc. To accommodate these transients, the shock train adjusts on both temporal and spatial scales. Sufficiently large downstream pressure increases in the combustor play a critical role during unstart, whereby strong compression waves propagate upstream in the isolator to ultimately force the shock train out the inlet, if it cannot stabilize to the pressure increase. This critically reduces the engine mass flow rate, potentially causing a catastrophic reduction in thrust and loss of the vehicle. A nearly instantaneous unstart can, in principle, be predicted by a ratio of the normal shock pressure jump in an inviscid isolator. However, lower-magnitude perturbations, below this bound, also induce unstarts. Understanding the transient behavior inherent in unsteady isolator operation is critical to reducing safety margin in isolator length and refining engine control algorithms to improve engine performance. The Isolator Dynamics Research Lab (IDRL) is a cold-flow, direct-connect isolator facility ideal for assessing the fundamental fluid dynamics of shock trains. The internal flow of a scramjet isolator is modeled using a 2D, symmetric Mach 2.5 nozzle followed by a 695 mm long, 25.4x50.8 mm rectangular isolator test section. A blockage cone on a high-speed ball screw actuator downstream of the isolator is pulsed to simulate the throttling behavior in an operational combustor. Wall flush high-speed static pressure transducers detect the streamwise shock train within the isolator. Time delays between the onset of the throttling and the response motion of the shock train are analyzed to assess the controlling scales and parameters. From this analysis, the response time of the isolator shock train system is found to depend primarily on the isolator back pressure time-rate-of-change but quickly saturates to a maximum propagation speed along the isolator. These results inform dual-mode scramjet engine designers of the critical parameters and limits of an unstart resistant engine.

scramjet↗

Hybridization of Bogoliubov Quasiparticles between Adjacent CuO 2 Layers in the Triple-Layer Cuprate Bi 2 Sr 2 Ca 2 Cu 3 O 10 + δ Studied by Angle-Resolved Photoemission Spectroscopy

Hybridization of Bogoliubov quasiparticles (BQPs) between the CuO 2 layers in the triple-layer cuprate high-temperature superconductor Bi 2 Sr 2 Cu 2 Cu 3 O 10+δ is studied by angle-resolved photoemission spectroscopy (ARPES). In the superconducting state, an anticrossing gap opens between the outer- and inner-BQP bands, which we attribute primarily to interlayer single-particle hopping with possible contributions from interlayer Cooper pairing. We find that the d-wave superconducting gap of both BQP bands smoothly develops with momentum without an abrupt jump in contrast to a previous ARPES study. Hybridization between the BQPs also gradually increases in going from the off nodal to the antinodal region, which is explained by the momentum dependence of the interlayer single-particle hopping. As possible mechanisms for the enhancement of the superconducting transition temperature, the hybridization between the BQPs as well as the combination of phonon modes of the triple CuO 2 layers and spin fluctuations represented by a four-well model are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗

Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.

Max, C. E.↗

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using as a model system ZnO, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

36 MATERIALS SCIENCE↗

Tuning the critical Li intercalation concentrations for MoX 2 (X = S, Se, or Te) bilayer phase transitions

Transition metal dichalcogenides (TMDs), such as MoS 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. In this work we study the phase stability and electronic structure of H and T’ Li-intercalated MoX 2 bilayers with X=S, Se, or Te. Here, using first-principles calculations in combination with classical and machine learning modeling approaches, we find that the H phase is more stable at low Li concentration for all X, and the critical Li concentration at which the T’→H transition occurs decreases with increasing mass of X. Furthermore the relative free energy of the two phases becomes less sensitive to Li insertion with increasing atomic mass of the chalcogen atom X. While the electronic conductivity increases with increasing ion concentration at low concentrations, we do not observe a (positive) conductivity jump at the phase transition from H to T’.

36 MATERIALS SCIENCE↗

Spin Berry points as crucial for ultrafast demagnetization

Laser-induced ultrafast demagnetization has puzzled researchers around the world for over two decades. Intrinsic complexity in electronic, magnetic and phononic subsystems is difficult to understand microscopically. So far, it is not possible to explain demagnetization using a single mechanism, which suggests a crucial piece of information still missing. In this paper, we return to a fundamental aspect of physics: spin and its change within each band in the entire Brillouin zone. We employ face-centered cubic (fcc) Ni as an example and use an extremely dense k mesh to map out spin changes for every band close to the Fermi level along all the high symmetry lines. To our surprise, spin angular momentum at some special k points abruptly changes from [Formula: see text] to [Formula: see text] simply by moving from one crystal momentum point to the next. This explains why intraband transitions, which the spin superdiffusion model is based upon, can induce a sharp spin moment reduction, and why electric current can change spin orientation in spintronics. These special k points, which are called spin Berry points [M. V. Berry, Proc. R. Soc. Lond. A 393 (1984) 45], are not random and appear when several bands are close to each other, so the Berry potential of spin majority states is different from that of spin minority states. Although within a single band, spin Berry points jump, when we group several neighboring bands together, they form distinctive smooth spin Berry lines. It is the band structure that disrupts those lines. Spin Berry points are crucial to laser-induced ultrafast demagnetization and spintronics.

Physics↗

Replica wormhole and information retrieval in the SYK model coupled to Majorana chains

Motivated by recent studies of the information paradox in (1+1)-D anti-de Sitter spacetime with a bath described by a (1+1)-D conformal field theory, we study the dynamics of second Ŕenyi entropy of the Sachdev-Ye-Kitaev (SYK) model (χ) coupled to a Majorana chain bath (ψ). The system is prepared in the thermofield double (TFD) state and then evolved by HL + HR. For small system-bath coupling, we find that the second Rényi entropy $S^(2)_{χL,χR}$ of the SYK model undergoes a first order transition during the evolution. In the sense of holographic duality, the long-time solution corresponds to a “replica wormhole”. The transition time corresponds to the Page time of a black hole coupled to a thermal bath. We further study the information scrambling and retrieval by introducing a classical control bit, which controls whether or not we add a perturbation in the SYK system. The mutual information between the bath and the control bit shows a positive jump at the Page time, indicating that the entanglement wedge of the bath includes an island in the holographic bulk.

1/N Expansion↗