SEARCH · Engineering Papers
Results for “linearized distribution”
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.
Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma
Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.
Electric Field Analysis in a Cold-Ion Source Using Stark Spectroscopy of Rydberg Atoms
Here, we analyze electric fields in ion sources generated by quasicontinuous photoionization of cold Rb atoms trapped in the focal spot of a near-concentric, in-vacuum cavity for 1064-nm laser light. Ion streams are extracted with an external electric field, F. Stark effects of Rb 57F and of nearby high-angular-momentum Rydberg levels, which exhibit large, linear Stark shifts, are employed to study the net electric-field probability distribution within the ion-source region over an extraction-field range of 0 < F < 0.35 V/cm. For F = 0, we also investigate ion-field-induced Stark spectra of the 60P 1/2 state, which exhibits a (lesser) quadratic electric-field response that affords a simplified electric-field analysis. Experimental Rydberg spectra are compared with theoretical Stark spectra, which are weighted with net electric-field distributions obtained from classical ion-trajectory simulations that include Coulomb interactions. Experiments and models are in good agreement. At small F and high ion-source rates, the field approximately follows a Holtsmark distribution, and the ion streams are degraded by the Coulomb microfields. With increasing F and at lower ion-source rates, the fields become narrowly distributed around F, resulting in directional ion streams that are less degraded by microfields. Our results are of interest for monitoring cold-ion sources for focused-ion-beam applications, where Coulomb interactions are of concern, and for studies of electric fields in cold plasmas.
Average nonlinear dynamics of particles in gravitational pulses: Effective Hamiltonian, secular acceleration, and gravitational susceptibility
Particles interacting with a prescribed quasimonochromatic gravitational wave (GW) exhibit secular (average) nonlinear dynamics that can be described by Hamilton’s equations. We derive the Hamiltonian of this “ponderomotive” dynamics to the second order in the GW amplitude for a general background metric. For the special case of vacuum GWs, we show that our Hamiltonian is equivalent to that of a free particle in an effective metric, which we calculate explicitly. We also show that already a linear GW pulse displaces a particle from its unperturbed trajectory by a finite distance that is independent of the GW phase and proportional to the integral of the pulse intensity.This effect is independent from the nonlinear memory effects that has been known. We calculate the particle displacement analytically and show that our result is in agreement with numerical simulations. We also show how the Hamiltonian of the nonlinear averaged dynamics naturally leads to the concept of the linear gravitational susceptibility of a particle gas with an arbitrary phase-space distribution. We calculate this susceptibility explicitly to apply it, in a follow-up paper, toward studying self-consistent GWs in inhomogeneous media within the geometrical-optics approximation
Tracing cosmic gas in filaments and halos: Low-redshift insights from the kinematic Sunyaev-Zel’dovich effect
In this work, we leverage cosmic microwave background (CMB) data from the Atacama Cosmology Telescope (ACT) and LSS data from the imaging survey conducted by the Dark Energy Spectroscopic Instrument (DESI) to study the distribution of gas around galaxy groups at low redshift, z ≈ 0.3, via the kinematic Sunyaev-Zel'dovich (kSZ) effect. In particular, we perform velocity-weighted stacking on the photometric galaxies from the Bright Galaxy Survey (BGS) to isolate the monopole and quadrupole of the kSZ signal, orienting the stacked images along 2D filaments identified using the Hessian of the projected gravitational potential. We find a 7.2σ detection in the monopole of the signal (i.e., the gas density profile) and a 4σ detection in the quadrupole (m = 2), constituting the first measurement of the alignment between gas distribution and the cosmic web through the kSZ effect. As it is a linear probe of the local gas density, the kSZ has heightened sensitivity to the warm-hot intergalactic medium (WHIM), which is believed to house the majority of the "missing baryons."Mapping out the gas density at low redshifts, as enabled by our measurements, is crucial for weak lensing surveys, for which the impact of baryons on small scales is a major impediment. We compare the anisotropic signal against two hydrodynamical simulations, TNG300-1 and Illustris, which have very different baryonic feedback prescriptions. We find that the anisotropic signal measured in the data is comparable but slightly larger and more extended compared with the simulations. Further, this suggests that there is excess accretion and feedback taking place through the filaments, hinting at the possible presence of spin-filament alignment of the BGS objects.
Revisiting the gluon density from the Balitsky-Kovchegov equation with kinematical constraints and large 𝑥 terms
We perform analysis of the small 𝑥 nonlinear evolution equation formulated in momentum space supplemented by higher order terms. The equation is defined in the wide range of transverse momentum and longitudinal momentum fractions extending previous studies performed in Kutak and Kwiecinski [Screening effects in the ultrahigh-energy neutrino interactions, Eur. Phys. J. C 29, 521 (2003).] and in Kutak and Stasto [Unintegrated gluon distribution from modified BK equation, Eur. Phys. J. C 41, 343 (2005).]. The linear part of the equation is motivated by the renormalization group improved small 𝑥 approach, which accounts for resummation of higher orders, and includes collinear splitting functions and kinematical constraints. The solution to the equation is then used to perform the fit to deep inelastic scattering reduced cross section data.
Scrambling Dynamics with Imperfections in a Solvable Model
We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.
Results of the second Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+)
Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf–ocean modelling, the second Ice Shelf–Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0–2 (Asay-Davis et al., 2016), which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice–ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf–ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice–ocean boundary layer and ice shelf cavity circulation.
Stability of a steady, large amplitude whistler wave.
Study of the behavior of weak electrostatic waves in a collisionless magnetoplasma supporting a steady large amplitude whistler wave. All waves are assumed to propagate parallel to a uniform background magnetic field B sub zero. In the presence of the whistler wave fields each particle executes an oscillatory motion parallel to B sub zero, in addition to a translation along B sub zero and transverse motions. This oscillation causes the Landau resonance to be replaced by a series of new resonances between particles and the electrostatic modes. A distribution function for the perturbed plasma is constructed by solving the Vlasov equation, linearized in the electrostatic wave amplitudes. A dispersion relation is obtained and solved approximately for the growth/damping rate of the perturbations. Growing electrostatic modes are found to be approximately uncoupled. Trapped particles have a strong influence on the stability of the system.
Infrared sky noise study
The hardware and techniques to measure and compare sky noise at several sites were studied, and a device was developed that would maximize its output and minimize its output for modulation. The instrument and its functions are described. The nature of sky emissions and the fluctuation, gaseous sources of sky noise, and aerosol sources are discussed. It is concluded that sky noise really exists, and the spatial distribution of the sky noise sources are such that observed noise values are linear functions of chopping stroke.
Apollo 16 far ultraviolet imagery of the polar auroras, tropical airglow belts, and general airglow
Far-ultraviolet imagery of the earth in the wavelength ranges from 1050 to 1600 A and from 1250 to 1600 A was obtained from the lunar surface during the Apollo 16 mission on Apr. 21, 1972. The images have an angular resolution of about 2 arcmin (230-km linear resolution) and have been quantitatively analyzed to obtain absolute intensities and spatial distributions of the polar auroras (both wavelength ranges), the day and night airglow, and tropical airglow belts (1250-A to 1600-A wavelength range). The observations are consistent with previous results obtained from the OGO-4 spacecraft, but they have also provided details on the spatial distributions of the various emissions over an entire hemisphere at a single time. A general night airglow, at least in the Northern Hemisphere, is indicated.
The kinetic effects of Alfven wave pressure in the solar wind
The results of a kinetic model for the radial evolution of the proton component of the solar wind in the presence of Alfven waves are presented. The calculation is based on general quasi-linear equations developed to describe the temporal and spatial evolution of the ion distribution functions of a multispecies plasma in presence of waves, using a short wavelength expansion. These equations include new wave-particle interaction terms arising from temporal and spatial inhomogeneities in the plasma. Numerical solutions are obtained of these equations specialized to the case of Alfven waves in a spherically symmetric solar wind. The Alfven wave effects on the proton distribution function vary strongly in velocity space. Protons with small transverse velocities are primarily decelerated with respect to the wave rest frame. This deceleration becomes less important with increasing transverse velocity, as wave induced diffusion to larger transverse velocity becomes the dominant effect. The competition of these effects results in interesting distortions of evolving proton distribution functions which give rise the wave acceleration well known from fluid theory.
A 1.3 centimeter wavelength interferometric study of Jupiter
Maps of the brightness distribution of Jupiter at 1.3 cm wavelength are reported for the two orthogonal linear polarizations; one approximately parallel to the equator and one approximately parallel to the pole. The synthesized beam is 6.0 arc sec by 12.0 arc sec, which provides about 55 picture elements over the visible disk. The two maps do not differ significantly. There are no small real features of the size of one resolution element with brightness greater than 6 K. The disk brightness is approximately uniform, which is consistent with saturated ammonia vapor as the source of opacity in the atmosphere. Uniformly mixed ammonia, molecular hydrogen, and methane are strongly excluded as major sources of the opacity. Forcing the maps to have symmetry about the pole, as expected for the rapidly rotating planet, reduces the residual noise level in the maps so that there appears to be a small amount of limb darkening along the equator and a small amount of brightening at the poles.
The application of the least squares finite element method to Abel's integral equation
Abel's integral equation is the governing equation for certain problems in physics and engineering, such as radiation from distributed sources. The finite element method for the solution of this non-linear equation is presented for problems with cylindrical symmetry and the extension to more general integral equations is indicated. The technique was applied to an axisymmetric glow discharge problem and the results show excellent agreement with previously obtained solutions
Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques
Refinements to the modeling procedure developed to compute the magnetic fields at satellite orbit due to current distributions in the ionosphere and magnetosphere are described. The modeling technique utilizes a linear current element representation of the large scale space current system. A model polar current system is presented and magnetic field perturbations resulting from this system are computed along two hypothetical satellite orbits.
Investigation of the effects of external current systems on the MAGSAT data utilizing grid cell modeling techniques
The status of the initial testing of the modeling procedure developed to compute the magnetic fields at satellite orbit due to current distributions in the ionosphere and magnetosphere is reported. The modeling technique utilizes a linear current element representation of the large scale space-current system.
Bifurcation and stability of low-order steady flows in horizontally and vertically forced convection
A nonlinear spectral model of two-dimensional, shallow Boussinesq convection which responds to heating in both the horizontal and vertical directions is examined. The governing partial differential system is converted to an infinite set of ordinary differential equations and truncated to a small set to permit detailed study of the number and types of transitions from one flow configuration to another. The Hadley number and the Rayleigh number are defined as the horizontal and vertical thermal forcing mechanisms, respectively, for inclusion in the nonlinear spectral model, which is composed of three equations. The model is then used to describe steady states, linearly stable solutions, and balancing factors in unstable stratification. The number and the distribution of the steady states are found to be qualitatively independent of the aspect ratio and the Prandtl number.
Dynamic response of shock waves in transonic diffuser and supersonic inlet - An analysis with the Navier-Stokes equations and adaptive grid
An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.