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.
A study is made in the physical plane of the flow properties near a conical cross flow sonic line. The shapes of the Mach lines at regular and singular sonic points are presented. Through a study of the jump conditions at a singular sonic point, the possibility of a double sonic line is discussed. The properties of bubbles of supersonic cross flow are examined. Numerical results are presented to support some of the theoretical observations.
We have developed and implemented a retrospective data assimilation system (RDAS) as an upgrade to the operational DAO/Terra data assimilation system. This formulation aims at improving analysis over filter analysis by the dynamically consistent incorporation of observation information past a given analysis time. The current implementation of the RDAS uses the adjoint of the tangent linear model of a simplified version of the Terra general circulation model and extensions to the physical-space statistical analysis system to propagate observation information back in time. The RDAS adopts the same assumptions of the regular data assimilation system, particularly, no explicit propagation of error covariances are involved therefore rendering a procedure that is computationally affordable. In this study, we show results of experiments conducted to investigate the performance of the 6-hour (lag-1) RDAS. Statistical results obtained over one month during a winter season indicate that the RDAS represents considerable improvement over the regular assimilation. Plans for implementation of the RDAS capability in our new finite-volume data assimilation system will also be presented at the time of the conference.
Degradation of lithium trivanadate ( Li x V 3 O 8 ) cathodes has been widely reported in the literature, but studies have offered little insight towards developing a detailed understanding of the evolution of the active material, and have been inconclusive as to the root cause of degradation. Here, we refit a validated physics-based model to discharge curves over the course of cycling at C/5, and use the evolution of the model parameters to track evolution of the cathode. A regularization penalty for adjusting model parameters from the validated model is introduced as a framework to identify which model parameters can explain a significant portion of the observed change in the voltage profile over the course of cycling. SEM reveals that lithium trivandate particles fracture upon cycling at C/5, consistent with the results of the parameter estimation, deactivation of lithium trivanadate and faster diffusion of lithium within the active particles. The physics-based model is then used to design modified cycling protocols which identify the phase transformation during discharge of lithium trivanadate as the driver of the particle fracture and capacity fade.
Abstract Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v ⊥ / v A , is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas ‘might as well’ be occurring in a cylinder.
High-temperature superconductor (HTS) tapes have shown promising characteristics of high critical current, which are prerequisites for applications in high-field magnets. Due to the unstable growth conditions in the HTS manufacturing process, however, the frequent occurrences of dropouts in the critical current impede the consistent performance of HTS tapes. To manufacture HTS tapes with large scale, high yield, and uniform performance, it is essential to develop novel data analysis approaches for modeling the dropouts and identifying the related important process parameters. Conventional methods for modeling recurrent events, such as the point process, require the extraction of events from quality measurements. As the critical current is a continuous process, it may not comprehensively represent the drop patterns by transforming the time-series measurements into a set of events. Here, to solve this issue, we develop a novel quantile regression-enriched event modeling (QREM) framework that integrates the non-homogeneous Poisson process for modeling the occurrence of dropouts and the quantile regression for capturing the drop patterns. By incorporating the feature selection and regularization, the proposed framework identifies a set of significant process parameters that can potentially cause the dropouts of HTS tapes. The proposed method is tested on real HTS tapes produced using an advanced manufacturing process, successfully identifying important parameters that influence dropout events including the substrate temperature and voltage. The results demonstrate that the proposed QREM method outperforms the standard point process in predicting the occurrence of dropouts.
G2Aero is a flexible and practical tool for design and deformation of 2D airfoils and 3D blades using data-driven approaches. G2Aero utilizes the geometry of matrix manifolds—specifically the Grassmannian—to build a novel framework for representing physics-based separable deformations of shapes. G2Aero offers the flexibility to generate perturbations in a customizable way over any portion of the blade. The G2Aero framework utilizes data-driven methods based on a curated database of physically relevant airfoils. Specific tools include: - principal geodesic analysis over normal coordinate neighborhoods of matrix manifolds; - a variety of data-regularized deformations to nominal 2D airfoil shapes; - Riemannian interpolation connecting a sequence of airfoil cross-sections to build 3D blades from 2D data; - consistent perturbations over the span of interpolated 3D blades based on dominant modes from the data-driven analysis.
The global constants associated with orthotropic slab/plate equations are discussed, and the rotational behavior of the modulus/compliance components associated with orthotropic slabs/plates are addressed. It is concluded that one cluster constant is less than or equal to unity for all physically possible materials. Rotationally anomalous behavior is found in two materials, and a simple inequality which can be used to identify regular or anomalous behavior is presented and discussed in detail.
Loop corrections to finite-time correlation functions in quantum field theories away from equilibrium can be calculated using the in-in path integral approach. In this paper, we calculate the unequal-time two-point correlator for different massless self-interacting scalar quantum field theories on a Minkowski background, starting the field evolution at an arbitrary initial time. We find the counterterms that need to be added to UV-renormalize the result, including usual in-out counterterms in the dynamics and additional initial state counterterms that are required to cancel all UV divergences. We find that the late-time limit of the renormalized correlation function exhibits a linear or logarithmic growth in time, depending on whether the interaction strength is dimension-one or dimensionless, respectively. The late-time correlations match those obtained in our companion paper and, as shown there, the divergences do not indicate a real IR issue, consistent with what one would expect in Minkowski.
We study fractional-derivative Maxwell theory, as appears in effective descriptions of, for example, large N f QED 3 , graphene, and some types of surface defects. We argue that when the theory is realized on a lattice, monopole condensation leads to a confining phase via the Polyakov confinement mechanism.
There is an ambiguity in choosing field-strength renormalization factors in the $ \overline{\mathrm{MS}} $ scheme starting from the 3-loop order in perturbation theory. More concerning, trivially choosing Hermitian factors has been shown to produce divergent renormalization group functions, which are commonly understood to be finite quantities. We demonstrate that the divergences of the RG functions are such that they vanish in the RG equation due to the Ward identity associated with the flavor symmetry. It turns out that any such divergences can be removed using the renormalization ambiguity and that the use of the flavor-improved β-function is preferred. We show how our observations resolve the issue of divergences appearing in previous calculations of the 3-loop SM Yukawa β-functions and provide the first calculation of the flavor-improved 3-loop SM β-functions in the gaugeless limit.
A measurable chaos indicator is used as the online optimization objective in tuning a complicated nonlinear system—the National Synchrotron Light Source-II storage ring. Through analyzing the Shannon entropy in measured Poincaré maps, not only can the commonly used nonlinear characterizations be extracted, but more importantly, the chaos can be quantified and then used for an online regularization of these maps. The method itself is general and applicable to other tunable nonlinear systems as well. Published by the American Physical Society 2025
Swimming bacteria of the species Bacillus subtilis require and consume oxygen. In static liquid cultures the cells' swimming behaviour leads them to accumulate up oxygen concentration gradients generated by consumption and supply. Since the density of bacterial cells exceeds that of the fluid in which they live, fluid regions where cells have accumulated are denser than depleted regions. These density variations cause convection. The fluid motion is dynamically maintained by the swimming of the cells toward regions of attraction: the air-fluid interface and the fluctuating advecting attractors, gradients of oxygen concentration that are embedded in the convecting fluid. Because of the fluid dynamical conservation laws, these complex physical and biological factors generate patterns ordered over distances > 10000 bacterial cell diameters. The convection enhances long-range transport and mixing of oxygen, cells and extracellular products by orders of magnitude. Thus, through the interplay of physical and biological factors, a population of undifferentiated selfish cells creates functional dynamic patterns. Populations of bacteria that have organised themselves into regularly patterned regions of vigorous convection and varying cell concentration interact with their environment as if they were one purposeful, coherent multicellular individual. The mathematical and experimental ingredients of these remarkable phenomena are presented here.
NSLS-II accelerator physicists provided operation support for user operations and machine start-ups, carried out routine lattice characterization, and produced regular reports on beam dynamics. The commissioning of a superconducting wiggler for the HEX beamline is complete, and the wiggler effects on beam dynamics have been studied. A Facility Improvement Project for absolute lattice characterization has been approved and funded. We implemented and tested a new orbit-based technique of lattice characterization and correction. We developed a theory to calculate the beam-induced power in ceramic vacuum chambers and a new convergence map technique for nonlinear lattice characterization. We continued working on the NSLS-II upgrade including lattice development, commissioning a low-energy complex bend prototype, analyzing beam sensitivity to various errors, implementing higher-harmonic RF cavities, studying collective effects, and assessing FEL options. We also supported lattice optimization and beam diagnostics development for the EIC project. Main research accomplishments achieved in 2023, are summarized in this report.
An attempt is made to provide a perspective on the problem of spatial structuring on scales smaller than can presently be directly and regularly observed from the ground or with which current space-based instrumentation can be anticipated. There is abundant evidence from both observations and theory that such spatial structuring of the solar outer atmosphere is ubiquitous not only on the observed scales, but also on spatial scales down to (at least) the subarcsecond range. This is not to say that the results to be obtained from observations on these small scales can be anticipated: quite the opposite. What is clear instead is that many of the classic problems of coronal and chromospheric activity - involving the basic dissipative nature of magnetized plasmas - will be seen from a novel perspective at these scales, and that there are reasons for believing that dynamical processes of importance to activity on presently-resolved scales will themselves begin to be resolved on the sub-arcsecond level. Since the Sun is the only astrophysical laboratory for which there is any hope of studying these processes in any detail, this observatioinal opportunity is an exciting prospect for any student of magnetic activity in astrophysics.
Characterizing the entanglement of matrix degrees of freedom is essential for understanding the holographic emergence of spacetime. The Quantum Hall Matrix Model is a gauged U(N) matrix quantum mechanics with two matrices whose ground state is known exactly and describes an emergent spatial disk with incompressible bulk dynamics. We define and compute an entanglement entropy in the ground state associated to a cut through the disk. There are two contributions. A collective field describing the eigenvalues of one of the matrices gives a gauge-invariant chiral boundary mode leading to an expected logarithmic entanglement entropy. Further, the cut through the bulk splits certain ‘off-diagonal’ matrix elements that must be duplicated and associated to both sides of the cut. Sewing these duplicated modes together in a gauge-invariant way leads to a bulk ‘area law’ contribution to the entanglement entropy. All of these entropies are regularized by finite N.
Advection-dispersion equations (ADEs) are commonly used to describe transport phenomena in porous media. Even though mature discretization-based numerical methods for ADEs exist, some challenges still remain, especially when it comes to solving advection-dominated forward ADEs and diffusion-dominated backward ADEs. The latter problem usually arises in the source identification context and leads to numerically unstable grid-based solutions that require a form of regularization or should be treated as an inverse problem that is computationally more expensive because it requires solving the forward problem multiple times. In this study, we propose a discretization-free approach based on the physics-informed neural network (PINN) method for solving coupled ADE and Darcy flow equations with space-dependent hydraulic conductivity. In this approach, the hydraulic conductivity, hydraulic head, and concentration fields are approximated with deep neural networks (DNNs). We assume that the conductivity field is given by its values on a grid, and we use these values to train the conductivity DNN. The head and concentration DNNs are trained by minimizing the residuals of the flow equation and ADE and using the initial and boundary conditions as additional constraints. The PINN method is applied to one- and two-dimensional forward ADE problems, where its performance for various P\'{e}clet numbers ($Pe$) is compared with the analytical and numerical solutions. We find that the PINN method is accurate with errors of less than 1\% and outperforms some conventional discretization-based methods for $Pe$ larger that 100. Next, we demonstrate that the PINN method remains accurate for the backward ADEs, with the relative errors in most cases staying under 5\% compared to the reference concentration field. Finally, we show that when available, the concentration measurements can be easily incorporated in the PINN method and significantly improve (by more than 50\% in the considered cases) the accuracy of the PINN solution of the backward ADE.
Not Available