Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 1,369 records · Page 76

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗

Single classical field description of interacting scalar fields

In this report we test the degree to which interacting Bosonic systems can be approximated by a classical field as total occupation number is increased. This is done with our publicly available code repository, QIBS, a new massively parallel solver for these systems. We use a number of toy models well studied in the literature and track when the classical field description admits quantum corrections, called the quantum breaktime. This allows us to test claims in the literature regarding the rate of convergence of these systems to the classical evolution. We test a number of initial conditions, including coherent states, number eigenstates, and field number states. We find that of these initial conditions, only number eigenstates do not converge to the classical evolution as occupation number is increased. We find that systems most similar to scalar field dark matter exhibit a logarithmic enhancement in the quantum breaktime with total occupation number. Systems with contact interactions or with field number state initial conditions, and linear dispersions, exhibit a power law enhancement. Finally, we find that the breaktime scaling depends on both model interactions and initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

A Bayesian approach to evaluation of soil biogeochemical models

Abstract. To make predictions about the carbon cycling consequences of rising global surface temperatures, Earth system scientists rely on mathematical soil biogeochemical models (SBMs). However, it is not clear which models have better predictive accuracy, and a rigorous quantitative approach for comparing and validating the predictions has yet to be established. In this study, we present a Bayesian approach to SBM comparison that can be incorporated into a statistical model selection framework. We compared the fits of linear and nonlinear SBMs to soil respiration data compiled in a recent meta-analysis of soil warming field experiments. Fit quality was quantified using Bayesian goodness-of-fit metrics, including the widely applicable information criterion (WAIC) and leave-one-out cross validation (LOO). We found that the linear model generally outperformed the nonlinear model at fitting the meta-analysis data set. Both WAIC and LOO computed higher overfitting risk and effective numbers of parameters for the nonlinear model compared to the linear model, conditional on the data set. Goodness of fit for both models generally improved when they were initialized with lower and more realistic steady-state soil organic carbon densities. Still, testing whether linear models offer definitively superior predictive performance over nonlinear models on a global scale will require comparisons with additional site-specific data sets of suitable size and dimensionality. Such comparisons can build upon the approach defined in this study to make more rigorous statistical determinations about model accuracy while leveraging emerging data sets, such as those from long-term ecological research experiments.

54 ENVIRONMENTAL SCIENCES↗

Inventory of forest and rangeland and detection of forest stress

The author has identified the following significant results. At the Atlanta site (226B) it was found that bulk color composites for October 15, 1972, and April 13, 1973, can be interpreted together to disclose the location of the perennial Kudzu vine (Pyeraria lobata). Land managers concerned with Kudzu eradication could use ERTS-1 to inventory locations over 200 meters (660 feet) square. Microdensitometer data collected on ERTS-1 Bulk photographic products for the Manitou test site (226C) have shown that the 15-step gray-scale tablets are not of systematic equal values corresponding to 1/14 the maximum radiant energy incident on the MSS sensor. The gray-scale values present a third-order polynomial function rather than a direct linear relationship. Although data collected on step tablets for precision photographic products appear more discrete, the density variation within blocks in almost as great as variations between blocks. These system errors will cause problems when attempting to analyze radiometric variances among vegetation and land use classes.

Heller, R. C.↗

Forestry, geology and hydrological investigations from ERTS-1 imagery in two areas of Ecuador, South America

The author has identified the following significant results. In the Oriente area, well-drained forests containing commercially valuable hardwoods can be recognized confidently and delineated quickly on the ERTS imagery. In the tropical rainforest, ERTS can provide an abundance of inferential information about large scale geologic structures. ERTS imagery is better than normal aerial photography for recognizing linears. The imagery is particularly useful for updating maps of the distributary system of the Guagas River Basin and of any other river with a similarly rapid changing channel pattern.

Moreno, N. V.↗

Evidence for a non-random magnetization of the moon

An approximate inversion, using a specific model, of Apollo subsatellite vector magnetometer data has provided estimates for the locations, scale sizes, and magnetization characteristics of strongly magnetized lunar crustal materials across a portion of the lunar far side. The nearly model-independent local bulk directions of magnetization, while otherwise random, exhibit a tendency to lie in or near planes perpendicular to the present lunar spin axis. The apparent existence of such a preferred plane of magnetization must be regarded as unfavorable for proposed local, presumably random sources of large-scale, coherent magnetization and implies that any assumed internal magnetizing dipole or external linear magnetizing field must also have been oriented mainly perpendicular to the present spin axis. The otherwise random orientations of these vectors when projected onto the equatorial plane suggests a simple explanation for the negligible lunar global dipole moment.

Hood, L.↗

An evaluation of linear acoustic theory for a hovering rotor

Linear acoustic calculations are compared with previously reported data for a small-scale hovering rotor operated at high tip Mach numbers. A detailed calculated description of the distributions of blade surface pressure and shear stress due to skin friction is presented. The noise due to skin friction and loading, in the rotor disk plane, is small compared to thickness noise. The basic conclusions of Boxwell et al about the importance of nonlinear effects are upheld. Some approximations involved in the current theories for the inclusion of nonlinear effects are discussed. Using a model nonlinear problem, it is shown that to use the acoustic analogy, good knowledge of the flowfield is required.

Morris, C. E. K., Jr.↗

Linear modes of convection in the solar envelope

The structure of the solar convective envelope calculated according to mixing-length theory is not consistent with the transport of the convective heat flux by any superposition of the fundamental linear convective modes. If the amplitudes of these modes can be estimated by consideration of the most important nonlinear terms in the equations of motion, it becomes possible to find solar models which are consistent with mode transport of the heat through the convective zone. These models may be characterized by a mixing length which varies from point to point through the envelope. The linear modes in such models show peaks in velocity at certain preferred length scales, but these scales do not necessarily correspond to those of the solar motions, such as the supergranulation.

Bogart, R. S.↗

On the effect of boundary layer growth on the stability of compressible flows

The method of multiple scales is used to describe a formally correct method based on the nonparallel linear stability theory, that examines the two and three dimensional stability of compressible boundary layer flows. The method is applied to the supersonic flat plate layer at Mach number 4.5. The theoretical growth rates are in good agreement with experimental results. The method is also applied to the infinite-span swept wing transonic boundary layer with suction to evaluate the effect of the nonparallel flow on the development of crossflow disturbances.

El-Hady, N. M.↗

Three-dimensional stability of growing boundary layers

A theory is developed for the linear stability of three-dimensional growing boundary layers. The method of multiple scales is used to derive partial-differential equations describing the temporal and spatial evolution of the complex amplitudes and wavenumbers of the disturbances. In general, these equations are elliptic unless certain conditions are satisfied. For a monochromatic disturbance, these conditions demand that the ratio of the components of the complex group velocity be real and thereby relate the direction of growth of the disturbance to the disturbance wave angle. For a nongrowing boundary layer, this condition reduces to d-alpha/d-beta being real, in agreement with the result obtained by using the saddle-point method. For a wavepacket, these conditions demand that the components of the group velocity be real.

Nayfeh, A. H.↗

Systematic variation of magnetic-field penetration depth in high-T(c) superconductors studied by muon spin relaxation

The muon relaxation rate (sigma) was measured in the high critical temperature superconductors YBa2Cu3O(x) for x = 6.66, 6.95, 7.0, and La1.85 Sr0.15 CuO4 in transverse external magnetic fields 1 is approximately 4kG. A simple relation is found which connects the transition temperature T(c), the magnetic field penetration depth lambda(L), the carrier concentration n(s) and the effective mass m* as T(c) varies as sigma which varies as 1/lambda(L) squared which varies as n(s)/m*. The linear dependence T(c) varies as n(s)/m* suggests a high energy scale for the coupling between superconducting carriers.

Uemura, Y. J.↗

Numerical simulation of homogeneous, isotropic, compressible turbulence

A new numerical algorithm is developed and applied to the direct simulation of compressible, homogeneous turbulent flows at low Mach numbers. A split method in time first solves a subset of the equations explicitly, followed by an implicit treatment of the pressure terms. The definition of an average sound speed over the entire field allows the implicit equations to be solved analytically, while relaxing the severe time limit imposed by the large sound speeds. Results from direct simulations on 96(3) grids provide a data base against which a new subgrid-scale model for compressible homogeneous turbulence is tested. This model reduces to the linear combination model in the absence of compressibility.

Erlebacher, G.↗

Systematic variation of magnetic-field penetration depth in high-Tc superconductors studied by muon-spin relaxation

The muon relaxation rate (sigma) was measured in the high critical temperature superconductors YBa2Cu3O(x) for x = 6.66, 6.95, 7.0, and La1.85 SrO.15 CuO4 in transverse external magnetic fields 1 is approximately 4 kG. A simple relation is found which connects the transition temperature T(c), the magnetic field penetration depth lambda(L), the carrier concentration n(s) and the effective mass m* as T(c) varies as sigma which varies as 1/lambda(L) squared which varies as n(s)/m*. The linear dependence T(c) varies as n(s)/m* suggests a high energy scale for the coupling between superconducting carriers.

Uemura, Y. J.↗

Instability of a supersonic shock free elliptic jet

This paper presents a comparison of the measured and the computed spatial stability properties of an aspect ratio 2 supersonic shock free elliptic jet. The shock free nature of the elliptic jet provides an ideal test of validity of modeling the large scale coherent structures in the initial mixing region of noncircular supersonic jets with linear hydrodynamic stability theory. Both aerodynamic and acoustic data were measured. The data are used to compute the mean velocity profiles and to provide a description of the spatial composition of pressure waves in the elliptic jet. A hybrid numerical scheme is applied to solve the Rayleigh problem governing the inviscid linear spatial stability of the jet. The measured mean velocity profiles are used to provide a qualitative model for the cross sectional geometry and the smooth velocity profiles used in the stability analysis. Computational results are presented for several modes of instability at two jet cross sections. The acoustic measurements show that a varicose instability is the jet's perferred mode of motion. The stability analysis predicts that the Strouhal number varies linearly as a function of axial distance in the jet's initial mixing region, which is in good qualitative agreement with previous measurements.

Baty, Roy S.↗

An improved reconstruction method for cosmological density fields

This paper proposes some improvements to existing reconstruction methods for recovering the initial linear density and velocity fields of the universe from the present large-scale density distribution. We derive the Eulerian continuity equation in the Zel'dovich approximation and show that, by applying this equation, we can trace the evolution of the gravitational potential of the universe more exactly than is possible with previous approaches based on the Zel'dovich-Bernoulli equation. The improved reconstruction method is tested using N-body simulations. When the Zel'dovich-Bernoulli equation describes the formation of filaments, then the Zel'dovich continuity equation also follows the clustering of clumps inside the filaments. Our reconstruction method recovers the true initial gravitational potential with an rms error about 3 times smaller than previous methods. We examine the recovery of the initial distribution of Fourier components and find the scale at which the recovered phases are scrambled with respect their true initial values. Integrating the Zel'dovich continuity equation back in time, we can improve the spatial resolution of the reconstruction by a factor of about 2.

Gramann, Mirt↗

Risk Reduction on X-33/RLV Engines

Risk management has received considerable attention in the X-33 and Reusable Launch Vehicle (RLV) program due to aggressive schedules, limited funding. and planned private investment to develop the commercial VentureStar vehicle. As an X-33 and RLV team member and main propulsion supplier, Boeing Rocketdyn Propulsion and Power has addressed risk through a methodical application of systems engineering in identifying, assessing, and mitigating risks. The methods employed involve rigorous risk mitigation planning early in development, continuous risk monitoring and assessment during the course of development, and the systematic verification of compliance with technical requirements prior to delivery. In addition, an engine system reliability analysis was conducted to reduce risk. In July 1996, NASA selected Lockheed Martin's "Skunk Works" (LMSW) as the lead contractor for the X-33 and RLV program. The X-33 vehicle is a half-scale pathfinder for the full-scale RLV. The LMSW RLV design is a lifting body shaped vehicle employing linear aerospike engine provided propulsion. The initial X-33 flight is planned for the summer of 2000, and the initial VentureStar flight is planned for between 2005 and 2007.

Crowley, Tim↗

Extending unbiased stereology of brain ultrastructure to three-dimensional volumes

OBJECTIVE: Analysis of brain ultrastructure is needed to reveal how neurons communicate with one another via synapses and how disease processes alter this communication. In the past, such analyses have usually been based on single or paired sections obtained by electron microscopy. Reconstruction from multiple serial sections provides a much needed, richer representation of the three-dimensional organization of the brain. This paper introduces a new reconstruction system and new methods for analyzing in three dimensions the location and ultrastructure of neuronal components, such as synapses, which are distributed non-randomly throughout the brain. DESIGN AND MEASUREMENTS: Volumes are reconstructed by defining transformations that align the entire area of adjacent sections. Whole-field alignment requires rotation, translation, skew, scaling, and second-order nonlinear deformations. Such transformations are implemented by a linear combination of bivariate polynomials. Computer software for generating transformations based on user input is described. Stereological techniques for assessing structural distributions in reconstructed volumes are the unbiased bricking, disector, unbiased ratio, and per-length counting techniques. A new general method, the fractional counter, is also described. This unbiased technique relies on the counting of fractions of objects contained in a test volume. A volume of brain tissue from stratum radiatum of hippocampal area CA1 is reconstructed and analyzed for synaptic density to demonstrate and compare the techniques. RESULTS AND CONCLUSIONS: Reconstruction makes practicable volume-oriented analysis of ultrastructure using such techniques as the unbiased bricking and fractional counter methods. These analysis methods are less sensitive to the section-to-section variations in counts and section thickness, factors that contribute to the inaccuracy of other stereological methods. In addition, volume reconstruction facilitates visualization and modeling of structures and analysis of three-dimensional relationships such as synaptic connectivity.

NASA Discipline Neuroscience↗