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 109 records · Page 6

Hybrid Star Models in the Light of New Multimessenger Data

Abstract Recent astrophysical mass inferences of compact stars HESS J1731-347 and PSR J0952-0607, with extremely small and large masses respectively, as well as the measurement of the neutron skin of Ca in the CREX experiment challenge and constrain the models of dense matter. We examine the concept of hybrid stars—objects containing quark cores surrounded by nucleonic envelopes—as models that account for these new data along with other inferences. We employ a family of 81 nucleonic equations of state (EOSs) with variable skewness and slope of symmetry energy at saturation density and a constant speed-of-sound EOS for quark matter. For each nucleonic EOS, a family of hybrid EOSs is generated by varying the transition density, the energy jump, and the speed of sound. These models are tested against the data from GW170817 and J1731-347, which favor low-density soft EOS and J0592-0607 and J0740+6620, which require high-density stiff EOS. The addition of J0592-0607's mass measurement to the constraints has no significant impact on the parameter space of the admissible EOS, but allows us to explore the potential effect of pulsars more massive than J0740+6620, if such exists. We then examine the occurrence of twin configurations and quantify the ranges of masses and radii that they can possess. It is shown that including J1731-347 data favors EOSs that predict low-mass twins with M ≲ 1.3 M ⊙ that can be realized if the deconfinement transition density is low. If combined with large speed of sound in quark matter such models allow for maximum masses of hybrid stars in 2.0–2.6 M ⊙ .

Astronomy & Astrophysics↗

Modeling and Simulation of Rotary Sloshing in Launch Vehicles

The nonlinear dynamics of propellant sloshing during orbital ascent are usually neglected in the flight control analysis of large boost vehicles under the assumption that the viscous damping of the fluid is sufficient to suppress nonlinear phenomena and confine the fluid to small, planar free surface displacements. In this case, the sloshing dynamics can be modeled using a spring-mass-damper or linearized pendulum mechanical analog. However, large, smooth-wall tanks without significant internal hardware or ring baffles are still susceptible to nonlinear effects. In particular, rotary sloshing can present a risk to flight control as it involves the formation of a stable limit cycle which can lead to undesirable roll coupling. The underlying phenomena of jump resonance does not manifest in linear models, but can be reproduced using a nonlinear spherical pendulum or the Bauer paraboloid model developed during the Apollo/Saturn program. In this paper, a detailed analysis of the rotary sloshing dynamics of these mechanical analogs is presented, and discussed in the context of flight control stability. High-fidelity simulations of a representative boost vehicle are used to verify the semi-analytical predictions of the nonlinear dynamic response.

Modeling↗

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling of gap conductance for LWR fuel rods applied in the BISON code

A new gap conductance model is proposed in this study as a combination of Toptan’s model and the Ross-Stoute model. A variance-based sensitivity analysis is performed to understand how simulation results depend on all input parameters of the proposed model. Additionally, new modeling options (e.g. fill gas thermal conductivity, temperature jump distance, thermal accommodation coefficient, etc.) are added into the nuclear fuel performance code, BISON. The need for further investigation of the gap heat transfer between fuel and cladding in BISON motivated this study to evaluate its impact on the code’s predictions. New gap conductance modeling is proposed. A series of integral-effects validation tests is performed: (1) to demonstrate the impact of the proposed model on the code’s fuel temperature predictions at the beginning of life and through the reactor’s life; (2) to ensure that the proposed model is capable of accurately modeling gap heat transfer characteristics in real-world problems; and (3) to investigate the impact of the estimation of fission gas release on the fuel temperature predictions with the proposed model. Overall, the results indicate that the proposed gap conductance model improves BISON’s predictions.

42 ENGINEERING↗

Electrostatic potential jump across fast-mode collisionless shocks

The electrostatic potential jump across fast-mode collisionless shocks is examined by comparing published observations, hybrid simulations, and a simple model, in order to better characterize its dependence on the various shock parameters. In all three, it is assumed that the electrons can be described by an isotropic power-law equation of state. The observations show that the cross-shock potential jump correlates well with the shock strength but shows very little correlation with other shock parameters. Assuming that the electrons obey an isotropic power law equation of state, the correlation of the potential jump with the shock strength follows naturally from the increased shock compression and an apparent dependence of the power law exponent on the Mach number which the observations indicate. It is found that including a Mach number dependence for the power law exponent in the electron equation of state in the simple model produces a potential jump which better fits the observations. On the basis of the simulation results and theoretical estimates of the cross-shock potential, it is discussed how the cross-shock potential might be expected to depend on the other shock parameters.

Mandt, M. E.↗

Spatio-temporal electric field distributions in an atmospheric plasma jet impinging on a microchannel array surface

The electric field distribution in the ionization waves (IWs) propagating over a microchannel array dielectric surface, with the channels either empty or filled with distilled water, is measured by ps electric field induced second harmonic generation. The surface IW is initiated by the atmospheric pressure N 2 –Ar plasma jet impinging on the surface vertically and powered by ns pulse discharge bursts. The results show that the electric field inside the microchannels, specifically its horizontal component, is enhanced by up to a factor of 2. The field enhancement region is localized within the channels. The vertical electric field inside the channels lags in time compared to the field measured at the ridges, indicating the transient reversal of the IW propagation direction across the channels (toward the jet). This is consistent with the phase-locked plasma emission images and confirmed by the kinetic modeling predictions, which show that the IW ‘jumps’ over the empty channels and propagates into the channels only after the jump between the adjacent ridges. When the channels are filled with water, the wave speed increases by up to 50%, due to the higher effective dielectric constant of the surface. No evidence of a significant electric field enhancement near the dielectric surface (ceramic or water) has been detected, within the spatial resolution of the present diagnostic, ~100 μ m.

plasmas and liquids↗

The Hubble Space Telescope observations of X-ray nova Muscae 1991 and its spectral evolution

Hubble Space Telescope Faint Object Spectrograph (FOS) Faint Object Camera (FOC) and observations of Nova Muscae 1991 obtained on May 14-15, 1991 are presented and discussed. A nearly featureless continuum with broad 2200 A absorption feature and no Balmer jump is found. The FOS spectrum is modeled, including previous multiepoch IUE, optical, and X-ray data using a simple blackbody accretion disk model. It is found that the mass transfer rate decays exponentially with a characteristic time of about 43 days. The cooling front predicted by the disk instability model should have been observable, but was not seen.

Cheng, F. H.↗

Ultraviolet spectropolarimetry of the Be star PP Carinae with the Wisconsin Ultraviolet Photo-Polarimeter Experiment

We present the first ultraviolet spectropolarimetric observations of the Be star PP Car, obtained with the Wisconsin Ultraviolet Photo-Polarimeter Experiment (WUPPE) aboard the Astro 1 mission. Usable polarization data were obtained from 1400 to 2330 A, along with a good spectrum from 1400 to 3200 A. These data show a lower polarization shortward of the Balmer jump than had been predicted by standard models, and a broad UV polarization dip around 1900 A is seen. These results are in agreement with those found from the WUPPE observations of two other Be stars, Xi Tau and Pi Aqr, which were published earlier. All these observations are an important probe of the Be circumstellar envelopes and demonstrate the need for the inclusion of metal-line effects in circumstellar disk models of Be star UV polarization.

Bjorkman, K. S.↗

The Effect of Entrainment Through Atmospheric Boundary Layer Growth on Observed and Modeled Surface Ozone in the Colorado Front Range

Ozone concentrations at the Earth's surface are controlled by meteorological and chemical processes and are a function of advection, entrainment, deposition, and net chemical production/loss. The relative contributions of these processes vary in time and space. Understanding the relative importance of these processes controlling surface ozone concentrations is an essential component for designing effective regulatory strategies. Here we focus on the diurnal cycle of entrainment through atmospheric boundary layer (ABL) growth in the Colorado Front Range. Aircraft soundings and surface observations collected in July/August 2014 during the DISCOVER-AQ/FRAPPÉ (Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality/Front Range Air Pollution and Photochemistry Éxperiment) campaigns and equivalent data simulated by a regional chemical transport model are analyzed. Entrainment through ABL growth is most important in the early morning, fumigating the surface at a rate of ~5 ppbv/h. The fumigation effect weakens near noon and changes sign to become a small dilution effect in the afternoon on the order of -1 ppbv/h. The chemical transport model WRF-Chem (Weather Research and Forecasting Model with chemistry) underestimates ozone at all altitudes during this study on the order of 10-15 ppbv. The entrainment through ABL growth is overestimated by the model in the order of 0.6-0.8 ppbv/h. This results from differences in boundary layer growth in the morning and ozone concentration jump across the ABL top in the afternoon. This implicates stronger modeled fumigation in the morning and weaker modeled dilution after 11:00 LT.

Entrainment through ABL Growth↗

Mode Jumping of an Isogrid Panel Under Quasi-Static Compression

A wide column test of a composite isogrid panel subjected to quasi-static, axial compression is modeled with a hybrid-static dynamic computational method. The data from the test panel exhibited discontinuous responses in the compressive load for slowly increased end-shortening. The computational model was developed to corroborate these discontinuities with the phenomenon of mode jumping. Mode jumping refers to the transient response of the panel from an unstable bifurcation point on a postbuckled equilibrium path to a second stable equilibrium state on a new equilibrium path. On the new equilibrium path, both the analysis and test show that the panel can resist increased endshortening beyond that of the unstable critical point. Fair agreement is achieved between the analysis and test.

Muheim, Danniella M.↗

HybridSystemsSimulations.jl - Solving the Merchant Collocated Facilities with JuMP

The development of new clean-generation technologies also leads to new plant-level architectures that combine several generation and storage assets behind the point of connection. These co-located generation resources (Hybrid Systems) primarily operate as merchant assets that employ automated market bidding models and internal Energy Management Systems (EMS) to comply with the operator's signals. Formulating an optimal bidding model requires embedding the EMS control model into the bidding algorithm, resulting in a bi-level optimization problem. In this presentation, we first showcase using JuMP to formulate and solve this problem effectively for multiple merchant systems and the bidding outcomes considering different model formulations. Second, the bidding outcomes are later integrated into a PowerSimulations.jl (also built with JuMP) simulation to study the system-level effects of the various merchant bidding and the interactions between market-clearing models and the embedded EMS model. We will showcase simulations conducted in the RTS system considering different levels of merchant hybrid systems participation. The presentation provides the following specific insights on JuMP usage: 1) the Formulation of specialized bi-level problems with custom cuts to solve the merchant hybrid system bidding problem; 2) the integration of a modular model within a complex simulation workflow supported by JuMP in PowerSimulations.jl; 3) Accelerating the solution of power systems operations simulation that employ agent optimization problems using JuMP.

energy markets↗

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗

Origin of gaplike behaviors in URu 2 Si 2 : Combined study via quasiparticle scattering spectroscopy and resistivity measurements

In this work, we address two long-standing questions regarding the hidden order in URu 2 Si 2 : Is it associated with the hybridization process, and what are the distinct roles played by the localized and itinerant electrons? Our quasiparticle scattering spectroscopy reveals a hybridization gap ubiquitous in the entire phase space spanned by P and Fe substitutions in URu 2 Si 2 , including the no-order and antiferromagnetic regions, with minimal change upon crossing the phase boundary. This indicates its opening is not associated with the ordering, and thus localized electrons must be the major player. Towards a consistent understanding of all the other gaplike behaviors observed only below transition temperatures, we analyze the electrical resistivity using a model in which gapped bosonic excitations are the dominant scattering source. With their stiffness set to follow an unusual temperature dependence (decreasing with decreasing temperature), this model fits all of our resistivity data well including the jump at the transition. Remarkably, the extracted gap increases slowly with increasing Fe content, similarly to the gap detected by inelastic neutron scattering at $Q_1$ = (1.4, 0, 0), suggesting a common origin. Such a model can also naturally explain the Hall effect temperature dependence without invoking Fermi surface gapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiquanta flux jumps in superconducting fractal

We study the magnetic field response of millimeter scale fractal Sierpinski gaskets (SG) assembled of superconducting equilateral triangular patches. Directly imaged quantitative induction maps reveal hierarchical periodic filling of enclosed void areas with multiquanta magnetic flux, which jumps inside the voids in repeating bundles of individual flux quanta Φ 0 . The number N s of entering flux quanta in different triangular voids of the SG is proportional to the linear size s of the void, while the field periodicity of flux jumps varies as 1/s. We explain this behavior by modeling the triangular voids in the SG with effective superconducting rings and by calculating their response following the London analysis of persistent currents, J s , induced by the applied field H a and by the entering flux. With changing H a , J s reaches a critical value in the vertex joints that connect the triangular superconducting patches and allows the giant flux jumps into the SG voids through phase slips or multiple Abrikosov vortex transfer across the vertices. The unique flux behavior in superconducting SG patterns, may be used to design tunable low-loss resonators with multi-line high-frequency spectrum for microwave technologies.

36 MATERIALS SCIENCE↗

Secular rotational motions and the mechanical structure of a dynamical viscoelastic earth

A survey is presented of analytical methods for computing the linear responses of the rotational axis of a layered viscoelastic earth to surface loading. Theoretical research in this area is first summarized, and the differences between the mechanical boundary conditions to be applied at the interface separating the upper and lower mantles for an adiabatically and chemically stratified mantle are discussed. Some examples of polar wander and secular variation of the spin rate from glacial excitation are presented for various types of chemical and viscosity stratifications. The effects of an artificial density jump at the base of the lithosphere in models are examined, and certain issues concerning the fluid tidal Love number for different types of density stratification are addressed. The meaning of effective plate thickness over geological time scales for rotational dynamics is discussed.

Yuen, D. A.↗

Towards a Theory of Sampled-Data Piecewise-Deterministic Markov Processes

The analysis and design of practical control systems requires that stochastic models be employed. Analysis and design tools have been developed, for example, for Markovian jump linear continuous and discrete-time systems, piecewise-deterministic processes (PDP's), and general stochastic hybrid systems (GSHS's). These model classes have been used in many applications, including fault tolerant control and networked control systems. This paper presents initial results on the analysis of a sampled-data PDP representation of a nonlinear sampled-data system with a jump linear controller. In particular, it is shown that the state of the sampled-data PDP satisfies the strong Markov property. In addition, a relation between the invariant measures of a sampled-data system driven by a stochastic process and its associated discrete-time representation are presented. As an application, when the plant is linear with no external input, a sufficient testable condition for the convergence in distribution to the invariant delta Dirac measure is given.

Herencia-Zapana, Heber↗

Range image segmentation using Zernike moment-based generalized edge detector

The authors proposed a novel Zernike moment-based generalized step edge detection method which can be used for segmenting range and intensity images. A generalized step edge detector is developed to identify different kinds of edges in range images. These edge maps are thinned and linked to provide final segmentation. A generalized edge is modeled in terms of five parameters: orientation, two slopes, one step jump at the location of the edge, and the background gray level. Two complex and two real Zernike moment-based masks are required to determine all these parameters of the edge model. Theoretical noise analysis is performed to show that these operators are quite noise tolerant. Experimental results are included to demonstrate edge-based segmentation technique.

Ghosal, S.↗

A guide to the use of the pressure disk rotor model as implemented in INS3D-UP

This is a guide for the use of the pressure disk rotor model that has been placed in the incompressible Navier-Stokes code INS3D-UP. The pressure disk rotor model approximates a helicopter rotor or propeller in a time averaged manner and is intended to simulate the effect of a rotor in forward flight on the fuselage or the effect of a propeller on other aerodynamic components. The model uses a modified actuator disk that allows the pressure jump across the disk to vary with radius and azimuth. The cyclic and collective blade pitch angles needed to achieve a specified thrust coefficient and zero moment about the hub are predicted. The method has been validated with experimentally measured mean induced inflow velocities as well as surface pressures on a generic fuselage. Overset grids, sometimes referred to as Chimera grids, are used to simplify the grid generation process. The pressure disk model is applied to a cylindrical grid which is embedded in the grid or grids used for the rest of the configuration. This document will outline the development of the method, and present input and results for a sample case.

Chaffin, Mark S.↗