Engineering Papers⌕ Search

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.

At least 847 records · Page 47

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extended nonergodic regime and spin subdiffusion in disordered SU(2)-symmetric Floquet systems

In this work, we explore thermalization and quantum dynamics in a one-dimensional disordered SU(2)-symmetric Floquet model, where a many-body localized phase is prohibited by the non-Abelian symmetry. Despite the absence of localization, we find an extended nonergodic regime at strong disorder where the system exhibits nonthermal behaviors. In the strong disorder regime, the level spacing statistics exhibit neither a Wigner-Dyson nor a Poisson distribution, and the spectral form factor does not show a linear-in-time growth at early times characteristic of random matrix theory. The average entanglement entropy of the Floquet eigenstates is subthermal, although violating an area-law scaling with system sizes. We further compute the expectation value of local observables and find strong deviations from the eigenstate thermalization hypothesis. The infinite-temperature spin autocorrelation function decays at long times as $t^{–β}$ with $\textit{β}$ < 0.5, indicating subdiffusive transport at strong disorders.

36 MATERIALS SCIENCE↗

Analytic solutions of the DGLAP evolution and theoretical uncertainties

The energy dependence for the singlet sector of Parton Distributions Functions (PDFs) is described by an entangled pair of ordinary linear differential equations. Although there are no exact analytic solutions, it is possible to provide approximated results depending on the assumptions and the methodology adopted. These results differ in their sub-leading, neglected terms and ultimately they are associated with different treatments of the theoretical uncertainties. In this work, a novel analytic approach in Mellin space is presented and a new methodology for obtaining closed and exponentiated analytic solutions is devised. Different results for the DGLAP evolution at Next-Leading-Order are compared, discussing advantages and disadvantages for each solution. The generalizations to higher orders are addressed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Performance Analysis and Optimal Node-aware Communication for Enlarged Conjugate Gradient Methods

Krylov methods are a key way of solving large sparse linear systems of equations but suffer from poor strong scalability on distributed memory machines. Furthermore, this is due to high synchronization costs from large numbers of collective communication calls alongside a low computational workload. Enlarged Krylov methods address this issue by decreasing the total iterations to convergence, an artifact of splitting the initial residual and resulting in operations on block vectors. In this article, we present a performance study of an enlarged Krylov method, Enlarged Conjugate Gradients (ECG), noting the impact of block vectors on parallel performance at scale. Most notably, we observe the increased overhead of point-to-point communication as a result of denser messages in the sparse matrix-block vector multiplication kernel. Additionally, we present models to analyze expected performance of ECG, as well as motivate design decisions. Most importantly, we introduce a new point-to-point communication approach based on node-aware communication techniques that increases efficiency of the method at scale.

97 MATHEMATICS AND COMPUTING↗

Analysis of the doubly radiative decays $η^{(')}$ → $π^0γγ$ and $η^{(')}$ → $ηγγ$

We study the doubly radiative decays $η^{(')}$ → $π^0γγ$ and $η^{(')}$ → $ηγγ$ in terms of vector and scalar meson exchange contributions using, respectively, the frameworks of Vector Meson Dominance and the Linear Sigma model. Predictions for the ${m}_{γγ}^{2}$ invariant mass distribution and branching ratios are given, and a discussion of how they compare to experimental data is carried out.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron pitch-angle diffusion driven by oblique whistler-mode turbulence.

A general description of cyclotron harmonic resonant pitch-angle scattering is presented. Quasi-linear diffusion coefficients are prescribed in terms of the wave normal distribution of plasma wave energy. Numerical computations are performed for the specific case of relativistic electrons interacting with a band of low frequency whistler-mode turbulence. A parametric treatment of the wave energy distribution permits normalized diffusion coefficients to be presented graphically solely as a function of the electron pitch-angle. The diffusion coefficients generally decrease with increasing cyclotron harmonic number. Higher harmonic diffusion is insignificant at very small electron pitch-angles, but becomes increasingly important as the pitch-angle increases. One thus expected the rate of pitch-angle scattering to decrease with increasing electron energy, since the resonant value of the latter varies proportionately with harmonic number. This indicates that, in mirror-type magnet field geometrics, such as the earth's radiation belts, the diffusion losses of high energy electrons are likely to be appreciably slower than those at low energy.

Lyons, L. R.↗

Numerical solution of laminar jet mixing with and without free stream.

Systematic numerical solutions of two-dimensional and axisymmetrical laminar jet of an incompressible fluid with and without free stream have been obtained. For the case without free stream, the exact numerical solutions have been obtained for different initial velocity profiles and compared with experimental results and similarity solutions. For the case with free stream, the numerical solutions have been compared with the linearized analytical solutions. An approximate numerical solution for the axial velocity distribution of the three-dimensional laminar jet is proposed. The accuracy of the method has been determined from the exact solutions of the two limiting cases of three-dimensional jets, i.e., the two-dimensional and axisymmetric cases.

Pai, S. I.↗

Combustion instability with partial length acoustic liners

The development of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors is discussed. The studies involved nonlinear wave propagation effects and nonlinear combustion zone, liner, and nozzle responses. A second primary extension of the basic theory was directed at the inclusion of the effects of injector face baffles. Attention was restricted to the linear (small amplitude) problem, finite Mach numbers, three dimensional oscillations, and distributed combustion effects.

Mitchell, C. E.↗

Estimation of leading-edge thrust for supersonic wings of arbitrary planform

A numerical method for the estimation of leading edge thrust for supersonic wings of arbitrary planform was developed and was programmed as an extension to an existing high speed digital computer method for prediction of wing pressure distributions. The accuracy of the method was assessed by comparison with linearized theory results for a series of flat delta wings. Application of the method to wings of arbitrary planform, both flat and cambered, is illustrated in several examples.

Carlson, H. W.↗

VLA observations of stellar planetary nebulae

Coordinates, dimensions, 4885-MHz flux densities, and brightness temperatures of K3-2, NGC 6833, Ps 1, II 5117, Me 2-2, Hb 12, Vy 1-1, and M1-5 are reported. In two other cases, H3-29 and H3-75, confused extended structure was detected in which the nebula could not be identified with certainty. He 2-467, M1-2, and Peterson's H-alpha object in M15 were also included in the observations but not detected with an upper limit of less than 10 mJy. The observations are compared with some of the previous optical and radio data, such as log S(H-beta). Distances are computed from the present data with standard assumptions. Corresponding linear radii range below 0.1 pc, among the smallest in previous distributions of radius.

Johnson, H. M.↗

The statistical character of Pc 4 magnetic pulsations at synchronous orbit

A statistical study of 215 magnetic pulsation events identified on the basis of waveform and period as Pc 4 (45- to 150-s period) observed at synchronous orbit by ATS 6 is reported. The study was designed to be similar to that done with Pc 3 (10- to 45-s period) by Arthur et al. (1977). These Pc 4 are found to occur most often near 1800 local time but are observed at all local times. The dominant period is approximately 100 s, but a secondary peak in the distribution occurred at approximately 53 s. They tend to be nearly linearly polarized, principally transverse to the ambient magnetic field. The Pc 4 are found to divide naturally into two classes (radial and azimuthal) on the basis of the azimuth of the major axis of polarization. These classes divide further into low-frequency (less than or equal to 0.015 Hz) and high-frequency (greater than or equal to 0.015 Hz) classes. These four classes organize the remaining wave characteristics quite well. Division of the Pc 4 events into these classes does not lead to identification of a single source mechanism as it did with Pc 3. However, some correspondence between the classes and various proposed mechanisms does exist.

Arthur, C. W.↗

The dynamics and control of large flexible space structures

The dynamics and attitude and shape control of very large, inherently flexible spacecraft systems were investigated. Increasingly more complex examples were examined, beginning with a uniform free-free beam, next a free-free uniform plate/platform and finally by considering a thin shallow spherical shell structure in orbit. The effects devices were modeled. For given sets of assumed actuator locations, the controllability of these systems was first established. Control laws for each of the actuators were developed based on decoupling techniques (including distributed modal control) pole placement algorithms and a application of the linear regulator problem for optical control theory.

Bainum, P. M.↗

Subsonic panel method for designing wing surfaces from pressure distribution

An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.

Bristow, D. R.↗

The cosmic-ray source composition

The cosmic-ray source composition is shown to be related to the observed cosmic-ray abundances by a set of linear equations which are valid for either ad hoc path length distributions or quantitative solutions of galactic propagation models. In the limit of no ionization energy loss, this method is exact for all path length distributions and can replace more detailed numerical integrations; it provides a simple way to calculate directly the source composition implied by a particular set of cross sections, observed abundances, and propagation model. This method is used to determine the source abundances of the galactic cosmic rays in a model-independent manner. Compared to a simple exponential, the path length distribution required by observations is deficient in the contributions from short (not greater than 0.5 g/sq cm) paths at energies of a few GeV per nucleon.

Margolis, S. H.↗

A probabilistic Hu-Washizu variational principle

A Probabilistic Hu-Washizu Variational Principle (PHWVP) for the Probabilistic Finite Element Method (PFEM) is presented. This formulation is developed for both linear and nonlinear elasticity. The PHWVP allows incorporation of the probabilistic distributions for the constitutive law, compatibility condition, equilibrium, domain and boundary conditions into the PFEM. Thus, a complete probabilistic analysis can be performed where all aspects of the problem are treated as random variables and/or fields. The Hu-Washizu variational formulation is available in many conventional finite element codes thereby enabling the straightforward inclusion of the probabilistic features into present codes.

Liu, W. K.↗

The behavior of planetary wave 2 in preconditioned zonal flows

The response of linear planetary wave 2 to changes in the isentropic zonally symmetric distribution of potential vorticity (PV) is investigated numerically. Wave 2 is sensitive to the width and position of a region where the meridional derivative of PV is weak, denoted the 'surf zone', in the middle stratosphere. A narrow surf zone leads to an amplification of wave 2, and confines some of its Eliassen-Palm wave activity in high latitudes. For wider surf zones the wave activity is concentrated near the associated critical surface, and the amplitude of the wave decreases. Large changes in the wave amplitude and in the distribution of its activity are associated with the subtle changes in the zonal winds produced by modest modifications in the distribution of PV. Basic states that include regions of reversed meridional gradients of PV lead to wave overreflection and strong poleward focusing of wave activity. The amplitude of wave 2 is enhanced in the presence of negative gradients, with large responses occurring for eastward traveling waves.

Robinson, Walter A.↗