Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONSERVATION”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Path-Integrated Stereo X-Ray Digital Image Correlation: Resolving the Violation of Conservation of Intensity

In this study, x-ray imaging addresses many challenges with visible light imaging in extreme environments, where optical diagnostics such as digital image correlation (DIC) and particle image velocimetry (PIV) suffer biases from index of refraction changes and/or cannot penetrate visibly occluded objects. However, conservation of intensity—the fundamental principle of optical image correlation algorithms—may be violated if ancillary components in the X-ray path besides the surface or fluid of interest move during the test. The test series treated in this work sought to characterize the safe use of fiber-epoxy composites in aerospace and aviation industries during fire accident scenarios. Stereo X-ray DIC was employed to measure test unit deformation in an extreme thermal environment involving a visibly occluded test unit, incident surface heating to temperatures above 600°C, and flames and soot from combusting epoxy decomposition gasses. The objective of the current work is to evaluate two solutions to resolve the violation of conservation of intensity that resulted from both the test unit and the thermal chamber deforming during the test. The first solution recovered conservation of intensity by pre-processing the path-integrated X-ray images to isolate the DIC pattern of the test unit from the thermal chamber components. These images were then correlated with standard, optical DIC software. The second solution, called Path-Integrated Digital Image Correlation (PI-DIC), modified the fundamental matching criterion of DIC to account for multiple, independently-moving components contributing to the final image intensity. The PI-DIC algorithm was extended from a 2D framework to a stereo framework and implemented in a custom DIC software. Both solutions accurately measured the cylindrical shape of the undeformed test unit, recovering radii values of R = 76.20±0.12 mm compared to the theoretical radius of R theor = 76.20 mm. During the test, the test unit bulged asymmetrically as decomposition gasses pressurized the interior and eventually burned in a localized jet. Both solutions measured the heterogeneous radius of this bulge, which reached a maximum of approximately R = 91 mm, with a discrepancy of 2–3% between the two solutions. Two solutions that resolve the violation of conservation of intensity for path-integrated X-ray images were developed. Both were successfully employed in a stereo X-ray DIC configuration to measure deformation of an aluminum-skinned, fiber-epoxy composite test unit in a fire accident scenario.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

High-order, stable, and conservative boundary schemes for central and compact finite differences

Stable and conservative numerical boundary schemes are constructed such that they do not diminish the overall accuracy of the method for interior schemes of orders 4, 6, and 8 using both explicit (central) and compact finite differences. Previous attempts to develop stable numerical boundary schemes for non-linear problems have resulted in schemes which significantly reduced the global accuracy and/or required some form of artificial dissipation. Thus, the schemes developed in this paper are the first to not require this tradeoff, while also ensuring discrete conservation and allowing for direct boundary condition enforcement. After outlining a general procedure for the construction of conservative boundary schemes of any order, a simple, yet novel, optimization strategy which focuses directly on the compressible Euler equations is presented. Furthermore, the result of this non-linear optimization process is a set of high-order, stable, and conservative numerical boundary schemes which demonstrate excellent stability and convergence properties on an array of linear and non-linear hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A Legendre–Fourier spectral method with exact conservation laws for the Vlasov–Poisson system

In this study, we present the design and implementation of an L 2 -stable spectral method for the discretization of the Vlasov–Poisson model of a collisionless plasma in one space and velocity dimension. The velocity and space dependence of the Vlasov equation are resolved through a truncated spectral expansion based on Legendre and Fourier basis functions, respectively. The Poisson equation, which is coupled to the Vlasov equation, is also resolved through a Fourier expansion. The resulting system of ordinary differential equation is discretized by the implicit second-order accurate Crank–Nicolson time discretization. The non-linear dependence between the Vlasov and Poisson equations is iteratively solved at any time cycle by a Jacobian-Free Newton–Krylov method. In this work we analyze the structure of the main conservation laws of the resulting Legendre–Fourier model, e.g., mass, momentum, and energy, and prove that they are exactly satisfied in the semi-discrete and discrete setting. The L 2 -stability of the method is ensured by discretizing the boundary conditions of the distribution function at the boundaries of the velocity domain by a suitable penalty term. The impact of the penalty term on the conservation properties is investigated theoretically and numerically. An implementation of the penalty term that does not affect the conservation of mass, momentum and energy, is also proposed and studied. A collisional term is introduced in the discrete model to control the filamentation effect, but does not affect the conservation properties of the system. Numerical results on a set of standard test problems illustrate the performance of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A fully implicit, conservative, non-linear, electromagnetic hybrid particle-ion/fluid-electron algorithm

The quasi-neutral hybrid model with kinetic ions and fluid electrons is a promising approach for bridging the inherent multi-scale nature of many problems in space and laboratory plasmas. A novel, implicit, particle-in-cell based scheme for the hybrid model is derived for fully 3D electromagnetic problems with multiple ion species, which features global mass, momentum and energy conservation. The scheme includes sub-cycling and orbit-averaging for the ions, with cell-centered finite differences and implicit midpoint time advance. To reduce discrete particle noise, the scheme allows arbitrary-order shape functions for the particle-mesh interpolations and the application of conservative binomial smoothing. The algorithm is verified for a number of test problems to demonstrate the correctness of the implementation, the unique conservation properties, and the favorable stability properties of the new scheme. In particular, there is no indication of unstable growth of the finite-grid instability for a population of cold ions drifting through a uniform spatial mesh, in a set-up where several commonly used non-conservative schemes are highly unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Conservation of the separase regulatory domain

We report a protein sequence analysis of the cell cycle regulatory protease, separase. The sequence and structural conservation of the C-terminal protease domain has long been recognized, whereas the N-terminal regulatory domain of separase was reported to lack detectable sequence similarity. Here we reveal significant sequence conservation of the separase regulatory domain and report a discovery of a cysteine motif (CxCxxC) conserved in major lineages of Metazoa including nematodes and vertebrates. This motif is found in a solvent exposed linker region connecting two TPR-like helical motifs. Mutation of this motif in Caenorhabditis elegans separase leads to a temperature sensitive hypomorphic protein. Conservation of this motif in organisms ranging from C. elegans to humans suggests its functional importance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Machine Learning for Conservative-to-Primitive in Relativistic Hydrodynamics

The numerical solution of relativistic hydrodynamics equations in conservative form requires root-finding algorithms that invert the conservative-to-primitive variables map. These algorithms employ the equation of state of the fluid and can be computationally demanding for applications involving sophisticated microphysics models, such as those required to calculate accurate gravitational wave signals in numerical relativity simulations of binary neutron stars. This work explores the use of machine learning methods to speed up the recovery of primitives in relativistic hydrodynamics. Artificial neural networks are trained to replace either the interpolations of a tabulated equation of state or directly the conservative-to-primitive map. The application of these neural networks to simple benchmark problems shows that both approaches improve over traditional root finders with tabular equation-of-state and multi-dimensional interpolations. In particular, the neural networks for the conservative-to-primitive map accelerate the variable recovery by more than an order of magnitude over standard methods while maintaining accuracy. Neural networks are thus an interesting option to improve the speed and robustness of relativistic hydrodynamics algorithms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A pseudospectral implicit particle-in-cell method with exact energy and charge conservation

The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

What is a particle-conserving Topological Superfluid? The fate of Majorana modes beyond mean-field theory

We investigate Majorana modes of number-conserving fermionic superfluids from both basic physics principles, and concrete models perspectives. After reviewing a criterion for establishing topological superfluidity in interacting systems, based on many-body fermionic parity switches, we reveal the emergence of zero-energy modes anticommuting with fermionic parity. Those many-body Majorana modes are constructed as coherent superpositions of states with different number of fermions. While realization of Majorana modes beyond mean field is plausible, we show that the challenge to quantum-control them is compounded by particle-conservation, and more realistic protocols will have to balance engineering needs with astringent constraints coming from superselection rules. Majorana modes in number-conserving systems are the result of a peculiar interplay between quantum statistics, fermionic parity, and an unusual form of spontaneous symmetry breaking. We test these ideas on the Richardson–Gaudin–Kitaev chain, a number-conserving model solvable by way of the algebraic Bethe ansatz, and equivalent in mean field to a long-range Kitaev chain.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The conserved Wdr8-hMsd1/SSX2IP complex localises to the centrosome and ensures proper spindle length and orientation

The centrosome plays a pivotal role in a wide range of cellular processes and its dysfunction is causally linked to many human diseases including cancer and developmental and neurological disorders. This organelle contains more than one hundred components, and yet many of them remain uncharacterised. Here we identified a novel centrosome protein Wdr8, based upon the structural conservation of the fission yeast counterpart. We showed that Wdr8 constitutively localises to the centrosome and super resolution microscopy uncovered that this protein is enriched at the proximal end of the mother centriole. Furthermore, we identified hMsd1/SSX2IP, a conserved spindle anchoring protein, as one of Wdr8 interactors by mass spectrometry. Wdr8 formed a complex and partially colocalised with hMsd1/SSX2IP. Intriguingly, knockdown of Wdr8 or hMsd1/SSX2IP displayed very similar mitotic defects, in which spindle microtubules became shortened and misoriented. Indeed, Wdr8 depletion resulted in the reduced recruitment of hMsd1/SSX2IP to the mitotic centrosome, though the converse is not true. Together, we propose that the conserved Wdr8-hMsd1/SSX2IP complex plays a critical role in controlling proper spindle length and orientation. - Highlights: • Human Wdr8 is a centrosomal protein enriched in the proximal end of the centriole. • Wdr8 and hMsd1/SSX2IP form a complex conserved in fungi. • Depletion of Wdr8 results in shorter, tilted spindle microtubules. • Depletion phenotypes of Wdr8 are very similar to those of hMsd1/SSX2IP knockdown.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Unexpected role of a conserved domain in the first extracellular loop in G protein-coupled receptor trafficking

Highlights: • Most Class A GPCRs contain a conserved WxFG motif. • Substitution of the W with a non-aromatic amino acid disrupts GPCR functioning. • Disrupted GPCR functioning is the due to aberrant plasma membrane localization. G protein-coupled receptors are the largest superfamily of cell surface receptors in the Metazoa and play critical roles in transducing extracellular signals into intracellular responses. This action is mediated through conformational changes in the receptor following ligand binding. A number of conserved motifs have critical roles in GPCR function, and here we focus on a highly conserved motif (WxFG) in extracellular loop one (EL1). A phylogenetic analysis documents the presence of the WxFG motif in ∼90% of Class A GPCRs and the motif is represented in 17 of the 19 Class A GPCR subfamilies. Using site-directed mutagenesis, we mutagenized the conserved tryptophan residue in eight receptors which are members of disparate class A GPCR subfamilies from different taxa. The modification of the Drosophila leucokinin receptor shows that substitution of any non-aromatic amino acid for the tryptophan leads to a loss of receptor function. Additionally, leucine substitutions at this position caused similar signaling defects in the follicle-stimulating hormone receptor (FSHR), Galanin receptor (GALR1), AKH receptor (AKHR), corazonin receptor (CRZR), and muscarinic acetylcholine receptor (mACHR1). Visualization of modified receptors through the incorporation of a fluorescent tag revealed a severe reduction in plasma membrane expression, indicating aberrant trafficking of these modified receptors. Taken together, these results suggest a novel role for the WxFG motif in GPCR trafficking and receptor function.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A heuristic tool to assess regional impacts of renewable energy infrastructure on conservation areas

Wind and hydropower are important renewable components of national energy portfolios, but their infrastructure negatively affects biodiversity. Regional development requires identification of scenarios that minimize the cumulative impacts of multiple facilities. We introduce the cumulative impact plot (CIP) to quantify cumulative impacts of renewable energy development within a region. Summed impacts of facilities are plotted as a function of the number of facilities, with facilities ranked by increasing (the best-case scenario) or decreasing (worst-case scenario) individual impact. These curves represent lower and upper bounds to which alternative development scenarios (e.g., facilities ranked by generating capacity) can be compared. We used CIPs to assess overlap of potential wind and hydropower facilities with two types of conservation area in the United States: (1) federally protected lands and (2) critical habitats of federally threatened and endangered species. Here, we evaluated two alternative scenarios: facilities ranked by (A) decreasing generating capacity and (B) increasing distance from urban centers. Differences between the best-case and worst-case scenarios were large, thus revealing opportunities to develop facilities with limited impact on conservation areas. Alternative scenarios maximizing energy density generally resulted in conservation area overlap intermediate to best-case and worst-case scenarios. CIPs can also identify the proportion of wind versus hydropower favored under alternative scenarios. Build-out scenarios aimed at minimizing conservation area overlap favor hydropower, whereas alternative scenarios favor wind power. We conclude that CIPs can (1) complement—but should not replace—project-level environmental impact assessment, (2) integrate strategic, cumulative, and scenario-based assessments, and (3) harness rapidly growing geospatial data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Evaluating impacts of development and conservation projects using sustainability indicators: Opportunities and challenges

There has been an increased interest in using sustainability indicators for evaluating the impacts of development and conservation projects. Past and recent experiences have shown that sustainability indicators can be powerful tools for measuring the outcomes of various interventions, when used appropriately and adequately. Currently, there is a range of methods for applying sustainability indicators for project impact evaluation at the environment–development interface. At the same time, a number of challenges persist which have implication for impact evaluation processes especially in developing countries. We highlight some key and recurrent challenges, using three cases from Kenya, Indonesia and Brazil. In this study, we have conducted a comparative analysis across multiple projects from the three countries, which aimed to conserve biodiversity and improve livelihoods. The assessments of these projects were designed to evaluate their positive, negative, short-term, long term, direct and indirect impacts. We have identified a set of commonly used sustainability indicators to evaluate the projects and have discussed opportunities and challenges associated with their application. Our analysis shows that impact evaluation processes present good opportunities for applying sustainability indicators. On the other hand, we find that project proponents (e.g. managers, evaluators, donors/funders) face challenges with establishing full impacts of interventions and that these are rooted in monitoring and evaluation processes, lack of evidence-based impacts, difficulties of measuring certain outcomes and concerns over scale of a range of impacts. We outline key lessons learnt from the multiple cases and propose ways to overcome common problems. Results from our analysis demonstrate practical experiences of applying sustainability indicators in developing countries context where there are different prevailing socio-economic, cultural and environmental conditions. The knowledge derived from this study may therefore be useful to a wider range of audience who are concerned with sustainable integration of development and environmental conservation. - Highlights: • Sustainability indicators are increasingly used for evaluating project impacts. • Lessons learnt are based on case studies from Africa, Asia and South America. • Similar challenges when assessing impacts of development and conservation projects • Need for pragmatic solutions to overcome challenges when assessing project impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Surrogate modeling for efficiently, accurately and conservatively estimating measures of risk

We present a surrogate modeling framework for conservatively estimating measures of risk from limited realizations of an expensive physical experiment or computational simulation. Risk measures combine objective probabilities with the subjective values of a decision maker to quantify anticipated outcomes. Given a set of samples, we construct a surrogate model that produces estimates of risk measures that are always greater than their empirical approximations obtained from the training data. These surrogate models limit over-confidence in reliability and safety assessments and produce estimates of risk measures that converge much faster to the true value than purely sample-based estimates. We first detail the construction of conservative surrogate models that can be tailored to a stakeholder’s risk preferences and then present an approach, based on stochastic orders, for constructing surrogate models that are conservative with respect to families of risk measures. Our surrogate models include biases that permit them to conservatively estimate the target risk measures. We provide theoretical results that show that these biases decay at the same rate as the L 2 error in the surrogate model. Numerical demonstrations confirm that risk-adapted surrogate models do indeed overestimate the target risk measures while converging at the expected rate.

42 ENGINEERING↗