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 253 records · Page 14

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: The Two Dimensional Time Marching Schemes - 1

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Chang, Sin-Chung↗

Standing shocks in a two-fluid solar wind

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Habbal, Shadia R.↗

An Integrated 0-1 Hour First-Flash Lightning Nowcasting, Lightning Amount and Lightning Jump Warning Capability

Using satellite-based methods that provide accurate 0-1 hour convective initiation (CI) nowcasts, and rely on proven success coupling satellite and radar fields in the Corridor Integrated Weather System (CIWS; operated and developed at MIT-Lincoln Laboratory), to subsequently monitor for first-flash lightning initiation (LI) and later period lightning trends as storms evolve. Enhance IR-based methods within the GOES-R CI Algorithm (that must meet specific thresholds for a given cumulus cloud before the cloud is considered to have an increased likelihood of producing lightning next 90 min) that forecast LI. Integrate GOES-R CI and LI fields with radar thresholds (e.g., first greater than or equal to 40 dBZ echo at the -10 C altitude) and NWP model data within the WDSS-II system for LI-events from new convective storms. Track ongoing lightning using Lightning Mapping Array (LMA) and pseudo-Geostationary Lightning Mapper (GLM) data to assess per-storm lightning trends (e.g., as tied to lightning jumps) and outline threat regions. Evaluate the ability to produce LI nowcasts through a "lightning threat" product, and obtain feedback from National Weather Service forecasters on its value as a decision support tool.

remote sensing↗

Non-conformal interface-cohesive modeling with the shifted boundary method

The accurate simulation of boundary- and interface-dominated problems on complex geometries remains challenging when boundary- or interface-fitted meshes are difficult to generate, particularly for curved boundaries, polycrystalline microstructures, and dense interface networks. The Shifted Boundary Method (SBM) alleviates this meshing burden by shifting the enforcement of boundary conditions from the true boundary to a nearby surrogate boundary and recovering the effect of the true boundary through geometric correction terms, thereby enabling standard finite element spaces on non-boundary-fitted meshes. In this report, we develop a general shiftedboundary and shifted-interface framework within the open-source MOOSE framework. We first present a general SBM implementation for complex geometries on non-boundary-fitted meshes. We then adopt the Shifted Interface Method (SIM) for internal interfaces and develop a unified shifted-interface treatment in which the interface law is enforced on a surrogate interface and the effect of the true interface is recovered through shifted jumps, fluxes, and tractions. This perspective brings scalar thermal-contact and vector-valued cohesive-zone mechanics into a single framework, the latter realized as the Shifted Cohesive Zone Method (SCZM) and coupled with history-dependent constitutive models from NEML2. We further extend the MOOSE mesh infrastructure to support cohesive-zone calculations on distributed meshes. The framework is verified and demonstrated through three progressive studies: Poisson’s equation on a smoothed starshaped domain, a manufactured thermal-contact problem on a non-interface-fitted mesh, and a two-dimensional polycrystalline representative volume element combining crystal plasticity with cohesive grain-boundary interfaces. Across these studies, the shifted formulations reproduce boundary- and interface-fitted reference solutions with high fidelity, indicating that the proposed framework provides an accurate and efficient route to boundary- and interface-dominated simulations on arbitrary geometries without requiring fitted meshes.

Yang, Cheng-Hau↗

The effects of 8-week complex training on lower-limb strength and power of Chinese elite female modern pentathlon athletes

Complex training (CT) is a combination training method that alternates between performing high-load resistance training (RT) and plyometric training within one single session. The study aimed to examine the effects of CT on lower-limb strength and power of elite female modern pentathlon athletes under the new modern pentathlon format and competition rules. Ten female participants (age: 23.55 ± 2.22 years, weight: 60.59 ± 3.87 kg, height: 169.44 ± 4.57 cm, and training experience: 6.90 ± 2.08 years) of the national modern pentathlon team completed 8 weeks of RT as followed by 8 weeks of CT, with 2 weeks of break. Then, the participants conducted 8 weeks of CT, which included RT combined with plyometric training (e.g., drop jump and continuous jump). All stages of training were designed by the linear strength training period theories, requiring participants to train twice for the first 4 weeks and three times for the second 4 weeks. The one-repetition maximum (1RM) of squat, isometric mid-thigh pull (IMTP), counter-movement jump (CMJ), squat jump (SJ), pre-stretch augmentation percentage (PSAP), and reaction strength index (RSI) were assessed before and after both RT and CT training. One-way repeated-measure ANOVA models revealed that the 1RM of squat was significantly improved ( p < 0.001) after RT as compared to pre-RT. No significant improvement in IMTP ( p = 0.055), CMJ ( p = 0.194), SJ ( p = 0.692), PSAP ( p = 0.087), and RSI ( p = 0.238) was not observed. After CT, 1RM of squat ( p < 0.001), IMTP ( p < 0.035), CMJ ( p < 0.001), SJ ( p < 0.008), RSI ( p < 0.006) were significant improved as compared to pre-RT, post-RT and pre-CT, while significant improvements in PSAP were observed as compared to pre-RT ( p = 0.003) and pre-CT ( p = 0.027), but not to post-RT ( p = 0.156). This pilot study showed the promise of CT following RT to improve lower-limb strength and power in elite female modern pentathlon athletes. The findings are worthwhile to be confirmed in future studies with larger sample size and randomized design.

Qiao, Zining↗

The effects of an ion-thruster exhaust plume on S-band carrier transmission

The study reported here was undertaken (1) to develop models of the effects of an ion-thruster exhaust plume on S-band signals, and (2) to measure the effects. The results show that an S-band signal passing through an ion-thruster plume is reduced in amplitude and advanced in phase. The mathematical models gave reasonable estimates of the average signal attenuation and phase shift. Negligible fluctuations in the signal amplitude and phase were measured during steady-state thruster operation. However, large jumps in phase occurred when changes were made in the thruster operating state. This study confirms that the thruster plume can have a significant effect on S-band communication link performance; hence the plume effects must be considered in S-band link calculations when electric thrusters are used for spacecraft propulsion.

Ackerknecht, W. E.↗

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗

Coronado Ecological Conservation: Assessing Vegetation Change Due to Border Wall Construction and Shifting Social Trails

Species monitoring is essential for mitigating the impacts of plant invasion, such as radical changes in an area’s ecosystem, degraded soil health, increased wildfire severity, landslides, and increased flooding. For this project, NASA DEVELOP partnered with the National Park Service (NPS) to investigate invasive species in disturbed lands: specifically, areas affected by off-trail travel and U.S.-Mexico border construction activities. The team assessed how construction has impacted the distribution of Lehmann’s lovegrass and Russian thistle invasives throughout Coronado National Memorial, AZ from 1986-2022. Using data from Landsat 5 and 8, Sentinel-2, NAIP, and PlanetScope, the team computed NDVI, NDMI, MSAVI2, EVI, and Tasseled Cap Wetness, Brightness, and Greenness transformations as vegetation health indicators to input into various machine learning algorithms. To minimize noise, the team conducted Principal Component Analysis on vegetation indices and spectral bands before running k-means clustering and random forest classification algorithms. Between all datasets, the team found that the median area fully overtaken by invasive plants was 5.37% of the park’s total area in 2022. The NPS will use end products to help increase restoration efforts in disturbed areas with high concentrations of invasive plants, and this project can serve as a jumping off point for future invasive species monitoring. The NPS’s collection of ground data for 2022-2023, in conjunction with future data collection, will notably improve the accuracy of classification models, leading to more precise monitoring of invasive species spread over time.

Coronado National Memorial↗

[Results of Simulation Studies]

Lattice Monte Carlo and off-lattice molecular dynamics simulations of h(sub 1)t(sub 4) and h(sub 4)t(sub l) (head/tail) amphiphile solutions have been performed as a function of surfactant concentration and temperature. The lattice and off-lattice systems exhibit quite different self-assembly behavior at equivalent thermodynamic conditions. We found that in the weakly aggregating regime (no preferred-size micelles), all models yield similar micelle size distributions at the same average aggregation number, albeit at different thermodynamic conditions (temperatures). In the strongly aggregating regime, this mapping between models (through temperature adjustment) fails, and the models exhibit qualitatively different micellization behavior. Incipient micellization in a model self-associating telechelic polymer solution results in a network with a transient elastic response that decays by a two-step relaxation: the first is due to a heterogeneous jump-diffusion process involving entrapment of end-groups within well-defined clusters and this is followed by rapid diffusion to neighboring clusters and a decay (terminal relaxation) due to cluster disintegration. The viscoelastic response of the solution manifests characteristics of a glass transition and entangled polymer network.

Source record↗

The essential character of the neutrino mechanism of core-collapse supernova explosions

Abstract Calibrating with detailed 2D core-collapse supernova (CCSN) simulations, we derive a simple CCSN explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behaviour of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ∼90 per cent of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code F ornax adheres to Nature and issues concerning the neutrino–matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.

Astronomy & Astrophysics↗

Staircase ion signature in the polar cusp - A case study

On 15 October, 1981 Dynamics Explorer 2 crossed the polar cusp at 1015 MLT and observed three distinct ion populations as it was moving poleward. These three populations had peak-flux energy around 2.7 keV, 850 eV, and 360 eV. At the time of observation, the IMF was southward. The first step coincided with a rotation of the flow from westward to poleward and then eastward. The second and third steps showed a flow directed principally poleward. Furthermore, the magnetic and electric perturbations in the first step are well fitted by an elongated flux tube footprint model. These results suggest that three consecutive Flux Transfer Events (FTEs) have injected solar wind plasma into the ionosphere forming the polar cusp. The individual FTE signatures can only be identified by the jumps in the precipitation pattern. The newest reconnected FTE footprint was crossed near the edge of the event while the two oldest ones were crossed around the event center. The small latitudinal size of these FTE footprints (40 km) and their short recurrence rate (3, 6 min) is consistent with an intermittent reconnection taking place at the subsolar point on a short time scale.

Escoubet, C. P.↗

Atomistic Investigation of Plastic Deformation and Dislocation Motion in Uranium Mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. However, many aspects of its mechanical behavior, particularly at reactor-relevant conditions, remain unclear. In this study, molecular dynamics (MD) simulations were employed to investigate the deformation behavior and dislocation motion in UN. We found that the Kocevski potential predicts the principal slip system as $\frac{1}{2}$ $\langle110\rangle${110}, aligning with experimental data. On the other hand, the Tseplyaev potential predicts slip to primarily occur on $\frac{1}{2}$ $\langle110\rangle${111}. MD simulations of stress–strain behavior were used to estimate the nanoindentation hardness, revealing that the Kocevski potential accurately predicts hardness even though it fails to model dynamic plasticity. Complete dislocation mobility functions have been fitted for the edge and screw dislocations in both the thermally activated and phonon-drag regimes. The 300 K linear mobility of the edge dislocation using the Tseplyaev potential was found to be 817 Pa -1 ·s -1 , whereas that of the screw dislocation using the Kocevski potential was found to be 4546 Pa -1 ·s -1 . At intermediate stresses, we observed that the subsonic steady-state motion of the edge dislocation in UN is intermittently interrupted by velocity jumps, reaching the average sound velocity. Finally, the threshold Schmid stress is calculated as 179–197 MPa, which gives an upper-limit estimate of the uniaxial yield stress of polycrystalline UN of 548–603 MPa. These findings, including the fitted dislocation mobility function, provide essential input for future plasticity and dislocation dynamics models of nuclear fuels.

36 MATERIALS SCIENCE↗

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

GTD analysis of airborne antennas radiating in the presence of lossy dielectric layers

The patterns of monopole or aperture antennas mounted on a perfectly conducting convex surface radiating in the presence of a dielectric or metal plate are computed. The geometrical theory of diffraction is used to analyze the radiating system and extended here to include diffraction by flat dielectric slabs. Modified edge diffraction coefficients valid for wedges whose walls are lossy or lossless thin dielectric or perfectly conducting plates are developed. The width of the dielectric plates cannot exceed a quarter of a wavelength in free space, and the interior angle of the wedge is assumed to be close to 0 deg or 180 deg. Systematic methods for computing the individual components of the total high frequency field are discussed. The accuracy of the solutions is demonstrated by comparisons with measured results, where a 2 lambda by 4 lambda prolate spheroid is used as the convex surface. A jump or kink appears in the calculated pattern when higher order terms that are important are not included in the final solution. The most immediate application of the results presented here is in the modelling of structures such as aircraft which are composed of nonmetallic parts that play a significant role in the pattern.

Rojas-Teran, R. G.↗

Corequake and shock heating model of the 5 March 1979 gamma ray burst

Ramatry, et al. proposed a model to account for the 5 March 1979 gamma ray burst in terms of a neutron star corequake and subsequent shock heating of the neutron star atmosphere. This model is extended by examining the overall energetics and characteristics of these shocks, taking into account the e(+)-e(-) pair production behind the shock. The effects of a dipole magnetic field in the shock jump conditions are also examined and it is concluded that the uneven heating produced by such a field can account for the temperature difference between pole and equator implied by the pulsating phase of the burst. The overall energetics and distribution of energy between e(+)-(-) pairs and photons appears to be in agreement with observations if this event is at a distance of 55 kpc as implied by its association with the Large Magellanic Cloud.

Ellison, D. C.↗

Corequake and shock heating model of the 5 March 1979 gamma ray burst

Ramatry, et al. proposed a model to account for the 5 March 1979 gamma ray burst in terms of a neutron star corequake and subsequent shock heating of the neutron star atmosphere. This model is extended by examining the overall energetics and characteristics of these shocks, taking into account the e(+)-e(-) pair production behind the shock. The effects of a dipole magnetic field in the shock jump conditions are also examined and it is concluded that the uneven heating produced by such a field can account for the temperature difference between pole and equator implied by the pulsating phase of the burst. The overall energetics and distribution of energy between e(+)-e(-) pairs and photons appear to be in agreement with observations if this event is at a distance of 55 kpc as implied by its association with the Large Magellanic Cloud. Previously announced in STAR as N83-31568

Ellison, D. C.↗

Vortex creep and the internal temperature of neutron stars Timing noise in pulsars

Vortex creep theory is used to construct model noise power spectra for three physically distinct types of events which might give rise to pulsar timing noise. These are 'pure' events, in which vortex unpinning is the source of the initial frequency jump; 'mixed' events, in which the initial frequency jump is produced by some physical process other than vortex unpinning but leads to the unpinning of some vortices; and 'external' events, in which the initial frequency jumps responsible for noise do not involve any vortex unpinning. For the first two types of events, it is found that relaxation processes in the region responsible for the noise will give rise to structure in the observed power spectra, while for external events, the resulting noise spectra will not be influenced by vortex creep. The theoretical results are compared with observed power spectra for 25 pulsars. The absence of structure in the observed power spectra of the Crab and Vela pulsars within the range of time scales which characterize their postglitch behavior indicates that the pinning regions which play a role in postglitch behavior do not experience the small unpinning events leading to timing noise.

Alpar, M. Ali↗

The energetic ion signature of an O-type neutral line in the geomagnetic tail

An energetic ion signature is presented which has the potential for remote sensing of an O-type neutral line embedded in a current sheet. A source plasma with a tailward flowing Kappa distribution yields a strongly non-Kappa distribution after interacting with the neutral line: sharp jumps, or ridges, occur in the velocity space distribution function f(nu-perpendicular, nu-parallel) associated with both increases and decreases in f. The jumps occur when orbits are reversed in the x-direction: a reversal causing initially earthward particles (low probability in the source distribution) to be observed results in a decrease in f, while a reversal causing initially tailward particles to be observed produces an increase in f. The reversals, and hence the jumps, occur at approximately constant values of perpendicular velocity in both the positive nu parallel and negative nu parallel half planes. The results were obtained using single particle simulations in a fixed magnetic field model.

Martin, R. F., Jr.↗