Design and analysis of a rate augmented digital-to-analog converter
Rate augmented digital-to-analog converter for computed time-dependent data
COMPUTER SIMULATION↗
SEARCH · Engineering Papers
Results for “TIME DEPENDENCE”
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.
Predicting Saturated Hydraulic Conductivity over Time for Degrading Saltstone Vault Concrete - 20376
At the Savannah River Site (SRS), low-level liquid tank waste is mixed with cementitious materials to produce a solid, monolithic waste form referred to as saltstone. The liquid saltstone mixture is pumped into massive concrete structures referred to as saltstone disposal units (SDUs) for curing and long-term storage. SDUs are vaults made of high performance concrete forming a barrier between saltstone and the surrounding environment. Predicting the long-term performance of SDU vault concrete, as well as the saltstone material itself, is critical to understanding the potential releases of radionuclides to the surrounding environment as described in the Saltstone Disposal Facility (SDF) performance assessment (PA). The hydraulic conductivity profile for SDU vault concrete was conservatively assumed to degrade linearly from the initial saturated hydraulic conductivity of the intact concrete over a long period of time. The initial hydraulic properties of saturated SDU vault concrete have been measured, and the final hydraulic properties of degraded vault concrete are assumed to be the same as the porous backfill material that will surround the vault after closure. The assumed final condition of the vault concrete is also conservative in that it provides a greater estimate of water penetration through the vault than would generally be expected. Therefore, insights that provide more realistic but conservative estimates as to how degradation processes may change the hydraulic properties of the SDU vault concrete, the time dependent relationship of the vault concrete hydraulic properties to degradation (and thus time), and the final hydraulic properties of the vault concrete provides useful additions to the SDF PA. This paper focuses on the assumption that the saltstone vault concrete saturated hydraulic conductivity would decrease linearly as a function of degradation time; this assumption resulted from a worst-case, weighted (by relative layer thickness), arithmetic average of intact and backfill material hydraulic conductivities representing flow parallel to a layered system. This research indicates that the effective saturated hydraulic conductivity of the degraded SDU vault concrete should be based on the relative thickness-weighted geometric mean of the intact and degraded layers hydraulic conductivities. This alternative, which is supported extensively in the relevant literature and evaluated in this paper, leads to a more realistic and defensible estimate of the effective saturated hydraulic conductivity in degraded vault concrete while remaining conservative. The impact of using the more realistic estimate of the effective saturated hydraulic conductivity for the degraded concrete decreases the effective saturated hydraulic conductivity by up to several orders of magnitude for early times (i.e., when the SDU vault concrete is assumed degraded via sulfate attack and carbonation). This revised estimate of the effective saturated hydraulic conductivity has been implemented in the recently revised SDF PA. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Constituent-Based Life Models Being Developed for SiC/SiC Composites
For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.
DiCarlo, James A.↗
Numerical solution of the compressible Navier-Stokes equations using density gradients as additional dependent variables
Numerical solution of two dimensional, time dependent, compressible viscous Navier-Stokes equations about arbitrary bodies was treated using density gradients as additional dependent variables. Thus, six dependent variables were computed with the SOR iteration method. Besides formulation for pressure gradient terms, a formulation for computing the body density was presented. To approximate the governing equations, an implicit finite difference method was employed. In computing the solution for the flow about a circular cylinder, a problem arose near the wall at both stagnation points. Thus, computations with various conditions were tried to examine the problem. Also, computations with and without formulations are compared. The flow variables were computed on 37 by 40 field first, then on an 81 by 40 field.
Kwon, J. H.↗
The Coordination Behaviour of Cu I Photosensitizers Bearing Multidentate Ligands Investigated by X-ray Absorption Spectroscopy
A systematic series of four novel homo- and heteroleptic CuI photosensitizers based on tetradentate 1,10-phenanthroline ligands of the type X^N^N^X containing two additional donor moieties in the 2,9-position (X=SMe or OMe) were designed. Their solid-state structures were assessed by X-ray diffraction. Cyclic voltammetry, UV-vis absorption, emission and X-ray absorption spectroscopy were then used to determine their electrochemical, photophysical and structural features in solution. Following, time-resolved X-ray absorption spectroscopy in the picosecond time scale, coupled with time-dependent density functional theory calculations, provided in-depth information on the excited state electron configurations. For the first time, a significant shortening of the Cu–X distance and a change in the coordination mode to a pentacoordinated geometry is shown in the excited states of the two homoleptic complexes. These findings are important with respect to a precise understanding of the excited state structures and a further stabilization of this type of photosensitizers.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities
We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Statistical learning for predicting density–matrix-based electron dynamics
In this study, we consider the problem of learning density-dependent molecular Hamiltonian matrices from time series of electron density matrices, all in the context of Hartree–Fock theory. Prior work developed a solution to this problem for small molecular systems with density and Hamiltonian matrices of size at most 6 × 6. Here, using a battery of techniques, we scale prior methods to larger molecular systems with, for example, 29 × 29 matrices. This includes systems that either have more electrons or are expressed in large basis sets such as 6-311++G**. Scaling the method to larger systems enhances its relevance for realistic applications in chemistry and physics. To achieve this scaling, we apply dimensionality reduction, ridge regression and analytic computation of Hessians. Through the combination of these techniques, we are able to learn Hamiltonians by minimizing an objective function that encodes local propagation error. Importantly, these learned Hamiltonians can then be used to predict electron dynamics for thousands of steps: When we use our learned Hamiltonians to numerically solve the time-dependent Hartree–Fock equation, we obtain predicted dynamics that are in close quantitative agreement with ground truth dynamics. This includes field-off trajectories similar to the training data and field-on trajectories outside of the training data.
59 BASIC BIOLOGICAL SCIENCES↗
Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities
Here we propose a high dimensional Bayesian inference framework for learning heterogeneous dynamics of a COVID-19 model, with a specific application to the dynamics and severity of COVID-19 inside and outside long-term care (LTC) facilities. We develop a heterogeneous compartmental model that accounts for the heterogeneity of the time-varying spread and severity of COVID-19 inside and outside LTC facilities, which is characterized by time-dependent stochastic processes and time-independent parameters in ~ 1500 dimensions after discretization. To infer these parameters, we use reported data on the number of confirmed, hospitalized, and deceased cases with suitable post-processing in both a deterministic inversion approach with appropriate regularization as a first step, followed by Bayesian inversion with proper prior distributions. To address the curse of dimensionality and the ill-posedness of the high-dimensional inference problem, we propose use of a dimension-independent projected Stein variational gradient descent method, and demonstrate the intrinsic low-dimensionality of the inverse problem. We present inference results with quantified uncertainties for both New Jersey and Texas, which experienced different epidemic phases and patterns. Moreover, we also present forecasting and validation results based on the empirical posterior samples of our inference for the future trajectory of COVID-19.
42 ENGINEERING↗
Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations
A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.
97 MATHEMATICS AND COMPUTING↗
Volumetric Rendering on Wavelet-Based Adaptive Grid
Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.
Vezolainen, Alexei V.↗
Neural Network Approaches for Mobile Spectroscopic Gamma-Ray Source Detection
Artificial neural networks (ANNs) for performing spectroscopic gamma-ray source identification have been previously introduced, primarily for applications in controlled laboratory settings. To understand the utility of these methods in scenarios and environments more relevant to nuclear safety and security, this work examines the use of ANNs for mobile detection, which involves highly variable gamma-ray background, low signal-to-noise ratio measurements, and low false alarm rates. Simulated data from a 2” × 4” × 16” NaI(Tl) detector are used in this work for demonstrating these concepts, and the minimum detectable activity (MDA) is used as a performance metric in assessing model performance.In addition to examining simultaneous detection and identification, binary spectral anomaly detection using autoencoders is introduced in this work, and benchmarked using detection methods based on Non-negative Matrix Factorization (NMF) and Principal Component Analysis (PCA). On average, the autoencoder provides a 12% and 23% improvement over NMF- and PCA-based detection methods, respectively. Additionally, source identification using ANNs is extended to leverage temporal dynamics by means of recurrent neural networks, and these time-dependent models outperform their time-independent counterparts by 17% for the analysis examined here. The paper concludes with a discussion on tradeoffs between the ANN-based approaches and the benchmark methods examined here.
Bilton, Kyle J. (ORCID:0000000184553689)↗
Evolutionary clock - Nonconstancy of rate in different species.
By using various methods for comparing polypeptide sequences we find that the evolutionary divergence of rattlesnake cytochrome c from cytochromes c of species in other classes has been more rapid than that of cytochrome c of another reptile, the snapping turtle. This suggests that the evolutionary rate of change of cytochromes c is species-dependent as well as time-dependent.
Jukes, T. H.↗
Pull-out fibers from composite materials at high rate of loading
Numerical and experimental results are presented on the pullout phenomenon in composite materials at a high rate of loading. The finite element method was used, taking into account the existence of a virtual shear deformation layer as the interface between fiber and matrix. Experimental results agree well with those obtained by the finite element method. Numerical results show that the interlaminar shear stress is time dependent, in addition, it is shown to depend on the applied load time history. Under step pulse loading, the interlaminar shear stress fluctuates, finally decaying to its value under static loading.
Amijima, S.↗
Equivalent source modeling of the main field using Magsat data
Magsat dipole solution models were obtained with both 32 deg and 21 deg resolution based on a data set extending over 4 months. Time dependence was modeled using first time derivatives for the dipole magnetization vector components. This doubles the n total number of parameters in the solution. The solutions displayed a very slow convergence in time derivatives although at each iteration the conversion of the dipole parameters to spherical harmonic coefficients showed close agreement with the MGST (12/80) spherical harmonic model to degree 13 in the constant terms and 8 in the secular variation terms. The program error in the option to simultaneously estimate observatory anomaly biases was found and corrected. A 32 deg resolution dipole model was generated using the selected magnetic observatory data from 1960-1977 used in the GSFC (9/80) spherical harmonic model.
Source record↗
A Lagrangian method for a two layer shallow water model
One of the biggest drawbacks of Eulerian methods is the presence of the nonlinear convective terms which leads to inaccurate representations of advection and fronts. An alternative and conceptually simpler approach is the Lagrangian approach in which the fluid particles themselves are tracked and equations are derived based on local spatial interaction. In this formulation the nonlinear convective terms do not appear. The main drawback with the Lagrangian formulation, however, is that the local spatial interactions are time dependent. That is, at each time level, one must know the current neighbors of a given particle to accurately compute the forces acting on the particle. The Voroni Mesh which is continuously deforming and always linking nearest neighbors is needed. The effectiveness of a simple Lagrangian method for simulating the onset of a Kelvin-Helmholtz instability in a 2 layer fluid is demonstrated.
Augenbaum, J. M.↗
Plume formation in the D-double prime-layer and the roughness of the core-mantle boundary
A series of very-high-resolution finite element calculations of plume formation in the D-double prime-layer has been performed for several plausible rheologies and boundary conditions in order to study both the early and later stages of boundary layer development. The results show that plumes are initiated by coalescence of small-scale convective instabilities within the low-viscosity region immediately above the core-mantle boundary (CMB). These instabilities support topographic roughness on the CMB having horizontal scales of 20-50 km and provide a source for scattered P-waves seen as precursors to the phases PKIKP and PKKP. The calculated structure of fully developed plumes emerging from the D-double prime-layer consists of 5-50 cm/yr flow confined to 50-100 km thick vertical conduits. With strongly temperature-dependent viscosity, plumes exhibit time-dependent behavior, including upward propagating solitary conduit waves, which may contribute to episodicity in hotspot volcanism.
Olson, Peter↗
Temperature Fluctuations During Crystal Growth
Technique developed to deconvolve period and relative amplitude of fluctuations of heat flow inBridgman crystal growth. Temperature-measuring device with enough sensitivity and frequency response to make desired measurements inserted close as possible to substance monitored. Time-dependent temperature response recorded, and time domain response converted to frequency domain for further analysis. Fast fourier transform (FFT) of data on temperature oscillations shows particular behavior at some specific frequencies corresponding with striations, or defects observed on crystal. Useful procedure to help determine sources of some growth-induced crystalline defects. Other processes sensitive to small temperature fluctuations, such as diffusion, precipitation, and corrsion, benefit from technique.
Fripp, Archibald L., Jr.↗
Turbulent supersonic convection in three dimensions
Previous numerical calculations of two-dimensional, compressible convection are extended to three dimensions, using a higher order Godunov scheme. The results show that the flow readily becomes supersonic in the upper boundary layer, where shock structures form intermittently in the vicinity of the strong downflow lanes. The convection as a whole is strongly time-dependent and evolves on a time scale comparable to the sound crossing time. The motions in the upper layers are characterized by the rapid expansion of the upward-moving fluid elements. In the interior, most of the heat is carried by a small fraction of the fluid residing in strong, highly coherent downflows. The remaining fluid is dominated by small-scale, disorganized turbulent motions.
Malagoli, Andrea↗