Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space methods”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Machine-learning-based dynamic-importance sampling for adaptive multiscale simulations

Multiscale simulations are a well-accepted way to bridge the length and time scales required for scientific studies with the solution accuracy achievable through available computational resources. Traditional approaches either solve a coarse model with selective refinement or coerce a detailed model into faster sampling, both of which have limitations. Here, we present a paradigm of adaptive, multiscale simulations that couple different scales using a dynamic-importance sampling approach. Our method uses machine learning to dynamically and exhaustively sample the phase space explored by a macro model using microscale simulations and enables an automatic feedback from the micro to the macro scale, leading to a self-healing multiscale simulation. As a result, our approach delivers macro length and time scales, but with the effective precision of the micro scale. Our approach is arbitrarily scalable as well as transferable to many different types of simulations. Overall, our method made possible a multiscale scientific campaign of unprecedented scale to understand the interactions of RAS proteins with a plasma membrane in the context of cancer research running over several days on Sierra, which is currently the second-most-powerful supercomputer in the world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Semiclassical theory and the Koopman-van Hove equation *

Abstract The phase space Koopman-van Hove (KvH) equation can be derived from the asymptotic semiclassical analysis of partial differential equations. Semiclassical theory yields the Hamilton–Jacobi equation for the complex phase factor and the transport equation for the amplitude. These two equations can be combined to form a nonlinear semiclassical version of the KvH equation in configuration space. There is a natural injection of configuration space solutions into phase space and a natural projection of phase space solutions onto configuration space. Hence, every solution of the configuration space KvH equation satisfies both the semiclassical phase space KvH equation and the Hamilton–Jacobi constraint. For configuration space solutions, this constraint resolves the paradox that there are two different conserved densities in phase space. For integrable systems, the KvH spectrum is the Cartesian product of a classical and a semiclassical spectrum. If the classical spectrum is eliminated, then, with the correct choice of Jeffreys–Wentzel–Kramers–Brillouin (JWKB) matching conditions, the semiclassical spectrum satisfies the Einstein–Brillouin–Keller quantization conditions which include the correction due to the Maslov index. However, semiclassical analysis uses different choices for boundary conditions, continuity requirements, and the domain of definition. For example, use of the complex JWKB method allows for the treatment of tunneling through the complexification of phase space. Finally, although KvH wavefunctions include the possibility of interference effects, interference is not observable when all observables are approximated as local operators on phase space. Observing interference effects requires consideration of nonlocal operations, e.g. through higher orders in the asymptotic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol–disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol–disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. Furthermore, as an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol–disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.

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

Uniform-in-phase-space data selection with iterative normalizing flows

Improvements in computational and experimental capabilities are rapidly increasing the amount of scientific data that are routinely generated. In applications that are constrained by memory and computational intensity, excessively large datasets may hinder scientific discovery, making data reduction a critical component of data-driven methods. Datasets are growing in two directions: the number of data points and their dimensionality. Whereas dimension reduction typically aims at describing each data sample on lower-dimensional space, the focus here is on reducing the number of data points. A strategy is proposed to select data points such that they uniformly span the phase-space of the data. The algorithm proposed relies on estimating the probability map of the data and using it to construct an acceptance probability. An iterative method is used to accurately estimate the probability of the rare data points when only a small subset of the dataset is used to construct the probability map. Instead of binning the phase-space to estimate the probability map, its functional form is approximated with a normalizing flow. Therefore, the method naturally extends to high-dimensional datasets. The proposed framework is demonstrated as a viable pathway to enable data-efficient machine learning when abundant data are available.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1962

VISUAL CONTROL OF RENDEZVOUS

Space rendezvous involving human occupants will be required in many phases of space flight. In order that rendezvous capability be concurrent with other phases of manned space flight to insure that such tasks as the rotation and addition of crews for manned space station, supply of stations and launching complexes for deep-space probes, and the emergency rescue of personnel may be carried out, it is desirable to determine the ability of a human pilot to control the rendezvous maneuver. The National Aeronautics and Space Administration has conducted a study program to determine vehicle requirements and control techniques for a number of proposed methods for performing piloted space rendezvous. The final step in each phase of the rendezvous study program has been a simulation of the man-machine combination to test its feasibility. Initially, an analytical study was conducted to determine a range of conditions that would exist at the beginning of the terminal phase of rendezvous (ref. 1). Piloted systems were investigated for rendezvous control using these initial conditions. One such study of piloted rendezvous in which all the required control information was displayed to the simulator pilot on instruments is presented in reference 2. While the control of the rendezvous maneuver by use of the instrument display was shown to be feasible, the sufficiently accurate measurement of the necessary information to operate the instrument was realized to be a difficult problem. Consideration was therefore given to the possibility of maneuvering range and range rate required for the control of the rendezvous from visual sightings made during the initial control application. To provide this capability, an analytical study was made. This analytical study, presented in reference 3, derived three techniques for obtaining range and range rate visually and conducted a simplified preliminary situation to test the three analytical techniques. The study reported herein is a continuation of the simulation work begun in reference 3, and is aimed at determining the ability of a pilot to use visual information to reduce or eliminate the need for some or all of the data that would normally be displayed on instruments.

Rendezvous↗
DOE OSTI2024

Topological and magnetic properties of the interacting Bernevig-Hughes-Zhang model

We investigate the effects of electronic correlations on the Bernevig-Hughes-Zhang model using the real-space density matrix renormalization group (DMRG) algorithm. We introduce a method to probe topological phase transitions in systems with strong correlations using DMRG, substantiated by an unsupervised machine learning methodology that analyzes the orbital structure of the real-space edges. Including the full multi-orbital Hubbard interaction term, we construct a phase diagram as a function of a gap parameter (m) and the Hubbard interaction strength (U) via exact DMRG simulations on N×4 cylinders. Our analysis confirms that the topological phase persists in the presence of interactions, consistent with previous studies, but it also reveals an intriguing phase transition from a paramagnetic to a stripey antiferromagnetic topological insulator. The combination of the magnetic structure factor, strength of magnetic moments, and the orbitally resolved density, provides real-space information on both topology and magnetism in a strongly correlated system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2010

Braze Development of Graphite Fiber for Use in Phase Change Material Heat Sinks

Hamilton Sundstrand (HS), together with NASA Johnson Space Center, developed methods to metallurgically join graphite fiber to aluminum. The goal of the effort was to demonstrate improved thermal conductance, tensile strength and manufacturability compared to existing epoxy bonded techniques. These improvements have the potential to increase the performance and robustness of phase change material heat sinks that use graphite fibers as an interstitial material. Initial work focused on evaluating joining techniques from 4 suppliers, each consisting of a metallization step followed by brazing or soldering of one inch square blocks of Fibercore graphite fiber material to aluminum end sheets. Results matched the strength and thermal conductance of the epoxy bonded control samples, so two suppliers were down-selected for a second round of braze development. The second round of braze samples had up to a 300% increase in strength and up to a 132% increase in thermal conductance over the bonded samples. However, scalability and repeatability proved to be significant hurdles with the metallization approach. An alternative approach was pursued which used nickel and active braze allows to prepare the carbon fibers for joining with aluminum. This approach was repeatable and scalable with improved strength and thermal conductance when compared with epoxy bonding.

Quinn, Gregory↗
NASA NTRS2011

Braze Development of Graphite Fiber for Use in Phase Change Material Heat Sinks

Hamilton Sundstrand (HS), together with NASA Johnson Space Center, developed methods to metallurgically join graphite fiber to aluminum. The goal of the effort was to demonstrate improved thermal conductance, tensile strength and manufacturability compared to existing epoxy bonded techniques. These improvements have the potential to increase the performance and robustness of phase change material heat sinks that use graphite fibers as an interstitial material. Initial work focused on evaluating joining techniques from four suppliers, each consisting of a metallization step followed by brazing or soldering of one inch square blocks of Fibercore graphite fiber material to aluminum end sheets. Results matched the strength and thermal conductance of the epoxy bonded control samples, so two suppliers were down-selected for a second round of braze development. The second round of braze samples had up to a 300% increase in strength and up to a 132% increase in thermal conductance over the bonded samples. However, scalability and repeatability proved to be significant hurdles with the metallization approach. An alternative approach was pursued which used a nickel braze allow to prepare the carbon fibers for joining with aluminum. Initial results on sample blocks indicate that this approach should be repeatable and scalable with good strength and thermal conductance when compared with epoxy bonding.

Quinn, Gregory↗
DOE OSTI2024

Retrofittable Thermal Switches for Dynamic Building Envelopes Integrated with Thermal Energy Storage: Preprint

Buildings in the United States consume about 40 quadrillion BTU of primary energy annually, which accounts for the nation's 40% of total energy use, 75% of all electricity use, and 35% of the net carbon emissions. Deploying thermal energy storage in the form of phase change material (PCM) in building envelopes is an effective method to reduce space heating/cooling loads, provide load shedding, and shift demand to periods of lower energy cost. However, the full potential of PCM-integrated envelopes can only be realized if the PCM undergoes complete phase change using free ambient heating/cooling, and the stored energy is effectively transferred between the exterior and the interior environments. Conventional thermal insulation (with a fixed thermal resistance) limits PCM utilization, particularly with the increasing emphasis on higher R-value in building envelopes, which negatively affects the energy-saving potential of a PCM-integrated envelope. In contrast, dynamic building envelopes integrated with PCMs provide the option of varying the thermal resistance based on the indoor and outdoor conditions, thereby enhancing utilization of free ambient cooling and heating to charge/discharge the PCM thermal storage, reducing the buildings' heating and cooling load, and shifting the peak energy demand. In this study, we demonstrate innovative retrofittable thermal switches in the form of the insertable plugs inside an insulation to provide variable thermal resistance depending on the operating temperature and direction of temperature gradient, thus allowing preferential directional heat flow. Notably, they are passive in nature, requiring no external power, and work solely based on the ambient temperature.

buildings↗
DOE OSTI2026

Vapor-phase pillarization of MXenes for engineering hierarchical interlayer porosity

MXenes, a family of two-dimensional (2D) multilamellar materials, possess excellent thermal and electronic properties for a range of applications. Their use in heterogeneous catalysis, however, is limited by the low surface area resulting from stacked layers. Pillarization with inorganic oxides can create more open, mesoporous MXene structures, improving accessibility for guest species to diffuse, reside or react in the space between 2D layers. A previous liquid-phase pillarization method, however, involves excessive use of solvent-based precursors and multiple processing steps. Here, we report a vapor-phase pillarization (VPP) strategy to introduce pillars, exemplified by silica pillars, with high pillar precursor usage efficiency and a simplified processing workflow. The resulting silica-pillared mesoporous MXene exhibits significantly increased surface area and porosity. These textural properties can be easily tuned by the VPP synthesis conditions. When applied as a ruthenium (Ru) catalyst support for the hydrogenolysis of low-density polyethylene (LDPE), the silica-pillared MXene enabled high Ru dispersion and catalytic activity. This study highlights the potential of the VPP method for engineering mesoporous, 2D MXene materials and demonstrates the effectiveness of mesoporous MXene as a catalyst support in overcoming mass transport and active-site accessibility challenges in heterogeneous catalysis involving bulky substances, such as plastics upcycling.

Luo, Song [University of Delaware, Newark, DE (Uni↗
DOE OSTI2026

Enhanced coercivity in Fe5C2/SiO2 core/shell nanocrystals

Rod-shaped Fe5C2 and core/shell Fe5C2/SiO2 nanocrystals were synthesized via a solution-based chemical method. Structural analysis confirmed the monoclinic phase of Fe5C2 with space group C2/c. Zero-field-cooling (ZFC) and field-cooling (FC) magnetization curves revealed distinct magnetic behaviors: uncoated Fe5C2 exhibited a low-temperature FC plateau indicative of strong dipolar interactions, while Fe5C2/SiO2 showed a monotonic increase in FC magnetization, suggesting reduced dipolar interactions due to SiO2 surface passivation. Isothermal remanent magnetization (IRM) and DC demagnetization (DCD) measurements supported this trend, with δM plots confirming weaker dipolar interactions in the coated sample. Bloch’s law fitting of temperature-dependent saturation magnetization showed a smaller Bloch’s constant for pure Fe5C2 and a larger value for Fe5C2/SiO2, reflecting enhanced surface disorder and reduced exchange coupling in the latter. Notably, Fe5C2/SiO2 demonstrated increased coercivity, attributed to decreased dipolar interaction and elevated surface anisotropy. Kneller’s law fitting yielded higher blocking temperatures for Fe5C2 (476 K) than Fe5C2/SiO2 (456 K), highlighting the impact of dipolar interactions on magnetic relaxation. These findings illustrate how SiO2 coatings effectively modulate dipolar interactions and enhance coercivity in Fe5C2 nanocrystals.

Joshi, Pramanand [Department of Physics, Universit↗
DOE OSTI2021

Topological Phase Transition and Phonon-Space Dirac Topology Surfaces in ZrTe5

We use first-principles methods to demonstrate that, in ZrTe 5 , a layered van der Waals material like graphite, atomic displacements corresponding to five of the six zone-center A g (symmetry-preserving) phonon modes can drive a topological transition from a strong to a weak topological insulator with a Dirac semimetal state emerging at the transition, giving rise to a Dirac topology surface in the multidimensional space formed by the A g phonon modes. This implies that the topological transition in ZrTe 5 , can be realized with many different settings of external stimuli capable of penetrating through the phonon-space Dirac surface without breaking the crystallographic symmetry. Furthermore, we predict that domains with effective mass of opposite signs can be created by laser pumping and will host Weyl modes of opposite chirality propagating along the domain boundaries. Studying phonon-space topology surfaces provides a new route to understanding and utilizing the exotic physical properties of ZrTe 5 , and related quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

Radial diffusion in Saturn's radiation belts - A modeling analysis assuming satellite and ring E absorption

A modeling analysis is carried out of six experimental phase space density profiles for nearly equatorially mirroring protons using methods based on the approach of Thomsen et al. (1977). The form of the time-averaged radial diffusion coefficient D(L) that gives an optimal fit to the experimental profiles is determined under the assumption that simple satellite plus Ring E absorption of inwardly diffusing particles and steady-state radial diffusion are the dominant physical processes affecting the proton data in the L range that is modeled. An extension of the single-satellite model employed by Thomsen et al. to a model that includes multisatellite and ring absorption is described, and the procedures adopted for estimating characteristic satellite and ring absorption times are defined. The results obtained in applying three representative solid-body absorption models to evaluate D(L) in the range where L is between 4 and 16 are reported, and a study is made of the sensitivity of the preferred amplitude and L dependence for D(L) to the assumed model parameters. The inferred form of D(L) is then compared with that which would be predicted if various proposed physical mechanisms for driving magnetospheric radial diffusion are operative at Saturn.

Hood, L. L.↗
NASA NTRS1984

Geometric accuracy assessment of Landsat-4 Thematic Mapper p-tapes

An analysis was performed on the P-format computer compatible tapes for the Landsat-4 Thematic Mapper (TM) scene 40109-15140 acquired November 2, 1982 over Washington, D.C. and environs. Three tests of sensor geometry were undertaken: (1) band-to-band registration, (2) line-to-line registration between swaths, and (3) geometric correction of sensor and spacecraft/ephemeris characteristics. The band-to-band and line-to-line registration was measured at one hundred pixel spacings along a line using the phase correlation image alignment method developed by Kuglin and Hines and adapted to a one-dimensional FFT correlation technique. The sensor and spacecraft geometric calibration analysis was checked by identifying 75 ground control points in the TM scene and on 1:24,000 topographic maps. A least squares fit of the gaps was computed and the vector offsets of the residuals was tabulated and plotted.

Bryant, N. A.↗
DOE OSTI2022

Electron Cooling Experiment for Proton Beams with Intense Space-Charge in IOTA

Electron cooling as a method of creating intense ion beams has a practical upper limit when it comes to the peak phase space density of ion beams which can be achieved in practice. We describe a new experiment to study electron cooling of 2.5 MeV protons at the intensity limit using the Integrable Optics Test Accelerator (IOTA), which is a storage ring dedicated to beam physics research at Fermilab. This system will enable the study of magnetized electron cooling of a proton beam with transverse incoherent tune shifts approaching -0.5 due to the presence of intense space-charge forces. We present an overview of the hardware design, simulations and specific experiments planned for this project.

43 PARTICLE ACCELERATORS↗
NASA NTRS2009

Estimating Transmitted-Signal Phase Variations for Uplink Array Antennas

A method of estimating phase drifts of microwave signals distributed to, and transmitted by, antennas in an array involves the use of the signals themselves as phase references. The method was conceived as part of the solution of the problem of maintaining precise phase calibration required for proper operation of an array of Deep Space Network (DSN) antennas on Earth used for communicating with distant spacecraft at frequencies between 7 and 8 GHz. The method could also be applied to purely terrestrial phased-array radar and other radio antenna array systems. In the DSN application, the electrical lengths (effective signal-propagation path lengths) of the various branches of the system for distributing the transmitted signals to the antennas are not precisely known, and they vary with time. The variations are attributable mostly to thermal expansion and contraction of fiber-optic and electrical signal cables and to a variety of causes associated with aging of signal-handling components. The variations are large enough to introduce large phase drifts at the signal frequency. It is necessary to measure and correct for these phase drifts in order to maintain phase calibration of the antennas. A prior method of measuring phase drifts involves the use of reference-frequency signals separate from the transmitted signals. A major impediment to accurate measurement of phase drifts over time by the prior method is the fact that although DSN reference-frequency sources separate from the transmitting signal sources are stable and accurate enough for most DSN purposes, they are not stable enough for use in maintaining phase calibrations, as required, to within a few degrees over times as long as days or possibly even weeks. By eliminating reliance on the reference-frequency subsystem, the present method overcomes this impediment. In a DSN array to which the present method applies (see figure), the microwave signals to be transmitted are generated by exciters in a signal-processing center, then distributed to the antennas via optical fibers. At each antenna, the signals are used to drive a microwave power-amplifier train, the output of which is coupled to the antenna for transmission. A small fraction of the power-amplifier-train output is sent back to the signal-processing center along another optical fiber that is part of the same fiber-optic cable used to distribute the transmitted signal to the antenna. In the signal-processing center, the signal thus returned from each antenna is detected and its phase is compared with the phase of the signal sampled directly from the corresponding exciter. It is known, from other measurements, that the signal-propagation path length from the power-amplifier-train output port to the phase center of each antenna is sufficiently stable and, hence, that sampling the signal at the power-amplifier-train output port suffices for the purpose of characterizing the phase drift of the transmitted signal at the phase center of the antenna

Paal, Leslie↗
DOE OSTI2023

Evaluation protocol for revealing magnonic contrast in TR-STXM measurements

We present a statistically motivated method to extract magnonic contrast from time-resolved scanning transmission x-ray microscopy (TR-STXM) measurements. TR-STXM is an element-specific method for resolving spin-dynamics in space and time. It offers nanometer spatial resolution and picosecond temporal resolution. The presented method makes it possible to obtain phase and amplitude profiles of spin-waves from STXM measurements. Furthermore, it allows for a rigorous transformation to reciprocal magnon k⃗-space, revealing k⃗-dependent magnon properties such as the magnon dispersion in three dimensions and for all directions of the magnetic anisotropy. We demonstrate our method using X-band ferromagnetic resonance on a micrometer-sized permalloy assembly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗