Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
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 OSTI2013

A discontinuous Galerkin conservative level set scheme for interface capturing in multiphase flows

The accurate conservative level set (ACLS) method of Desjardins et al. [O. Desjardins, V. Moureau, H. Pitsch, An accurate conservative level set/ghost fluid method for simulating turbulent atomization, J. Comput. Phys. 227 (18) (2008) 8395–8416] is extended by using a discontinuous Galerkin (DG) discretization. DG allows for the scheme to have an arbitrarily high order of accuracy with the smallest possible computational stencil resulting in an accurate method with good parallel scaling. This work includes a DG implementation of the level set transport equation, which moves the level set with the flow field velocity, and a DG implementation of the reinitialization equation, which is used to maintain the shape of the level set profile to promote good mass conservation. A near second order converging interface curvature is obtained by following a height function methodology (common amongst volume of fluid schemes) in the context of the conservative level set. Various numerical experiments are conducted to test the properties of the method and show excellent results, even on coarse meshes. The tests include Zalesak’s disk, two-dimensional deformation of a circle, time evolution of a standing wave, and a study of the Kelvin–Helmholtz instability. Finally, this novel methodology is employed to simulate the break-up of a turbulent liquid jet.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A consistent and conservative scheme for MHD flows with complex boundaries on an unstructured Cartesian adaptive system

The numerical simulation of Magnetohydrodynamics (MHD) flows with complex boundaries has been a topic of great interest in the development of a fusion reactor blanket for the difficulty to accurately simulate the Hartmann layers and side layers along arbitrary geometries. An adaptive version of a consistent and conservative scheme has been developed for simulating the MHD flows. Besides, the present study forms the first attempt to apply the cut-cell approach for irregular wall-bounded MHD flows, which is more flexible and conveniently implemented under adaptive mesh refinement (AMR) technique. It employs a Volume-of-Fluid (VOF) approach to represent the fluid–conducting wall interface that makes it possible to solve the fluid–solid coupling magnetic problems, emphasizing at how electric field solver is implemented when conductivity is discontinuous in cut-cell. For the irregular cut-cells, the conservative interpolation technique is applied to calculate the Lorentz force at cell-center. On the other hand, it will be shown how consistent and conservative scheme is implemented on fine/coarse mesh boundaries when using AMR technique. Then, the applied numerical schemes are validated by five test simulations and excellent agreement was obtained for all the cases considered, simultaneously showed good consistency and conservative properties.

42 ENGINEERING↗
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↗
DOE OSTI2025

Carbon-conserving bioproduction of malate in an E. coli -based cell-free system

Formate, a biologically accessible form of CO 2 , has attracted interest as a renewable feedstock for bioproduction. However, approaches are needed to investigate efficient routes for biological formate assimilation due to its toxicity and limited utilization by microorganisms. Cell-free systems hold promise due to their potential for efficient use of carbon and energy sources and compatibility with diverse feedstocks. However, bioproduction using purified cell-free systems is limited by costly enzyme purification, whereas lysate-based systems must overcome loss of flux to background reactions in the cell extract. Here, we engineer an E. coli-based system for an eight-enzyme pathway from DNA and incorporate strategies to regenerate cofactors and minimize loss of flux through background reactions. We produce the industrial di-acid malate from glycine, bicarbonate, and formate by engineering the carbon-conserving reductive TCA and formate assimilation pathways. We show that in situ regeneration of NADH drives metabolic flux towards malate, improving titer by 15-fold. Background reactions can also be reduced 6-fold by diluting the lysate following expression and introducing chemical inhibitors of competing reactions. Together, these results establish a carbon-conserving, lysate-based cell-free platform for malate production, producing 64 μM malate after 8 h. This system conserves 43 % of carbon otherwise lost as CO 2 through the TCA cycle and incorporates 0.13 mol CO 2 equivalents/mol glycine fed. Finally, techno-economic analysis of cell-free malate production from formate revealed that the high cost of lysate is a key challenge to the economic feasibility of the process, even assuming efficient cofactor recycling. This work demonstrates the capabilities of cell-free expression systems for both the prototyping of carbon-conserving pathways and the sustainable bioproduction of platform chemicals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Principles of regulatory information conservation between mouse and human

To broaden our understanding of the evolution of gene regulation mechanisms, we generated occupancy profiles for 34 orthologous transcription factors (TFs) in human–mouse erythroid progenitor, lymphoblast and embryonic stem-cell lines. By combining the genome-wide transcription factor occupancy repertoires, associated epigenetic signals, and co-association patterns, here we deduce several evolutionary principles of gene regulatory features operating since the mouse and human lineages diverged. The genomic distribution profiles, primary binding motifs, chromatin states, and DNA methylation preferences are well conserved for TF-occupied sequences. However, the extent to which orthologous DNA segments are bound by orthologous TFs varies both among TFs and with genomic location: binding at promoters is more highly conserved than binding at distal elements. Notably, occupancy-conserved TF-occupied sequences tend to be pleiotropic; they function in several tissues and also co-associate with many TFs. Lastly, single nucleotide variants at sites with potential regulatory functions are enriched in occupancy-conserved TF-occupied sequences.

60 APPLIED LIFE SCIENCES↗