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 289 records · Page 16

Nucleation of Grain Boundary Phases

Here, we derive a theory that describes homogeneous nucleation of grain boundary (GB) phases. Our analysis takes account of the energy resulting from the GB phase junction, the line defect separating two different GB structures, which is necessarily a dislocation as well as an elastic line force due to the jump in GB stresses. The theory provides analytic forms for the elastic interactions and the core energy of the GB phase junction that, along with the change in GB energy, determines the nucleation barrier. We apply the resulting nucleation model to simulations of GB phase transformations in tungsten. Our theory explains why under certain conditions GBs cannot spontaneously change their structure even to a lower energy state.

36 MATERIALS SCIENCE↗

Predicted continuum spectra of type II supernovae - LTE results

The continuum spectral energy distribution of the flux emerging from type II supernovae is calculated from quasi-static radiative transfer through a power-law density gradient, assuming radiative equilibrium and LTE. It is found that the Balmer jump disappears at high effective temperatures and low densities, while the spectrum resembles that of a dilute blackbody but is flatter with a sharper cutoff at the short-wavelength end. A significant UV excess is found in all models calculated. The calculation should be considered exploratory because of significant effects which are anticipated to arise from departure from LTE.

Shaviv, G.↗

MHD intermediate shocks in coronal mass ejections

A simplified model of coronal mass ejections is considered in which at least a portion of the interaction with the background corona involves a shock wave, and the allowable shock solutions and their compressive signatures are examined. The MHD shock-jump equations have a maximum of three possible types of solutions with an entropy rise for fixed values of the physical variables (slow, intermediate, and fast shocks). However, one of the three solution classes (the intermediate shock) is widely believed to not occur in nature and is regarded as nonevolutionary or extraneous. Without the intermediate shock, there is no multiplicity of solutions in that only one shock (or none) can occur for given physical values. All three potential shock types are considered, and it is shown solely on the basis of the shock-jump equations, that intermediate shocks must exist along some segment of the shock front for certain parametric regimes and for conditions that probably occur in some coronal mass ejections.

Steinolfson, R. S.↗

Model atmospheres for the early stages of novae in outburst

For the photospheres of classical novae during the early stages of their outbursts, continuum and line blanketed models are presented. The expanding envelopes are characterized by small density gradients that lead to very large geometrical extensions and large temperature differences between inner and outer parts. The spectra show large excesses in the IR and small Balmer jumps which may either be in absorption or emission. For the parameters considered most lines are in absorption if powerlaw density distributions and radiative equilibrium are assumed. The effects of modifications in the temperature structure (e.g., heating from shock fronts) and in the density stratifications as well as changes in the abundances are discussed.

Wehrse, R.↗

Weak shock propagation through a turbulent atmosphere

Consideration is given to the propagation through turbulence of transient pressure waveforms whose initial onset at any given point is an abrupt shock. The work is motivated by the desire to eventually develop a mathematical model for predicting statistical features, such as peak overpressures and spike widths, of sonic booms generated by supersonic aircraft. It is argued that the transient waveform received at points where x greater than 0 will begin with a pressure jump and a formulation is developed for predicting the amount of this jump and the time derivatives of the pressure waveform immediately following the jump.

Pierce, Allan D.↗

On the Stability of Jump-Linear Systems Driven by Finite-State Machines with Markovian Inputs

This paper presents two mean-square stability tests for a jump-linear system driven by a finite-state machine with a first-order Markovian input process. The first test is based on conventional Markov jump-linear theory and avoids the use of any higher-order statistics. The second test is developed directly using the higher-order statistics of the machine s output process. The two approaches are illustrated with a simple model for a recoverable computer control system.

Patilkulkarni, Sudarshan↗

First lithium experiments in HIDRA and evidence of helium retention during quasi-steady-state stellarator plasma operations

Recent experiments in Hybrid Illinois Device for Research and Applications (HIDRA) have had operational discharges between t discharge = 60 and 1000 s using electron cyclotron resonant heating (ECRH) of the plasma. This means that quasi-steady-state plasma discharges reach conditions to study long-pulse plasma material interactions (PMIs). The newly commissioned HIDRA-Material Analysis Test-stand PMI diagnostic is used to place a drop of lithium onto a heated tungsten surface, transfer the sample in-vacuo and expose it in a helium plasma. Helium is of interest as there is an open question to whether lithium will be able to remove helium ash in real fusion devices. The introduction of the W-Li sample in HIDRA resulted in evaporation of lithium into the helium plasma during a 600 s pulse and caused a reduction of over 90% in neutral pressure during the discharge. It was also observed that the plasma density and temperature increased by over 2.5 times. Using spectroscopy and a helium collisional radiative model, the peak temperature and density of the helium plasma can be monitored during the discharge. During lithium evaporation, as significant lithium ionization occurs, there is a 85% drop in the HIDRA vessel neutral pressure, despite a constant flow rate of He gas. This reduction in neutral pressure is supported by spectroscopy data with corresponding reductions in He I line intensities (587 nm, 667 nm, 706 nm, and 728 nm), as well as those of other impurities. At one point in the discharge a lithium plasma is created, as indicated by an increase in Li + emission and a complete reduction in He + emission, but the electron density jumps from ne = 3 × 10 18 m –3 to over ne = 8 × 10 18 m –3 while the core temperature stays relatively constant between Te = 16 eV and 20 eV. Once lithium has completely evaporated from the sample and the majority of the ionized lithium has diffused from the plasma to the vessel walls, pressure and spectroscopy data paired with He collisional radiative model calculations shows a re-establishment of a helium plasma in a low recycling regime. In this regime, the density drops down to ne = 2 × 10 18 m –3 and the electron temperature increases from T e = 20 eV to over T e = 50 eV indicating an increase in helium heating efficiency. This is also indicated by the He+ emission re-establishing and having a higher intensity. Here, we show the results from the first lithium campaign in HIDRA. In the presence of lithium, and in particular when lithium ions are present, the helium disappears from the plasma via an as of yet unknown complex relationship that needs to be further studied. The most likely explanation is that the lithium ions are distributed around the vessel and able to trap helium to the surface turning HIDRA into a large gettering surface. These results have potential implications on future plasma facing component design using liquid lithium for impurity and recycling control using limiters and divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rotational Dynamics and Transition Mechanisms of Surface-Adsorbed Proteins

Assembly of biomolecules at solid-water interfaces requires molecules to traverse complex orientation-dependent energy landscapes through processes that are poorly understood, largely due to the dearth of in-situ single molecule measurements and statistical analyses of the rotational dynamics that define directional selection. Emerging capabilities in high-speed atomic force microscopy and machine learning have allowed us to directly determine the orientational energy landscape and observe and quantify the rotational dynamics for protein nanorods on the surface of muscovite mica under a variety of conditions. Comparisons with kinetic Monte Carlo simulations show that the transition rates between adjacent orientation-specific energetic minima can largely be understood through traditional models of in-plane Brownian rotation across a biased energy landscape, with resulting transition rates that are exponential in the energy-barriers between states. However, transitions between more distant angular states are decoupled from barrier height, with jump-size distributions showing a power-law decay that is characteristic of a non-classical Levy-flight random walk, indicating that large jumps are enabled by alternative modes of motion via activated states. The findings provide new insights into the dynamics of biomolecules at solid-liquid interfaces that lead to self-assembly, epitaxial matching and other orientationally anisotropic outcomes and define a general procedure for exploring such dynamics with implications for hybrid biomolecular-inorganic materials design.

Orientational energy landscapes, Rotational dynami↗

Observation of Electromagnetic Transients in a Nb3Sn 4-layer Dipole Mirror Magnet

During testing of the Stress Managed Cosine-Theta dipole mirror magnet SMCTM1, magnet quenches were observed following large voltage spikes in the half-coil voltage taps which preceded normal quench initiation. Following some recent work at CERN measuring power current converter transients in relation to flux jumps, an additional differential current probe was added to the magnet instrumentation. Measurements of this probe, as well as the half coil voltage taps show a fast and substantial current change preceding the quench, hinting at coil motion. A simple model based on relative shift of the inner and outer layer coils was developed which suggests the magnitude of coil shift events is on the estimated relative coil shift from the magnet constraint conditions. The behavior is shown to be reversible and characteristic of stick-slip.

Krave, Steven [Fermilab] (ORCID:0000000217521457)↗

Software For Design Of Life-Support Systems

Design Assistant Workstation (DAWN) computer program is prototype of expert software system for analysis and design of regenerative, physical/chemical life-support systems that revitalize air, reclaim water, produce food, and treat waste. Incorporates both conventional software for quantitative mathematical modeling of physical, chemical, and biological processes and expert system offering user stored knowledge about materials and processes. Constructs task tree as it leads user through simulated process, offers alternatives, and indicates where alternative not feasible. Also enables user to jump from one design level to another.

Rudokas, Mary R.↗

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗

The Prelude to and Aftermath of the Giant Flare of 2004 December 27: Persistent and Pulsed X-ray Properties of SGR 1806-20 from 1993 to 2005

We report on the evolution of key spectral and temporal parameters of SGR 1806-20 prior to and following the highly energetic giant flare of 2004 December 27. Using RXTE, we track the pulse frequency of the SGR and find that the spin-down rate varied erratically in the months before and after the flare. Contrary to the giant flare in SGR 1900+14, we find no evidence for a discrete jump in spin frequency at the time of the December 27th flare (Absolute value of (Delta * upsilon/upsilon) < 5 10 (exp -6)). In the months surrounding the flare, we find a strong correlation between pulsed flux and torque consistent with the model for magnetar magnetosphere electrodynamics proposed. As with the flare in SGR 1900+14, the pulse morphology of SGR 1806-20 changes drastically following the flare. Using Chandra and other publicly available imaging X-ray detector observations, we construct a spectral history of SGR 1806-20 from 1993 to 2005. The usual magnetar persistent emission spectral model of a power law plus a blackbody provides an excellent fit to the data. We confirm the earlier finding of increasing spectral hardness of SGR 1806-20 between 1993 and 2004. However, our results indicate significant differences in the temporal evolution of the spectral hardening. Rather than a direct correlation between torque and spectral hardness, we find evidence for a sudden torque change that preceded a gradual hardening of the energy spectrum on a timescale of years. Interestingly, the spectral hardness, spin-down rate, phase-averaged flux, and pulsed flux of SGR 1806-20 all peak months before the flare epoch.

Woods, Peter M.↗

Radar-anomalous, high-altitude features on Venus

Over nearly all of the surface of Venus the reflectivity and emissivity at centimeter wavelengths are about 0.15 and 0.85 respectively. These values are consistent with moderately dense soils and rock populations, but the mean reflectivity is about a factor of 2 greater than that for the Moon and other terrestrial planets. Pettingill and Ford, using Pioneer Venus reflectivities and emissivities, found a number of anomalous features on Venus that showed much higher reflectivities and much lower emissivities with both values approaching 0.5. These include Maxwell Montes, a number of high regions in Aphrodite Terra and Beta Regio, and several isolated mountain peaks. Most of the features are at altitudes above the mean radius by 2 to 3 km or more. However, such features have been found in the Magellan data at low altitudes and the anomalies do not exist on all high structures, Maat Mons being the most outstanding example. A number of papers have been written that attempt to explain the phenomena in terms of the geochemistry balance of weathering effects on likely surface minerals. The geochemists have shown that the fundamentally basaltic surface would be stable at the temperatures and pressures of the mean radius in the form of magnetite, but would evolve to pyrite and/or pyrrhotite in the presence of sulfur-bearing compounds such as SO2. Pyrite will be stable at altitudes above 4 or 5 km on Venus. Although the geochemical arguments are rather compelling, it is vitally important to rationally look at other explanations for radar and radio emission measurements such as that presented by Tryka and Muhleman. The radar reflectivity values are retrieved from the raw Magellan backscatter measurements by fitting the Hagfors' radar scattering model in which a surface roughness parameters and a normal incidence electrical reflectivity are estimated. The assumptions of the theory behind the model must be considered carefully before the results can be believed. These include that the surface roughness exists only at horizontal scales large compared to the wavelength, the vertical deviations are gaussianly distributed, there is no shadowing, and that the reflection occurs at the interface of two homogeneous dielectric half-spaces. Probably all these conditions are violated at the anomalous features under discussion. The most important of these is the homogeneity of the near surface of Venus, particularly in highlands. Under the assumptions of the theory, all of the radio energy is reflected by the impedance jump at the very boundary. However, in heterogeneous soil some fraction of the illuminating energy is propagated into the soil and then scattered back out by impedance discontinuities such as rock, voids, and cracks. In light soils, the latter effect can overwhelm the scattering effects of the true surface and greatly enhance the backscatter power, suggesting a much higher value of an effective dielectric constant that would be estimated from Hagfors' model.

Muhleman, Duane O.↗

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗

Quantum algorithm to simulate Lindblad master equations

We present a quantum algorithm for simulating a family of Markovian master equations that can be realized through a probabilistic application of unitary channels and state preparation. Our approach employs a second-order product formula for the Lindblad master equation, achieved by decomposing the dynamics into dissipative and Hamiltonian components and replacing the dissipative segments with randomly compiled, easily implementable elements. The sampling approach eliminates the need for ancillary qubits to simulate the dissipation process and reduces the gate complexity in terms of the number of jump operators. We provide a rigorous performance analysis of the algorithm. We also extend the algorithm to time-dependent Lindblad equations, generalize the family of Markovian master equations it can be applied to, and explore applications beyond the Markovian noise model. A new error bound, in terms of the diamond norm, for second-order product formulas for time-dependent Liouvillians is provided that might be of independent interest. Published by the American Physical Society 2025

Borras, Evan (ORCID:000900017709037X)↗

Tracker Fits for shot 1s1674

Gas gun shock-to-detonation experiments provide data for calibrating HE burn models. At LANL the principal diagnostic is a magnetic gauge package that provides data on the lead shock trajectory, which is used to determine the run-distance to detonation, and the lead shock speed with run distance, which together with jump-off velocity of the velocity gauges determines Hugoniot locus points for the reactants.

42 ENGINEERING↗

A numerical solution of the axisymmetric jet counterflow problem

A numerical computation of a normal shock case in which the shock waves and shear layers are modeled as normal and tangential discontinuities bounding the regions of an inviscid flow is discussed. Within these regions the unsteady gasdynamic equations are solved in finite-volume form on a moving computational mesh. The positions of discontinuities are determined from the jump conditions. The results obtained for values of the ratio of jet total pressure to free-stream total pressure ranging from 10.0 to 100.0 are found to be in agreement with the experimental results obtained by Fleeman and Nelson (1974).

Schiff, L. B.↗

A vector-based failure detection and isolation algorithm for a dual fail-operational redundant strapdown inertial measurement unit

A vector-based failure detection and isolation technique for a skewed array of two degree-of-freedom inertial sensors is developed. Failure detection is based on comparison of parity equations with a threshold, and isolation is based on comparison of logic variables which are keyed to pass/fail results of the parity test. A multi-level approach to failure detection is used to ensure adequate coverage for the flight control, display, and navigation avionics functions. Sensor error models are introduced to expose the susceptibility of the parity equations to sensor errors and physical separation effects. The algorithm is evaluated in a simulation of a commercial transport operating in a range of light to severe turbulence environments. A bias-jump failure level of 0.2 deg/hr was detected and isolated properly in the light and moderate turbulence environments, but not detected in the extreme turbulence environment. An accelerometer bias-jump failure level of 1.5 milli-g was detected over all turbulence environments. For both types of inertial sensor, hard-over, and null type failures were detected in all environments without incident. The algorithm functioned without false alarm or isolation over all turbulence environments for the runs tested.

Morrell, Frederick R.↗