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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗

DOE OSTI2021

On the non-existence of stepped-pressure equilibria far from symmetry

The Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502) has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. Here we found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov–Arnold–Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.

3D equilibrium↗

NASA NTRS1994

HST eclipse mapping of dwarf nova OY Carinae in quiescence: An 'Fe II curtain' with Mach approx. = 6 velocity dispersion veils the white dwarf

Hubble Space Telescope (HST) observations of the eclipsing dwarf nova OY Car in its quiescent state are used to isolate the ultraviolet spectrum (1150-2500 A at 9.2 A Full Width at Half Maximum (FWHM) resolution) of the white dwarf, the accretion disk, and the bright spot. The white dwarf spectrum has a Stark-broadened photospheric L(alpha) absorption, but is veiled by a forest of blended Fe II features that we attribute to absorption by intervening disk material. A fit gives T(sub w) approx. = 16.5 x 10(exp 3) K for the white dwarf with a solar-abundance, log g = 8 model atmosphere, and T approx. = 10(exp 4) K, n(sub e) approx. = 10(exp 13)/cu cm, N(sub H) approx. = 10(exp 22) sq cm, and velocity dispersion delta V approx. = 60 km/s for the veil of homogeneous solar-abundance local thermodynamic equilibrium (LTE) gas. The veil parameters probably measure characteristic physical conditions in the quiescent accretion disk or its chromosphere. The large velocity dispersion is essential for a good fit; it lowers (chi square)/778 from 22 to 4. Keplerian shear can produce the velocity dispersion if the veiling gas is located at R approx. = 5 R(sub W) with (delta R)/R approx. = 0.3, but this model leaves an unobscured view to the upper hemisphere of the white dwarf, incompatible with absorptions that are up to 80% deep. The veiling gas may be in the upper atmosphere of the disk near its outer rim, but we then require supersonic (Mach approx. = 6) but sub-Keplerian (delta V/V(sub Kep) approx. = 0.07) velocity disturbances in this region to produce both the observed radial velocity dispersion and vertical motions sufficient to elevate the gas to z/R = cos i = 0.12. Such motions might be driven by the gas stream, since it may take several Kepler periods to reestablish the disk's vertical hydrostatic equilibrium. The temperature and column density of the gas we see as Fe II absorption in the ultraviolet are similar to what is required to produce the strong Balmer jump and line emissions seen in optical spectra of OY Car and similar quiescent dwarf novae. The outer accretion disk is detected at mid-eclipse with a spectrum that rises from 0.05 to 0.3 mJy between 2000 and 2500 A, consistent with combinations of cool blackbodies, blended Fe II emission lines, and Balmer continuum emission. The total disk flux density is 0.5 mJy at 2500 A, and this shallow disk eclipse implies a roughly flat surface brightness distribution. The bright spot, somewhat bluer than the disk, has a flux density rising from 0.05 to 0.15 mJy between 1600 and 2500 A. The C IV emission line has a broad shallow eclipse, but the radial velocity variations observed during the eclipse do not clearly distinguish between a disk or wind origin. The only possible indications of boundary layer emission are fast UV flares that appear to arise from near the central object -- not from the bright spot.

Horne, Keith↗

NASA NTRS2020

Fisher's Peak Ecological Forecasting: Mapping Biomass to Inform Conservation Planning of a Future State Park in Southern Colorado

Fisher’s Peak is a 77.5 km2 property southeast of Trinidad, Colorado that is planned to become Colorado’s newest state park. The area has experienced limited anthropogenic disturbance and is home to an abundance of unique habitats and species. A rapid, approximately 900 m change in elevation over the extent of the area nurtures a variety of plants and animals, including the endangered New Mexico meadow jumping mouse. In 2019, the State of Colorado obtained Fisher’s Peak with plans to make it Colorado’s second largest state park. A diverse group of collaborators, including the Colorado State Forest Service and The Nature Conservancy, worked closely to design the state park to maximize recreation opportunity while conserving the property’s rich habitats and biodiversity. The Fisher’s Peak Ecological Forecasting Team utilized Light Detection and Ranging (LiDAR) surveys, in situ forest inventory data, and Earth observations from Landsat8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Sentinel-2 Multispectral Instrument (MSI), Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) and the Shuttle Radar Topography Mission (SRTM) to quantify and map biomass over the extent of the study area. The results from modeling biomass had an out-of-bag root mean square error of 55 Mg/ha and an R2 of 12. The resulting map indicates areas where carbon storage on the property is high, informing decision-making processes for future park development. While more in situ training data may improve modeling capacity for biomass in the Fisher’s Peak area, this work represents a feasible attempt to better understand biomass distribution using earth observation.

DEVELOP Tech Paper↗

NASA NTRS2001

Interaction of the body, head, and eyes during walking and turning

Body, head, and eye movements were measured in five subjects during straight walking and while turning corners. The purpose was to determine how well the head and eyes followed the linear trajectory of the body in space and whether head orientation followed changes in the gravito-inertial acceleration vector (GIA). Head and body movements were measured with a video-based motion analysis system and horizontal, vertical, and torsional eye movements with video-oculography. During straight walking, there was lateral body motion at the stride frequency, which was at half the frequency of stepping. The GIA oscillated about the direction of heading, according to the acceleration and deceleration associated with heel strike and toe flexion, and the body yawed in concert with stepping. Despite the linear and rotatory motions of the head and body, the head pointed along the forward motion of the body during straight walking. The head pitch/roll component appeared to compensate for vertical and horizontal acceleration of the head rather than orienting to the tilt of the GIA or anticipating it. When turning corners, subjects walked on a 50-cm radius over two steps or on a 200-cm radius in five to seven steps. Maximum centripetal accelerations in sharp turns were ca.0.4 g, which tilted the GIA ca.21 degrees with regard to the heading. This was anticipated by a roll tilt of the head of up to 8 degrees. The eyes rolled 1-1.5 degrees and moved down into the direction of linear acceleration during the tilts of the GIA. Yaw head deviations moved smoothly through the turn, anticipating the shift in lateral body trajectory by as much as 25 degrees. The trunk did not anticipate the change in trajectory. Thus, in contrast to straight walking, the tilt axes of the head and the GIA tended to align during turns. Gaze was stable in space during the slow phases and jumped forward in saccades along the trajectory, leading it by larger angles when the angular velocity of turning was greater. The anticipatory roll head movements during turning are likely to be utilized to overcome inertial forces that would destabilize balance during turning. The data show that compensatory eye, head, and body movements stabilize gaze during straight walking, while orienting mechanisms direct the eyes, head, and body to tilts of the GIA in space during turning.

Non-NASA Center↗

DOE OSTI2023

Strain-tuned quantum criticality in electronic Potts-nematic systems

Motivated by recent observations of threefold rotational symmetry breaking in twisted moiré systems, cold-atom optical lattices, quantum Hall systems, and triangular antiferromagnets, we phenomenologically investigate the strain-temperature phase diagram of the electronic 3-state Potts-nematic order. While in the absence of strain the quantum Potts-nematic transition is first-order, quantum critical points (QCP) emerge when uniaxial strain is applied, whose nature depends on whether the strain is compressive or tensile. In one case, the nematic amplitude jumps between two non-zero values while the nematic director remains pinned, leading to a symmetry-preserving metanematic transition that terminates at a quantum critical end-point. For the other type of strain, the nematic director unlocks from the strain direction and spontaneously breaks an in-plane twofold rotational symmetry, which in twisted moiré superlattices triggers an electric polarization. Such a piezoelectric transition changes from first to second-order upon increasing strain, resulting in a quantum tricritical point. Using a Hertz-Millis approach, we show that these QCPs share interesting similarities with the widely studied Ising-nematic QCP. The existence of three minima in the nematic action also leaves fingerprints in the strain-nematic hysteresis curves, which display multiple loops. At non-zero temperatures, because the upper critical dimension of the 3-state Potts model is smaller than three, the Potts-nematic transition is expected to remain first-order in 3D, but to change to second-order in 2D. As a result, the 2D strain-temperature phase diagram displays two first-order transition wings bounded by lines of critical end-points or tricritical points, reminiscent of the phase diagram of metallic ferromagnets. Furthermore, we discuss how our results can be used to unambiguously identify spontaneous Potts-nematic order.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2020

Fisher's Peak Ecological Forecasting - Mapping Biomass to Inform Conservation Planning of a Future State Park in Southern Colorado

Fisher’s Peak is a 77.5 km2 property southeast of Trinidad, Colorado that is planned to become Colorado’s newest state park. The area has experienced limited anthropogenic disturbance and is home to an abundance of unique habitats and species. A rapid, approximately 900 m change in elevation over the extent of the area nurtures a variety of plants and animals, including the endangered New Mexico meadow jumping mouse. In 2019, the State of Colorado obtained Fisher’s Peak with plans to make it Colorado’s second largest state park. A diverse group of collaborators, including the Colorado State Forest Service and The Nature Conservancy, worked closely to design the state park to maximize recreation opportunity while conserving the property’s rich habitats and biodiversity. The Fisher’s Peak Ecological Forecasting Team utilized Light Detection and Ranging (LiDAR) surveys, in situ forest inventory data, and Earth observations from Landsat 8 Operational Land Imager (OLI), Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), Sentinel-2 Multispectral Instrument (MSI), Advanced Land Observing Satellite 2 (ALOS-2) Phased Array type L-band Synthetic Aperture Radar (PALSAR-2) and the Shuttle Radar Topography Mission (SRTM) to quantify and map biomass over the extent of the study area. The results from modeling biomass had an out-of-bag root mean square error of 55 Mg/ha and an R2 of 12. The resulting map indicates areas where carbon storage on the property is high, informing decision-making processes for future park development. While more in situ training data may improve modeling capacity for biomass in the Fisher’s Peak area, this work represents a feasible attempt to better understand biomass distribution using earth observation

Lauren Lad↗

DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗

DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗

DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗

DOE OSTI2021

Criticality of Two-Dimensional Disordered Dirac Fermions in the Unitary Class and Universality of the Integer Quantum Hall Transition

Two-dimensional (2D) Dirac fermions are a central paradigm of modern condensed matter physics, describing low-energy excitations in graphene, in certain classes of superconductors, and on surfaces of 3D topological insulators. At zero energy E=0, Dirac fermions with mass m are band insulators, with the Chern number jumping by unity at m=0. This observation lead Ludwig et al. [Phys. Rev. B 50, 7526 (1994)PRBMDO0163-182910.1103/PhysRevB.50.7526] to conjecture that the transition in 2D disordered Dirac fermions (DDF) and the integer quantum Hall transition (IQHT) are controlled by the same fixed point and possess the same universal critical properties. Given the far-reaching implications for the emerging field of the quantum anomalous Hall effect, modern condensed matter physics, and our general understanding of disordered critical points, it is surprising that this conjecture has never been tested numerically. In this work, we report the results of extensive numerics on the phase diagram and criticality of 2D DDF in the unitary class. We find a critical line at m=0, with an energy-dependent localization length exponent. At large energies, our results for the DDF are consistent with state-of-the-art numerical results ν IQH =2.56-2.62 from models of the IQHT. At E=0, however, we obtain ν 0 =2.30-2.36 incompatible with ν IQH . This result challenges conjectured relations between different models of the IQHT, and several interpretations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Employing Interacting Qubits for Distributed Microgrid Control

To empower flexible and scalable operations, distributed control of multi-inverter microgrids, based on classical communication networks among distributed energy resources, has attracted considerable attention as it can guarantee synchronization and provide suitable remedies to the problem of improper power sharing. Notwithstanding this, resilience of the current schemes on classical communication makes microgrids vulnerable to cyber attacks. Inspired by recent revolutionary breakthroughs in quantum communication, in this paper, we devise a novel synchronization mechanism. We extend the synchronization framework utilized in distributed control algorithms to networks of quantum systems by generating pinning terms and coupling mechanism for the new synchronization rule via exploiting proper quantum jump operators and observables, and show that the quantum system will converge to a time-variant target state. Our devised quantum distributed controller (QDC) gives rise to a novel quantum communication scheme for distributed control of microgrids and enables microgrids to exploit the state-of-the-art quantum communication frameworks as communication infrastructure. Finally, test results on two representative AC and DC networked microgrids validate the efficacy and universality of the quantum distributed control.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRSDate not supplied

Anomaly Detection for the Roman Space Telescope Wide Field Instrument’s Science Data Processing Pipeline

The Roman Space Telescope (RST) Wide Field Instrument (WFI) will be utilizing a preliminary Science Data Processing (SDP) pipeline during its Integration and Test, and to some extent during Operations, to track basic statistics and identify known features such as cosmic rays, snowballs as well as possible anomalies in raw detector data. In our detectors, these anomalies appear as jumps in the ramp of a readout and are classified as cosmic rays if they appear as a streak or snowballs if they’re more circular. The WFI employs an array of 18 H4RG-10 detectors that collect image samples. Each set of raw frames within a non-destructive exposure is packaged by the SDP pipeline into image cubes for each detector. Each cube is a time series of 4096 × 4096 accumulating pixel frames. The preliminary analysis pipeline is used to locate anomalies in these time-series accumulation frames and identify the type of anomaly, either natural phenomena or detector characteristic. To compare different methods, we’ve implemented both heuristic-based and data-driven methods to identify anomalies. For the heuristic-based approach, we identify snowballs and cosmic rays by the size and shape of outlier pixel clusters between consecutive frames. For data driven methods, we evaluated a Convolutional Neural Network (CNN) model, and more traditional methods like Principal Component Analysis (PCA). CNN is a supervised learning/classification method. Thus, we used a labeled dataset of anomalies to perform segmentation of the image and identify anomalies. We used previously identified cosmic rays and snowballs to measure the accuracy and efficiency of the mentioned approaches. In evaluating these methods, we aim to pick the best fit for the SDP pipeline’s anomaly detection in terms of both performance and runtime.

Paul Horton↗

DOE OSTI2020

Hydrogen Diffusivity Measurements of YH 1.87 Moderator Material with Incoherent Quasielastic Neutron Scattering

Yttrium hydride is an optimal choice of moderator material for thermal neutron spectrum reactors requiring small core volumes. The Transformation Challenge Reactor (TCR) program at Oak Ridge National Laboratory (ORNL) aims to develop an additively manufactured microreactor core by leveraging recent advances in materials, manufacturing, data analytics, and high-fidelity modeling and simulation. Yttrium hydride was selected as the moderator material for TCR due to its high moderating power, which enables the spatially efficient reactor design and thermal stability that is superior when compared with zirconium hydride. Hydrogen desorption from yttrium hydride is expected at elevated temperatures, posing safety and performance concerns. As a lifetime component in advanced reactors, appropriate strategies to mitigate hydrogen release from yttrium hydride over long-term reactor operation are required. Understanding hydrogen mobility in yttrium hydride provides the scientific basis for developing a hydrogen barrier for hydride moderator and predicting the hydrogen redistribution within the material during various operational conditions. This study investigated the hydrogen motion in YH 1.87 as a function of temperature ranging from 200 to 1,173 K using incoherent quasielastic neutron scattering (IQNS) at the Spallation Neutron Source (SNS). The results are presented and discussed in terms of hydrogen self-diffusion coefficients, activation energy for translational hydrogen motion, and hydrogen jump distances. Based on the IQNS data, YH 1.87 is found to be stable at high temperatures with no translational H motion below 1,023 K within the IQNS detection limits. In the temperature range of 1,073–1,173 K, hydrogen diffusivity in YH 1.87 is $D = 4.57 X 10^{⁻3} [m^2/s] exp (⁻\frac{1.73[eV]}{RT}$). The hydrogen concentration’s dependence on hydrogen diffusivity is also discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

A dynamical mechanism for the Page curve from quantum chaos

If the evaporation of a black hole formed from a pure state is unitary, the entanglement entropy of the Hawking radiation should follow the Page curve, increasing from zero until near the halfway point of the evaporation, and then decreasing back to zero. The general argument for the Page curve is based on the assumption that the quantum state of the black hole plus radiation during the evaporation process is typical. In this paper, we show that the Page curve can result from a simple dynamical input in the evolution of the black hole, based on a recently proposed signature of quantum chaos, without resorting to typicality. Our argument is based on what we refer to as the “operator gas” approach, which allows one to understand the evolution of the microstate of the black hole from generic features of the Heisenberg evolution of operators. One key feature which leads to the Page curve is the possibility of dynamical processes where operators in the “gas” can “jump” outside the black hole, which we refer to as void formation processes. Such processes are initially exponentially suppressed, but dominate after a certain time scale, which can be used as a dynamical definition of the Page time. In the Hayden-Preskill protocol for young and old black holes, we show that void formation is also responsible for the transfer of information from the black hole to the radiation. We conjecture that void formation may provide a microscopic explanation for the recent semi-classical prescription of including islands in the calculation of the entanglement entropy of the radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Negative string tension of a higher-charge Schwinger model via digital quantum simulation

We study some properties of generalized global symmetry for the charge- q Schwinger model in the Hamiltonian formalism, which is the (1 + 1)D quantum electrodynamics with a charge- q Dirac fermion. This model has the |$\mathbb {Z}_q\, 1$| -form symmetry, which is a remnant of the electric |$U(1)\, 1$| -form symmetry in the pure Maxwell theory. It is known that, if we put the theory on closed space, then the Hilbert space is decomposed into q distinct sectors, called universes, and some states with higher energy density do not decay to the ground state due to the selection rule of the 1-form symmetry. Even with open boundaries, we can observe the stability of such states by seeing a negative string tension behavior, meaning that opposite charges repel each other. In order to see negative string tensions, the vacuum angle θ has to be large enough and the standard path-integral Monte Carlo method suffers from the sign problem. We develop a method based on the adiabatic state preparation to see this feature with digital quantum simulation and confirm it using a classical simulator of quantum devices. In particular, we measure the local energy density and see how it jumps between the inside and outside of the insertion of the probe charges. We explicitly see that the energy density inside is lower than that outside. This is a clear signature of the negative string tension.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2019

Evolution and Impact of Saturn V on Space Launch System from a Guidance, Navigation, and Mission Analysis Perspective

The Saturn V launch vehicle represented a jump in capability for heavy lift launch vehicles, enabling the Lunar Orbit Rendezvous approach to planetary exploration employed by the Apollo program 50 years ago. Following Apollo, and the development of the Space Transportation System, the NASA space exploration program shifted focus from lunar exploration to long-term, sustained, re-usable access to Low Earth Orbit. With the recent focus of NASA on the Artemis program and continued exploration of cislunar space as a precursor to Martian exploration, the shift has swung back to heavy lift capability. To meet this need, NASA has developed the Space Launch System. While the vehicle is a new design, it is heavily influenced by the engineering solutions and approach used on the Saturn V while taking advantage of the state of the art of launch vehicle design. The approach to abort, for example, shares many familiarities with the triggers and concept of operations used on Saturn V. Analysis approaches to dispersed trajectory performance are also very similar, but advances in computing technology have enabled a much more expanded set of inputs that can be modelled and assessed in a rapid manner. Additionally, guided flight algorithms share similar first principles but have expanded to include day of launch wind information. Trajectory optimization has also advanced significantly due to the availability of computing resources, but similar maneuvers and profiles are flown across both vehicles. Also, while the approach of onboard inertial navigation has been maintained between the two programs, the shift from platform to strapdown systems enables reduced complexity in the system design while maintaining required performance. As described, the Space Launch System is the evolution of NASA launch vehicle designs, owing a large heritage to the Saturn vehicle program and incorporating advances in propulsion systems, avionics, computing, and sensor technology over the past 50 years.

Ahmad, Naeem↗

NASA NTRS2018

Effects of Multiscale Phase-Mixing and Interior Conductance in the Lunar-like Pickup Ion Plasma Wake. First Results from 3-D Hybrid Kinetic Modeling

The study of multiscale pickup ion phase-mixing in the lunar plasma wake with a hybrid model is the main subject of our investigation in this paper. Photoionization and charge exchange of protons with the lunar exosphere are the ionization processes included in our model. The computational model includes the self-consistent dynamics of the light (hydrogen ions or hydrogen molecule ion or helium ion), and heavy (sodium ion ) pickup ions. The electrons are considered as a fluid. The lunar interior is considered as a weakly conducting body. In this paper we considered for the first time the cumulative effect of heavy neutrals in the lunar exosphere (e.g., Aluminum, Argon), an effect which was simulated with one species of sodium ion but with a tenfold increase in total production rates. We find that various species produce various types of plasma tail in the lunar plasma wake. Specifically, sodium ion and helium ion pickup ions form a cycloid-like tail, whereas the hydrogen ion or hydrogen molecule ion pickup ions form a tail with a high density core and saw-like periodic structures in the flank region. The length of these structures varies from 1:5 R(sub M) to 3:3 R(sub M) depending on the value of gyro radius for hydrogen ion or hydrogen molecule ion pickup ions. The light pickup ions produce more symmetrical jump in the density and magnetic field at the Mach cone which is mainly controlled by the conductivity of the interior, an effect previously unappreciated. Although other pickup ion species had little effect on the nature of the interaction of the Moon with the solar wind, the global structure of the lunar tail in these simulations appeared quite different when the hydrogen molecule ion production rate was high.

Lipatov, A.S.↗