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,027 records · Page 57

Motion tracking in the liver: Validation of a method based on 4D ultrasound using a nonrigid registration technique

Purpose: Treatments like radiotherapy and focused ultrasound in the abdomen require accurate motion tracking, in order to optimize dosage delivery to the target and minimize damage to critical structures and healthy tissues around the target. 4D ultrasound is a promising modality for motion tracking during such treatments. In this study, the authors evaluate the accuracy of motion tracking in the liver based on deformable registration of 4D ultrasound images. Methods: The offline analysis was performed using a nonrigid registration algorithm that was specifically designed for motion estimation from dynamic imaging data. The method registers the entire 4D image data sequence in a groupwise optimization fashion, thus avoiding a bias toward a specifically chosen reference time point. Three healthy volunteers were scanned over several breathing cycles (12 s) from three different positions and angles on the abdomen; a total of nine 4D scans for the three volunteers. Well-defined anatomic landmarks were manually annotated in all 96 time frames for assessment of the automatic algorithm. The error of the automatic motion estimation method was compared with interobserver variability. The authors also performed experiments to investigate the influence of parameters defining the deformation field flexibility and evaluated how well the method performed with a lower temporal resolution in order to establish the minimum frame rate required for accurate motion estimation. Results: The registration method estimated liver motion with an error of 1 mm (75% percentile over all datasets), which was lower than the interobserver variability of 1.4 mm. The results were only slightly dependent on the degrees of freedom of the deformation model. The registration error increased to 2.8 mm with an eight times lower temporal resolution. Conclusions: The authors conclude that the methodology was able to accurately track the motion of the liver in the 4D ultrasound data. The authors believe that the method has potential in interventions on moving abdominal organs such as MR or ultrasound guided focused ultrasound therapy and radiotherapy, pending the method is enabled to run in real-time. The data and the annotations used for this study are made publicly available for those who would like to test other methods on 4D liver ultrasound data.

60 APPLIED LIFE SCIENCES↗

Differential Credibility of Climate Modes in CMIP6

This work quantitatively evaluates the fidelity with which the northern annular mode (NAM), southern annular mode (SAM), Pacific–North American pattern (PNA), El Niño–Southern Oscillation (ENSO), Pacific decadal oscillation (PDO), Atlantic multidecadal oscillation (AMO), and the first-order mode interactions are represented in Earth system model (ESM) output from the CMIP6 archive. Several skill metrics are used as part of a differential credibility assessment (DCA) of both spatial and temporal characteristics of the modes across ESMs, ESM families, and specific ESM realizations relative to ERA5. The spatial patterns and probability distributions are generally well represented but skill scores that measure the degree to which the frequencies of maximum variance are captured are consistently lower for most ESMs and climate modes. Substantial variability in skill scores manifests across realizations from individual ESMs for the PNA and oceanic modes. Further, the ESMs consistently overestimate the strength of the NAM–PNA first-order interaction and underestimate the NAM–AMO connection. Furthermore, these results suggest that the choice of ESM and ESM realizations will continue to play a critical role in determining climate projections at the global and regional scale at least in the near term.

54 ENVIRONMENTAL SCIENCES↗

Implementation of a feature selection algorithm in FARM to identify important state variables and time-invariant matrices

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM aids the HERON software module in the evaluation of the optimal dispatch for the different IES components. Set-point trajectories are required to meet limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To evaluate the feasibility of HERON generated set-points and to do so in an acceptable time, FARM employs reduced order models to represent the dynamic behavior of the systems to be dispatched. These surrogate models take the form of a linear dynamic system with sets of Linear Parameter Varying (LPV) matrices that are mapped to the system operating space. These matrices are derived from the trajectories of system state variables and system output variables during transients. The accuracy of LPV matrices depends on the selection of state variables. In previous reports, state variables were selected by adopting a complicated workflow requiring multiple software licenses and an advanced level of user expertise. In this report, a new workflow that automates the state variable selection process is presented. It significantly reduces the frequency of user interventions and does not require multiple software licenses. Each module in the new workflow is described in detail, and the input / output examples in each step of the workflow are provided. It was demonstrated that this workflow can greatly reduce the complexity of the state variable selection process, and that the updated FARM-Gamma and FARM-Delta validators can benefit from this workflow when solving the power dispatch problem of a representative IES test case. Finally, some code improvements that can further enhance the efficiency are suggested.

42 ENGINEERING↗

Toward more realistic projections of soil carbon dynamics by Earth system models

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool-and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. Furthermore, we recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.

54 ENVIRONMENTAL SCIENCES↗

Differential Organic Carbon Mineralization Responses to Soil Moisture in Three Different Soil Orders Under Mixed Forested System

Soil microbial respiration is one of the largest sources of carbon (C) emissions to the atmosphere in terrestrial ecosystems, which is strongly dependent on multiple environmental variables including soil moisture. Soil moisture content is strongly dependent on soil texture, and the combined effects of texture and moisture on microbial respiration are complex and less explored. Therefore, this study examines the effects of soil moisture on the mineralization of soil organic C Soil organic carbon in three different soils, Ultisol, Alfisol and Vertisol, collected from mixed forests of Georgia, Missouri, and Texas, United States , respectively. A laboratory microcosm experiment was conducted for 90 days under different moisture regimes. Soil respiration was measured weekly, and destructive harvests were conducted at 1, 15, 60, and 90 days after incubation to determine extractable organic C (EOC), phospholipid fatty acid based microbial community, and C-acquiring hydrolytic extracellular enzyme activities (EEA). The highest cumulative respiration in Ultisol was observed at 50% water holding capacity (WHC), in Alfisol at 100% water holding capacity, and in Vertisol at 175% WHC. The trends in Extractable Organic Carbon were opposite to that of cumulative microbial respiration as the moisture levels showing the highest respiration showed the lowest EOC concentration in all soil types. Also, extracellular enzyme activities increased with increase in soil moisture in all soils, however, respiration and EEA showed a decoupled relationship in Ultisol and Alfisol soils. Soil moisture differences did not influence microbial community composition.

54 ENVIRONMENTAL SCIENCES↗

STRONG UV AND X-RAY VARIABILITY OF THE NARROW LINE SEYFERT 1 GALAXY WPVS 007-ON THE NATURE OF THE X-RAY LOW STATE

We report on multi-wavelength observations of the X-ray transient Narrow Line Seyfert 1 (NLS1) galaxy WPVS 007. The galaxy was monitored with Swift between 2005 October and 2013 July, after it had previously undergone a dramatic drop in its X-ray flux. For the first time, we are able to repeatedly detect this NLS1 in X-rays again. This increased number of detections in the last couple of years may suggest that the strong absorber that has been found in this active galactic nucleus (AGN) is starting to become leaky and may eventually disappear. The X-ray spectra obtained for WPVS 007 are all consistent with a partial covering absorber model. A spectrum based on the data during the extreme low X-ray flux states shows that the absorption column density is of the order of 4 Multiplication-Sign 10{sup 23} cm{sup -2} with a covering fraction of 95%. WPVS 007 also displays one of the strongest UV variabilities seen in NLS1s. The UV continuum variability anti-correlates with the optical/UV slope {alpha}{sub UV}, which suggests that the variability may be primarily due to reddening. The UV variability timescales are consistent with moving dust ''clouds'' located beyond the dust sublimation radius of R{sub sub} Almost-Equal-To 20 lt-days. We present for the first time near-infrared JHK data of WPVS 007, which reveal a rich emission-line spectrum. Recent optical spectroscopy does not indicate significant variability in the broad permitted and Fe II emission lines, implying that the ionizing continuum seen by those gas clouds has not significantly changed over the last decades. All X-ray and UV observations are consistent with a scenario in which an evolving broad absorption line (BAL) flow obscures the continuum emission. As such, WPVS 007 is an important target for our understanding of BAL flows in low-mass AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗

Relaxation approximations to second-order traffic flow models by high-resolution schemes

A relaxation-type approximation of second-order non-equilibrium traffic models, written in conservation or balance law form, is considered. Using the relaxation approximation, the nonlinear equations are transformed to a semi-linear diagonilizable problem with linear characteristic variables and stiff source terms with the attractive feature that neither Riemann solvers nor characteristic decompositions are in need. In particular, it is only necessary to provide the flux and source term functions and an estimate of the characteristic speeds. To discretize the resulting relaxation system, high-resolution reconstructions in space are considered. Emphasis is given on a fifth-order WENO scheme and its performance. The computations reported demonstrate the simplicity and versatility of relaxation schemes as numerical solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

How much carbon can be added to soil by sorption?

Abstract Quantifying the upper limit of stable soil carbon storage is essential for guiding policies to increase soil carbon storage. One pool of carbon considered particularly stable across climate zones and soil types is formed when dissolved organic carbon sorbs to minerals. We quantified, for the first time, the potential of mineral soils to sorb additional dissolved organic carbon (DOC) for six soil orders. We compiled 402 laboratory sorption experiments to estimate the additional DOC sorption potential, that is the potential of excess DOC sorption in addition to the existing background level already sorbed in each soil sample. We estimated this potential using gridded climate and soil geochemical variables within a machine learning model. We find that mid- and low-latitude soils and subsoils have a greater capacity to store DOC by sorption compared to high-latitude soils and topsoils. The global additional DOC sorption potential for six soil orders is estimated to be 107 $$\pm$$ ± 13 Pg C to 1 m depth. If this potential was realized, it would represent a 7% increase in the existing total carbon stock.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using large eddy simulations to reveal the size, strength, and phase of updraft and downdraft cores of an Arctic mixed-phase stratocumulus cloud

Three-dimensional large eddy simulations (LES) are used to analyze a springtime Arctic mixed-phase stratocumulus observed on 26 April 2008 during the Indirect and Semi-Direct Aerosol Campaign. Two subgrid-scale turbulence parameterizations are compared. The first scheme is a 1.5-order turbulent kinetic energy (1.5-TKE) parameterization that has been previously applied to boundary layer cloud simulations. The second scheme, Cloud Layers Unified By Binormals (CLUBB), provides higher-order turbulent closure with scale awareness. The simulations, in comparisons with observations, show that both schemes produce the liquid profiles within measurement variability but underpredict ice water mass and overpredict ice number concentration. The simulation using CLUBB underpredicted liquid water path more than the simulation using the 1.5-TKE scheme, so the turbulent length scale and horizontal grid box size were increased to increase liquid water path and reduce dissipative energy. The LES simulations show this stratocumulus cloud to maintain a closed cellular structure, similar to observations. The updraft and downdraft cores self-organize into a larger meso-γ-scale convective pattern with the 1.5-TKE scheme, but the cores remain more isotropic with the CLUBB scheme. Additionally, the cores are often composed of liquid and ice instead of exclusively containing one or the other. Furthermore, these results provide insight into traditionally unresolved and unmeasurable aspects of an Arctic mixed-phase cloud. From analysis, this cloud's updraft and downdraft cores appear smaller than other closed-cell stratocumulus such as midlatitude stratocumulus and Arctic autumnal mixed-phase stratocumulus due to the weaker downdrafts and lower precipitation rates.

54 ENVIRONMENTAL SCIENCES↗

Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations

In shallow water flow and transport modeling, the monotonic upstream-centered scheme for conservation laws (MUSCL) is widely used to extend the original Godunov scheme to second-order accuracy. The most important step in MUSCL-type schemes is MUSCL reconstruction, which calculate-extrapolates the values of independent variables from the cell center to the edge. The monotonicity of the scheme is preserved with the help of slope limiters that prevent the occurrence of new extrema during reconstruction. On structured grids, the calculation of the slope is straightforward and usually based on a 2-point stencil that uses the cell centers of the neighbor cell and the so-called far-neighbor cell of the edge under consideration. On unstructured grids, the correct choice for the upwind slope becomes nontrivial. 2 novel total variation diminishing schemes are developed based on different techniques for calculating the upwind slope and the downwind slope. An additional treatment that stabilizes the scheme is discussed. The proposed techniques are compared to 2 existing MUSCL reconstruction techniques, and a detailed discussion of the results is given. Finally, it is shown that the proposed MUSCL reconstruction schemes obtain more accurate results with less numerical diffusion and higher efficiency.

54 ENVIRONMENTAL SCIENCES↗

Tracking of Nuclear Cable Insulation Polymer Structural Changes using the Gel Fraction and Uptake Factor Method

Cross-linked polyethylene (XLPE) cable insulation samples were exposed to heat and gamma radiation at a series of temperatures, dose rates, and exposure times to evaluate the effects of these variables on material degradation. The samples were tested using the solvent incubation method to collect gel fraction and uptake factor data in order to assess the crosslinking and chain scission occurring in polymer samples with aging. Consistent with previous reports, gel fraction values were observed to increase and uptake factor values to decrease with radiation and thermal exposure. The trends seen were also more prominent as exposure time increased, suggesting this to be a viable method of tracking structural changes in the XLPE-insulated cable material over extended periods. For the conditions explored, the cable insulation material evaluated did not indicate signs of anomalous aging such as inverse temperature effect in which radiation-induced aging is more severe at lower temperature. Ongoing aging under identical radiation conditions and at lower temperature will further inform conclusions regarding the importance of inverse temperature effects for this material under these conditions.

cable aging↗

Hadron mass corrections in semi-inclusive deep-inelastic scattering

We found that the spin-dependent cross sections for semi-inclusive lepton-nucleon scattering are derived in the framework of collinear factorization, including the effects of masses of the target and produced hadron at finite Q 2 . At leading order the cross sections factorize into products of parton distribution and fragmentation functions evaluated in terms of new, mass-dependent scaling variables. Furthermore, the size of the hadron mass corrections is estimated at kinematics relevant for current and future experiments, and the implications for the extraction of parton distributions from semi-inclusive measurements are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

Comparative estrogenicity of endogenous, environmental and dietary estrogens in pregnant women II: Total estrogenicity calculations accounting for competitive protein and receptor binding and potency

Evaluating the biological significance of human-relevant exposures to environmental estrogens involves assessing the individual and total estrogenicity of endogenous and exogenous estrogens found in serum, for example from biomonitoring studies. We developed a method for this assessment by integrating approaches for (i) measuring total hormone concentrations by mass spectrometry (Fleck et al., 2017), (ii) calculating hormone bioavailable concentrations in serum and, (iii) solving multiple equilibria between estrogenic ligands and receptors, and (iv) quantitatively describing key elements of estrogen potency. The approach was applied to endogenous (E1, E2,E3, E4), environmental (BPA), and dietary Genistein (GEN), Daidzein (DDZ) estrogens measured in the serum of thirty pregnant women. Fractional receptor occupancy (FRO) based estrogenicity was dominated by E1, E2 and E3 (ER-a, 94.4–99.2% (median: 97.3%), ER-ß, 82.7–97.7% (median: 92.8%), as was the total response (TR), which included ligand specific differences in recruitment of co-activator proteins (RCA). The median FRO for BPA was at least fve orders of magnitude lower than E1, E2 and E3, and three orders of magnitude lower than the fetal derived E4 and GEN and DDZ. BPA contributed less than 1/1000th of the normal daily variability in total serum estrogenicity in this cohort of pregnant women.

60 APPLIED LIFE SCIENCES↗

Plasmonic trace sensing below the photon shot noise limit

Plasmonic sensors are important detectors of biochemical trace compounds, but those that utilize optical readout are approaching their absolute limits of detection as defined by the Heisenberg uncertainty principle in both differential intensity and phase readout. However, the use of more general minimum uncertainty states in the form of squeezed light can push the noise floor in these sensors below the shot noise limit (SNL) in one analysis variable at the expense of another. Here, we demonstrate a quantum plasmonic sensor whose noise floor is reduced below the SNL in order to perform index of refraction measurements with sensitivities unobtainable with classical plasmonic sensors. The increased signal-to-noise ratio can result in faster detection of analyte concentrations that were previously lost in the noise. As a result, these benefits are the hallmarks of a sensor exploiting quantum readout fields in order to manipulate the limits of the Heisenberg uncertainty principle.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Evolution of Structure-Reactivity Correlations for the Hydrogen Abstraction Reaction by Chlorine Atom

Empirical structure-reactivity correlations are developed for log k298, the gas-phase rate constants for the reaction (Cl + HCR3 ClH + CR3 ). It has long been recognized that correlation with rH is weak. The poor performance of the linear Evans-Polanyi formulation s illustrated and was little improved by adding a quadratic term, e.g., by making its slope smoothly dependent on rH [ ( rH rHmin) / ( rHmax rHmin)]. The polar effect ( -Cl---H---CR3 +) has also been long discussed but there is no formalization of this dependence based on widely available independent variable(s). Using the sum of Hammett constants for the R substituents also gave at best modest correlations, either for para or for its dissection into F (field/inductive) and R (resonance) effects. Much greater success was achieved by combining these approaches with the preferred independent variable set being either [( rH)2, rH, F, and R] or [ , rH, F, and R]. For 64 rate constants which span 7 orders of magnitude, these correlation formulations give r2 > 0.87 and a mean unsigned deviation of <0.5 log k units, with even better performance if primary secondary, and tertiary reaction centers are treated separately.

chlorine atom↗

Multilevel Monte Carlo Predictions of First Passage Times in Three-Dimensional Discrete Fracture Networks: A Graph-Based Approach

In this paper, we present a method combining multilevel Monte Carlo (MLMC) and a graph-based primary subnetwork identification algorithm to provide estimates of the mean and variance of the distribution of first passage times in fracture media at significantly lower computational cost than standard Monte Carlo (MC) methods. Simulations of solute transport are performed using a discrete fracture network (DFN), and instead of using various grid resolutions for levels in the MLMC, which is standard practice in MLMC, we identify a hierarchy of subnetworks in the DFN based on the shortest topological paths through the network using a graph-based method. While the mean of these ensembles is of critical importance, the variance is also essential in fractured media where uncertainty is an overarching theme, and understanding variability across an ensemble is a requirement for safety assessments. The method provides good estimates of the mean and variance at two orders of magnitude lower computational cost than MC.

54 ENVIRONMENTAL SCIENCES↗