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 379 records · Page 21

DRIPS: A framework for dimension reduction and interpolation in parameter space

Reduced-order models are often used to describe the behavior of complex systems, whose simulation with a full model is too expensive, or to extract salient features from the full model’s output. We introduce a new model-reduction framework DRIPS (dimension reduction and interpolation in parameter space) that combines the offline local model reduction with the online parameter interpolation of reduced-order bases (ROBs). The offline step of this framework relies on dynamic mode decomposition (DMD) to build a low-rank linear surrogate model, equipped with a local ROB, for quantities of interest derived from the training data generated by repeatedly solving the (nonlinear) high-fidelity model for multiple parameter points. The online step consists of the construction of a parametric reduced-order model for each target/test point in the parameter space, with the interpolation of ROBs done on a Grassman manifold and the interpolation of reduced-order operators done on a matrix manifold. The DMD component enables DRIPS to model (typically low-dimensional) quantities of interest directly, without having to access the (typically high-dimensional and possibly nonlinear) operators in a high-fidelity model that governs the dynamics of the underlying high-dimensional state variables, as required in projection-based reduced-order modeling. A series of numerical experiments suggests that DRIPS yields a model reduction, which is computationally more efficient than the commonly used projection-based proper orthogonal decomposition; it does so without requiring a prior knowledge of the governing equation for quantities of interest. Furthermore, for the nonlinear systems considered, DRIPS is more accurate than Gaussian-process interpolation (Kriging).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deep multiscale model learning

The objective of this article is to design novel multi-layer neural networks for multiscale simulations of flows taking into account the observed fine data and physical modeling concepts. Our approaches use deep learning techniques combined with local multiscale model reduction methodologies to predict flow dynamics. Using reduced-order model concepts is important for constructing robust deep learning architectures since the reduced-order models provide fewer degrees of freedom. We consider flow dynamics in porous media as multi-layer networks in this work. More precisely, the solution (e.g., pressures and saturation) at the time instant n+1 depends on the solution at the time instant n and input parameters, such as permeability fields, forcing terms, and initial conditions. One can regard the solution as a multi-layer network, where each layer, in general, is a nonlinear forward map and the number of layers relates to the internal time steps. We will rely on rigorous model reduction concepts to define unknowns and connections between layers. It is critical to use reduced-order models for this purpose, which will identify the regions of influence and the appropriate number of variables. Furthermore, due to the lack of available observed fine data, the reduced-order model can provide us sufficient inexpensive data as needed. The designed deep neural network will be trained using both coarse simulation data which is obtained from the reduced-order model and observed fine data. We will present the main ingredients of our approach and numerical examples. Numerical results show that using deep learning with data generated from multiscale models as well as available observed fine data, we can obtain an improved forward map which can better approximate the fine scale model.

97 MATHEMATICS AND COMPUTING↗

Resistance controllability and variability improvement in a TaO{sub x}-based resistive memory for multilevel storage application

In order to obtain reliable multilevel cell (MLC) characteristics, resistance controllability between the different resistance levels is required especially in resistive random access memory (RRAM), which is prone to resistance variability mainly due to its intrinsic random nature of defect generation and filament formation. In this study, we have thoroughly investigated the multilevel resistance variability in a TaO{sub x}-based nanoscale (<30 nm) RRAM operated in MLC mode. It is found that the resistance variability not only depends on the conductive filament size but also is a strong function of oxygen vacancy concentration in it. Based on the gained insights through experimental observations and simulation, it is suggested that forming thinner but denser conductive filament may greatly improve the temporal resistance variability even at low operation current despite the inherent stochastic nature of resistance switching process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Monitoring the variability of intrinsic absorption lines in quasar spectra , ,

We have monitored 12 intrinsic narrow absorption lines (NALs) in five quasars and seven mini-broad absorption lines (mini-BALs) in six quasars for a period of 4-12 yr (1-3.5 yr in the quasar rest-frame). We present the observational data and the conclusions that follow immediately from them, as a prelude to a more detailed analysis. We found clear variability in the equivalent widths (EWs) of the mini-BAL systems but no easily discernible changes in their profiles. We did not detect any variability in the NAL systems or in narrow components that are often located at the center of mini-BAL profiles. Variations in mini-BAL EWs are larger at longer time intervals, reminiscent of the trend seen in variable BALs. If we assume that the observed variations result from changes in the ionization state of the mini-BAL gas, we infer lower limits to the gas density ∼10{sup 3}-10{sup 5} cm{sup –3} and upper limits on the distance of the absorbers from the central engine of the order of a few kiloparsecs. Motivated by the observed variability properties, we suggest that mini-BALs can vary because of fluctuations of the ionizing continuum or changes in partial coverage while NALs can vary primarily because of changes in partial coverage.

79 ASTRONOMY AND ASTROPHYSICS↗

Natural variability and antioxidant properties of commercially cultivated switchgrass extractives

Optimum utilization of lignocellulosic bioenergy feedstock, such as switchgrass, via recovery of value-added phytochemicals and reduction of waste streams will increase the economic feasibility as well as the sustainability of integrated biorefineries. Switchgrass contains significant amount of overlooked non-structural components, i.e. extractives, which encompasses phenolic derivatives with potential antimicrobial and antioxidant properties. As a perennial crop, it is essential to assess the natural variability of switchgrass extractives content and composition in order to determine the feasibility of its extraction and utilization. In this study, three switchgrass cultivars were harvested from four commercial farms in their second and third growing seasons. Their respective extractive fractions were characterized for the presence of free sugars, total phenolics (TPC), hydroxycinnamic acid (HCAs) derivatives, total flavonols, and other important phytochemicals. Biomass harvested in the third growing season displayed greater amounts of TPC, HCAs, and total flavonols. Environmental conditions such as total precipitation, and cultural practices including harvest time, played important roles in the accumulation of TPC. As a result, the switchgrass extractives also exhibited DPPH (2,2-diphenyl-1-picrylhydrazyl) radical scavenging activity with a positive correlation between inhibiting concentration (IC 50 ) value and total flavonols content.

59 BASIC BIOLOGICAL SCIENCES↗

An artificial neural network framework for reduced order modeling of transient flows

Here, this paper proposes a supervised machine learning framework for the non-intrusive model order reduction of unsteady fluid flows to provide accurate predictions of non-stationary state variables when the control parameter values vary. Our approach utilizes a training process from full-order scale direct numerical simulation data projected on proper orthogonal decomposition (POD) modes to achieve an artificial neural network (ANN) model with reduced memory requirements. This data-driven ANN framework allows for a nonlinear time evolution of the modal coefficients without performing a Galerkin projection. Our POD-ANN framework can thus be considered an equation-free approach for latent space dynamics evolution of nonlinear transient systems and can be applied to a wide range of physical and engineering applications. Within this framework we introduce two architectures, namely sequential network (SN) and residual network (RN), to train the trajectory of modal coefficients. We perform a systematic analysis of the performance of the proposed reduced order modeling approaches on prediction of a nonlinear wave-propagation problem governed by the viscous Burgers equation, a simplified prototype setting for transient flows. We find that the POD-ANN-RN yields stable and accurate results for test problems assessed both within inside and outside of the database range and performs significantly better than the standard intrusive Galerkin projection model. Finally, our results show that the proposed framework provides a non-intrusive alternative to the evolution of transient physics in a POD basis spanned space, and can be used as a robust predictive model order reduction tool for nonlinear dynamical systems.

42 ENGINEERING↗

Radiative effects of interannually varying vs. interannually invariant aerosol emissions from fires

Abstract. Open-burning fires play an important role in the earth's climate system. In addition to contributing a substantial fraction of global emissions of carbon dioxide, they are a major source of atmospheric aerosols containing organic carbon, black carbon, and sulfate. These “fire aerosols” can influence the climate via direct and indirect radiative effects. In this study, we investigate these radiative effects and the hydrological fast response using the Community Atmosphere Model version 5 (CAM5). Emissions of fire aerosols exert a global mean net radiative effect of −1.0 W m−2, dominated by the cloud shortwave response to organic carbon aerosol. The net radiative effect is particularly strong over boreal regions. Conventionally, many climate modelling studies have used an interannually invariant monthly climatology of emissions of fire aerosols. However, by comparing simulations using interannually varying emissions vs. interannually invariant emissions, we find that ignoring the interannual variability of the emissions can lead to systematic overestimation of the strength of the net radiative effect of the fire aerosols. Globally, the overestimation is +23 % (−0.2 W m−2). Regionally, the overestimation can be substantially larger. For example, over Australia and New Zealand the overestimation is +58 % (−1.2 W m−2), while over Boreal Asia the overestimation is +43 % (−1.9 W m−2). The systematic overestimation of the net radiative effect of the fire aerosols is likely due to the non-linear influence of aerosols on clouds. However, ignoring interannual variability in the emissions does not appear to significantly impact the hydrological fast response. In order to improve understanding of the climate system, we need to take into account the interannual variability of aerosol emissions.

54 ENVIRONMENTAL SCIENCES↗

Effectiveness of the food recovery at the retailing stage under shelf life uncertainty: An application to Italian food chains

Highlights: • The food recovery is seen as suitable way to manage food near to its expiry date. • The variability of the products shelf life must be taken into account. • The paper addresses the mathematic modeling of the profit related to food recovery. • The optimal time to withdraw the products is determinant for food recovery. - Abstract: Food losses represent a significant issue affecting food supply chains. The possibility of recovering such products can be seen as an effective way to reduce such a phenomenon, improve supply chain performances and ameliorate the conditions of undernourished people. The topic has been already investigated by a previous paper enforcing the hypothesis of deterministic and constant Shelf Life (SL) of products. However, such a model cannot be properly extended to products affected by uncertainties of the SL as it does not take into account the deterioration costs and loss of profits due to the overcoming of the SL within the cycle time. Thus the present paper presents an extension of the previous one under stochastic conditions of the food quality. Differently from the previous publication, this work represents a general model applicable to all supply chains, especially to those managing fresh products characterized by uncertain SL such as fruits and vegetables. The deterioration costs and loss of profits are included in the model and the optimal time at which to withdraw the products from the shelves as well as the quantities to be shipped at each alternative destination have been determined. A comparison of the proposed model with that reported in the previous publication has been carried out in order to underline the impact of the SL variability on the optimality conditions. The results show that the food recovery strategy in the presence of uncertainty of the food quality is rewarding, even if the optimal profit is lower than that of the deterministic case.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

VariableR Reclustering in Multiple Top Quark and W Boson Events

VariableR jet reclustering is an innovative technique that allows for the reconstruction of boosted object over a wide range of kinematic regimes. Such capability enables the efficient identification of events with multiple boosted top quarks which is a typical signature for new physics processes such as the production of the supersymmetric partner of the gluon. In order to evaluate the performance of the algorithm, the VariableR reclustered jets are compared with fixed radius reclustered jets. The flexibility of the algorithm is tested by reconstructing both boosted top quarks and boosted W bosons. The VariableR reclustering method is found to be more efficient than the fixed radius algorithm at identifying top quarks and W bosons in events with four top quarks, therefore enhancing the sensitivity for gluino searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

VariableR reclustering in multiple top quark events - Oral Presentation

VariableR jet reclustering is an innovative technique that allows for the reconstruction of boosted object over a wide range of kinematic regimes. Such capability enables the efficient identification of events with multiple boosted top quarks which is a typical signature for new physics processes such as the production of the supersymmetric partner of the gluon. In order to evaluate the performance of the algorithm, the VariableR reclustered jets are compared with fixed radius reclustered jets. The flexibility of the algorithm is tested by reconstructing both boosted top quarks and boosted W bosons. The VariableR reclustering method is found to be more efficient than the fixed radius algorithm at identifying top quarks and W bosons in events with four top quarks, therefore enhancing the sensitivity for gluino searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-E-J-113: Effects of Deformable Registration On First-Order Texture Maps Calculated From Thoracic Lung CT Scans

Purpose: To determine the stability of eight first-order texture features following the deformable registration of serial computed tomography (CT) scans. Methods: CT scans at two different time points from 10 patients deemed to have no lung abnormalities by a radiologist were collected. Following lung segmentation using an in-house program, texture maps were calculated from 32×32-pixel regions of interest centered at every pixel in the lungs. The texture feature value of the ROI was assigned to the center pixel of the ROI in the corresponding location of the texture map. Pixels in the square ROI not contained within the segmented lung were not included in the calculation. To quantify the agreement between ROI texture features in corresponding pixels of the baseline and follow-up texture maps, the Fraunhofer MEVIS EMPIRE10 deformable registration algorithm was used to register the baseline and follow-up scans. Bland-Altman analysis was used to compare registered scan pairs by computing normalized bias (nBias), defined as the feature value change normalized to the mean feature value, and normalized range of agreement (nRoA), defined as the range spanned by the 95% limits of agreement normalized to the mean feature value. Results: Each patient’s scans contained between 6.8–15.4 million ROIs. All of the first-order features investigated were found to have an nBias value less than 0.04% and an nRoA less than 19%, indicating that the variability introduced by deformable registration was low. Conclusion: The eight first-order features investigated were found to be registration stable. Changes in CT texture maps could allow for temporal-spatial evaluation of the evolution of lung abnormalities relating to a variety of diseases on a patient-by-patient basis. SGA and HA receives royalties and licensing fees through the University of Chicago for computer-aided diagnosis technology. Research reported in this publication was supported by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number R25GM109439.

60 APPLIED LIFE SCIENCES↗

Improved accuracy in degenerate variational integrators for guiding centre and magnetic field line flow

First-order-accurate degenerate variational integration (DVI) was introduced in Ellison et al. ( Phys. Plasmas , vol. 25, 2018, 052502) for systems with a degenerate Lagrangian, i.e. one in which the velocity-space Hessian is singular. In this paper we introduce second-order-accurate DVI schemes, both with and without non-uniform time stepping. We show that it is not in general possible to construct a second-order scheme with a preserved two-form by composing a first-order scheme with its adjoint, and discuss the conditions under which such a composition is possible. We build two classes of second-order-accurate DVI schemes. We test these second-order schemes numerically on two systems having non-canonical variables, namely the magnetic field line and guiding centre systems. Variational integration for Hamiltonian systems with non-uniform time steps, in terms of an extended phase space Hamiltonian, is generalized to non-canonical variables. It is shown that preservation of proper degeneracy leads to single-step (one-step) methods without parasitic modes, i.e. to non-uniform time step DVIs. This extension applies to second-order-accurate as well as first-order schemes, and can be applied to adapt the time stepping to an error estimate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Frankenstein’s ROMster: Avoiding pitfalls of reduced-order model development

Reduced-order models (ROMs) are a widely used and powerful approach to reducing the complexity of predictive physics-based numerical simulations for a wide range of applications, including electronics and fluid mechanics such as geologic CO 2 sequestration (GCS). ROMs are critical for optimization, sensitivity analysis, model calibration and uncertainty quantification where full-order models cannot be feasibly executed many times. Traditional approaches generate a single ROM for each simulated response (e.g., CO 2 injection rates, pH changes) based on a set of training simulations. Here, we demonstrate that a single ROM can display excellent overall predictive statistics, but have predictions that dramatically and unacceptably deviate from simulator responses especially when the response variable has a large range (i.e., vary over multiple orders of magnitude). For example, we show that a traditional statistically-high-performing GCS ROM (coefficient of determination R 2 of 0.99) can have average absolute relative errors of over 200%. To address this, we propose a new and novel approach where a set of sub-ROMs are generated to overcome the potential pitfalls in traditional single ROM development. The effectiveness of the proposed approach—the ROMster framework—is demonstrated using a case study of predicted CO 2 injection rates for GCS. We find our approach is a robust and general framework for ROM development, reducing the average “error” from 200% to only 4% in our case study.

54 ENVIRONMENTAL SCIENCES↗

THE BERLIN EXOPLANET SEARCH TELESCOPE II CATALOG OF VARIABLE STARS. I. CHARACTERIZATION OF THREE SOUTHERN TARGET FIELDS

A photometric survey of three southern target fields with BEST II yielded the detection of 2406 previously unknown variable stars and an additional 617 stars with suspected variability. This study presents a catalog including their coordinates, magnitudes, light curves, ephemerides, amplitudes, and type of variability. In addition, the variability of 17 known objects is confirmed, thus validating the results. The catalog contains a number of known and new variables that are of interest for further astrophysical investigations, in order to, e.g., search for additional bodies in eclipsing binary systems, or to test stellar interior models. Altogether, 209,070 stars were monitored with BEST II during a total of 128 nights in 2009/2010. The overall variability fraction of 1.2%-1.5% in these target fields is well comparable to similar ground-based photometric surveys. Within the main magnitude range of R in [11, 17], we identify 0.67(3)% of all stars to be eclipsing binaries, which indicates a completeness of about one third for this particular type in comparison to space surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving the accuracy of discretisations of the vector transport equation on the lowest-order quadrilateral Raviart-Thomas finite elements

Within finite element models of fluids, vector-valued fields such as velocity or momentum variables are commonly discretised using the Raviart-Thomas elements. However, when using the lowest-order quadrilateral Raviart-Thomas elements, standard finite element discretisations of the vector transport equation typically have a low order of spatial accuracy. This paper describes two schemes that improve the accuracy of transporting such vector-valued fields on two-dimensional curved manifolds. The first scheme that is presented reconstructs the transported field in a higher-order function space, where the transport equation is then solved. The second scheme applies a mixed finite element formulation to the vector transport equation, simultaneously solving for the transported field and its vorticity. In this work, an approach to stabilising this mixed vector-vorticity formulation is presented that uses a Streamline Upwind Petrov-Galerkin (SUPG) method. These schemes are then demonstrated, along with their accuracy properties, through some numerical tests. Two new test cases are used to assess the transport of vector-valued fields on curved manifolds, solving the vector transport equation in isolation. The improvement of the schemes is also shown through two standard test cases for rotating shallow-water models.

97 MATHEMATICS AND COMPUTING↗

The arbitrary order mixed mimetic finite difference method for the diffusion equation

Here, we propose an arbitrary-order accurate mimetic finite difference (MFD) method for the approximation of diffusion problems in mixed form on unstructured polygonal and polyhedral meshes. As usual in the mimetic numerical technology, the method satisfies local consistency and stability conditions, which determines the accuracy and the well-posedness of the resulting approximation. The method also requires the definition of a high-order discrete divergence operator that is the discrete analog of the divergence operator and is acting on the degrees of freedom. The new family of mimetic methods is proved theoretically to be convergent and optimal error estimates for flux and scalar variable are derived from the convergence analysis. A numerical experiment confirms the high-order accuracy of the method in solving diffusion problems with variable diffusion tensor. It is worth mentioning that the approximation of the scalar variable presents a superconvergence effect.

97 MATHEMATICS AND COMPUTING↗

Towards an improved understanding of the influence of subpixel vegetation structure on pixel-level spectra: a simulation approach

The planned NASA Hyperspectral Infrared Imager (HyspIRI) mission, equipped with an imaging spectrometer that has the capability of monitoring ecosystems globally, will provide an unprecedented opportunity to address scientific challenges related to ecosystem function and change. However, uncertainty remains around the impact of subpixel vegetation structure, in combination with the point spread function, on pixel-level imaging spectroscopy data. We estimated structural parameters, e.g., leaf area index (LAI), canopy cover, and tree location, from HyspIRI spectral data, with the goal of assessing how subpixel variation in these parameters impact pixel-level imaging spectroscopy data. The fine-scale variability of real vegetation structure makes this a challenging endeavor. Therefore, we utilized a simulation-based approach to counter the time-consuming and often destructive sampling needs of vegetation structural analysis and to simultaneously generate synthetic HyspIRI data pre-launch. Three virtual scenes were constructed, corresponding to the actual vegetation structure in the National Ecological Observatory Network’s (NEON) Pacific Southwest Domain (Fresno, CA). These included an oak savanna, a dense coniferous forest, and a conifer-manzanita-mixed forest. Simulated spectroscopy data for these scenes were then generated using the Digital Imaging and Remote Sensing Image Generation (DIRSIG) model. Simulations first were used to verify the physical model, virtual scene geometrical information, and simulation parameters. This was followed by simulations of HyspIRI data, where within-pixel structural variability was introduced, e.g., by iteratively changing per-pixel canopy cover and tree placement, tree clustering, leaf area index (LAI), etc., between simulation runs for the virtual scenes. Finally, narrow-band vegetation indices (VIs) were extracted from the data in an attempt to describe the variability of the subpixel structural parameters; this was done in order to assess VI robustness to changes in structural “levels”, as well as placement of trees/canopies within the instrument’s instantaneous field-of-view (IFOV). Our ultimate goal is not only to better understand how such subpixel variability influence imaging spectroscopy outputs, but also to better estimate vegetation structural parameters using spectra. We constructed regression models for LAI (R2 = 0.92) and canopy cover (R2 = 0.97) with narrow-band VIs via this simulation approach. Our models ultimately are intended to improve the HyspIRI mission’s ability to monitor global vegetation structure.

Van Aardt, Jan↗

A robust and efficient finite volume method for compressible inviscid and viscous two-phase flows [A robust and efficient finite volume method for compressible viscous two-phase flows]

A robust and efficient density-based finite volume method is developed for solving the six-equation single pressure system of two-phase flows at all speeds on hybrid unstructured grids. Unlike conventional approaches where an expensive exact Riemann solver is normally required for computing numerical fluxes at the two-phase interfaces in addition to AUSM-type fluxes for single-phase interfaces in order to maintain stability and robustness in cases involving interactions of strong pressure and void-fraction discontinuities, a volume-fraction coupling term for the AUSM + -up fluxes is introduced in this work to impart the required robustness without the need of the exact Riemann solver. The resulting method is significantly less expensive in regions where otherwise the Riemann solver would be invoked. A transformation from conservative variables to primitive variables is presented and the primitive variables are then solved in the implicit method in order for the current finite volume method to be able to solve, effectively and efficiently, low Mach number flows in traditional multiphase applications, which otherwise is a great challenge for the standard density-based algorithms. Here, a number of benchmark test cases are presented to assess the performance and robustness of the developed finite volume method for both inviscid and viscous two-phase flow problems. The numerical results indicate that the current density-based method provides an attractive and viable alternative to its pressure-based counterpart for compressible two-phase flows at all speeds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗