Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space physics”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Novel Spatio-Temporal Regime Tracking Method for Impact Simulations

In this proposal, we present a novel method of tracking the rheological regimes activated during impact cratering events that will allow researchers to gain new insights into cratering mechanics. Rheology describes the stress-strain response of rocks to different conditions. Planetary impact cratering events often occur on too large of a scale to be feasibly captured in controlled experiments. Instead, these dynamic events are primarily studied using multi-physics codes equipped with complex material models that enable calculations of the impact event at scale. However, determining which physical processes are required for the problem of interest is challenging. Because the dominant rheological regimes change with space and time during crater formation, it is difficult to link numerical simulations with observable features of craters at the end of the event. The basis of this work is the implementation of numerical flags that track the activation of each rheological regime throughout impact simulations. We demonstrate this with the ‘Rock Model’ implemented in the CTH shock-physics code. We will use this rheology tracking method to ’zoom in’ on a specific region within an event and track the conditions the rock experiences over time. This work will develop community benchmarks to validate and distribute our implemented rheological models. We will focus on improving the melt models used by the planetary impact modeling community by developing an EOS-aware rheological transition from solid to melt. Through analysis of cell and tracer-particle based tracking data, we will study the effects of different rheologies on modeled outcomes, particularly on the volume and distributions of impacts melts. We will also use this method to link observable features with the rheological mechanisms responsible for them. The deliverables (peer-reviewed papers) from the proposed work are (i) tests of the implemented rheologic processes and demonstrations of the tracking flags; (ii) the first calculations of the spatio-temporal evolution of the dominant rheologies during impact cratering events; and (iii) application to delivery of impactor iron during basin-scale impacts.

58 GEOSCIENCES↗
DOE OSTI2023

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Design and construction of the CMS Outer Tracker for the Phase-2 Upgrade

The High Luminosity LHC (HL-LHC) is expected to deliver an integrated luminosity of $3000-4000\;$fb$^{-1}$ after 10 years of operation with peak instantaneous luminosity reaching about $5-7.5\times10^{34}$cm$^{-2}$s$^{-1}$. During Long Shutdown 3, several components of the CMS detector will undergo major changes, called Phase-2 upgrades, to be able to operate in the challenging environment of the HL-LHC. The current CMS silicon strip tracker has to be replaced with a new detector. The Phase-2 Outer Tracker (OT) will have high radiation tolerance, higher granularity and the capability to handle higher data rates. Another key feature of the OT will be to provide tracking information to the Level-1 trigger, allowing trigger rates to be kept at a sustainable level without sacrificing physics potential. For this, the OT will be made out of modules which have two closely spaced sensors read out by front-end ASICs, which can correlate hits in the two sensors creating short track segments, called stubs. The stubs will be used for tracking in the L1 track finder. In this contribution, the design of the CMS Phase-2 OT, the technological choices and first results with prototype devices are reported.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Learning PDFs through interpretable latent representations in Mellin space

Representing the parton distribution functions (PDFs) of the proton and other hadrons through flexible, high-fidelity parametrizations has been a long-standing goal of particle physics phenomenology. This is particularly true since the chosen parametrization methodology can play an influential role in the ultimate PDF uncertainties as extracted in QCD global analyses; these, in turn, are often determinative of the reach of experiments at the LHC and other facilities to nonstandard physics, including at large 𝑥, where parametrization effects can be significant. In this study, we explore a series of encoder-decoder machine-learning (ML) models with various neural-network topologies as efficient means of reconstructing PDFs from meaningful information stored in an interpretable latent space. Given recent effort to pioneer synergies between QCD analyses and lattice-gauge calculations, we formulate a latent representation based on the behavior of PDFs in Mellin space, i.e., their integrated moments, and test the ability of various models to decode PDFs from this information faithfully. We introduce a numerical package, PDFdecoder, which implements several encoder-decoder models to reconstruct PDFs with high fidelity and use this end-to-end tool to explore how such neural-network-based models might connect PDF parametrizations to underlying properties like their Mellin moments. We additionally dissect patterns of learned correlations between encoded Mellin moments and reconstructed PDFs that suggest opportunities for further improvements to ML-based approaches to PDF parametrizations and uncertainty quantification.

Machine learning↗
DOE OSTI2026

GP Cosmology Surrogate v1.0

GP Cosmology Surrogate is a Python library for building and training a generalized multi-output Gaussian process (GP) framework of @takhtaganov2021cosmic. In this approach, the surrogate is constructed sequentially, guided by a Bayesian optimization acquisition function that targets reduction of emulation error in the regions most consistent with the observational data. This adaptive design concentrates computational resources where they have the greatest impact on inference accuracy. The library supports efficient training for separable GP kernels, which allows the use of Kronecker algebra to handle high-dimensional input spaces and large numbers of correlated outputs. This makes it well suited for applications such as modeling cosmological power spectra, large-scale physical simulations, and multi-output hyperparameter tuning. By combining scalable multi-output GP modeling with data-driven adaptive sampling, GPsurrogate enables parameter inference and optimization with substantially fewer simulations than conventional space-filling designs.

Lukic, Zarija [Lawrence Berkeley National Laborato↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Magnetic Reconnection in High Energy Density Laboratory Plasmas

Understanding the physics of collisionless, magnetized plasmas is key to unlocking a number of important problems in space and astrophysics. Magnetic reconnection is a fundamental mechanism which can explosively release stored magnetic energy, convert it to plasma heat and flows, and accelerate particles. In this project we developed technique and conducted experiments and simulations of magnetic reconnection in laser produces plasmas. Through this project, we also developed experimental techniques to generate magnetized collisionless shocks, and we therefore expanded the scope of the project to develop these experiments. These projects make fundamental progress on understanding the dynamics of magnetized laser-produced plasmas. In additional to plasma astrophysics, these results have broad interest for laboratory high energy density plasmas, since magnetic fields can be used to confine and control plasmas in many contexts and may offer new capabilities toward laboratory ignited plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Physics basis for design of 3D coils in tokamaks

Recent progress in 3D tokamak modeling is now leveraged to create a conceptual design of new external 3D field coils for the DIII-D tokamak. In this work generalized perturbed equilibrium code is used to determine optimally efficient spectrum for driving total, core, and edge neoclassical toroidal viscosity torque. These fundamental modes of 3D control are shown to have consistent outboard structures across a wide variety of plasma scenarios and machines. Given these target spectra, the currents and 3D geometry of multiple coils can be optimized to increase efficient drive for the physics of interest without undesired secondary effects. Here, this nonlinear optimization is demonstrated using the flexible optimized coils using space-curves code. The optimized coils are individually distorted in space, creating toroidal 'arrays' containing a variety of shapes that often wrap around a significant poloidal extent of the machine. Importantly, efficient coupling can be maintained even when enforcing large distances between coils and the plasma during the geometric optimization of coil designs. Finally, the physics-driven optimization presented here thus provides a practical path to utilizing coils built on the exterior of the vacuum chamber in future reactors to obtain the powerful 3D field benefits demonstrated on current machines with close, internal coils.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Possible routes to a strong science case for a joint CLEoPATRA-Roman microlensing survey

The Contemporaneous LEnsing Parallax and Autonomous TRansient Array (CLEoPATRA) is a mis sion concept for a space-based microlensing survey, proposed to enhance science for the Roman Space Telescope (RST) microlensing survey. From CLEoPATRA’s planned vantage point in geosynchronous equatorial orbit (GEO), it will observe microlensing events simultaneously and from a different point of view to RST at L2. This joint L2-GEO observation baseline gives a projected angular separation of order the size of the characteristic angular lensing scale of Free-Floating-Planets (FFPs) in the Galaxy, making it powerful for constraining FFP properties. We find that a strong science case will demonstrate which parts of the FFP parameter space a joint RST-CLEoPATRA survey will give lead ing lens constraints compared with competing ground- and space- based RST observatory companions (e.g., PRIME, Euclid). Both the simulations and framework that connect telescope capabilities to physical constraints already exist for RST FFPs and can be used to demonstrate the power of joint RST-CLEoPATRA survey, with the addition of a simulated CLEoPATRA event blending distribution. Population-level FFP constraints could also be calculated by adapting ongoing black hole microlensing research at Lawrence Livermore National Laboratory. We also find that considering CLEoPATRA as a third space-based observatory (in addition to RST and Euclid/PRIME), which will break microlensing degeneracies leading to unambiguous FFP constraints, is worth investigating as part of a science case. Finally, we highlight possible auxiliary CLEoPATRA science cases including relieving tension between the RST Core Community Surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Parameter inference from event ensembles and the top-quark mass

One of the key tasks of any particle collider is measurement. In practice, this is often done by fitting data to a simulation, which depends on many parameters. Sometimes, when the effects of varying different parameters are highly correlated, a large ensemble of data may be needed to resolve parameter-space degeneracies. An important example is measuring the top-quark mass, where other physical and unphysical parameters in the simulation must be profiled when fitting the top-quark mass parameter. We compare four different methodologies for top-quark mass measurement: a classical histogram fit similar to one commonly used in experiment augmented by soft-drop jet grooming; a 2D profile likelihood fit with a nuisance parameter; a machine-learning method called DCTR; and a linear regression approach, either using a least-squares fit or with a dense linearly-activated neural network. Despite the fact that individual events are totally uncorrelated, we find that the linear regression methods work most effectively when we input an ensemble of events sorted by mass, rather than training them on individual events. Although all methods provide robust extraction of the top-quark mass parameter, the linear network does marginally best and is remarkably simple. For the top study, we conclude that the Monte-Carlo-based uncertainty on current extractions of the top-quark mass from LHC data can be reduced significantly (by perhaps a factor of 2) using networks trained on sorted event ensembles. More generally, machine learning from ensembles for parameter estimation has broad potential for collider physics measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Finding wombling boundaries in LHC data with Voronoi and Delaunay tessellations

We address the problem of finding a wombling boundary in point data generated by a general Poisson point process, a specific example of which is an LHC event sample distributed in the phase space of a final state signature, with the wombling boundary created by some new physics. We discuss the use of Voronoi and Delaunay tessellations of the point data for estimating the local gradients and investigate methods for sharpening the boundaries by reducing the statistical noise. The outcome from traditional wombling algorithms is a set of boundary cell candidates with relatively large gradients, whose spatial properties must then be scrutinized in order to construct the boundary and evaluate its significance. Here we propose an alternative approach where we simultaneously form and evaluate the significance of all possible boundaries in terms of the total gradient flux. We illustrate our method with several toy examples of both straight and curved boundaries with varying amounts of signal present in the data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Traversable wormhole solution with a background Kalb–Ramond field

Highlights: • Wormhole solutions of gravity in a Lorentz breaking symmetry scenario are found. • Light deflection by the wormhole in a Lorentz breaking symmetry scenario is calculated. • It is shown the influence of the vacuum expectation value of a Kalb–Ramond field in deflection of light. We obtain a static spherically symmetric wormhole solution due to the vacuum expectation value (vev) of a Kalb–Ramond field. The Kalb–Ramond (vev) is a background tensor field which produces a local Lorentz symmetry breaking (LSB) of spacetime. Considering a non-minimal coupling between the Kalb–Ramond (vev) and the Ricci tensor, we found an exact traversable wormhole solution sustained by matter sources with a negative isotropic pressure. The matter source satisfies the energy conditions at the throat for particular choices of the LSB parameter. Moreover, we employ the Gauss–Bonnet method to analyze the deflection angle of light in the weak limit approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Switching from batch to continuous granulation: A case study of metoprolol succinate ER tablets

Continuous manufacturing (CM) has been used to produce several immediate release drug products. No extended-release (ER) product manufactured employing CM technology has been approved yet. Herein this study investigated the critical aspects of switching from the batch mode of high shear granulation to the continuous operation of twin-screw granulation for extended-release tablets. Metoprolol succinate ER tablets was used as a model ER formulation for this purpose. A central composite design (CCD) was employed to determine the effects of high shear granulator (HSG) parameters, namely impeller speed, granulation time, and binder liquid feeding rate, on the critical granulation characteristics important for product performance. These critical granulation characteristics served as a guide for switching from the batch processing to the continuous operation for achieving the same breaking strength and dissolution for this ER metoprolol tablets. The granulation time was the most critical factor affecting the bulk properties of granules which contributed to tablet dissolution. The higher density and lower compressibility of granules were attained at the longest granulation time of 5.4 min with the fastest liquid feeding rate of 75 g/min. The granules’ density was the primary factor negatively affecting the dissolution of metoprolol tablets. However, the breaking strength of tablets confounded the effect of granules density on metoprolol dissolution. Switching the processing parameters of high shear granulation to twin-screw granulation achieved similar dissolution profiles (F2 greater than 50). The screw speed was not found to affect bulk properties of granules. The root cause of granulation failures in twin-screw granulation, such as premature consolidation, excessive swelling, poor cohesion, inconsistent shearing effects, and formation of deformed agglomerates, were identified. In conclusion, the use of critical granulation characteristics through a performance-based approach of ER tablets facilitated the switching of manufacturing of an ER formulation form batch to continuous operation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Universal corner symmetry and the orbit method for gravity

A universal symmetry algebra organizing the gravitational phase space has been recently found. It corresponds to the subset of diffeomorphisms that become physical at corners – codimension-2 surfaces supporting Noether charges. It applies to both finite distance and asymptotic corners. In this paper, we study this algebra and its representations, via the coadjoint orbit method. We show that generic orbits of the universal algebra split into sub-orbits spanned by finite distance and asymptotic corner symmetries, such that the full universal symmetry algebra gives rise to a unified treatment of corners in a manifold. We then identify the geometric structure that captures these algebraic properties on corners, which is the Atiyah Lie algebroid associated to a principal GL(2,R) $\ltimes$ R 2 -bundle. This structure is suggestive of the existence of a novel quantum gravitational theory which would unitarily glue such geometric structures, with spacetime geometries appearing as semi-classical configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Effect of Ice Number Concentration on the Evolution of Boundary Layer Clouds During Arctic Marine Cold‐Air Outbreaks

Abstract Marine cold‐air outbreaks (MCAOs) are crucial for Arctic Ocean heat loss, featuring convective cloud rolls that transition into convection cells downstream. Understanding factors controlling this transformation is the key for improving MCAO cloud representation in climate models. This study employs large‐eddy simulations to investigate how cloud ice number concentrations () affect cloud evolution using a case from the Cold‐Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) campaign. The simulations, performed in a Lagrangian framework following an air mass trajectory, are driven by ERA5 reanalysis data. Initially, all simulations produce similar cloud patterns, but higher leads to earlier breakup of cloud rolls. Between 4 and 10 hr, surface precipitation rates are similar across simulations, but precipitation initiates earlier, and the cloud‐base precipitation rates are higher when is higher. The stronger precipitation evaporation leads to increased stability of the boundary layer and reduced intensity of vertical mixing between the surface and cloud layer. An increased sink of cloud layer moisture via precipitation and decreased source through diminished vertical transport result in earlier cloud breakup in higher conditions. Simulations with different sea surface temperatures (SST) indicate that this cloud breakup mechanism remains valid for MCAOs of different strengths, although the cloud organization is more sensitive to SST changes in low environments. This work highlights the importance of accurate representations of ice processes in simulating MCAO clouds and suggests the need for observational constraints of ice nucleating particles and over the mixed‐phase cloud regimes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Neural network-based order parameter for phase transitions and its applications in high-entropy alloys

Phase transition is one of the most important phenomena in nature and plays a central role in materials design. All phase transitions are characterized by suitable order parameters, including the order–disorder phase transition. However, finding a representative order parameter for complex systems is non-trivial, such as for high-entropy alloys. Given the strength of dimensionality reduction of a variational autoencoder (VAE), we introduce a VAE-based order parameter. Here, we propose that the Manhattan distance in the VAE latent space can serve as a generic order parameter for order–disorder phase transitions. The physical properties of our order parameter are quantitatively interpreted and demonstrated by multiple refractory high-entropy alloys. Using this order parameter, a generally applicable alloy design concept is proposed by mimicking the natural mixing process of elements. Our physically interpretable VAE-based order parameter provides a computational technique for understanding chemical ordering in alloys, which can facilitate the development of rational alloy design strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Theoretical prediction of magnetic exchange coupling constants from broken-symmetry coupled cluster calculations

Exchange coupling constants ( J ) are fundamental to the understanding of spin spectra of magnetic systems. In this study, we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. J values calculated from CC in this fashion converge smoothly toward the full configuration interaction result with increasing level of CC excitation. We compare this BS-CC scheme to the complementary equation-of-motion CC approach on a selection of bridged molecular cases and give results from a few other methodologies for context.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Recent developments in the PySCF program package

PySCF is a Python-based general-purpose electronic structure platform that supports first-principles simulations of molecules and solids as well as accelerates the development of new methodology and complex computational workflows. Here, we explain the design and philosophy behind PySCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PySCF as a development environment. We then summarize the capabilities of PySCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PySCF across the domains of quantum chemistry, materials science, machine learning, and quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗