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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the intial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given, in order to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS2000

Supporting Development of Satellite's Guidance Navigation and Control Software: A Product Line Approach

The NASA Goddard Space Flight Center Flight Software Branch (FSB) is developing a Guidance, Navigation, and Control (GNC) Flight Software (FSW) product line. The demand for increasingly more complex flight software in less time while maintaining the same level of quality has motivated us to look for better FSW development strategies. The GNC FSW product line has been planned to address the core GNC FSW functionality very similar on many recent low/near Earth missions in the last ten years. Unfortunately these missions have not accomplished significant drops in development cost since a systematic approach towards reuse has not been adopted. In addition, new demands are continually being placed upon the FSW which means the FSB must become more adept at providing GNC FSW functionality's core so it can accommodate additional requirements. These domain features together with engineering concepts are influencing the specification, description and evaluation of FSW product line. Domain engineering is the foundation for emerging product line software development approaches. A product line is 'A family of products designed to take advantage of their common aspects and predicted variabilities'. In our product line approach, domain engineering includes the engineering activities needed to produce reusable artifacts for a domain. Application engineering refers to developing an application in the domain starting from reusable artifacts. The focus of this paper is regarding the software process, lessons learned and on how the GNC FSW product line manages variability. Existing domain engineering approaches do not enforce any specific notation for domain analysis or commonality and variability analysis. Usually, natural language text is the preferred tool. The advantage is the flexibility and adapt ability of natural language. However, one has to be ready to accept also its well-known drawbacks, such as ambiguity, inconsistency, and contradictions. While most domain analysis approaches are functionally oriented, the idea of applying the object-oriented approach in domain analysis is not new. Some authors propose to use UML as the notation underlying domain analysis. Our work is based on the same idea of merging UML and domain analysis. Further, we propose a few extensions to UML in order to express variability, and we define precisely their semantics so that a tool can support them. The extensions are designed to be implemented on the API of a popular industrial CASE tool, with obvious advantages in cost and availability of tool support. The paper outlines the product line processes and identifies where variability must be addressed. Then it describes the product line products with respect to how they accommodate variability. The Celestial Body subdomain is used as a working example. Our results to date are summarized and plans for the future are described.

McComas, David↗
DOE OSTI2018

Compact finite difference schemes for the time fractional diffusion equation with nonlocal boundary conditions

Fourth-order compact finite difference schemes for solving the time fractional diffusion equation with nonlocal boundary conditions are considered in this paper. The second-order derivative in one space variable is approximated by a fourth-order compact finite difference. Two different time derivatives, the Caputo derivative and Riemann–Liouville derivative are considered in the equations that we discuss. The Caputo time derivative is discretized by the L1-approximation and a second-order approximation is used for the Riemann–Liouville derivative. Stability analysis using the discrete energy method for Scheme I is given, and error estimation of this scheme is derived by conducting the local truncation error analysis, using the stability result. Numerical results are provided to show that both two fully discrete schemes have fourth-order accuracy in space, and 2−γ order of accuracy in time for Scheme I with γ being the order of the fractional derivative, while Scheme II has second-order accuracy in time when the theoretical solutions are smooth with zero initial values.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Coupled Multimodal Dynamics of Hydrogen-Containing Ion Networks in Water-Deficient, Sodium Hydroxide-Aluminate Solutions

The (meta)stability of low water activity sodium hydroxide/aluminate (Na + OH–/Al(OH) 4 –) electrolytes dictates kinetics in the Bayer process for aluminum refining and high-level nuclear waste processing. We utilized quasi-elastic neutron scattering (QENS) and proton nuclear magnetic resonance spectroscopy ( 1 H NMR) in extremely concentrated sodium aluminate solutions to investigate the picosecond (ps) to microsecond (ms) timescale motions of H-bearing species (Al(OH) 4 – monomers/clusters, OH– and H 2 O). In the QENS data, in contrast to typical liquids, no short-time translational diffusion was observed at 293 K, but two types of localized motions were found: (i) local backbone tumbling or a formation of large hydrated ion clusters on the order of 40–60 ps; and (ii) much slower, complex, and collective dynamics of the ensemble of H-bearing species on the order of 350–750 ps. Variable temperature, pulsed field gradient, diffusion-ordered 1H NMR was used to determine the ensemble translational motion along with relaxometry to calculate rotational correlation coefficients. The ensemble rotational correlation times were on the order of 184–300 ps from 1 H NMR, which is consistent with the timescale of the QENS components. Complementary molecular dynamics simulation of NaOH solutions exhibit extensive ion networks potentially responsible for the observed dynamical coupling of water with the motion of large hydrated ion clusters. Understanding these collective motions will aid in predicting the behavior of complex solutions during aluminum production and during nuclear waste processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Using Cosmic Ray Muons to Assess Geological Characteristics in the Subsurface

Cosmic rays are energetic nuclei and elementary particles that originate from stars and intergalactic events. The interaction of these particles with the upper atmosphere produces a wide range of secondary particles that reach the surface of the earth, of which muons are the most prominent. With enough energy, muons can travel up to a few kilometers beneath the surface of the earth before being stopped completely. The terrestrial muon flux profile and associated zenith angle can be utilized to determine geological characteristics of a location (e.g., rock overburden and density) without having to use conventional methods such as boreholes. This work uses a low-power plastic scintillator-based muon detection system as a prototype for this non-destructive geological assay methodology. Four custom designed 102 cm x 51 cm x 5 cm plastic scintillation panels are used to realize two orthogonal detection planes. Optical photons from each scintillation panel are read using OnSemi J-Series 4x4 silicon photomultiplier (SiPM) arrays in conjunction with preamplifiers. Simultaneous triggers between detectors from two planes indicate a coincidence event which is recorded using the QuarkNet data acquisition system (DAQ) from Fermi National Accelerator Laboratory. A custom detector holder was designed to securely mount the detection system and rotate the panels along the zenith to collect data at variable angles. In order to quantify the systematic uncertainties associated with the detector, such as energy depositions and angular resolution of the detector design, a Monte Carlo (MC) simulation using Geant4 is being developed. Cosmic ray flux prediction will be included in the project by adding the CORSIKA MC code to the simulation toolchain. Simulated and experimental data will drive the development and validation of a reconstruction algorithm that, upon completion, is expected to predict average overburden and rock density. Extended detector exposure to muons can be used as a means to understand changes in the surrounding environment like rock porosity. On the experimental front, muons will initially be measured at the surface, establishing the baseline flux. This is followed by recording the muon flux at variable depths and zenith angles, where the data will be used by the reconstruction algorithm to predict the overburden. The result will be benchmarked against geological surveys. The measured flux data will also be used to benchmark independent and established models. Successful proof-of-concept demonstration of this technology can open doors for long term non-invasive geological monitoring. The detector design, experimental methodology, and the benchmarking efforts are detailed in this work.

Gadey, Harish Reddy↗
DOE OSTI2019

Boundary-Value Problems with Impulsive Conditions for Parabolic Equations with Degenerations

We consider the first boundary-value problem and the unilateral boundary-value problem with impulsive conditions with respect to the time variable for a linear parabolic equation with power singularities of any order in the space variables. With the help of the maximum principle and a priori estimates, we establish the existence and uniqueness of the solutions of the posed problems in the Hölder spaces with power weights.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Sensitivity of Regional Climate to Deforestation in the Amazon Basin

The deforestation results in several adverse effect on the natural environment. The focus of this paper is on the effects of deforestation on land-surface processes and regional climate of the Amazon basin. In general, the effect of deforestation on climate are likely to depend on the scale of the defrosted area. In this study, we are interested in the effects due to deforestation of areas with a scale of about 250 km. Hence, a meso-scale climate model is used in performing numerical experiments on the sensitivity of regional climate to deforestation of areas with that size. It is found that deforestation results in less net surface radiation, less evaporation, less rainfall, and warmer surface temperature. The magnitude of the of the change in temperature is of the order 0.5 C, the magnitudes of the changes in the other variables are of the order of IO%. In order to verify some of he results of the numerical experiments, the model simulations of net surface radiation are compared to recent observations of net radiation over cleared and undisturbed forest in the Amazon. The results of the model and the observations agree in the following conclusion: the difference in net surface radiation between cleared and undisturbed forest is, almost, equally partioned between net solar radiation and net long-wave radiation. This finding contributes to our understanding of the basic physics in the deforestation problem.

Eltahir, Elfatih A. B.↗
DOE OSTI2023

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗
NASA NTRS1976

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitational and rotational terms in the equations are of first order in the space variables, the pressure-gradient terms are of second order, and the turbulent-viscosity term is of third order. The presence of turbulent viscosity ensures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial flow is always inward and allows collapse to occur (axially) even when the rotation is large. An approximate solution of the governing partial differential equations is also given in order to study the spatial distributions of the density and velocity.

Deissler, R. G.↗
DOE OSTI2018

Mass-preserving time second-order explicit–implicit domain decomposition schemes for solving parabolic equations with variable coefficients

In this paper, the new mass-preserving time second-order explicit–implicit domain decomposition (DD) schemes for solving parabolic equations with various coefficients are proposed. In the schemes, first, the interface values on the interfaces of sub-domains are explicitly calculated, in which two kinds of explicit difference schemes are used. Second, the interior values in sub-domains are determined by the implicit schemes. Third, the interface values are recalculated explicitly to preserve mass. The important feature is that the developed schemes are more than first order in time step and the correction operators at the third step ensure the schemes’ mass conservative. The stability and convergence are analyzed for the two schemes. The optimal error estimates are obtained in the H{sup 1} semi-norm under a stability condition. Numerical experiments confirm the theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A fast and accurate method for perturbative resummation of transverse momentum-dependent observables

Here, we propose a novel strategy for the perturbative resummation of transverse momentum-dependent (TMD) observables, using the q T spectra of gauge bosons (γ*, Higgs) in pp collisions in the regime of low (but perturbative) transverse momentum q T as a specific example. First we introduce a scheme to choose the factorization scale for virtuality in momentum space instead of in impact parameter space, allowing us to avoid integrating over (or cutting off) a Landau pole in the inverse Fourier transform of the latter to the former. The factorization scale for rapidity is still chosen as a function of impact parameter b, but in such a way designed to obtain a Gaussian form (in ln b) for the exponentiated rapidity evolution kernel, guaranteeing convergence of the b integral. We then apply this scheme to obtain the q T spectra for Drell-Yan and Higgs production at NNLL accuracy. In addition, using this scheme we are able to obtain a fast semi-analytic formula for the perturbative resummed cross sections in momentum space: analytic in its dependence on all physical variables at each order of logarithmic accuracy, up to a numerical expansion for the pure mathematical Bessel function in the inverse Fourier transform that needs to be performed just once for all observables and kinematics, to any desired accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A variational multiscale immersed meshfree method for heterogeneous materials

Abstract We introduce an immersed meshfree formulation for modeling heterogeneous materials with flexible non-body-fitted discretizations, approximations, and quadrature rules. The interfacial compatibility condition is imposed by a volumetric constraint, which avoids a tedious contour integral for complex material geometry. The proposed immersed approach is formulated under a variational multiscale based formulation, termed the variational multiscale immersed method (VMIM). Under this framework, the solution approximation on either the foreground or the background can be decoupled into coarse-scale and fine-scale in the variational equations, where the fine-scale approximation represents a correction to the residual of the coarse-scale equations. The resulting fine-scale solution leads to a residual-based stabilization in the VMIM discrete equations. The employment of reproducing kernel (RK) approximation for the coarse- and fine-scale variables allows arbitrary order of continuity in the approximation, which is particularly advantageous for modeling heterogeneous materials. The effectiveness of VMIM is demonstrated with several numerical examples, showing accuracy, stability, and discretization efficiency of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A molecular dynamics survey study of impact-induced amorphization in yttria nanoparticles

The impact of a single yttria (Y 2 O 3 ) nanoparticle onto a yttria substrate was studied as a function of particle velocity (300–1200 m/s) and diameter (5–12 nm) for three particle orientations and two substrate orientations using molecular dynamics simulations. Impacts resulted in minimal dislocation motion but extensive disordering of the crystal lattice without approaching the melting temperature. The resulting amorphous microstructure was both precisely coincident with plastic strain and retained until the end of the simulation, suggesting a deformation mechanism of amorphization-facilitated viscous flow. Further analysis revealed that the degree of amorphization depended primarily on the impact velocity and secondarily on the orientations of the particle and substrate, with significant second-order effects from variable interactions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Relationships among diurnal variations of polar night cloud, precipitation, surface temperatures, and the fair-weather return current of the global electric circuit (GEC)

Over the past five years (June 2017-current), the vertical electric field (E z ) as well as numerous cloud, precipitation and radiation properties have been monitored at the Department of Energy-Atmospheric Radiation Measurement (DOE ARM) North Slope of Alaska (NSA) field site. Comparisons between the composite diurnal averaged fair-weather E z , and composite ceilometer derived cloud base height during the polar night, reveal a significant correlation between the parameters (r = 0.62), supporting previous studies that there is high correlation between local electric field and cloud properties. With the use of extensive instrumentation at the site, such as the Micro Pulse Lidar (MPL), Ka-band Zenith Radar (KAZR), ceilometer, SKYRAD, Precipitation Imaging Package (PIP), among others, this study provides a more comprehensive examination of the diurnal cycle of cloud and precipitation properties along with the localized fair-weather return current of the larger Global Electric Circuit (GEC) system. Comparisons between the composite diurnal averaged fair-weather E z , and cloud thickness, maximum column backscatter, and precipitation particle counts all show similar diurnal variability during the polar night, indicating that during the largest magnitude fair weather E z time periods, clouds bases tend to be higher, clouds are thicker, have a larger column backscatter, and display more precipitating particles. Furthermore, a slight diurnal variability in the polar night surface temperature was found to be highly correlated (r = 0.87) to the longwave downwelling irradiance measured by SKYRAD, indicating that the variations in the physical properties of local clouds could modulate the diurnal polar night surface temperature variability on the order of 0.5 °C/day. In conclusion, these findings emphasize the importance and global nature of the GEC system, with the global aggregate of thunderstorms and electrified clouds potentially influencing polar night cloud properties as well as diurnal wintertime polar surface temperature variation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spectral Fidelity of Earth's Terrestrial and Aquatic Ecosystems

The Surface Biology and Geology (SBG) investigation will create global maps of spectral surface reflectance and emissivity at a cadence of 16 days or better, with coverage to address global questions about Earth's geology, cryosphere and ecosystems. The revolutionary potential poses a commensurate challenge: creating contiguous maps free from regional biases induced by atmosphere, observation geometry, or inversion error. This will require an accurate calibration with precise knowledge of each channel's spectral response. In this work, we quantify the impact of spectral calibration on SBG's aquatic and terrestrial ecosystem objectives. We find that contemporary algorithms for ecosystem trait retrieval demand more accurate spectral calibration than historical missions. Errors due to drift or spatial nonuniformity in the wavelength calibration that have previously been considered acceptable can cause systematic errors larger than the instrument noise and of the same order as the variability SBG aims to measure. Moreover, their impact on atmospheric correction can induce climate-dependent systematic errors that thwart comparisons between ecosystems. These results underscore the importance of spectral response accuracy in SBG mission design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Probing the influence of ion-pairing on ligand-field excited-state dynamics

Exploration of the photophysical and photochemical properties of transition metal complexes has driven ground-breaking advancements in solar energy conversion technologies, including photoredox catalysis. While significant research has been devoted to understanding excited state properties of second- and third-row transition metal complexes, earth-abundant first-row metal complexes have received comparatively little attention in this context until very recently. In particular, the role of ion-pairing – which has been identified as a potentially significant factor for Ir(III)-based photosensitizers – has not been examined with regard to its influence on the ligand-field excited states that dominate much of first-row photophysics. A key challenge in studying ion-pair interactions lies in quantifying the extent and nature of ion-pairing, particularly in non-aqueous media where the vast majority of photophysical studies are performed. Cobalt(III) polypyridyl complexes provide an attractive platform to address such questions due to their demonstrated potential for applications in photoredox catalysis involving ligand-field excited states. In the present study, we prepared a cobalt(III) polypyridyl complex, [Co(4,4′-OMebpy) 3 ](BAr F 4 ) 3 (where 4,4′-OMebpy is 4,4′-dimethoxy-2,2′-bipyridine and BAr F 4 is tetrakis(3,5-bis(trifluoromethyl)-phenyl)borate) to probe ion-pairing in non-aqueous solutions. Specifically, analysis of data acquired from both variable-temperature diffusion ordered spectroscopy (DOSY) NMR and 1-D rotating-frame nuclear Overhauser effect (ROE) experiments allowed us to identify and differentiate between solvent-separated ion pairs in high-dielectric media and contact ion pairs in a low-dielectric solvent. Time-resolved absorption spectroscopy was then used to measure ground-state recovery dynamics under these varying conditions of ion-pairing, the results of which revealed an increase in excited-state lifetime for contact ion-pairs that we suggest arises from a reduction in outer-sphere reorganization energy relative to conditions which favored solvent-separated ion pairs. We believe this study demonstrates that one can leverage broadly available NMR-based methods to understand ion-pairing in non-aqueous solutions, which in turn can provide a microscopic picture of intermolecular interactions that can impact the photophysical properties of transition metal-based chromophores.

Ghosh, Atanu [Michigan State University, East Lans↗
DOE OSTI2019

Modeling memoryless degradation under variable stress

Accelerated degradation tests can be used as the basis for predicting the performance or state of health of products and materials at use conditions over time. Measurements acquired at accelerated levels of stress are used to develop models that relate to the degradation of one or more performance measures. Frequently, products/materials of interest are subjected to variable stress levels during their lifetimes. However, testing is usually performed only at a few fixed stress levels. In such cases, cumulative degradation models are developed and assessed by using data acquired under those fixed stress conditions. The degradation rate at any stress condition within the range of the model can be estimated by the derivative of the cumulative model at that stress condition. It follows that, to predict cumulative degradation over variable use conditions, one might integrate the fluctuating degradation rate over time. Existing approaches for doing this consider degradation rates that depend only on the current stress level. Here, we propose to allow the degradation rate to also depend on the current state of health as indicated by the associated performance measure(s). The resulting modeling approach is capable of portraying a broader range of degradation behavior than existing approaches. The assertion of memoryless degradation by using this or any other approach should be assessed experimentally with data acquired under variable stress in order to increase confidence that the integrated rate model is accurate. In this article, we demonstrate the additional capability of the proposed approach by developing empirical memoryless rate-based degradation models to predict resistance increase and capacity decrease in lithium-ion cells that are being evaluated for use in electric vehicles. Here, we then assess the plausibility of these models.

42 ENGINEERING↗