Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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,243 records · Page 69

Seismic modeling with radial basis function-generated finite differences (RBF-FD) – a simplified treatment of interfaces

In a previous study of seismic modeling with radial basis function-generated finite differences (RBF-FD), we outlined a numerical method for solving 2-D wave equations in domains with material interfaces between different regions. The method was applicable on a mesh-free set of data nodes. It included all information about interfaces within the weights of the stencils (allowing the use of traditional time integrators), and was shown to solve problems of the 2-D elastic wave equation to 3rd-order accuracy. In the present paper, we discuss a refinement of that method that makes it simpler to implement. It can also improve accuracy for the case of smoothly-variable model parameter values near interfaces. We give several test cases that demonstrate the method solving 2-D elastic wave equation problems to 4th-order accuracy, even in the presence of smoothly-curved interfaces with jump discontinuities in the model parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗

Correlation of macro and micro cardiovascular function during weightlessness and simulated weightlessness

The investigation of cardiovascular function necessarily involves a consideration of the exchange of substances at the capillary. If cardiovascular function is compromised or in any way altered during exposure to zero gravity in space, then it stands to reason that microvascular function is also modified. We have shown that an increase in cardiac output similar to that reported during simulated weightlessness is associated with a doubling of the number of post-capillary venules and a reduction in the number of arterioles by 35%. If the weightlessness of space travel produces similar changes in cardiopulmonary volume and cardiac output, a reasonable expectation is that astronauts will undergo venous neovascularization. We have developed an animal model in which to correlate microvascular and systemic cardiovascular function. The microcirculatory preparation consists of a lightweight, thermo-neutral chamber implanted around intact skeletal muscle on the back of a rat. Using this technique, the performed microvasculature of the cutaneous maximus muscle may be observed in the conscious, unanesthetized animal. Microcirculatory variables which may be obtained include venular and arteriolar numbers, lengths and diameters, single vessel flow velocities, vasomotion, capillary hematocrit anastomoses and orders of branching. Systemic hemodynamic monitoring of cardiac output by electromagnetic flowmetry, and arterial and venous pressures allows correlation of macro- and microcirculatory changes at the same time, in the same animal. Observed and calculated hemodynamic variables also include pulse pressure, heart rate, stroke volume, total peripheral resistance, aortic compliance, minute work, peak aortic flow velocity and systolic time interval. In this manner, an integrated assessment of total cardiovascular function may be obtained in the same animal without the complicating influence of anesthetics.

NASA Discipline Number 00-00↗

Analytical Model for the Higher Order Moments of Midlatitude Atmospheric Temperature Distributions

Observed distributions of atmospheric temperature are non-Gaussian. Therefore, moments beyond variance are necessary in determining the frequency of extreme temperature events. Here we propose a simple kinematic model for atmospheric mid-latitude temperature variability based on symmetric advection from a non-symmetric background temperature profile. We then use this model to derive analytical expressions for the higher order moments of temperature distributions. Our results show that nonzero skewness and kurtosis arise due to the nonlinearity of the time-mean meridional temperature profile. The analytical model matches an idealized Held-Suarez atmospheric model, indicating nonlinearity of time-mean temperature in latitude is the dominant contribution to nonzero skewness and kurtosis in synoptic temperature variations. Model analysis further shows decrease in higher order moments due to climate change come roughly equally from changes in mixing length and changes in the background temperature profiles.

54 ENVIRONMENTAL SCIENCES↗

Cloud Giovanni: Reining in Costs and Improving Performance with Analytical Data Stores Using Scalable Serverless Architecture

Giovanni is the Geospatial Interactive Online Visualization ANd aNalysis Infrastructure developed at NASA GES DISC which provides a simple and intuitive way to visualize, analyze, and access vast amounts of Earth science data. It receives large number of user requests each day for a variety of analysis and visualization services, which leads to the big data challenge of serving gradually increasing large data volumes with diverse statistical algorithms. We hereby propose a multi-dimensional accumulation method which provides fast and cost-efficient cloud analysis for diverse services including both area averaging and time averaging. This method involves the weighted volume integration over multiple variable dimensions (time and space), and is implemented in AWS using Athena providing serverless and highly scalable data analysis. Compared to the standard method, this approach dramatically reduced the computational time by order of magnitude with a minimal AWS cost incurred. For example, for a benchmark of 10-year area averaging over the 1x1 degree daily variable, the computational time was reduced from minutes to seconds, and the Athena cost is only $5 for 100,000 requests.

Zhang, Hailiang↗

Statistical Analysis of Risk Factors in the Prebreathe Reduction Protocol

The 165 exposures from four 2-hour protocols were analyzed for correlations or trends between decompression sickness (DCS) or venous gas emboli (VGE), and variables that affect risk in the subject and astronaut populations. The assumption in this global survey is that the distributions of gender, age, body mass index, etc., are equally represented in all four tested procedures. We used Student t-test for comparisons between means and chi-square test between comparisons of proportions with p<0.05 defining a significant level. The type and distribution of the 19 cases of DCS were similar to historical cases. There was no correlation of age, gender, body mass index or fitness level with greater incidence of DCS or VGE. However increased age was associated with more Grade IV VGE in males. The duration and quantity of exercise during prebreathe is inversely related to risk of DCS and VGE. The latency time for VGE was longer (103 min +/- 56 SD, n = 15) when the ergometry was done approximately 15 min into the prebreathe than when done at the start of the prebreathe (53 min +/- 31, n = 13). The order of the ergometry did not influence the overall DCS and VGE incidence. We identified variables other than those of the prebreathe procedures that influence the DCS and VGE outcome. The analysis suggests that males over 40 years have a high incidence of Grade IV VGE.

Gerth, Wayne A.↗

Measurements of the pp → WZ inclusive and differential production cross sections and constraints on charged anomalous triple gauge couplings at $\sqrt{s} $ = 13 TeV

The WZ production cross section is measured in proton-proton collisions at a centre-of-mass energy $\sqrt{s} =$ 13 TeV using data collected with the CMS detector, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The inclusive cross section is measured to be $\sigma_{\text{tot}}$(pp $\to$ WZ$)$ = 48.09 $^{+1.00}_{-0.96}$ (stat) $^{+0.44}_{-0.37}$ (theo) $^{+2.39}_{-2.17}$ (syst) $\pm$ 1.39 (lumi) pb, resulting in a total uncertainty of $-$2.78/$+$2.98 pb. Fiducial cross section and ratios of charge-dependent cross section measurements are provided. Differential cross section measurements are also presented with respect to three variables: the Z boson transverse momentum $p_\mathrm{T}$, the leading jet $p_\mathrm{T}$, and the $m$(WZ) variable, defined as the invariant mass of the system composed of the three leptons and the missing transverse momentum. Differential measurements with respect to the W boson $p_\mathrm{T}$, separated by charge, are also shown. Results are consistent with standard model predictions, favouring next-to-next-to-leading-order predictions over those at next-to-leading order. Constraints on anomalous triple gauge couplings are derived via a binned maximum likelihood fit to the $m$(WZ) variable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗

Progress Towards a Neutral Current $\pi^0$ Cross Section Analysis in the NOvA Near Detector

The NOvA neutrino experiment is attempting to measure properties of neutrinos in order to figure out information about the universe. To detect the signal neutrino interactions, we must determine methods to identify and isolate background events. Research focused on a specific background interaction called a single prong neutral current $\pi^0$ interaction. To do this, a basic cuts based analysis was performed, followed by feeding data into a multi-variate analysis package using a boosted decision tree (BDT) algorithm. Using the BDT, a a new variable was generated which separates signal and background very efficiently. Further work must still be done in order to continue improving the performance of the BDT. This research is valuable to the field of studying neutrino cross sections as it is a background which will always be present in this type of analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Regional variability in the environmental controls of precipitation regimes in the tropics

The environmental factors that control precipitation regimes in Boreal winter rainfall in the tropics and their regional variabilities are examined using a simple statistical analysis and a machine learning model. Radar-derived precipitation from field campaigns at Darwin Australia, Manaus Brazil, and the Equatorial Indian Ocean, along with the corresponding large-scale environmental variables from ERA5 reanalysis, are used. The dependence of marginal distributions of frequencies of these regimes on five environmental variables are calculated. The variables are hourly column integrated precipitable water (PW), convective available potential energy (CAPE), convective inhibition (CIN), and lower and upper tropospheric wind shear. The simple machine learning model that predicts the probability of transition from suppressed to an active regime as a function of the environmental variables is designed and optimized for a potential application as a trigger function for convection parameterizations. To the first order the analysis shows an abrupt increase in the probability of an active regime near PW > 60 mm and CIN of <100 J kg -1 . The key differences in the frequencies of active regimes among the regions are found to be related to the fact that over Darwin there is strong variability in PW while CIN is generally low. Over Amazon, on the other hand, both PW and CIN are quite variable. Over the DYNAMO domain the comparatively frequent low PW is compensated for by consistently low CIN resulting in a moderate frequency of an active regime.

54 ENVIRONMENTAL SCIENCES↗

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗

Improving Land-Surface Model Hydrology: Is an Explicit Aquifer Model Better than a Deeper Soil Profile?

Land surface models (LSMs) are computer programs, similar to weather and climate prediction models, which simulate the storage and movement of water (including soil moisture, snow, evaporation, and runoff) after it falls to the ground as precipitation. It is not currently possible to measure all of the variables of interest everywhere on Earth with sufficient accuracy. Hence LSMs have been developed to integrate the available information, including satellite observations, using powerful computers, in order to track water storage and redistribution. The maps are used to improve weather forecasts, support water resources and agricultural applications, and study the Earth's water cycle and climate variability. Recently, the models have begun to simulate groundwater storage. In this paper, we compare several possible approaches, and examine the pitfalls associated with trying to estimate aquifer parameters (such as porosity) that are required by the models. We find that explicit representation of groundwater, as opposed to the addition of deeper soil layers, considerably decreases the sensitivity of modeled terrestrial water storage to aquifer parameter choices. We also show that approximate knowledge of parameter values is not sufficient to guarantee realistic model performance: because interaction among parameters is significant, they must be prescribed as a harmonious set.

Gulden, L. E.↗

Polarizabilities of two-electron positive ions and Rydberg levels of lithium

Second-order sums are computed for the two-electron positive ions Li(+), Be(2+), and B(3+) by representing the intermediate states as discrete pseudostates. The Hamiltonian for the system is given and translated to semi-Jacobi variables for its derivation in reduced Rydberg units. The interaction potential is developed in its multipole form and employed in second-order perturbation theory. The perturbation is always a multipole-potential term, and the initial and intermediate states are constructed from two-particle Hylleraas basis sets. The resulting pseudostates are shown to represent the continuum with good convergence, and the asymptotic-optical-potential method is applied to the two-body quantities to determine the energy levels of three-electron systems in high Rydberg states. The method is shown to give accurate values for the polarizabilities of the two-electron isoelectronic systems.

Bhatia, A. K.↗

Estimation of surface heat and moisture fluxes over a prairie grassland. II - Two-dimensional time filtering and site variability

The behavior of in situ measurements of surface fluxes obtained during FIFE 1987 is examined by using correlative and spectral techniques in order to assess the significance of fluctuations on various time scales, from subdiurnal up to synoptic, intraseasonal, and annual scales. The objectives of this analysis are: (1) to determine which temporal scales have a significant impact on areal averaged fluxes and (2) to design a procedure for filtering an extended flux time series that preserves the basic diurnal features and longer time scales while removing high frequency noise that cannot be attributed to site-induced variation. These objectives are accomplished through the use of a two-dimensional cross-time Fourier transform, which serves to separate processes inherently related to diurnal and subdiurnal variability from those which impact flux variations on the longer time scales. A filtering procedure is desirable before the measurements are utilized as input with an experimental biosphere model, to insure that model based intercomparisons at multiple sites are uncontaminated by input variance not related to true site behavior. Analysis of the spectral decomposition indicates that subdiurnal time scales having periods shorter than 6 hours have little site-to-site consistency and therefore little impact on areal integrated fluxes.

Crosson, William L.↗

Prediction of Probabilistic Detonation Threshold via Millimeter-Scale Microstructure-Explicit and Void-Explicit Simulations

We present an approach and relevant models for predicting the probabilistic shoc-to-detonation transition (SDT) behavior and Pop plot (PP) of heterogeneous energetic materials (HEM) via mesoscopic microstructure-explicit (ME) and void explicit (VE) simulations at the millimeter (mm) sample size scale. Although the framework here is general, the particular material considered in this paper is pressed Octahydro-1,3,5,7-tetranitro-1,2,3,5-tetrazocine (HMX). To systematically delineate the effects of material heterogeneities, four material cases are considered. These cases are homogeneous material, material with granular microstructure but no voids, homogeneous material with voids, and material with both granular microstructure and voids. Statistically equivalent microstructure sample sets (SEMSS) are generated and used. Eulerian hydrocode simulations explicitly resolve the material heterogeneities, voids, and the coupled mechanical-thermal-chemical processes. In particular, it is found that both microstructure and voids strongly influence the SDT behavior and PP. The effects of different combinations of microstructure heterogeneity and voids on the SDT process and PP are quantified and rank-ordered. The overall framework uses the Mie–Grüneisen equation of state and a history variable reactive burn model (HVRB). A novel probabilistic representation for quantifying the PP is developed, allowing the calculation of (1) the probability of observing SDT at a given combination of shock pressure and run distance, (2) the run-distance to detonation under a given combination of shock pressure and prescribed probability, and (3) the shock pressure required for achieving SDT at a given run distance with a prescribed probability. The results are in agreement with general trends in experimental data in the literature.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Superstratum symbiosis

Superstrata are smooth horizonless microstate geometries for the supersymmetric D1-D5-P black hole in type IIB supergravity. In the CFT, “superstratum states” are defined to be the component of the supergraviton gas that is obtained by breaking the CFT into “|00&rang;-strands” and acting on each strand with the “small,” anomaly-free superconformal generators. We show that the recently-constructed supercharged superstrata represent a final and crucial component for the construction of the supergravity dual of a generic superstratum state and how the supergravity solution faithfully represents all the coherent superstratum states of the CFT. For the supergravity alone, this shows that generic superstrata do indeed fluctuate as functions of three independent variables. Smoothness of the complete supergravity solution also involves “coiffuring constraints” at second-order in the fluctuations and we describe how these lead to new predictions for three-point functions in the dual CFT. We use a hybrid of the original and supercharged superstrata to construct families of single-mode superstrata that still have free moduli after one has fixed the asymptotic charges of the system. We also study scalar wave perturbations in a particular family of such solutions and show that the mass gap depends on the free moduli. This can have interesting implications for superstrata at non-zero temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupling Approaches with Non-matching Grids for Classical Linear Elasticity and Bond-based Peridynamic Models in 1D

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. To facilitate the coupling of the two models, non-matching grids are often desirable as nonlocal grids usually require a finer resolution than local grids. In that case, it is often convenient to resort to interpolation operators so that models can exchange information in the overlap regions when nodes from the two grids do not coincide. This paper studies three existing coupling approaches, namely 1) a method that enforces matching displacements in an overlap region, 2) a variant that enforces a constraint on the stresses instead, and 3) a method that considers a variable horizon in the vicinity of the interfaces. Further, the effect of the interpolation order and of the grid ratio on the performance of the three coupling methods with non-matching grids is carefully studied on one-dimensional examples using polynomial manufactured solutions. The numerical results show that the degree of the interpolants should be chosen with care to avoid introducing additional modeling errors, or simply minimize these errors, in the coupling approach.

97 MATHEMATICS AND COMPUTING↗

The quantum mechanics of high-order kinematic values

In this paper, within the framework of a unified mathematical model, the new formulation of quantum mechanics – quantum mechanics of higher order kinematic values – is proposed. In contrast to the well-known formulations of quantum mechanics in the new formulation, the wave function has not only a coordinate or momentum representation, but also representations through acceleration and accelerations of higher orders. The representations of the wave function in terms of higher order kinematic values make it possible to consider the dissipative systems problem and the variable entropy systems one within a single mathematical apparatus. The new formulation is not built phenomenologically, but from first principles and in a particular case contains the formulation of the de Broglie–Bohm &lang;wave-pilot&rang; quantum mechanics. The new formulation is based on the Vlasov equations chain. This approach allows us to obtain a natural connection between classical and quantum systems. Examples of such systems are considered in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗