Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic structure factor”

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 739 records · Page 41

ELSI — An open infrastructure for electronic structure solvers

Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn–Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. As a result, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.

97 MATHEMATICS AND COMPUTING↗

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion and Adsorption of 2-Methylpentane and 3-Methylpentane in MFI-Type Zeolite Crystals

Due to the intimate contact of guest molecules with the inner surface of microporous adsorbents, slight changes in the structure of the guest molecules may lead to remarkable changes in their microdynamics. As an impressive example, in this work, diffusion measurements via IR microimaging (IRM) show that the transport diffusivity of 2-methylpentane (2MP) exceeds the value attained for 3-methylpentane (3MP) by a factor of 2−4.5 in silicalite- 1, an MFI-type zeolite crystal and a technically relevant porous host system. This experimental finding is rationalized by comparison with the outcome of dynamically corrected transition state theory (dcTST) simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ThO 2 and Th 1– x U x O 2 Nanoscale Materials and Thin Films for Nuclear Science Applications

This study investigates the dynamics and mechanisms of solution combustion synthesis (SCS) for the preparation of nanoscale ThO 2 and Th 1–x U x O 2 materials, utilizing metal nitrates (Th(NO 3 ) 4 and UO 2 (NO 3 ) 2 ) and acetylacetone (C 5 H 8 O 2 ) as reactants dissolved in a 2-methoxyethanol (C 3 H 8 O 2 ) solvent. By combining thermodynamic calculations, dynamic time–temperature profile measurements with differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA), this research reveals how variations in acetylacetone concentration and uranium content influence the structural parameters of the synthesized oxides. The time–temperature measurements show that the heating rate and maximum combustion temperatures are sensitive to acetylacetone concentration. DSC-TGA results indicate shifts in exothermic peak temperatures as the uranium content changes. The complexation between thorium and acetylacetone emerges as a critical factor, impacting combustion parameters and the structural characteristics of the final products. The uniform distribution of Th and U in the Th 1–x U x O 2 solid solution and the formation of nanoscale particles with strained crystallites are considered essential for the low-temperature densification of these materials for nuclear fuel pellet applications. Additionally, high-quality ThO 2 and Th 1–x U x O 2 thin (100–150 nm) films are successfully synthesized via electrospray deposition of combustible solutions followed by a brief period of heat treatment. Furthermore, these films exhibit excellent structural and morphological uniformity, making them ideal candidates for nuclear measurements, irradiation damage studies, and investigations into the physical properties of both pure and mixed oxides.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of a Near-Threshold Proton Resonance in 11 B

A near-threshold proton resonance in 11 B at 𝐸 ex = 11.44 ± 0.04 MeV is observed via the reaction 10 Be⁢(𝑑,𝑛)⁢ 11 B → 10 Be + 𝑝 in inverse kinematics, measured with a beam of the radioactive isotope 10 Be. The resonance energy at 𝐸 res = 211⁢(40) keV is consistent with a proton signal observed by Ayyad et al. in the 𝛽-delayed proton decay of 11 Be. By comparison to a distorted wave Born approximation calculation, a 0.27(6) spectroscopic factor is extracted and a tentative (ℓ = 0) character is assigned for this resonance. The significant cross section in the proton-transfer (𝑑,𝑛) reaction, as well as the observation of its proton-decay signal, point to the threshold-resonance character of this state. The position of this state, its structure, and strong coupling to the 𝑠-wave continuum represent an ideal case to study quantum near-threshold many-body dynamics of unstable states. In conclusion, the presence of this state is an important step toward understanding the excessively large beta-delayed proton-decay branch of 11 Be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nanoklystron: A Monolithic Tube Approach to THz Power Generation

The authors propose a new approach to THz power generation: the nanoklystron. Utilizing silicon micromachining techniques, the design and fabrication concept of a monolithic THz vacuum-tube reflex-klystron source is described. The nanoklystron employs a separately fabricated cathode structure composed of densely packed carbon nanotube field emitters and an add-in repeller. The nanotube cathode is expected to increase the current density, extend the cathode life and decrease the required oscillation voltage to values below 100 V. The excitation cavity is based on ridged-waveguide and differs from the conventional cylindrical re-entrant structures found in lower frequency klystrons. A quasi-static field analysis of the cavity and output coupling structure show excellent control of the quality factor and desired field distribution. Output power is expected to occur through an iris coupled matched rectangular waveguide and integrated pyramidal feed horn. The entire circuit is designed so as to be formed monolithically from two thermocompression bonded silicon wafers processed using deep reactive ion etching (DRIE) techniques. To expedite prototyping, a 600 GHz mechanically machined structure has been designed and is in fabrication. A complete numeric analysis of the nanoklystron circuit, including the electron beam dynamics has just gotten underway. Separate evaluation of the nanotube cathodes is also ongoing. The authors will describe the progress to date as well as plans for the immediate implementation and testing of nanoklystron prototypes at 640 and 1250 GHz.

Siegel, Peter H.↗

An MSK Radar Waveform

The minimum-shift-keying (MSK) radar waveform is formed by periodically extending a waveform that separately modulates the in-phase and quadrature- phase components of the carrier with offset pulse-shaped pseudo noise (PN) sequences. To generate this waveform, a pair of periodic PN sequences is each passed through a pulse-shaping filter with a half sinusoid impulse response. These shaped PN waveforms are then offset by half a chip time and are separately modulated on the in-phase and quadrature phase components of an RF carrier. This new radar waveform allows an increase in radar resolution without the need for additional spectrum. In addition, it provides self-interference suppression and configurable peak sidelobes. Compared strictly on the basis of the expressions for delay resolution, main-lobe bandwidth, effective Doppler bandwidth, and peak ambiguity sidelobe, it appears that bi-phase coded (BPC) outperforms the new MSK waveform. However, a radar waveform must meet certain constraints imposed by the transmission and reception of the modulation, as well as criteria dictated by the observation. In particular, the phase discontinuity of the BPC waveform presents a significant impediment to the achievement of finer resolutions in radar measurements a limitation that is overcome by using the continuous phase MSK waveform. The phase continuity, and the lower fractional out-of-band power of MSK, increases the allowable bandwidth compared with BPC, resulting in a factor of two increase in the range resolution of the radar. The MSK waveform also has been demonstrated to have an ambiguity sidelobe structure very similar to BPC, where the sidelobe levels can be decreased by increasing the length of the m-sequence used in its generation. This ability to set the peak sidelobe level is advantageous as it allows the system to be configured to a variety of targets, including those with a larger dynamic range. Other conventionally used waveforms that possess an even greater spectral efficiency than the MSK waveform, such as linear frequency modulation (LFM) and Costas frequency hopping, have a fixed peak sidelobe level that is therefore not configurable, and can be exceeded by high contrast targets. Furthermore, in the case of a multistatic experiment observing a target in motion, self-interference from the transmitter to the receiver is mitigated by the MSK waveform. Waveforms that have delay Doppler coupling, such as LFM, provide no such protection.

Quirk, Kevin J.↗

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗

Tuning microwave losses in superconducting resonators

Abstract The performance of superconducting resonators, particularly cavities for particle accelerators and micro cavities and thin film resonators for quantum computations and photon detectors, has been improved substantially by recent material treatments and technological advances. As a result, the niobium cavities have reached the quality factors Q ∼ 10 11 at 1–2 GHz and 1.5 K and the breakdown radio-frequency (rf) fields H close to the dc superheating field of the Meissner state. These advances raise the questions of whether the state-of-the-art cavities are close to the fundamental limits, what these limits actually are, and to what extent the Q and H limits can be pushed by the materials nano structuring and impurity management. These issues are also relevant for many applications using high-Q thin film resonators, including single-photon detectors and quantum circuits. This topical review outlines basic physical mechanisms of the rf nonlinear surface impedance controlled by quasiparticles, dielectric losses and trapped vortices, as well as the dynamic field limit of the Meissner state. Sections cover methods of engineering an optimum quasiparticle density of states and superfluid density to reduce rf losses and kinetic inductance by pairbreaking mechanisms related to magnetic impurities, rf currents, and proximity-coupled metallic layers at the surface. A section focuses on mechanisms of residual surface resistance, which dominates rf losses at ultra low temperatures. Microwave losses of trapped vortices and their reduction by optimizing the concentration of impurities and pinning potential are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multiple Pathways of Influence for Tightly and Loosely Structured Organizations: Implications for Systems Resilience

Organizations play a key role in supporting various societal functions, ranging from environmental governance to the manufacturing of goods. Here, the behaviors of organization are impacted by various influences, including information, technology, authority, economic leverage, historical experiences, and external factors, such as regulations. This paper introduces a generalized framework, focused on the relative structure of an organization (tight vs. loose), that can be used to understand how different influence pathways can impact decision-making within differently structured organizations. This generalized framework is then translated into a modeling and simulation platform to support and assess implications of these structural differences in resilience to disinformation (measured by organizational behaviors of timeliness and inclusion of quality information) using a systems dynamics approach Preliminary results indicate that a tightly structured organization may be less timely at processing information but could be more resilient against using poor quality information in organizational decisions compared to a loosely structured organization. Ongoing work is underway to understand the robustness of these findings and to validate current model design activities with empirical insights.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Servo-aero-gravo-elastic (SAGE) scaling and its application to a 13-MW downwind turbine

Reduced scale wind turbines can be extremely cost-effective to test new rotor concepts since prototype costs are heavily dependent on the rotor diameter. Ideally, the scaled model would have the same non-dimensional deflections, dynamics, and control behavior as the full-scale model. This would provide a high-fidelity demonstration of the full-scale performance, which is ideal if the full-scale turbine has significant aeroelastic interactions. To this end, servo-aero-gravo-elastic (SAGE) scaling is developed and applied to a 13-MW turbine that is scaled to a 20% scale model. The scaling preserves the tip-speed ratio, the rotor speed normalized by the flapping frequency, and the tip deflections normalized by the blade length. In addition, the controller employs the same control structure (gain-scheduled pitch control and variable speed torque control) and is scaled dynamically (e.g., matching non-dimensional time constant of the pitch angle, etc.). Furthermore, the thrust, gravity, and centrifugal moments are scaled such that the load angles are preserved as a function of a non-dimensional wind speed. However, the environmental scaling must consider differences in Reynolds number (since this parameter cannot be held constant) and subsequent changes in the axial induction factor. While the presented results showcase these differences during operational conditions, the non-dimensional tip deflections remain comparable through all wind speed ranges, indicating the viability of the SAGE scaling method in matching full-scale aeroelastic responses.

Kaminski, Meghan (ORCID:0000000322334784)↗

A three-dimensional model of co-rotating streams in the solar wind. 2: Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are explored by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single-fluid, three-dimensional (3-D) hydrodynamic formulation. A simple, hypothetical 3-D stream structure is defined on a source surface located at 35 solar radius and carefully documents its evolution to 1 AU under the influence of solar rotation. By manipulating the structure of this prototype configuration at the source surface, it is possible to elucidate the factors most strongly affecting stream evolution: (1) the intrinsic correlations among density, temperature, and velocity existing near the source; (2) the amplitude of the stream; (3) the longitudinal breadth of the stream; (4) the latitudinal breadth of the stream; and (5) the heliographic latitude of the centroid of the stream.

Pizzo, V. J.↗

The Atmospheric Dynamics of alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for alpha Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from alpha Tau can be consistently understood via a model of upward-traveling Alfven waves in a gravitationally stratified atmosphere. These waves cause non-thermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

The Atmospheric Dynamics of alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for alpha Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from alpha Tau can be consistently understood via a model of upward-traveling Alfven waves in a gravitationally stratified atmosphere. These waves cause non-thermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

The Atmospheric Dynamics of Alpha Tau (K5 III) - Clues to Understanding the Magnetic Dynamo in Late-Type Giant Stars

Using HST/GHRS, HST/STIS and FUSE archival data for a Tau and the CHIANTI spectroscopic code, we have derived line shifts, volumetric emission measures, and plasma density estimates, and calculated filling factors for a number of UV lines forming between 10,000 K and 300,000 K in the outer atmosphere of this red giant star. The data suggest the presence of low-temperature extended regions and high-temperature compact regions, associated with magnetically open and closed structures in the stellar atmosphere, respectively. The signatures of UV lines from a Tau can be consistently understood via a model of upward-traveling Alfv6n waves in a gravitationally stratified atmosphere. These waves cause nonthermal broadening in UV lines due to unresolved wave motions and downward plasma motions in compact magnetic loops heated by resonant Alfven wave heating.

Carpenter, Kenneth G.↗

Flow field topology of transient mixing driven by buoyancy

Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.

Nonlinear Dynamics↗

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES↗

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗