SEARCH · Search NASA

Results for “fully linear 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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Mitigation of the Brighter-Fatter Effect in the LSST Camera

Thick, fully depleted charge-coupled devices (CCDs) are known to exhibit non-linear behavior at high signal levels due to the dynamic behavior of charges collecting in the potential wells of pixels, called the brighter-fatter effect (BFE). This particularly impacts bright calibration stars, which appear larger than their intrinsic shape, creating a flux-dependent point-spread function (PSF) that if left unmitigated, could make up a large fraction of the error budget in Stage IV weak-lensing (WL) surveys such as the Legacy Survey of Space and Time (LSST). In this paper, we analyze image measurements of flat fields and artificial stars taken at different illumination levels with the LSST Camera (LSSTCam) at SLAC National Accelerator Laboratory in order to quantify this effect in the LSST Camera before and after a previously introduced correction technique. We observe that the BFE evolves anisotropically as a function of flux due to higher-order BFEs, which violates the fundamental assumption of this correction method. We then introduce a new sampling method based on a physically motivated model to account these higher-order terms in the correction, and then we test the modified correction on both datasets. We find that the new method corrects the effect in flat fields better than it corrects the effect in artificial stars which we conclude is the result of a unmodeled curl component of the deflection field by the correction. We use these results to define a new metric for the full-well capacity of our sensors and advise image processing strategies to further limit the impact of the effect on LSST WL science pathways.

47 OTHER INSTRUMENTATION↗

DOE OSTI2013

A Parallel, Fully Coupled, Fully Implicit Solution to Reactive Transport in Porous Media Using the Preconditioned Jacobian-Free Newton-Krylov Method

Modeling large multicomponent reactive transport systems in porous media is particularly challenging when the governing partial differential algebraic equations (PDAEs) are highly nonlinear and tightly coupled due to complex nonlinear reactions and strong solution-media interactions. Here we present a preconditioned Jacobian-Free Newton-Krylov (JFNK) solution approach to solve the governing PDAEs in a fully coupled and fully implicit manner. A well-known advantage of the JFNK method is that it does not require explicitly computing and storing the Jacobian matrix during Newton nonlinear iterations. Our approach further enhances the JFNK method by utilizing physics-based, block preconditioning and a multigrid algorithm for efficient inversion of the preconditioner. This preconditioning strategy accounts for self- and optionally, cross-coupling between primary variables using diagonal and off-diagonal blocks of an approximate Jacobian, respectively. Numerical results are presented demonstrating the efficiency and massive scalability of the solution strategy for reactive transport problems involving strong solution-mineral interactions and fast kinetics. We found that the physics-based, block preconditioner significantly decreases the number of linear iterations, directly reducing computational cost; and the strongly scalable algebraic multigrid algorithm for approximate inversion of the preconditioner leads to excellent parallel scaling performance.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

A stable partitioned FSI algorithm for incompressible flow and deforming beams

An added-mass partitioned (AMP) algorithm is described for solving fluid–structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier–Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler–Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. The key AMP interface condition is a generalized Robin (mixed) condition on the fluid pressure. This condition, which is derived at a continuous level, has no adjustable parameters and is applied at the discrete level to couple the partitioned domain solvers. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet–Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Quantified Gamow shell model interaction for p s d -shell nuclei

Background: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interaction approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow shell model (GSM) provides such a framework as it is capable of describing resonant and nonresonant many-body states on equal footing. Purpose: To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the psdf-shell-model space. The resulting interaction is expected to describe nuclei with 5 ≤ A ≲ 12 at the p-sd-shell interface. Method: The one-body potential of the 4 He core is modeled by a Woods-Saxon + spin-orbit + Coulomb potential, and the finite-range nucleon-nucleon interaction between the valence nucleons consists of central, spin-orbit, tensor, and Coulomb terms. The GSM is used to compute key fit observables. The χ 2 optimization is performed using the Gauss-Newton algorithm augmented by the singular value decomposition technique. The resulting covariance matrix enables quantification of statistical errors within the linear regression approach. Results: The optimized one-body potential reproduces nucleon- 4 He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the helium, lithium, and beryllium isotopes up to A = 9 . A very good agreement with experimental results was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. In conclusion: The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the psd region of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2016

Quantifying PM2.5-Meteorology Sensitivities in a Global Climate Model

Climate change can influence fine particulate matter concentrations (PM2.5) through changes in air pollution meteorology. Knowledge of the extent to which climate change can exacerbate or alleviate air pollution in the future is needed for robust climate and air pollution policy decision-making. To examine the influence of climate on PM2.5, we use the Geophysical Fluid Dynamics Laboratory Coupled Model version 3 (GFDL CM3), a fully-coupled chemistry-climate model, combined with future emissions and concentrations provided by the four Representative Concentration Pathways (RCPs). For each of the RCPs, we conduct future simulations in which emissions of aerosols and their precursors are held at 2005 levels while other climate forcing agents evolve in time, such that only climate (and thus meteorology) can influence PM2.5 surface concentrations. We find a small increase in global, annual mean PM2.5 of about 0.21 micro-g/cu m3 (5%) for RCP8.5, a scenario with maximum warming. Changes in global mean PM2.5 are at a maximum in the fall and are mainly controlled by sulfate followed by organic aerosol with minimal influence of black carbon. RCP2.6 is the only scenario that projects a decrease in global PM2.5 with future climate changes, albeit only by -0.06 micro-g/cu m (1.5%) by the end of the 21st century. Regional and local changes in PM2.5 are larger, reaching upwards of 2 micro-g/cu m for polluted (eastern China) and dusty (western Africa) locations on an annually averaged basis in RCP8.5. Using multiple linear regression, we find that future PM2.5 concentrations are most sensitive to local temperature, followed by surface wind and precipitation. PM2.5 concentrations are robustly positively associated with temperature, while negatively related with precipitation and wind speed. Present-day (2006-2015) modeled sensitivities of PM2.5 to meteorological variables are evaluated against observations and found to agree reasonably well with observed sensitivities (within 10e50% over the eastern United States for several variables), although the modeled PM2.5 is less sensitive to precipitation than in the observations due to weaker convective scavenging. We conclude that the hypothesized "climate penalty" of future increases in PM2.5 is relatively minor on a global scale compared to the influence of emissions on PM2.5 concentrations.

PM2.5↗

DOE OSTI2014

Effect of driver impedance on dense plasma focus Z-pinch neutron yield

The Z-pinch phase of a dense plasma focus (DPF) heats the plasma by rapid compression and accelerates ions across its intense electric fields, producing neutrons through both thermonuclear and beam-target fusion. Driver characteristics have empirically been shown to affect performance, as measured by neutron yield per unit of stored energy. We are exploring the effect of driver characteristics on DPF performance using particle-in-cell (PIC) simulations of a kJ scale DPF. In this work, our PIC simulations are fluid for the run-down phase and transition to fully kinetic for the pinch phase, capturing kinetic instabilities, anomalous resistivity, and beam formation during the pinch. The anode-cathode boundary is driven by a circuit model of the capacitive driver, including system inductance, the load of the railgap switches, the guard resistors, and the coaxial transmission line parameters. It is known that the driver impedance plays an important role in the neutron yield: first, it sets the peak current achieved at pinch time; and second, it affects how much current continues to flow through the pinch when the pinch inductance and resistance suddenly increase. Here we show from fully kinetic simulations how total neutron yield depends on the impedance of the driver and the distributed parameters of the transmission circuit. Direct comparisons between the experiment and simulations enhance our understanding of these plasmas and provide predictive design capability for neutron source applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1999

Electrorheological Fluid Based Force Feedback Device

Parallel to the efforts to develop fully autonomous robots, it is increasingly being realized that there are applications where it is essential to have a fully controlled robot and "feel" its operating conditions, i.e. telepresence. This trend is a result of the increasing efforts to address tasks where humans can perform significantly better but, due to associated hazards, distance, physical limitations and other causes, only robots can be employed to perform these tasks. Such robots need to be assisted by a human that remotely controls the operation. To address the goal of operating robots as human surrogates, the authors launched a study of mechanisms that provide mechanical feedback. For this purpose, electrorheological fluids (ERF) are being investigated for the potential application as miniature haptic devices. This family of electroactive fluids has the property of changing the viscosity during electrical stimulation. Consequently, ERF can be used to produce force feedback haptic devices for tele-operated control of medical and space robotic systems. Forces applied at the robot end-effector due to a compliant environment are reflected to the user using an ERF device where a change in the system viscosity will occur proportionally to the transmitted force. Analytical model and control algorithms are being developed taking into account the non-linearities of these type of devices. This paper will describe the concept and the developed mechanism of ERF based force feedback. The test process and the physical properties of this device will be described and the results of preliminary tests will be presented.

Pfeiffer, Charles↗

DOE OSTI2017

Fast and accurate implementation of Fourier spectral approximations of nonlocal diffusion operators and its applications

This work is concerned with the Fourier spectral approximation of various integral differential equations associated with some linear nonlocal diffusion and peridynamic operators under periodic boundary conditions. For radially symmetric kernels, the nonlocal operators under consideration are diagonalizable in the Fourier space so that the main computational challenge is on the accurate and fast evaluation of their eigenvalues or Fourier symbols consisting of possibly singular and highly oscillatory integrals. For a large class of fractional power-like kernels, we propose a new approach based on reformulating the Fourier symbols both as coefficients of a series expansion and solutions of some simple ODE models. We then propose a hybrid algorithm that utilizes both truncated series expansions and high order Runge–Kutta ODE solvers to provide fast evaluation of Fourier symbols in both one and higher dimensional spaces. It is shown that this hybrid algorithm is robust, efficient and accurate. As applications, we combine this hybrid spectral discretization in the spatial variables and the fourth-order exponential time differencing Runge–Kutta for temporal discretization to offer high order approximations of some nonlocal gradient dynamics including nonlocal Allen–Cahn equations, nonlocal Cahn–Hilliard equations, and nonlocal phase-field crystal models. Numerical results show the accuracy and effectiveness of the fully discrete scheme and illustrate some interesting phenomena associated with the nonlocal models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Confinement and Processing Can Alter the Morphology and Periodicity of Bottlebrush Block Copolymers in Thin Films

Bottlebrush block copolymers (BBCPs) are intriguing architectural variations on linear BCPs with highly tunable structure. Confinement can have a significant impact on polymer assembly, giving rise to changes in morphology, assembly kinetics, and properties like the glass transition. Given that confinement leads to significant changes in the persistence length of bottlebrush homopolymers, it is reasonable to expect that BBCPs will see significant changes in their structure and periodicity relative to the bulk morphology. Understanding how confinement influences assembly will be important for designing BBCPs for thin film applications including membranes, integrated photonic structures, and potentially BCP lithography. In this work, in order to study the effects of confinement on BBCP conformation and morphology, a blade coating was used to prepare films with continuous variation in film thickness. Unlike thin films of linear BCPs, islands/holes were not observed, and instead mixtures of parallel and perpendicular morphologies emerge after annealing. The lamellar periodicity ($L_0$) of the morphologies is found to be thickness dependent, increasing $L_0$ with decreasing film thickness for blade coated films. Films coated out of tetrahydrofuran (THF) resulted in a single well-defined lamellar periodicity, verified through atomic force microscopy (AFM) and grazing incidence small-angle X-ray scattering (GISAXS), which increases dramatically from the bulk value (30.6 nm) and continues to increase as the film thickness decreases. The largest observed $L_0$ was 65.5 nm, and this closely approaches the estimated upper limit of 67 nm corresponding to a fully extended backbone in a bilayer arrangement. Films coated out of propylene glycol methyl ether acetate (PGMEA) resulted in a mixture of perpendicular lamellae and a smaller, likely cylindrical morphology. The lamellar portion of the film shows the same thickness dependence as the lamellae observed in the THF coated films. The scaling of the lamellar $L_0$ with respect to film thickness follows predictions for confined semiflexible polymers with weak excluded volume interactions and can be related to models for confinement of DNA. Spin coated films shows the same reduction in periodicity, although at very different film thicknesses. This result suggests that the material has shallow free-energy barriers to transitioning between different $L_0$ and morphologies, a property that could be taken advantage of for patterning diverse structures with a single material.

36 MATERIALS SCIENCE↗

NASA NTRS2001

Integral Twist Actuation of Helicopter Rotor Blades for Vibration Reduction

Active integral twist control for vibration reduction of helicopter rotors during forward flight is investigated. The twist deformation is obtained using embedded anisotropic piezocomposite actuators. An analytical framework is developed to examine integrally-twisted blades and their aeroelastic response during different flight conditions: frequency domain analysis for hover, and time domain analysis for forward flight. Both stem from the same three-dimensional electroelastic beam formulation with geometrical-exactness, and axe coupled with a finite-state dynamic inflow aerodynamics model. A prototype Active Twist Rotor blade was designed with this framework using Active Fiber Composites as the actuator. The ATR prototype blade was successfully tested under non-rotating conditions. Hover testing was conducted to evaluate structural integrity and dynamic response. In both conditions, a very good correlation was obtained against the analysis. Finally, a four-bladed ATR system is built and tested to demonstrate its concept in forward flight. This experiment was conducted at NASA Langley T~ansonic Dynamics Tunnel and represents the first-of-a-kind Mach-scaled fully-active-twist rotor system to undergo forward flight test. In parallel, the impact upon the fixed- and rotating-system loads is estimated by the analysis. While discrepancies are found in the amplitude of the loads under actuation, the predicted trend of load variation with respect to its control phase correlates well. It was also shown, both experimentally and numerically, that the ATR blade design has the potential for hub vibratory load reduction of up to 90% using individual blade control actuation. Using the numerical framework, system identification is performed to estimate the harmonic transfer functions. The linear time-periodic system can be represented by a linear time-invariant system under the three modes of blade actuation: collective, longitudinal cyclic, and lateral cyclic. A vibration minimizing controller is designed based on this result, which implements classical disturbance rejection algorithm with some modifications. The controller is simulated numerically, and more than 90% of the 4P hub vibratory load is eliminated. By accomplishing the experimental and analytical steps described in this thesis, the present concept is found to be a viable candidate for future generation low-vibration helicopters. Also, the analytical framework is shown to be very appropriate for exploring active blade designs, aeroelastic behavior prediction, and as simulation tool for closed-loop controllers.

Shin, SangJoon↗

DOE OSTI2013

On the spatial stability of a liquid jet in the presence of vapor cavities

A dispersion equation describing the effect of temperature differences on the stability of three-dimensional cylindrical liquid jets in the presence of vapor cavities is presented by the use of linear stability analysis. The mathematical model and its solving method are verified by comparing them with the data in the literature, and then the effect of temperature differences between jet and surrounding gas on the spatial stability of liquid jet is investigated. Some conclusions can be drawn from the results of this investigation: (1) the temperature difference destabilizes the liquid jet when the jet liquid is cooler than the surrounding gas, (2) the smallest atomized droplet without taking into account the effect of temperature differences is significantly larger than that when the effect of temperature differences is taken into account, (3) the effect of temperature differences on the stability of liquid jet has little relationship with azimuthal wave modes, (4) cavitation destabilizes the liquid jet when the value of the bubble volume fraction is not greater than 0.1 (0 ≤ α ≤ 0.1), and the temperature difference can weaken this effect of cavitation on the stability of liquid jet, and (5) cavitation is responsible for generating smaller droplets, the effect of cavitation on the critical wave number with and without taking into account the effect of temperature differences is quite different, and temperature difference is likely to fully restrain the effect of cavitation on the critical wave number; however, cavitation is again responsible for generating smaller droplets despite the effect of temperature differences when the bubble volume fraction α = 0.1. These findings may explain some observations of practical atomizer performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Performance analysis of fully explicit and fully implicit solvers within a spectral element shallow-water atmosphere model

In this study, explicit Runge–Kutta methods and implicit multistep methods utilizing a Newton–Krylov nonlinear solver are evaluated for a range of configurations of the shallow-water dynamical core of the spectral element community atmosphere model to evaluate their computational performance. These configurations are designed to explore the attributes of each method under different but relevant model usage scenarios including varied spectral order within an element, static regional refinement, and scaling to the largest problem sizes. This analysis is performed within the shallow-water dynamical core option of a full climate model code base to enable a wealth of simulations for study, with the aim of informing solver development within the more complete hydrostatic dynamical core used for climate research. The limitations and benefits to using explicit versus implicit methods, with different parameters and settings, are discussed in light of the trade-offs with Message Passing Interface (MPI) communication and memory and their inherent efficiency bottlenecks. Given the performance behavior across the configurations analyzed here, the recommendation for future work using the implicit solvers is conditional based on scale separation and the stiffness of the problem. For the regionally refined configurations, the implicit method has about the same efficiency as the explicit method, without considering efficiency gains from a preconditioner. The potential for improvement using a preconditioner is greatest for higher spectral order configurations, where more work is shifted to the linear solver. Finally, initial simulations with OpenACC directives to utilize a Graphics Processing Unit (GPU) when performing function evaluations show improvements locally, and that overall gains are possible with adjustments to data exchanges.

58 GEOSCIENCES↗

DOE OSTI2018

Comparative analysis of magnetic resonance in the polaron pair recombination and the triplet exciton-polaron quenching models

We present a comparative theoretical study of magnetic resonance within the polaron pair recombination (PPR) and the triplet exciton-polaron quenching (TPQ) models. Both models have been invoked to interpret the photoluminescence detected magnetic resonance (PLDMR) results in π-conjugated materials and devices. We show that resonance line shapes calculated within the two models differ dramatically in several regards. First, in the PPR model, the line shape exhibits unusual behavior upon increasing the microwave power: it evolves from fully positive at weak power to fully negative at strong power. In contrast, in the TPQ model, the PLDMR is completely positive, showing a monotonic saturation. Second, the two models predict different dependencies of the resonance signal on the photoexcitation power, P L . At low P L , the resonance amplitude Δ I/I is ∝ P L within the PPR model, while it is ∝ P 2 L crossing over to P 3 L within the TPQ model. On the physical level, the differences stem from different underlying spin dynamics. Most prominently, a negative resonance within the PPR model has its origin in the microwave-induced spin-Dicke effect, leading to the resonant quenching of photoluminescence. The spin-Dicke effect results from the spin-selective recombination, leading to a highly correlated precession of the on-resonance pair partners under the strong microwave power. This effect is not relevant for TPQ mechanism, where the strong zero-field splitting renders the majority of triplets off resonance. On the technical level, the analytical evaluation of the line shapes for the two models is enabled by the fact that these shapes can be expressed via the eigenvalues of a complex Hamiltonian. This bypasses the necessity of solving the much larger complex linear system of the stochastic Liouville equations. Lastly, our findings pave the way towards a reliable discrimination between the two mechanisms via cw PLDMR.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

In many recent applications, particularly in the field of atom-centered descriptors for interatomic potentials, tensor products of spherical harmonics have been used to characterize complex atomic environments. When coupled with a radial basis, the atomic cluster expansion (ACE) basis is obtained. However, symmetrization with respect to both rotation and permutation results in an overcomplete set of ACE descriptors with linear dependencies occurring within blocks of functions corresponding to particular generalized Wigner symbols. All practical applications of ACE employ semi-numerical constructions to generate a complete, fully independent basis. While computationally tractable, the resultant basis cannot be expressed analytically, is susceptible to numerical instability, and thus has limited reproducibility. Here we present a procedure for generating explicit analytic expressions for a complete and independent set of ACE descriptors. The procedure uses a coupling scheme that is maximally symmetric w.r.t. permutation of the atoms, exposing the permutational symmetries of the generalized Wigner symbols, and yields a permutation-adapted rotationally and permutationally invariant basis (PA-RPI ACE). Theoretical support for the approach is presented, as well as numerical evidence of completeness and independence. A summary of explicit enumeration of PA-RPI functions up to rank 6 and polynomial degree 32 is provided. The PA-RPI blocks corresponding to particular generalized Wigner symbols may be either larger or smaller than the corresponding blocks in the simpler rotationally invariant basis. Finally, we demonstrate that basis functions of high polynomial degree persist under strong regularization, indicating the importance of not restricting the maximum degree of basis functions in ACE models a priori.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2017

Analysis of Alfven eigenmodes destabilization by energetic particles in TJ-II using a Landau-closure model

Alfvén Eigenmodes (AE) can be destabilized by energetic particles in neutral beam injection (NBI) heated plasmas through inverse Landau damping and couplings with gap modes in the shear Alfvén continua. In this paper, we describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system using the reduced MHD equations, density and parallel velocity moments for the energetic particles as well as the geodesic acoustic wave dynamics. A closure relation adds the Landau damping and resonant destabilization effects in the model. We apply the model to study the Alfvén modes stability in TJ-II, performing a parametric analysis in a range of realistic values of energetic particle β ($\beta_{f}$ ), ratios of thermal/Alfvén velocities ($V_{\rm th}/V_{\rm A0}$ ), energetic particle density profiles and toroidal modes (n) including toroidal and helical couplings. The study predicts a large helical coupling between different toroidal modes and the destabilization of helical Alfvén eigenmodes (HAE) with frequencies similar to the AE activity measured in TJ-II, between 50–400 kHz. The analysis has also revealed the destabilization of GAE (global Alfvén eigenmodes), TAE (toroidal Alfvén eigenmodes) and EPM (energetic particle modes). For the modes considered here, optimized TJ-II operations require a $\rlap{-} \iota$ profile in the range of $[0.845, 0.979]$ to stabilize AEs in the inner and middle plasma. AEs in the plasma periphery cannot be fully stabilized, although for a configuration with $\rlap{-} \iota = [0.945, 1.079]$ , only $n=7, 11, 15$ AE are unstable with a growth rate 4 times smaller compared to the standard $\rlap{-} \iota = [1.54, 1.68]$ case and a frequency of 100 kHz. Finally, we reproduce the frequency sweeping evolution of the AE frequency observed in TJ-II as the $\rlap{-} \iota$ profile is varied. The AE frequency sweeping is caused by consecutive changes of the instability dominant modes between different helical families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Habemus superstratum! A constructive proof of the existence of superstrata

Here, we construct the first example of a superstratum: a class of smooth horizonless supergravity solutions that are parameterized by arbitrary continuous functions of (at least) two variables and have the same charges as the supersymmetric D1-D5-P black hole. We work in Type IIB string theory on T 4 or K 3 and our solutions involve a subset of fields that can be described by a six-dimensional supergravity with two tensor multiplets. The solutions can thus be constructed using a linear structure, and we give an explicit recipe to start from a superposition of modes specified by an arbitrary function of two variables and impose regularity to obtain the full horizonless solutions in closed form. We also give the precise CFT description of these solutions and show that they are not dual to descendants of chiral primaries. They are thus much more general than all the known solutions whose CFT dual is precisely understood. Hence our construction represents a substantial step toward the ultimate goal of constructing the fully generic superstratum that can account for a finite fraction of the entropy of the three-charge black hole in the regime of parameters where the classical black hole solution exists.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH 4 , C 2 H 4 , C 2 H 6 , C 3 H 6 , C 3 H 8 , and CO 2 across a temperature range of 250–400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure–temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption–desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300–550 K. Here, a temperature threshold of ~400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

The time evolution of gaps in tidal streams

We model the time evolution of gaps in tidal streams that are caused by the impact of a dark matter subhalo, while these orbit a spherical gravitational potential. To this end, we make use of the simple behavior of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential, the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale, and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10{sup 8} M{sub ⊙} directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kiloparsecs after a few gigayears. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericenter. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.

79 ASTRONOMY AND ASTROPHYSICS↗