Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservative”

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 OSTI2026

Conserved Macromolecular Architecture of Poplar Secondary Cell Walls Revealed by ssNMR and Atomistic Modeling

The macromolecular architecture of plant secondary cell walls governs wood's mechanical and biochemical properties, yet its natural intra-species variability remains poorly characterized. Here, we combined 13C solid-state NMR (ssNMR), multivariate statistical analysis, and molecular modeling to profile nanoscale structure across 13 genetically diverse Populus trichocarpa genotypes grown in 13C-enriched atmospheres. SsNMR-derived phenotypes spanning composition, structure, mobility, and inter-polymer proximities reveal a conserved architecture, with a subtle yet coordinated variation organizing into dominant structural and secondary mobility axes. A representative atomistic model captures these features and reproduces experimental metrics. Molecular dynamics simulations support a weak but consistent positive correlation between cellulose abundance and crystalline-like order, with interior cellulose chains enriched in tg (trans-gauche) conformations without expanding crystalline cores. Together, experiment and simulation reveal a genetically buffered, broadly conserved nanoscale architecture across genotypes, where subtle fine-tuning of cellulose bundling and matrix packing balances mechanical performance with biological function.

09 BIOMASS FUELS↗
DOE OSTI2021

Ultraconserved enhancer function does not require perfect sequence conservation

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Furthermore, our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Fluctuations in heavy ion collisions and global conservation effects

Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 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 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↗
DOE OSTI2023

Identification of conserved skeletal enhancers associated with craniosynostosis risk genes

Abstract Craniosynostosis, defined by premature fusion of one or multiple cranial sutures, is a common congenital defect affecting more than 1/2000 infants and results in restricted brain expansion. Single gene mutations account for 15%–20% of cases, largely as part of a syndrome, but the majority are nonsyndromic with complex underlying genetics. We hypothesized that the two noncoding genomic regions identified by a GWAS for craniosynostosis contain distal regulatory elements for the risk genes BMPER and BMP2. To identify such regulatory elements, we surveyed conserved noncoding sequences from both risk loci for enhancer activity in transgenic Danio rerio. We identified enhancers from both regions that direct expression to skeletal tissues, consistent with the endogenous expression of bmper and bmp2. For each locus, we also found a skeletal enhancer that also contains a sequence variant associated with craniosynostosis risk. We examined the activity of each enhancer during craniofacial development and found that the BMPER-associated enhancer is active in the restricted region of cartilage closely associated with frontal bone initiation. The same enhancer is active in mouse skeletal tissues, demonstrating evolutionarily conserved activity. Using enhanced yeast one-hybrid assays, we identified transcription factors that bind each enhancer and observed differential binding between alleles, implicating multiple signaling pathways. Our findings help unveil the genetic mechanism of the two craniosynostosis risk loci. More broadly, our combined in vivo approach is applicable to many complex genetic diseases to build a link between association studies and specific genetic mechanisms.

He(何璇), Xuan Anita (ORCID:0000000313221972)↗
DOE OSTI2022

Essential amino acids in the Plant-Conserved and Class-Specific Regions of cellulose synthases

The Plant-Conserved Region (P-CR) and the Class-Specific Region (CSR) are two plant-unique sequences in the catalytic core of cellulose synthases (CESAs) for which specific functions have not been established. Here, we used site-directed mutagenesis to replace amino acids and motifs within these sequences predicted to be essential for assembly and function of CESAs. We developed an in vivo method to determine the ability of mutated CesA1 transgenes to complement an Arabidopsis (Arabidopsis thaliana) temperature-sensitive root-swelling1 (rsw1) mutant. Replacement of a Cys residue in the CSR, which blocks dimerization in vitro, rendered the AtCesA1 transgene unable to complement the rsw1 mutation. Examination of the CSR sequences from 33 diverse angiosperm species showed domains of high-sequence conservation in a class-specific manner but with variation in the degrees of disorder, indicating a nonredundant role of the CSR structures in different CESA isoform classes. The Cys residue essential for dimerization was not always located in domains of intrinsic disorder. Expression of AtCesA1 transgene constructs, in which Pro417 and Arg453 were substituted for Ala or Lys in the coiled-coil of the P-CR, were also unable to complement the rsw1 mutation. Despite an expected role for Arg457 in trimerization of CESA proteins, AtCesA1 transgenes with Arg457Ala mutations were able to fully restore the wild-type phenotype in rsw1. Our data support that Cys662 within the CSR and Pro417 and Arg453 within the P-CR of Arabidopsis CESA1 are essential residues for functional synthase complex formation, but our data do not support a specific role for Arg457 in trimerization in native CESA complexes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Phylogenetic conservation of soil bacterial responses to simulated global changes.

Soil bacterial communities are altered by anthropogenic drivers such as climate change-related warming and fertilization. However, we lack a predictive understanding of how bacterial communities respond to such global changes. Here, we tested whether phylogenetic information might be more predictive of the response of bacterial taxa to some forms of global change than others. We analysed the composition of soil bacterial communities from perturbation experiments that simulated warming, drought, elevated CO2 concentration and phosphorus (P) addition. Bacterial responses were phylogenetically conserved to all perturbations. The phylogenetic depth of these responses varied minimally among the types of perturbations and was similar when merging data across locations, implying that the context of particular locations did not affect the phylogenetic pattern of response. We further identified taxonomic groups that responded consistently to each type of perturbation. These patterns revealed that, at the level of family and above, most groups responded consistently to only one or two types of perturbations, suggesting that traits with different patterns of phylogenetic conservation underlie the responses to different perturbations. We conclude that a phylogenetic approach may be useful in predicting how soil bacterial communities respond to a variety of global changes. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Filling constraints on translation invariant dipole conserving systems

Systems with conserved dipole moment have drawn considerable interest in light of their realization in recent experiments on tilted optical lattices. Here, an important question for such systems is delineating the conditions under which they admit a unique gapped ground state that is consistent with all symmetries. Here, we study one-dimensional translation-invariant lattices that conserve U(1) charge and Z L dipole moment, where discreteness of the dipole symmetry is enforced by periodic boundary conditions, with L the system size. We show that in these systems, a symmetric, gapped, and non-degenerate ground state requires not only integer charge filling, but also a fixed value of the dipole filling, while other fractional dipole fillings enforce either a gapless or symmetry-breaking ground state. In contrast with prior results in the literature, we find that the dipole filling constraint depends both on the charge filling as well as the system size, emphasizing the subtle interplay of dipole symmetry with boundary conditions. We support our results with numerical simulations and exact results.

1-dimensional systems↗
DOE OSTI2023

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗