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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Bottom-up derivation of conservative and dissipative interactions for coarse-grained molecular liquids with the conditional reversible work method

Molecular simulations of soft matter systems have been performed in recent years using a variety of systematically coarse-grained models. With these models, structural or thermodynamic properties can be quite accurately represented while the prediction of dynamic properties remains difficult, especially for multi-component systems. In this work, we use constraint molecular dynamics simulations for calculating dissipative pair forces which are used together with conditional reversible work (CRW) conservative forces in dissipative particle dynamics (DPD) simulations. The combined CRW-DPD approach aims to extend the representability of CRW models to dynamic properties and uses a bottom-up approach. Dissipative pair forces are derived from fluctuations of the direct atomistic forces between mapped groups. The conservative CRW potential is obtained from a similar series of constraint dynamics simulations and represents the reversible work performed to couple the direct atomistic interactions between the mapped atom groups. Neopentane, tetrachloromethane, cyclohexane, and n-hexane have been considered as model systems. These molecular liquids are simulated with atomistic molecular dynamics, coarse-grained molecular dynamics, and DPD. We find that the CRW-DPD models reproduce the liquid structure and diffusive dynamics of the liquid systems in reasonable agreement with the atomistic models when using single-site mapping schemes with beads containing five or six heavy atoms. For a two-site representation of n-hexane (3 carbons per bead), time scale separation can no longer be assumed and the DPD approach consequently fails to reproduce the atomistic dynamics.

2-2-DIMETHYLPROPANE↗
DOE OSTI2014

Momentum conserving Brownian dynamics propagator for complex soft matter fluids

We present a Galilean invariant, momentum conserving first order Brownian dynamics scheme for coarse-grained simulations of highly frictional soft matter systems. Friction forces are taken to be with respect to moving background material. The motion of the background material is described by locally averaged velocities in the neighborhood of the dissolved coarse coordinates. The velocity variables are updated by a momentum conserving scheme. The properties of the stochastic updates are derived through the Chapman-Kolmogorov and Fokker-Planck equations for the evolution of the probability distribution of coarse-grained position and velocity variables, by requiring the equilibrium distribution to be a stationary solution. We test our new scheme on concentrated star polymer solutions and find that the transverse current and velocity time auto-correlation functions behave as expected from hydrodynamics. In particular, the velocity auto-correlation functions display a long time tail in complete agreement with hydrodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Boltzmann-conserving classical dynamics in quantum time-correlation functions: “Matsubara dynamics”

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or “classical Wigner approximation”) results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N → ∞, such that the lowest normal-mode frequencies take their “Matsubara” values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of ħ{sup 2} at ħ{sup 0} (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting “Matsubara” dynamics is inherently classical (since all terms O(ħ{sup 2}) disappear from the Matsubara Liouvillian in the limit N → ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A conservative scheme of drift kinetic electrons for gyrokinetic simulation of kinetic-MHD processes in toroidal plasmas

A conservative scheme of drift kinetic electrons for gyrokinetic simulations of kinetic-magnetohydrodynamic processes in toroidal plasmas has been formulated and verified. Both vector potential and electron perturbed distribution function are decomposed into adiabatic part with analytic solution and non-adiabatic part solved numerically. The adiabatic parallel electric field is solved directly from the electron adiabatic response, resulting in a high degree of accuracy. Here, the consistency between electrostatic potential and parallel vector potential is enforced by using the electron continuity equation. Since particles are only used to calculate the non-adiabatic response, which is used to calculate the non-adiabatic vector potential through Ohm's law, the conservative scheme minimizes the electron particle noise and mitigates the cancellation problem. Linear dispersion relations of the kinetic Alfvén wave and the collisionless tearing mode in cylindrical geometry have been verified in gyrokinetic toroidal code simulations, which show that the perpendicular grid size can be larger than the electron collisionless skin depth when the mode wavelength is longer than the electron skin depth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy-conserving coupled trajectory mixed quantum–classical dynamics

The coupled-trajectory mixed quantum–classical method (CTMQC), derived from the exact factorization approach, has successfully predicted photo-chemical dynamics in a number of interesting molecules, capturing population transfer and decoherence from first principles. Furthermore, due to the approximations made, CTMQC does not guarantee energy conservation. We propose a modified algorithm, CTMQC-E, which redefines the integrated force in the coupled-trajectory term so to restore energy conservation, and demonstrate its accuracy on scattering in Tully’s extended coupling region model and photoisomerization in a retinal chromophore model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2016

Grasses suppress shoot-borne roots to conserve water during drought

Significance Grasses, whose members constitute key food and bioenergy crops worldwide, use unique developmental programs to establish the root system from the shoot. Shoot-borne crown roots originate near the soil surface and provide the main conduits through which the plant takes up water and nutrients. We show that crown root development is the major target of drought stress signaling. Water deficit-triggered crown root arrest provides an important mechanism to conserve water under drought, and this response is widely conserved across grass species. Substantial phenotypic variation exists in maize for this trait, which may be a useful target in breeding efforts to improve drought tolerance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A conserved molecular logic for neurogenesis to gliogenesis switch in the cerebral cortex

During development, neural stem cells in the cerebral cortex, also known as radial glial cells (RGCs), generate excitatory neurons, followed by production of cortical macroglia and inhibitory neurons that migrate to the olfactory bulb (OB). Understanding the mechanisms for this lineage switch is fundamental for unraveling how proper numbers of diverse neuronal and glial cell types are controlled. We and others recently showed that Sonic Hedgehog (Shh) signaling promotes the cortical RGC lineage switch to generate cortical oligodendrocytes and OB interneurons. During this process, cortical RGCs generate intermediate progenitor cells that express critical gliogenesis genes Ascl1, Egfr, and Olig2. The increased Ascl1 expression and appearance of Egfr + and Olig2 + cortical progenitors are concurrent with the switch from excitatory neurogenesis to gliogenesis and OB interneuron neurogenesis in the cortex. While Shh signaling promotes Olig2 expression in the developing spinal cord, the exact mechanism for this transcriptional regulation is not known. Furthermore, the transcriptional regulation of Olig2 and Egfr has not been explored. Here, we show that in cortical progenitor cells, multiple regulatory programs, including Pax6 and Gli3, prevent precocious expression of Olig2, a gene essential for production of cortical oligodendrocytes and astrocytes. We identify multiple enhancers that control Olig2 expression in cortical progenitors and show that the mechanisms for regulating Olig2 expression are conserved between the mouse and human. Our study reveals evolutionarily conserved regulatory logic controlling the lineage switch of cortical neural stem cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Biochemical Roles for Conserved Residues in the Bacterial Fatty Acid-binding Protein Family

Fatty acid kinase (Fak) is a ubiquitous Gram-positive bacterial enzyme consisting of an ATP-binding protein (FakA) that phosphorylates the fatty acid bound to FakB. In Staphylococcus aureus, Fak is a global regulator of virulence factor transcription and is essential for the activation of exogenous fatty acids for incorporation into phospholipids. The 1.2-Å x-ray structure of S. aureus FakB2, activity assays, solution studies, site-directed mutagenesis, and in vivo complementation were used to define the functions of the five conserved residues that define the FakB protein family (Pfam02645). The fatty acid tail is buried within the protein, and the exposed carboxyl group is bound by a Ser-93-fatty acid carboxyl-Thr-61-His-266 hydrogen bond network. The guanidinium of the invariant Arg-170 is positioned to potentially interact with a bound acylphosphate. The reduced thermal denaturation temperatures of the T61A, S93A, and H266A FakB2 mutants illustrate the importance of the hydrogen bond network in protein stability. The FakB2 T61A, S93A, and H266A mutants are 1000-fold less active in the Fak assay, and the R170A mutant is completely inactive. All FakB2 mutants form FakA(FakB2)2 complexes except FakB2(R202A), which is deficient in FakA binding. Allelic replacement shows that strains expressing FakB2 mutants are defective in fatty acid incorporation into phospholipids and virulence gene transcription. These conserved residues are likely to perform the same critical functions in all bacterial fatty acid-binding proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Domain alternation and active site remodeling are conserved structural features of ubiquitin E1

E1 enzymes for ubiquitin (Ub) and Ub-like modifiers (Ubls) harbor two catalytic activities that are required for Ub/Ubl activation: adenylation and thioester bond formation. Structural studies of the E1 for the Ubl small ubiquitin-like modifier (SUMO) revealed a single active site that is transformed by a conformational switch that toggles its competency for catalysis of these two distinct chemical reactions. Although the mechanisms of adenylation and thioester bond formation revealed by SUMO E1 structures are thought to be conserved in Ub E1, there is currently a lack of structural data supporting this hypothesis. Here, we present a structure of Schizosaccharomyces pombe Uba1 in which the second catalytic cysteine half-domain (SCCH domain) harboring the catalytic cysteine has undergone a 106° rotation that results in a completely different network of intramolecular interactions between the SCCH and adenylation domains and translocation of the catalytic cysteine 12 Å closer to the Ub C terminus compared with previous Uba1 structures. SCCH domain alternation is accompanied by conformational changes within the Uba1 adenylation domains that effectively disassemble the adenylation active site. Importantly, the structural and biochemical data suggest that domain alternation and remodeling of the adenylation active site are interconnected and are intrinsic structural features of Uba1 and that the overall structural basis for adenylation and thioester bond formation exhibited by SUMO E1 is indeed conserved in Ub E1. Finally, the mechanistic insights provided by the novel conformational snapshot of Uba1 presented in this study may guide efforts to develop small molecule inhibitors of this critically important enzyme that is an active target for anticancer therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Cell–cell adhesion in metazoans relies on evolutionarily conserved features of the α-catenin·β-catenin–binding interface

Stable tissue integrity during embryonic development relies on the function of the cadherin·catenin complex (CCC). The Caenorhabditis elegans CCC is a useful paradigm for analyzing in vivo requirements for specific interactions among the core components of the CCC, and it provides a unique opportunity to examine evolutionarily conserved mechanisms that govern the interaction between α- and β-catenin. HMP-1, unlike its mammalian homolog α-catenin, is constitutively monomeric, and its binding affinity for HMP-2/β-catenin is higher than that of α-catenin for β-catenin. A crystal structure shows that the HMP-1·HMP-2 complex forms a five-helical bundle structure distinct from the structure of the mammalian α-catenin·β-catenin complex. Deletion analysis based on the crystal structure shows that the first helix of HMP-1 is necessary for binding HMP-2 avidly in vitro and for efficient recruitment of HMP-1 to adherens junctions in embryos. HMP-2 Ser-47 and Tyr-69 flank its binding interface with HMP-1, and we show that phosphomimetic mutations at these two sites decrease binding affinity of HMP-1 to HMP-2 by 40–100-fold in vitro. In vivo experiments using HMP-2 S47E and Y69E mutants showed that they are unable to rescue hmp-2(zu364) mutants, suggesting that phosphorylation of HMP-2 on Ser-47 and Tyr-69 could be important for regulating CCC formation in C. elegans. Our data provide novel insights into how cadherin-dependent cell–cell adhesion is modulated in metazoans by conserved elements as well as features unique to specific organisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The GTPase-activating protein p120RasGAP has an evolutionarily conserved “FLVR-unique” SH2 domain

The Src homology 2 (SH2) domain has a highly conserved architecture that recognizes linear phosphotyrosine motifs and is present in a wide range of signaling pathways across different evolutionary taxa. A hallmark of SH2 domains is the arginine residue in the conserved FLVR motif that forms a direct salt bridge with bound phosphotyrosine. In this work, we solve the X-ray crystal structures of the C-terminal SH2 domain of p120RasGAP (RASA1) in its apo and peptide-bound form. We find that the arginine residue in the FLVR motif does not directly contact pTyr1087 of a bound phosphopeptide derived from p190RhoGAP; rather, it makes an intramolecular salt bridge to an aspartic acid. Unexpectedly, coordination of phosphotyrosine is achieved by a modified binding pocket that appears early in evolution. Using isothermal titration calorimetry, we find that substitution of the FLVR arginine R377A does not cause a significant loss of phosphopeptide binding, but rather a tandem substitution of R398A (SH2 position βD4) and K400A (SH2 position βD6) is required to disrupt the binding. These results indicate a hitherto unrecognized diversity in SH2 domain interactions with phosphotyrosine and classify the C-terminal SH2 domain of p120RasGAP as “FLVR-unique.”

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The effect of curvature induced broken potential vorticity conservation on drift wave turbulences

Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.

Physics↗