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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1936

Tank Tests of a Model of the NC Flying-boat Hull - N.A.C.A. Model 44

A 1/7.06 full-size model of the NC-type hull was tested in the N.A.C.A. tank by both the general method and the specific or free-to-trim method. The results of the tests are given in curves plotted as non dimensional coefficients and are compared with the test results of N.A.C.A. model 11-A. The NC model (N.A.C.A. model 44) shows higher resistance than model 11-A at hump speed but lower resistance at high speeds. Model 44 has a higher best trim angle at the jump and a lower maximum positive trimming moment than model 11-A. At high speeds the best trim angle and the trimming moments of the two models are approximately the same.

Bell, Joe W↗
DOE OSTI2025

STREAM: A technology planning and capacity expansion model for the industrial sector

The Strategic Technology Roadmapping and Energy, Environmental, and Economic Analysis Model—STREAM—is an optimization-based modeling tool and analysis framework to assist with strategic planning and technology investments of the industrial sector. This open-source framework is written in Julia using the JuMP package, which enables users to model future “pathways” for incumbent and future production technologies, costs, fuels and energy carriers, and energy and non-energy environmental impacts from industries as they transform in pursuit of a robust and competitive manufacturing sector. The model starts with an initial stock of industrial production technologies and assets at a facility level and then determines pathways that minimize cost, subject to an array of possible constraints on demand, market shares, environmental flows, and other exogenously specified operational considerations such as capacity utilization rates or regional energy costs. Key features of the framework include flexibility to model a wide range of industries and industrial technologies/processes at varying levels of granularity, ability to perform parametric sensitivity analyses, and ability to visualize model results using visualization objects.

capacity expansion↗
DOE OSTI2022

MPC4CLR (Model-Predictive-Control-for-Critical-Load-Restoration-in-Power-Distribution-Systems) [SWR-22-24]

Model predictive control (MPC) is a system or process control technique for making decisions under uncertainty via rolling look-ahead optimizations at each control step where only the current step decisions are applied, and the rest are discarded. In this work, we developed an MPC for a critical load restoration (CLR) in power distribution systems to recover system service (electricity delivery) following an extreme event-triggered substation outage. The method considers the problem of controlling distributed energy resources (DERs) of the distribution system with the objective of achieving maximum load pick up while satisfying distribution network flow and voltage constraints. A linearized optimal power flow (OPF) model is employed to represent the physics of the network. The problem formulation is augmented with a ramping (up) reserve product for the DERs to ensure improved and upward monotonic load restoration as time evolves. Simulation analysis and performance tests are performed using a modified IEEE 13-bus test feeder integrated with wind, solar, microturbine, and energy storage battery. The software is developed using various software packages in Julia and Python. The MPC model is implemented using the JuMP optimization language in Julia while the data analytics including renewable generation and load demand forecasts, running the MPC simulation and visualizations is performed in Python.

Eseye, Abinet Tesfaye↗
NASA NTRS1976

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗
NASA NTRS1990

Panel methods: An introduction

Panel methods are numerical schemes for solving (the Prandtl-Glauert equation) for linear, inviscid, irrotational flow about aircraft flying at subsonic or supersonic speeds. The tools at the panel-method user's disposal are (1) surface panels of source-doublet-vorticity distributions that can represent nearly arbitrary geometry, and (2) extremely versatile boundary condition capabilities that can frequently be used for creative modeling. Panel-method capabilities and limitations, basic concepts common to all panel-method codes, different choices that were made in the implementation of these concepts into working computer programs, and various modeling techniques involving boundary conditions, jump properties, and trailing wakes are discussed. An approach for extending the method to nonlinear transonic flow is also presented. Three appendices supplement the main test. In appendix 1, additional detail is provided on how the basic concepts are implemented into a specific computer program (PANAIR). In appendix 2, it is shown how to evaluate analytically the fundamental surface integral that arises in the expressions for influence-coefficients, and evaluate its jump property. In appendix 3, a simple example is used to illustrate the so-called finite part of the improper integrals.

Erickson, Larry L.↗
NASA NTRS2003

Closed-Loop Performance Measures for Flight Controllers Subject to Neutron-Induced Upsets

It has been observed that atmospheric neutrons can produce single event upsets in digital flight control hardware. The phenomenon has been studied extensively at the chip level, and now system level experiments are underway. In this paper analytical closed-loop performance measures for the tracking error are developed for a plant that is stabilized by a recoverable computer system subject to neutron induced upsets. The underlying model is a Markov jump-linear system with process noise. The steady-state tracking error is expressed in terms of a generalized observability Gramian.

Gray, W. Steven↗
DOE OSTI2025

Turbulent Processes that Influence Boundary-Layer Cloud Structure (Final Technical Report)

Many models of the atmosphere simulate vertical profiles of observable fields, such as temperature, moisture, and vertical velocity. To assess the realism of such model output, it is important to evaluate the modeled profiles against observations. However, many extant evaluations consist of comparing time averaged profiles from the model and the observations. Such time averages can smooth out unrealistic features that appear in instantaneous profiles. In this way, the time averages may obscure certain model errors from detection. For instance, many models misrepresent turbulent processes, leading to profiles that are too smoothed out or too noisy. E.g., if a model has coarse vertical resolution and overly strong diffusion, then the model may be unable to produce observed sharp peaks. On the other hand, a model may have artificial approximations that produce spurious jumps in the vertical. Much can be learned by browsing through instantaneous snapshots of model output. But the volume of output produced by a global model is far too great to rely purely on visual inspection. In addition, we need a way to statistically characterize the profiles, so that anomalies can be detected more conveniently.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Heat-pipe transient model for space applications

A two-dimensional model is developed for simulating heat pipes transient performance following changes in the input/rejection power or in the evaporator/condenser temperatures. The model employs the complete form of governing equations and momentum and energy jump conditions at the liquid-vapor interface. Although the model is capable of handling both cylindrical and rectangular geometries, the results reported are for a circular heat pipe with liquid lithium as the working fluid. The model incorporates a variety of other working fluids, such as water, ammonia, potassium, sodium and mercury, and offers combinations of isothermal, isoflux, convective and radiative heating/cooling conditions in the evaporator and condenser regions of the heat pipe. Results presented are for lithium heat pipes with exponential heating of the evaporator and isothermal cooling of the condenser.

Tournier, Jean-Michel↗
NASA NTRS1988

Calculation of unsteady Euler flows in turbomachinery using the linearized Euler equations

The model used by the present unsteady cascade flow calculation method is based on the linearized unsteady Euler equations and accounts for blade loading and geometry, shock motion, and wake motion. The steady flow that must be determined prior to the unsteady flow is ascertained by means of a Newton iteration procedure. A noteworthy feature of the procedure is the use of shock-fitting to determine steady and unsteady shock positions; the use of the Euler equations in conjunction with the Rankine-Hugoniot shock-jump conditions is found to directly model the generation of entropy and vorticity at shocks. Results are presented for both channel and cascade flows.

Hall, Kenneth C.↗
NASA NTRS1990

Linear evolution of current sheets in sheared force-free magnetic fields with discontinuous connectivity

Thin current sheets arising in tenuous, magnetized solar coronal plasmas may constitute an important mechanism for energy buildups and subsequent energy releases; they could arise from the continuous-and-random motion of magnetic footprints associated with photospheric velocity fields. A model is presented for study of the quasi-static evolution of current sheets due to shearing of the footpoints, in a highly idealized geometry that incorporates an abrupt jump in field-line connectivity. The model highlights that formation of thin current layers and allows large shearing motions prior to violation of the linear approximation. Excess energy comparable to that released by solar flares can be stored in a sheared field.

Wolfson, Richard↗
DOE OSTI2020

Probabilistic analysis of masked loads with aggregated photovoltaic production

In this paper we present a probabilistic analysis framework to estimate behind-the-meter photovoltaic generation in real time. We develop a forward model consisting of a spatiotemporal stochastic process that represents the photovoltaic generation and a stochastic differential equation with jumps that represents the demand. Here, we employ this model to disaggregate the behind-the-meter photovoltaic generation using net load and irradiance measurements.

14 SOLAR ENERGY↗
NASA NTRS1990

Interaction of minor ions with fast and slow shocks

The coronal slow shock was predicted to exist embedded in large coronal holes at 4 to 10 solar radii. A three-fluid model was used to study the jumps in minor ions propertes across the coronal slow shock. The jump conditions were formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. A simularity solution was also obtained for heavy ions. The results show that on the downstream side of the coronal slow shock the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.

Whang, Y. C.↗
NASA NTRS1990

Heating of minor ions by the coronal slow shock

The coronal slow shock has been predicted to exist embedded in large coronal holes at 4-10 solar radii. In this paper, a three-fluid model is used to study the jumps in minor ion properties across a slow shock such as the coronal slow shock. The jump conditions are formulated in the de Hoffmann-Teller frame of reference. The Rankine-Hugoniot solution determines the MHD flow and the magnetic field across the shocks. For each minor ion species, the fluid equations for the conservation of mass, momentum, and energy can be solved to determine the velocity and the temperature of the ions across the shock. Also obtained is a similarity solution for heavy ions. The results show that, on the downstream side of the slow shock, the ion temperatures are nearly proportional to the ion masses for He, O, Si, and Fe, in agreement with observed ion temperatures in the inner solar wind. This indicates that the possibly existing coronal slow shock can be responsible for the observed heating of minor ions in the solar wind.

Whang, Y. C.↗
DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗
NASA NTRS1988

Muon motion in titanium hydride

Motional narrowing of the transverse-field muon spin rotation signal was observed in gamma-TiH(x) for x = 1.83, 1.97, and 1.99. An analysis of the data for TiH1.99 near room temperature indicates that the mechanism responsible for the motion of the muon out of the octahedral site is thermally activated diffusion with an attempt frequency comparable to the optical vibrations of the lattice. Monte Carlo calculations to simulate the effect of muon and proton motion upon the muon field-correlation time were used to interpret the motional narrowing in TiH1.97 near 500 K. The interpretation is dependent upon whether the Bloembergen, Purcell, and Pound (BPP) theory or an independent spin-pair relaxation model is used to obtain the vacancy jump rate from proton NMR T1 measurements. Use of BPP theory shows that the field-correction time can be obtained if the rate of motion of the muon with respect to the rate of the motion for the protons is decreased. An independent spin-pair relaxation model indicates that the field-correlation time can be obtained if the rate of motion for the nearest-neighbor protons is decreased.

Kempton, J. R.↗
NASA NTRS1999

Condensation of Forced Convection Two-Phase Flow in a Miniature Tube

A physical/mathematical model of annular film condensation at the inlet of a miniature tube has been developed. In the model, the liquid flow is coupled with the vapor flow along the liquid-vapor interface through the interfacial temperature, heat flux, shear stress, and pressure jump conditions due to surface tension effects. The model predicts the shape of the liquid-vapor interface along the condenser and leads to the conclusion that there is complete condensation at a certain distance from the condenser inlet. The numerical results show that complete condensation of the incoming vapor is possible at comparatively low heat loads and that this is a special case of a more general condensation regime with two-phase bubbly flow downstream of the initial annular film condensation region. Observations from the flow visualization experiment confirm the existence and qualitative features of annular film condensation leading to the complete condensation phenomenon in a small diameter (3.25 mm) circular tube condenser.

Begg, E.↗
DOE OSTI2022

The Paschen Jump as a Diagnostic of the Diffuse Nebular Continuum Emission in Active Galactic Nuclei*

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 Å is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope Space Telescope Imaging Spectrograph spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multicomponent spectral models over the range 6800–9700 Å and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from ~10% to 50% at 8000 Å, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe ii emission, high-order Balmer line emission, and other components.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Computation of Direct Sensitivities of Spatial Multibody Systems With Joint Friction

Abstract Friction exists in most mechanical systems and may have a major influence on the dynamic performance of the system. The incorporation of friction in dynamic systems has been a subject of active research for several years owing to its high nonlinearity and its dependence on several parameters. Consequently, optimization of dynamic systems with friction becomes a challenging task. Gradient-based optimization of dynamical systems is a prominent technique for optimal design and requires the computation of model sensitivities with respect to the design parameters. The novel contribution of this paper is the derivation of the analytical methodology for the computation of direct sensitivities for smooth multibody systems with joint friction using the Lagrangian index-1 formulation. System dynamics have been computed using two different friction models; the Brown and McPhee, and the Gonthier et al. model. The methodology proposed to obtain model sensitivities has also been validated using the complex finite difference method. A case study has been conducted on a spatial multibody system to observe the effect of friction on the dynamics and model sensitivities, compare sensitivities with respect to different parameters and demonstrate the numerical and validation aspects. Since design parameters can have very different magnitudes and units, the sensitivities have been scaled with the parameters for comparison. Finally, a discussion has been presented on the interpretation of the case study results. Due to the incorporation of joint friction, ‘jumps’ or discontinuities are observed in the model sensitivities akin to those observed for hybrid dynamical systems.

Engineering↗