SEARCH · Engineering Papers
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.
Report on the evaluation of significant EV deployments on distribution feeders, mostly examining their load behavior during voltage sags from the transmission system. A particular emphasis on if a charger is grid-friendly or grid-unfriendly (especially in respect to inducing FIDVR behavior) is examined.
We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.
A computational method was developed for the fully-coupled solution of nongray, radiating gas flows with ablation product effects about blunt bodies during planetary entries. The treatment of radiation accounts for molecular band, continuum, and atomic line transitions with a detailed frequency dependence of the absorption coefficient. The ablation of the entry body was solved as part of the solution for a steady-state ablation process. The method was applied by results at typical conditions during entry to Venus. The radiative heating rates along the downstream region of the body can exceed the stagnation point value. The radiative heating to the body is attenuated in the boundary layer at the downstream region of the body and at the stagnation point of the body. A study of the radiation, inviscid flow about spherically capped, conical bodies during planetary entries shows that the nondimensional, radiative heating distributions are nonsimilar with entry conditions. Caution should be exercised in attempting to extrapolate results from known distributions to other entry conditions for which solutions have not yet been obtained.
The nonlinear space-charge effect plays a significant role in high-intensity accelerators and has been extensively studied using multi-particle tracking methods. In this paper, we present a novel 2.5- dimensional symplectic space-charge solver specifically designed for long beam bunches. We begin by detailing its application to a transverse Gaussian density distribution under open boundary conditions in a straight system, where a semi-analytical expression is derived. We then demonstrate the solver’s adaptation to arbitrary distributions in open space, as well as within rectangular and round conducting pipes. Finally, we discuss the extension of this solver to circular accelerator systems. This study shows that the fast 2.5-dimensional solver can be a good approximation to the fully three-dimensional solver for long bunches in large circular accelerators.
Interaction graphs on cascading outages of power systems provide valuable insights into how cascading outages evolve and propagate, and which components and links are critical to the propagation of cascading outages, enabling further development of mitigation strategies to support decision-making. However, the sensitivity of the interaction graph’s topology to the system’s loading condition has not been studied sufficiently. This paper compared interaction graphs under various loading conditions on the Northeastern Power Coordinating Council 140-bus system, and discovers the strong relationships between the graph topology, the cascade size distribution, and the load condition. Accordingly, three representative interaction graphs are constructed and illustrated.
Beam equation for elastically supported clamped- clamped beam solved for conditions of initial, boundary and distributed load, using Hermitian operator
Machine learning enables unbinned, highly-differential cross section measurements. A recent idea uses generative models to morph a starting simulation into the unfolded data. We show how to extend two morphing techniques, Schrödinger Bridges and Direct Diffusion, in order to ensure that the models learn the correct conditional probabilities. This brings distribution mapping (DM) to a similar level of accuracy as the state-of-the-art conditional generative unfolding methods. Numerical results are presented with a standard benchmark dataset of single jet substructure as well as for a new dataset describing a 22-dimensional phase space of Z+2 -jets.
Throughout the course of FY21, significant effort was put into investigating models within the LANL developed Hauser-Feshbach fission fragment decay code, CGMF, to understand and potentially solve the long-standing challenge of a too-soft prompt fission neutron spectrum, PFNS. Several inputs and models to CGMF were investigated, including the discrete nuclear levels, the optical model potential, level densities, and the fission fragment initial conditions. Some of the global models within CGMF led to a slight hardening of the neutron spectrum—particularly the likely incomplete discrete levels in through which γ-rays decay—but none of the changes where large enough for the tail of the PFNS to reproduce experimental data. A significant hardening of the spectrum tail was observed when the fission fragment initial conditions were optimized based on their sensitivities to the PFNS data for thermal incident neutrons. In this way, the parameters for the CGMF mass and total kinetic energy distributions, along with the spin cutoff factor were adjusted to better reproduce the experimental PFNS measurements. This optimization hardened the tail of the PFNS slightly but led to unphysical mass distributions for the fission fragments before neutron emission. It was clear from the above that we do not expect to be able to produce an evaluation-quality PFNS with CGMF in the near future. Challenges at thermal will persist–and possibly worsen–with increasing incident energy, where more models are needed to completely describe the fission. Basic-science research funding exceeding the amount available and scope of our NCSP funds would be needed to tackle this decade-long challenge impacting many fission-fragment event generator. And, in fact, Amy Lovell won LDRD ECR funding to do so over the next few years. Therefore, we focused in FY22 on extending evaluation capabilities beyond thermal incident neutrons, to take into account the incident energy dependence of the PFNS and fission fragment initial condition distributions in CGMF. We chose to set up the evaluation methodology to perform PFNS evaluations with CGMF across incident-neutron energies, in order to have it readily available for future NCSP evaluations when the PFNS from CGMF has improved. In this report, we outline the evaluation methodology, along with the results of the optimization, including full model calculations with CGMF using the evaluated parameters.
Abstract Gravitational lensing data is frequently collected at low resolution due to instrumental limitations and observing conditions. Machine learning-based super-resolution techniques offer a method to enhance the resolution of these images, enabling more precise measurements of lensing effects and a better understanding of the matter distribution in the lensing system. This enhancement can significantly improve our knowledge of the distribution of mass within the lensing galaxy and its environment, as well as the properties of the background source being lensed. Traditional super-resolution techniques typically learn a mapping function from lower-resolution to higher-resolution samples. However, these methods are often constrained by their dependence on optimizing a fixed distance function, which can result in the loss of intricate details crucial for astrophysical analysis. In this work, we introduce DiffLense , a novel super-resolution pipeline based on a conditional diffusion model specifically designed to enhance the resolution of gravitational lensing images obtained from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Our approach adopts a generative model, leveraging the detailed structural information present in Hubble space telescope (HST) counterparts. The diffusion model, trained to generate HST data, is conditioned on HSC data pre-processed with denoising techniques and thresholding to significantly reduce noise and background interference. This process leads to a more distinct and less overlapping conditional distribution during the model’s training phase. We demonstrate that DiffLense outperforms existing state-of-the-art single-image super-resolution techniques, particularly in retaining the fine details necessary for astrophysical analyses.
An intriguing regime of universal charge transport at high entropy density has been proposed for periodically driven interacting one-dimensional systems with Bloch bands separated by a large single-particle band gap. For weak interactions, a simple picture based on well-defined Floquet quasiparticles suggests that the system should host a quasisteady state current that depends only on the populations of the system's Floquet-Bloch bands and their associated quasienergy winding numbers. Here we show that such topological transport persists into the strongly interacting regime where the single-particle lifetime becomes shorter than the drive period. Analytically, we show that the value of the current is insensitive to interaction-induced band renormalizations and lifetime broadening when certain conditions are met by the system's nonequilibrium distribution function. We show that these conditions correspond to a quasisteady state. We support these predictions through numerical simulation of a system of strongly interacting fermions in a periodically modulated chain of Sachdev-Ye-Kitaev dots. Our paper establishes universal transport at high entropy density as a robust far from equilibrium topological phenomenon, which can be readily realized with cold atoms in optical lattices. Published by the American Physical Society 2024
A preliminary dynamic model of beaver spatial distribution and population growth was developed. The feasibility of locating beaver ponds on LANDSAT digital tapes, and of using this information to provide initial conditions of beaver spatial distribution for the model, and to validate model predictions is discussed. The techniques used to identify beaver ponds on LANDSAT are described.
The performance and noise of a high speed SR-7A model propeller under takeoff/landing conditions are considered. The blade loading distributions are obtained by solving the three-dimensional Euler equations and the sound pressure levels are computed using a time domain approach. At the nominal takeoff operating point, the blade sections near the hub are lightly or negatively loaded. The chordwise loading distributions are distinctly different from those of cruise conditions. The noise of the SR-7A model propeller at takeoff is dominated by the loading noise, similar to that at cruise conditions. The waveforms of the acoustic pressure signature are nearly sinusoidal in the plane of the propeller. The computed directivity of the blade passing frequency tone agrees fairly well with the data at nominal takeoff blade angle.
The performance and noise of a high speed SR-7A model propeller under takeoff/landing conditions are considered. The blade loading distributions are obtained by solving the three-dimensional Euler equations and the sound pressure levels are computed using a time domain approach. At the nominal takeoff operating point, the blade sections near the hub are lightly or negatively loaded. The chordwise loading distributions are distinctly different from those of cruise conditions. The noise of the SR-7A model propeller at takeoff is dominated by the loading noise, similar to that at cruise conditions. The waveforms of the acoustic pressure signature are nearly sinusoidal in the plane of the propeller. The computed directivity of the blade passing frequency tone agrees fairly well with the data at nominal takeoff blade angle.
The VSAERO low order panel method formulation is described for the calculation of subsonic aerodynamic characteristics of general configurations. The method is based on piecewise constant doublet and source singularities. Two forms of the internal Dirichlet boundary condition are discussed and the source distribution is determined by the external Neumann boundary condition. A number of basic test cases are examined. Calculations are compared with higher order solutions for a number of cases. It is demonstrated that for comparable density of control points where the boundary conditions are satisfied, the low order method gives comparable accuracy to the higher order solutions. It is also shown that problems associated with some earlier low order panel methods, e.g., leakage in internal flows and junctions and also poor trailing edge solutions, do not appear for the present method. Further, the application of the Kutta conditions is extremely simple; no extra equation or trailing edge velocity point is required. The method has very low computing costs and this has made it practical for application to nonlinear problems requiring iterative solutions for wake shape and surface boundary layer effects.
We compute the skewness of the galaxy distribution arising from the nonlinear evolution of arbitrary non-Gaussian intial conditions to second order in perturbation theory including the effects of nonlinear biasing. The result contains a term identical to that for a Gaussian initial distribution plus terms which depend on the skewness and kurtosis of the initial conditions. The results are model dependent; we present calculations for several toy models. At late times, the leading contribution from the initial skewness decays away relative to the other terms and becomes increasingly unimportant, but the contribution from initial kurtosis, previously overlooked, has the same time dependence as the Gaussian terms. Observations of a linear dependence of the normalized skewness on the rms density fluctuation therefore do not necessarily rule out initially non-Gaussian models. We also show that with non-Gaussian initial conditions the first correction to linear theory for the mean square density fluctuation is larger than for Gaussian models.
Comparison of STS-2 Shuttle flight heating data along the windward centerline has been made with two-dimensional nonequilibrium viscous shock-layer solutions obtained with shock and wall-slip conditions at an altitude range of 90 to 110 km. The shock slip condition used is the modified Rankine-Hugoniot relations of Cheng as used by Davis, and the wall-slip conditions are based on the first order consideration derived from kinetic theory as given by Scott and Hendricks. The results indicate that the calculated heating distributions with slip boundary conditions agree better with the flight data than those without slip conditions. The agreement improves when the accommodation coefficient or freestream density is decreased to one-half, suggesting the possibility of less than full accommodation for the tile surface and (or) an overestimate of freestream density using the Jacchia-Roberts model. Heating reduction due to the slip effect becomes very pronounced as the flow becomes more rarefied, and the effect is more significant for the stagnation region than the aft region of the vehicle.
The ability of recurrent neural networks (RNN) to model nonlinear dynamics of high dimensional process data has enabled data-driven RNN-based fault detection algorithms. Previous studies have focused on detecting faults by identifying the discrepancies in data distribution between the faulty and normal data, as reflected in prediction errors generated by RNN models. However, in industrial processes, variations in data distribution can also result from changes in normal control setpoints and compensatory control adjustments in response to disturbances, making it hard to differentiate between normal and faulty conditions. This paper proposes a fault detection method utilizing a long short-term memory (LSTM) encoder–decoder structure with regularized hidden dynamics and reversible instance normalization (RevIN) to compactly represent high-dimensional measurements for effective monitoring. During training, the hidden states of the model are regularized to form a low-dimensional latent space representation of the original multivariate time series data. As a result, the prediction errors of the latent states can be used to monitor the abnormal dynamic variations, while the reconstruction errors of the measured variables are used to monitor the abnormal static variations. Furthermore, the proposed indices can reflect operating conditions, even when the distribution of test data changes, which helps distinguish faults from normal adjustments and disturbances that controllers can settle. Here, data from numerical simulation and the Tennessee Eastman process are used to illustrate the effectiveness of the proposed fault detection method.
The dynamics of the satellite systems of Milky-Way-like galaxies offer a useful means by which to study the galaxy formation process in the cosmological context. It has been suggested that the currently observed anisotropic distribution of the satellites in such galaxy systems is inconsistent with the concordance ΛCDM cosmology model on the galactic scale if the observed satellites are random samples of the dark matter (DM) sub-halos that are nearly isotropically distributed around the central galaxy. In this study, we present original high-resolution zoom-in studies of central galaxies and satellite systems based upon initial conditions for the DM distribution from the Aquarius simulations but with substantial high-resolution baryon physics added. We find that the galaxy most like the Milky Way in this study does indeed contain a disk of satellites (DOS). Although one galaxy DOS system does not answer the question of how common such disks are, it does allow the opportunity to explore the properties and dynamics of the DOS system. Our investigation centers on the spatial arrangement (distances, angles, etc.) of satellites in this Milky-Way-like galaxy system with a specific emphasis on identifying and analyzing the disk-like structure along with its dynamical and morphological properties. Among the conclusions from this study, we find that the satellites and DM sub-halos in the galaxy simulations are anisotropically distributed. The dynamical properties of the satellites, however, indicate that the direction of the angular momentum vector of the whole satellite system is different from the normal direction of the fitted DOS and from the normal direction of the velocity dispersion of the system. Hence, the fitted DOS appears to be comprised of infalling sub-halos and is not a rotationally supported system.