Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservative variables”

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 163 records · Page 9

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗

A New Approach to Evaluate and Reduce Uncertainty of Model-Based Biodiversity Projections for Conservation Policy Formulation

Biodiversity projections with uncertainty estimates under different climate, land-use, and policy scenarios are essential to setting and achieving international targets to mitigate biodiversity loss. Evaluating and improving biodiversity predictions to better inform policy decisions remains a central conservation goal and challenge. A comprehensive strategy to evaluate and reduce uncertainty of model outputs against observed measurements and multiple models would help to produce more robust biodiversity predictions. We propose an approach that integrates biodiversity models and emerging remote sensing and in-situ data streams to evaluate and reduce uncertainty with the goal of improving policy-relevant biodiversity predictions. In this work, we describe a multivariate approach to directly and indirectly evaluate and constrain model uncertainty, demonstrate a proof of concept of this approach, embed the concept within the broader context of model evaluation and scenario analysis for conservation policy, and highlight lessons from other modeling communities.

essential biodiversity variables↗

Spatio-Temporal Distribution of Particulate and Dissolved Organic Matter in the Mississippi River Bight From Optical Measurements

The Mississippi River Bight is a highly dynamic region influenced by the seasonally variable outflow from the Mississippi River. In an effort to characterize the distribution of particulate and dissolved organic matter in the region, we conducted a two-year field program in the spring and fall (high and low flow river discharge) of 2000 and 2002. We collected a comprehensive set of bio-optical measurements consisting of vertical profiles (absorption, scattering, chlorophyll fluorescence and radiometry) and discrete measurements (pigment concentrations, particulate and CDOM absorption) that enabled us to obtain better insight into the seasonal and spatial variability of some important biogeochemical parameters. Our field measurements generally showed higher phytoplankton clorophyll concentrations in the plume waters (associated with lower surface salinities) and confirmed the high biological activity abserved in other studies. The seasonal flow of river discharge and advective currents due to wind forcing exerted a strong influence on the biological and optical properties of the region. An examination of absorption at 440 nm by the algal and non-algal fraction of the particulate pool and of CDOM revealed that at nearshore stations, contributions by the non-algal particles were high (about 40%) and decresed with increasing salinities. While CDOM absorption exhibited conservative mixing, its relative contribution to the total absorption was variable. Surface waters at most stations had lower salinities that generalliy increased with dept. Particulate matter and CDOM also decreased with depth as evidenced by absorption and scattering measurements. Good correlations in surface waters between concentrations of particulate and dissolved matter, the inherent optical properties of absorption and ackscattering and remote sensing reflectance values has allowed the development of robust empirical algorithms for phytoplankton chlorophyll and CDOM absorption.

DSa, E. J.↗

Variable Refrigerant Flow Performance: Humidity and Efficiency in Multifamily Buildings in New York, NY

Variable refrigerant flow (VRF) systems have become common in large multifamily buildings. The upfront cost is appealing to developers, as a single system can provide all heating and cooling. VRF systems also use no on-site fossil fuels which can help avoid greenhouse gas emissions. However, there have been very few data on how VRF systems actually perform in buildings. Researchers monitored three new, multifamily buildings with a focus on VRF system sizing, humidity control, and energy efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dendritic cells focus CTL responses toward highly conserved and topologically important HIV-1 epitopes

During early HIV-1 infection, immunodominant T cell responses to highly variable epitopes lead to the establishment of immune escape virus variants. Here we assessed a type 1-polarized monocyte-derived dendritic cell (MDC1)-based approach to selectively elicit cytotoxic T lymphocyte (CTL) responses against highly conserved and topologically important HIV-1 epitopes in HIV-1-infected individuals from the Thailand RV254/SEARCH 010 cohort who initiated antiretroviral therapy (ART) during early infection (Fiebig stages I-IV). Autologous MDC1 were used as antigen presenting cells to induce in vitro CTL responses against HIV-1 Gag, Pol, Env, and Nef as determined by flow cytometry and ELISpot assay. Ultra-conserved or topologically important antigens were respectively identified using the Epigraph tool and a structure-based network analysis approach and compared to overlapping peptides spanning the Gag proteome. MDC1 presenting either the overlapping Gag, Epigraph, or Network 14-21mer peptide pools consistently activated and expanded HIV-1-specific T cells to epitopes identified at the 9-13mer peptide level. Interestingly, some CTL responses occurred outside known or expected HLA associations, providing evidence of new HLA-associated CTL epitopes. Comparative analyses demonstrated more sequence conservation among Epigraph antigens but a higher magnitude of CTL responses to Network and Gag peptide groups. Importantly, CTL responses against topologically constrained Gag epitopes contained in both the Network and Gag peptide pools were selectively enhanced in the Network pool-initiated cultures. Our study supports the use of MDC1 as a therapeutic strategy to induce and focus CTL responses toward putative fitness-constrained regions of HIV-1 to prevent immune escape and control HIV-1 infection.

60 APPLIED LIFE SCIENCES↗

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗

An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables

We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.

97 MATHEMATICS AND COMPUTING↗

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗

A New Approach to Extreme Value Estimation Applicable to a Wide Variety of Random Variables

Designing reliable structures requires an estimate of the maximum and minimum values (i.e., strength and load) that may be encountered in service. Yet designs based on very extreme values (to insure safety) can result in extra material usage and hence, uneconomic systems. In aerospace applications, severe over-design cannot be tolerated making it almost mandatory to design closer to the assumed limits of the design random variables. The issue then is predicting extreme values that are practical, i.e. neither too conservative or non-conservative. Obtaining design values by employing safety factors is well known to often result in overly conservative designs and. Safety factor values have historically been selected rather arbitrarily, often lacking a sound rational basis. To answer the question of how safe a design needs to be has lead design theorists to probabilistic and statistical methods. The so-called three-sigma approach is one such method and has been described as the first step in utilizing information about the data dispersion. However, this method is based on the assumption that the random variable is dispersed symmetrically about the mean and is essentially limited to normally distributed random variables. Use of this method can therefore result in unsafe or overly conservative design allowables if the common assumption of normality is incorrect.

Holland, Frederic A., Jr.↗

End-Use Savings Shapes Measure Documentation: Variable Refrigerant Flow With 25% Upsizing Allowance for Heating

This documentation focuses on a single end-use savings shape upgrade - a variable refrigerant flow with heat recovery (VRF HR) heating and cooling system coupled with a dedicated outdoor air system (DOAS) for ventilation and where 25% upsizing (or 125% of the original size) is allowed for heating in colder climates (i.e., heating dominant regions). This document will primarily discuss the additional changes to the sizing algorithm and modeling approach, while a comprehensive overview of the fundamental modeling methodology and background of the VRF HR DOAS upgrade, including applicability and other key assumptions can be found in the original documentation: Variable Refrigerant Flow with Heat Recovery and Dedicated Outdoor Air System. To provide more context on the 25% upsizing algorithm, if the building is a cooling dominant (e.g., design cooling load higher than design heating load ), then the outdoor unit capacity of the VRF heat pump is sized based on the design cooling load. However, if the building is heating dominant, then this building is a candidate for the 25% upsizing allowance. Once the 25% upsizing is allowed, if the 25% upsized capacity (or 125% from the original size) represented with design condition exceeds the design heating load, then the design heating load is used for calculating the rated capacity of the outdoor unit. And if the 25% upsized capacity represented with the design condition does not exceed the design heating load, then the 25% upsized capacity represented with rated condition is used for the capacity of the outdoor unit while the remaining heating load is handled with the supplemental/backup electric resistance coil.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Population structure limits the use of genomic data for predicting phenotypes and managing genetic resources in forest trees

There is overwhelming evidence that forest trees are locally adapted to climate. Thus, genecological models based on population phenotypes have been used to measure local adaptation, infer genetic maladaptation to climate, and guide assisted migration. However, instead of phenotypes, there is increasing interest in using genomic data for gene resource management. We used whole-genome resequencing and common-garden experiments to understand the genetic architecture of adaptive traits in black cottonwood. We studied the potential of using genome-wide association studies (GWAS) and genomic prediction to detect causal loci, identify climate-adapted phenotypes, and inform gene resource management. We analyzed population structure by partitioning phenotypic and genomic (single-nucleotide polymorphism) variation among 840 genotypes collected from 91 stands along 16 rivers. Most phenotypic variation (60 to 81%) occurred among populations and was strongly associated with climate. Population phenotypes were predicted well using genomic data (e.g., predictive abilityr> 0.9) but almost as well using climate or geography (r> 0.8). In contrast, genomic prediction within populations was poor (r< 0.2). We identified many GWAS associations among populations, but most appeared to be spurious based on pooled within-population analyses. Hierarchical partitioning of linkage disequilibrium and haplotype sharing suggested that within-population genomic prediction and GWAS were poor because allele frequencies of causal loci and linked markers differed among populations. Given the urgent need to conserve natural populations and ecosystems, our results suggest that climate variables alone can be used to predict population phenotypes, delineate seed zones and deployment zones, and guide assisted migration.

Science & Technology - Other Topics↗

A 3-D upwind Euler solver for unstructured meshes

A three-dimensional finite-volume upwind Euler solver is developed for unstructured meshes. The finite-volume scheme solves for solution variables at vertices of the mesh and satisfies the integral conservation law on nonoverlapping polyhedral control volumes surrounding vertices of the mesh. The schene achieves improved solution accuracy by assuming a piecewise linear variation of the solution in each control volume. This improved spatial accuracy hinges heavily upon the calculation of the solution gradient in each control volume given pointwise values of the solution at vertices of the mesh. Several algorithms are discussed for obtaining these gradients. Details concerning implementation procedures and data structures are discussed. Sample calculations for inviscid Euler flow about isolated aircraft wings at subsonic and transonic speeds are compared with established Euler solvers as well as experiment.

Barth, Timothy J.↗

Combustion of liquid-fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both subcritical and supercritical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates variable thermophysical properties, finite-rate chemical kinetics, and a full treatment of liquid-vapor phase equilibrium at the drop surface. The governing equations and associated interfacial boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures in the range of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the critical pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure.

Shuen, J. S.↗

Kinetic simulation of plasma flows in the inner magnetosphere

A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.

Miller, Ronald H.↗

Application of the Space-Time Conservation Element and Solution Element Method to One-Dimensional Convection-Diffusion Problems

In the space-time conservation element and solution element (CE/SE) method, the independent marching variables used comprise not only the mesh value of the physical dependent variables but also, in contrast to it typical numerical method, the Mesh values of the spatial derivatives of the physical variables The use of the extra marching variables results from the need to construct the two-level explicit and nondissipative schemes which are at the core of the CE/SE development. It also results from the need to minimize the stencil while maintaining accuracy. In this paper using the 1D(sub (alpha)-mu) scheme as an example, the effect of this added complication on consistency, accuracy and operation count is assessed. As part of this effort, an equivalent yet more efficient form of the alpha-mu scheme in which the independent marching variables are the local fluxes tied to each mesh point is introduced. Also, the intriguing relations that exist among the alpha-mu. Leapfrog, and DuFort-Frankel schemes are further explored. In addition, the redundance of the Leapfrog, DUFort-Frankel, and Lax scheme and the remedy for this redundance are discussed. This paper is concluded with the construction and evaluation of a CE/SE solver for the inviscid Burger equation.

Chang, Sin-Chung↗

Multiplicative noise and the diffusion of conserved densities

Stochastic fluid dynamics governs the long time tails of hydrodynamic correlation functions, and the critical slowing down of relaxation phenomena in the vicinity of a critical point in the phase diagram. In this work we study the role of multiplicative noise in stochastic fluid dynamics. Multiplicative noise arises from the dependence of transport coefficients, such as the diffusion constants for charge and momentum, on fluctuating hydrodynamic variables. We study long time tails and relaxation in the diffusion of a conserved density (model B), and a conserved density coupled to the transverse momentum density (model H). Careful attention is paid to fluctuation-dissipation relations. We observe that multiplicative noise contributes at the same order as non-linear interactions in model B, but is a higher order correction to the relaxation of a scalar density and the tail of the stress tensor correlation function in model H.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new numerical framework for solving conservation laws: The method of space-time conservation element and solution element

A new numerical framework for solving conservation laws is being developed. It employs: (1) a nontraditional formulation of the conservation laws in which space and time are treated on the same footing, and (2) a nontraditional use of discrete variables such as numerical marching can be carried out by using a set of relations that represents both local and global flux conservation.

Chang, Sin-Chung↗

Africa Food Security & Agriculture: Predicting the Likelihood of Human-elephant Conflict and Assessing Elephant Habitat Conditions During Extreme Drought and Crop Deficit in the Kavango-Zambezi Area

Human-wildlife conflict is increasingly more common due to human population growth, habitat fragmentation, and changing climatic conditions. This conflict is particularly evident in the Kavango-Zambezi area, where over three million people share the landscape with an abundant megafauna population. As the changing climate continues to exacerbate drought severity and subsequent food availability, conflict between humans and wildlife has become more prevalent. In the Kavango-Zambezi area, conflict between elephants and humans has resulted in crop loss, property damage, and threats to public safety. In order to manage current and future conflict, The Ecoexist Project and Connected Conservation have been working to empower farmers and conserve the natural habitat. This DEVELOP project employed Earth observations to conduct a time series analysis of vegetation health change, elephant movement, and climate conditions, from 2017 to 2020. Data were aggregated into the wet (November through April) and dry (May through October) seasons. The resulting analysis demonstrated the potential to use Landsat 8 Operational Land Imager (OLI), Global Precipitation Measurement’s Integrated Multi-satellite Retrievals for GPM (GPM IMERG), and TerraClimate data to identify potential areas of conflict under increased seasonal variability. An improved understanding of conflict drivers will help support sustainable wildlife conservation and food security in the future.

DEVELOP Project Summary↗